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

Home » Public Forums » archive » Re: PSF Energy inside circle
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: PSF Energy inside circle [message #61576] Wed, 23 July 2008 13:52 Go to previous message
wlandsman is currently offline  wlandsman
Messages: 743
Registered: June 2000
Senior Member
On Jul 23, 2:03 pm, Michael Aye <kmichael....@googlemail.com> wrote:
> Dear all,

> So how could I find and integrate the next "ring" of pixels? How would
> I even calculate the ever growing circumference correctly, taking into
> account that I have to sum up ever more pixels?

The options mentioned so far seem fine for large radii. For small
radii, the use of integer pixel values can introduce "noise" since a
pixel must be either entirely within or outside the circle. The
program aper.pro
( http://idlastro.gsfc.nasa.gov/ftp/pro/idlphot/aper.pro )
sums using partial pixels, using the exact overlap area of a pixel
with a circle.

The program is somewhat awkward to use, in part because it is very
old, and in part because it propagates error sources. But if you
want to compute the flux within 10 concentric circles centered at the
position [500.2, 500.5] with radii of 3, 3.5, 4, 4.5, ... 7.5 then the
following is a start:

IDL> apr = 3. + findgen(10)/2.
IDL> aper,im, 500.2, 500.5, flux, eflux, sky,skyerr, 1, apr,/
flux,setsky=0
IDL> plot, apr, flux ;plot encircled flux vs. radius

In this case I force the "sky" (background) to have a value of
zero. You can have aper.pro compute the background by giving it an
inner and outer "sky" radii far from the central source. --Wayne
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: plot stddev bars for each data point
Next Topic: Re: Newbie question concerning summations/loops in IDL

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

Current Time: Thu Oct 09 14:27:15 PDT 2025

Total time taken to generate the page: 0.95619 seconds