X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=tools%2FclitkImageUncertainty.cxx;h=63b1f12173a161647d64043f6fa986fe5d1248bc;hb=d4377ec80fd4e6f0156bb0673bb661acbf5b8cec;hp=9240bd8278161746245374a99c79e47d32254ef6;hpb=a845c6a04ecfeedc8eced5f4051e06ee88c59e90;p=clitk.git diff --git a/tools/clitkImageUncertainty.cxx b/tools/clitkImageUncertainty.cxx index 9240bd8..63b1f12 100644 --- a/tools/clitkImageUncertainty.cxx +++ b/tools/clitkImageUncertainty.cxx @@ -34,13 +34,12 @@ #include "clitkImageUncertainty_ggo.h" #include "clitkImageCommon.h" #include "clitkCommon.h" +#include "clitkPortability.h" // itk include #include "itkImageRegionConstIterator.h" #include "itkImageRegionIterator.h" -#include - //==================================================================== int main(int argc, char * argv[]) { @@ -79,7 +78,7 @@ int main(int argc, char * argv[]) { double squared = pii.Get(); double mean = pi.Get(); double uncert = sqrt((NumberOfEvents*squared - mean*mean) / ((NumberOfEvents-1)*(mean*mean))); - if (uncert != uncert || uncert == std::numeric_limits::infinity() || uncert == -std::numeric_limits::infinity()) uncert = 1.; + if (!IsNormal(uncert)) uncert = 1.; po.Set(uncert); ++pi; ++pii;