Changeset 25433 for branches/eam_branches/20090715/psphot/src/psphotRadialProfileByAngle.c

Timestamp:

Sep 17, 2009, 2:26:32 PM (17 years ago)

Author:

eugene

Message:

change radius for extended sources to use footprint; clean up some of the visualizations; plug some leaks

File:

• branches/eam_branches/20090715/psphot/src/psphotRadialProfileByAngle.c (modified) (4 diffs)

branches/eam_branches/20090715/psphot/src/psphotRadialProfileByAngle.c

@@ -1 +1 @@
 # include "psphotInternal.h"
 
-// Given a source at (x,y), generate a collection of radial profiles at even angular
-// separations
-
+// Given a source at (x,y), generate a collection of radial profiles at even angular separations
+
+// These functions are used to calculate the stats in a rectangle at arbitrary orientation.
+// XXX Move these elsewhere (psLib?)
 float psphotMeanSectorValue (psImage *image, float x, float y, float dL, float dW, float theta);
 psVector *psphotBoxValues (psImage *image, float x0, float y0, float dL, float dW, float theta);
@@ -9 +10 @@
 psVector *psphotLineValuesBresen (psImage *image, int X1, int Y1, int X2, int Y2, int dW, int swapcoords);
 
-bool psphotRadialProfilesByAngles (pmPetrosian *petrosian, pmSource *source, int Nsec, float Rmax) {
+bool psphotRadialProfilesByAngles (pmSource *source, pmPetrosian *petrosian, int Nsec, float Rmax) {
 
     // we want to have an even number of sectors so we can do 180 deg symmetrizing
@@ -124 +125 @@
 
     psVector *values = psphotBoxValues (image, x, y, dL, dW, theta);
+    if (!values) goto escape;
+    if (!values->n) goto escape;
     
     psStats *stats = psStatsAlloc(PS_STAT_SAMPLE_MEDIAN);
@@ -134 +137 @@
     
     return value;
+
+escape:
+    psFree(values);
+    return NAN;    
 }