Changeset 25452 for branches/eam_branches/20090715/psphot/src/psphotRadialProfileByAngles.c

Timestamp:

Sep 19, 2009, 8:49:41 AM (17 years ago)

Author:

eugene

Message:

further work on the petrosian analysis & extended sources

File:

• branches/eam_branches/20090715/psphot/src/psphotRadialProfileByAngles.c (moved) (moved from branches/eam_branches/20090715/psphot/src/psphotRadialProfileByAngle.c ) (4 diffs)

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

@@ -10 +10 @@
 psVector *psphotLineValuesBresen (psImage *image, int X1, int Y1, int X2, int Y2, int dW, int swapcoords);
 
-bool psphotRadialProfilesByAngles (pmSource *source, pmPetrosian *petrosian, int Nsec, float Rmax) {
+bool psphotRadialProfilesByAngles (pmSource *source, int Nsec, float Rmax) {
 
     // we want to have an even number of sectors so we can do 180 deg symmetrizing
@@ -16 +16 @@
     float dtheta = 2.0*M_PI / Nsec;
 
-    psFree(petrosian->radii);
-    psFree(petrosian->fluxes);
-    psFree(petrosian->theta);
-
-    petrosian->radii = psArrayAllocEmpty(Nsec);
-    petrosian->fluxes = psArrayAllocEmpty(Nsec);
-    petrosian->theta = psVectorAllocEmpty(Nsec, PS_TYPE_F32);
+    pmSourceRadialProfile *profile = source->extpars->profile;
+    psFree(profile->radii);
+    psFree(profile->fluxes);
+    psFree(profile->theta);
+
+    profile->radii = psArrayAllocEmpty(Nsec);
+    profile->fluxes = psArrayAllocEmpty(Nsec);
+    profile->theta = psVectorAllocEmpty(Nsec, PS_TYPE_F32);
 
 
@@ -72 +73 @@
         }
 
-        psArrayAdd (petrosian->radii, 100, radius);
-        psArrayAdd (petrosian->fluxes, 100, flux);
-        psVectorAppend (petrosian->theta, theta);
+        psArrayAdd (profile->radii, 100, radius);
+        psArrayAdd (profile->fluxes, 100, flux);
+        psVectorAppend (profile->theta, theta);
 
         // psphotPetrosianVisualProfileByAngle (radius, flux);
@@ -84 +85 @@
     for (int i = 0; i < Nsec / 2; i++) {
 
-        psVector *r1 = petrosian->radii->data[i];
-        psVector *r2 = petrosian->radii->data[i+Nsec/2];
-
-        psVector *f1 = petrosian->fluxes->data[i];
-        psVector *f2 = petrosian->fluxes->data[i+Nsec/2];
+        psVector *r1 = profile->radii->data[i];
+        psVector *r2 = profile->radii->data[i+Nsec/2];
+
+        psVector *f1 = profile->fluxes->data[i];
+        psVector *f2 = profile->fluxes->data[i+Nsec/2];
 
         psAssert (r1->n == r2->n, "mis-matched vectors");