[Insight-users] need help image histogram

[Luis Ibanez]

Hi Hossein,


    Why are you choosing "4096" as the maximum value ?


Can you please run the MaximumMinimumImageCalculator in your
image and post back to the list the maximum and minimum values
that you get ?


    It is unlikely that 4096 is actually the maximum value.


You may not be considering the fact that the pixel intensities
may be linearly rescaled by the slope and intercept values
indicated in the DICOM header.


Please let us know what your find,


    Thanks


       Luis


-------------------
hossein nazem wrote:
> 
>  
> 
> Dear Luis
> 
> Hi,
> 
>  
> 
> Thank you for your consideration. It seems that a mapping take place 
> somewhere and convert data from uint16 to unit8.
> 
> These are my parameters:
> 
>  
> 
>   histogramGenerator->SetNumberOfBins( 4096 );
> 
>   histogramGenerator->SetMarginalScale( 1 );
> 
>  
> 
>   histogramGenerator->SetHistogramMin(  -0.5 );
> 
>   histogramGenerator->SetHistogramMax( 4096 );
> 
>  
> 
> I expect that I have frequency of nearly 1000 around bin no 1000 but its 
> 0 and I expect that some of frequency match the total sum of pixels but 
> its nearly half of the total sum.
> 
> I need the real histogram (0-4096).
> 
> And here is how I check the result for each bin (actually I use the 
> examples I mentioned before):
> 
>  
> 
> std::cout << histogram->GetFrequency( bin, 0 ) << std::endl;
> 
>  
> 
> And the total number of pixels must be equal to:
> 
>  
> 
> s=s+histogram->GetFrequency( bin, 0 );
> 
>  
> 
> Yours truly,
> 
> Hossein
> 
>  
>  
>  
>  
>  
>  
>  
>  
>  
> 
> --- On *Tue, 9/9/08, Luis Ibanez /<luis.ibanez at kitware.com>/* wrote:
> 
>     From: Luis Ibanez <luis.ibanez at kitware.com>
>     Subject: Re: [Insight-users] need help image histogram
>     To: hn_6112 at yahoo.com
>     Cc: insight-users at itk.org
>     Date: Tuesday, September 9, 2008, 9:22 PM
> 
> Hi Hossein,
> 
> 
> What values did you set up for the Max and Min of the Histogram ?
> 
> 
>    min.Fill(   -0.5 );  ??
>    max.Fill(  255.5 );  ??
> 
>    generator->SetHistogramMin( min );
>    generator->SetHistogramMax( max );
> 
> 
> If the total number of entries in the Histogram doesn't match the
> number of pixels in the image, then it is likely that your Min and
> Max setting in the histogram calculator are defining an intensity
> range smaller than the actual dynamic range of your DICOM image.
> 
> 
> Please let us know,
> 
> 
>    Thanks
> 
> 
>       Luis
> 
> 
> -------------------
> hossein nazem wrote:
>> Hi,
>> 
>> I’m trying to extract the histogram of a dicom image. The result is 
>> different to what I achieve from MATLAB. There is something wrong in itk 
>> because total pixels are fix and 512*512=262144 but in my code it’s 
>> variable when I change the bins,  and the shape of histogram is 
>> different comparing with MATLAB output. I work on a 512*512 dicom uint12 
>> image. I tried the Examples/Statistics/ImageHistogram1.cxx & 
>> Examples/Statistics/ImageHistogram2.cxx.
>> 
>>  
>> 
>> Yours truly,
>> 
>> Hossein
>> 
>> 
>> 
>> ------------------------------------------------------------------------
>> 
>> _______________________________________________
>> Insight-users mailing list
>> Insight-users at itk.org
>> http://www.itk.org/mailman/listinfo/insight-users
> 
>