![]() |
Institut für Astronomie und AstrophysikAbteilung AstronomieWaldhäuser Str. 64, D-72076 Tübingen, Germany |
![]() |
PIXWT
Circle-rectangle overlap area computation.
Compute the fraction of a unit pixel that is interior to a circle.
The circle has a radius r and is centered at (xc, yc). The center of
the unit pixel (length of sides = 1) is at (x, y).
CCD data processing
area = Pixwt( xc, yc, r, x, y )
xc, yc : Center of the circle, numeric scalars
r : Radius of the circle, numeric scalars
x, y : Center of the unit pixel, numeric scalar or vector
None.
None.
Function value: Computed overlap area.
What is the area of overlap of a circle with radius 3.44 units centered
on the point 3.23, 4.22 with the pixel centered at [5,7]
IDL> print,pixwt(3.23,4.22,3.44,5,7) ==> 0.6502
None.
Divides the circle and rectangle into a series of sectors and
triangles. Determines which of nine possible cases for the
overlap applies and sums the areas of the corresponding sectors
and triangles. Called by aper.pro
If improved speed is needed then a C version of this routines, with
notes on how to linkimage it to IDL is available at
ftp://ftp.lowell.edu/pub/buie/idl/custom/
Ported by Doug Loucks, Lowell Observatory, 1992 Sep, from the
routine pixwt.c, by Marc Buie.
[Home Page] [Software, Documentation] [IDL Documentation] [Quick Reference] [Feedback]