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 #70637 is a reply to message #70636] Mon, 26 April 2010 23:24 Go to previous messageGo to previous message
Aram Panasenco is currently offline  Aram Panasenco
Messages: 41
Registered: April 2010
Member
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
[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: Sat Oct 11 04:03:37 PDT 2025

Total time taken to generate the page: 0.48047 seconds