comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: How to find the pixel position
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: How to find the pixel position [message #70697 is a reply to message #70637] Wed, 28 April 2010 23:40 Go to previous message
sid is currently offline  sid
Messages: 50
Registered: January 1995
Member
On Apr 27, 11:24 am, Aram Panasenco <panasencoa...@gmail.com> wrote:
> sid wrote:
>> Hi,
>>      My data is in fits format. The is of 1024 * 1024 array. The counts
>> vary from 5000 to 6000 and I know that 5500 counts is there in my
>> data, but I need to know at which pixel this 5500 counts occur
>> exactly, without displaying the image, because I need to do this for
>> several files. So each time I can't display and check for the pixel
>> position. please helpout in this regard.
>> regards
>> sid
>
> I think what you are saying (correct me if I am wrong) is that you have
> a 1024x1024 array, and you want to find where the pixel values are equal
> to 5500.
>
> You can use the WHERE function:
>
> fitsData = readfits('filename.fits')
> countValue = 5500
>
> findIndices = where(fitsData eq countValue)
>
> Note that the WHERE function returns one-dimensional subscripts. You can
> convert them back to two-dimensional subscripts (if you need to) using
> the ARRAY_INDICES function:
>
> rectIndices = array_indices([1024,1024],findIndices,/dimensions)
>
> Cheers
> ~Aram Panasenco

Hi,
I did like this
raw=readfits('filename.fits')
b=where(raw eq 2832.90)
I know that it occurs at raw(5,5)
so now if I do
print,b
it should print 5, since where function returns one dimensional
subscripts.(am I right, correct me if it is wrong)
but instead it is printing -1. Please help me out.
regards
sid
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: IDL Page Downloads
Next Topic: Re: Regarding " WRITING DICOM IMAGES"

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 17 18:12:45 PDT 2025

Total time taken to generate the page: 0.53997 seconds