Changeset 26669 for branches/eam_branches/20091201/psModules/src/objects/pmPSF.c

Timestamp:

Jan 22, 2010, 11:30:32 AM (16 years ago)

Author:

Paul Price

Message:

Add function to return FWHM from PSF.

File:

• branches/eam_branches/20091201/psModules/src/objects/pmPSF.c (modified) (2 diffs)

branches/eam_branches/20091201/psModules/src/objects/pmPSF.c

@@ -42 +42 @@
 #include "pmErrorCodes.h"
 
+
+#define MAX_AXIS_RATIO 20.0             // Maximum axis ratio for PSF model
+
 /*****************************************************************************/
 /* FUNCTION IMPLEMENTATION - PUBLIC                                          */
@@ -405 +408 @@
     return psf;
 }
+
+
+float pmPSFtoFWHM(const pmPSF *psf, float x, float y)
+{
+    PS_ASSERT_PTR_NON_NULL(psf, NAN);
+
+    pmModel *model = pmModelFromPSFforXY(psf, x, y, 1.0); // Model of source
+    psF32 *params = model->params->data.F32; // Model parameters
+    psEllipseAxes axes = pmPSF_ModelToAxes(params, MAX_AXIS_RATIO); // Ellipse axes
+
+    // Curiously, the minor axis can be larger than the major axis, so need to check.
+    float fwhm = 2.355 * PS_MAX(axes.minor, axes.major); // FWHM, converted from sigma
+
+    psFree(model);
+
+    return fwhm;
+}