Changeset 25105 for branches/eam_branches/20090715/psphot/src/psphotPetrosianProfile.c

Timestamp:

Aug 18, 2009, 6:24:53 AM (17 years ago)

Author:

eugene

Message:

updates for petrosian analysis study

Location:

branches/eam_branches/20090715/psphot/src

Files:

• . (modified) (1 prop)
• psphotPetrosianProfile.c (modified) (3 diffs)

branches/eam_branches/20090715/psphot/src

@@ -18 +18 @@
 psphotVersionDefinitions.h
 psphotMomentsStudy
+psphotPetrosianStudy

@@ -18 +18 @@
 psphotVersionDefinitions.h
 psphotMomentsStudy
+psphotPetrosianStudy

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

@@ -1 +1 @@
 # include "psphotInternal.h"
+
+// generate the Petrosian radius and flux using elliptical contours
+
+// XXX much of this function is focused on generating the clean contours, which can be used by 
+// any number of aperture-like measurements.  probably will want to rename the pmPetrosian
+// structure to something the pmRadialProfile
 
 bool psphotPetrosianProfile (pmSource *source) {
 
+  // container to hold results from the radial profile analysis
   pmPetrosian *petrosian = pmPetrosianAlloc();
 
-  if (!psphotRadialProfilesByAngle (petrosian, source, Nsec, Rmax)) {
+  int Nsec = 24;
+  float Rmax = 64;
+  float fluxMin = 0.0;
+  float fluxMax = 1000.0;
+
+  // generate a series of radial profiles at Nsec evenly spaced angles.  the profile flux
+  // is measured by interpolation for small radii; for large radii, the pixels in a box
+  // are averaged to increase the S/N (XXX not yet done)
+  if (!psphotRadialProfilesByAngles (petrosian, source, Nsec, Rmax)) {
     psError (PS_ERR_UNKNOWN, false, "failed to measure radial profile for petrosian");
     return false;
   }
 
+  // use the radial profiles to determine the radius of a given isophote.  this isophote
+  // is used to determine the elliptical shape of the object, so it has a relatively high
+  // value (nominally 50% of the peak)
   if (!psphotRadiiFromProfiles (petrosian, fluxMin, fluxMax)) {
     psError (PS_ERR_UNKNOWN, false, "failed to measure isophotal radii from profiles");
@@ -15 +33 @@
   }
 
+  // convert the isophotal radius vs angle measurements to an elliptical contour
   if (!psphotEllipticalContour (petrosian)) {
     psError (PS_ERR_UNKNOWN, false, "failed to measure elliptical contour");
@@ -20 +39 @@
   }
   
+  // generate a single, normalized radial profile following the elliptical contours.
+  // the radius is normalized by the axis ratio so that on the major axis, 1 pixel = 1 pixel
+  if (!psphotEllipticalProfile (source, petrosian)) {
+    psError (PS_ERR_UNKNOWN, false, "failed to generate elliptical profile");
+    return false;
+  }
+  
+  // integrate the radial profile for radial bins defined for the petrosian measurement:
+  // SB_i (r_i) where \alpha r_i < r < \beta r_i
+  if (!psphotPetrosianRadialBins (source, petrosian, Rmax)) {
+    psError (PS_ERR_UNKNOWN, false, "failed to generate elliptical profile");
+    return false;
+  }
+  
+  // use the SB_i from above to calculate the petrosian radius and the flux within that radius
+  if (!psphotPetrosianStats (petrosian)) {
+    psError (PS_ERR_UNKNOWN, false, "failed to generate elliptical profile");
+    return false;
+  }
+  
+  psphotPetrosianFreeVectors(petrosian);
 
+  fprintf (stderr, "petrosian radius: %f, flux: %f, axis ratio: %f, angle: %f\n",
+           petrosian->petrosianRadius, petrosian->petrosianFlux, petrosian->axes.minor/petrosian->axes.major, PS_DEG_RAD*petrosian->axes.theta);
+
+  return true;
 }

@@ -18 +18 @@
 psphotVersionDefinitions.h
 psphotMomentsStudy
+psphotPetrosianStudy