aboutsummaryrefslogtreecommitdiff
path: root/noao/digiphot/apphot/center/aprefitcen.x
diff options
context:
space:
mode:
Diffstat (limited to 'noao/digiphot/apphot/center/aprefitcen.x')
-rw-r--r--noao/digiphot/apphot/center/aprefitcen.x193
1 files changed, 193 insertions, 0 deletions
diff --git a/noao/digiphot/apphot/center/aprefitcen.x b/noao/digiphot/apphot/center/aprefitcen.x
new file mode 100644
index 00000000..c50e2efb
--- /dev/null
+++ b/noao/digiphot/apphot/center/aprefitcen.x
@@ -0,0 +1,193 @@
+include <mach.h>
+include "../lib/apphotdef.h"
+include "../lib/apphot.h"
+include "../lib/noisedef.h"
+include "../lib/centerdef.h"
+include "../lib/center.h"
+
+define CONVERT .424660900 # conversion factor from fwhmpsf to sigma
+
+# APREFITCENTER -- Procedure to refit the centers assuming that the appropriate
+# pixel buffer is in memory. See apfitcenter for further information.
+
+int procedure aprefitcenter (ap, im, ier)
+
+pointer ap # pointer to the apphot structure
+pointer im # the input image descriptor
+int ier # previous error code
+
+int fier
+pointer cen, nse
+int ap_ctr1d(), ap_mctr1d(), ap_gctr1d(), ap_lgctr1d()
+
+begin
+ cen = AP_PCENTER(ap)
+ nse = AP_NOISE(ap)
+
+ # Initialize
+ AP_XCENTER(cen) = AP_CXCUR(cen)
+ AP_YCENTER(cen) = AP_CYCUR(cen)
+ AP_XSHIFT(cen) = 0.0
+ AP_YSHIFT(cen) = 0.0
+ AP_OXSHIFT(cen) = 0.0
+ AP_OYSHIFT(cen) = 0.0
+ AP_XERR(cen) = INDEFR
+ AP_YERR(cen) = INDEFR
+
+ # Return if the center is undefined.
+ if (IS_INDEFR(AP_CXCUR(cen)) || IS_INDEFR(AP_CYCUR(cen))) {
+ AP_OXINIT(cen) = INDEFR
+ AP_OYINIT(cen) = INDEFR
+ AP_OXCENTER(cen) = INDEFR
+ AP_OYCENTER(cen) = INDEFR
+ return (AP_CTR_NOAREA)
+ }
+
+ # Convert the coordinates.
+ switch (AP_WCSOUT(ap)) {
+ case WCS_WORLD, WCS_PHYSICAL:
+ call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
+ AP_OYINIT(cen), 1)
+ call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXCENTER(cen),
+ AP_OYCENTER(cen), 1)
+ case WCS_TV:
+ call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
+ AP_OYINIT(cen), 1)
+ call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXCENTER(cen),
+ AP_OYCENTER(cen), 1)
+ default:
+ AP_OXINIT(cen) = AP_CXCUR(cen)
+ AP_OYINIT(cen) = AP_CYCUR(cen)
+ AP_OXCENTER(cen) = AP_CXCUR(cen)
+ AP_OYCENTER(cen) = AP_CYCUR(cen)
+ }
+
+ # Return input coordinates if no center fitting.
+ if (AP_CENTERFUNCTION(cen) == AP_NONE)
+ return (AP_OK)
+
+ # Choose the centering algorithm.
+ switch (AP_CENTERFUNCTION(cen)) {
+
+ case AP_CENTROID1D:
+
+ # Compute the x and y centroids.
+ if (AP_CTHRESHOLD(cen) <= 0.0) {
+ fier = ap_mctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen),
+ AP_CNY(cen), AP_EPADU(nse), AP_XCENTER(cen),
+ AP_YCENTER(cen), AP_XERR(cen), AP_YERR(cen))
+ if (IS_INDEFR(AP_XERR(cen)))
+ AP_XCENTER(cen) = AP_CXC(cen)
+ if (IS_INDEFR(AP_YERR(cen)))
+ AP_YCENTER(cen) = AP_CYC(cen)
+ } else {
+ fier = ap_ctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen),
+ AP_CNY(cen), AP_EPADU(nse), AP_XCENTER(cen),
+ AP_YCENTER(cen), AP_XERR(cen), AP_YERR(cen))
+ if (IS_INDEFR(AP_XERR(cen)))
+ AP_XCENTER(cen) = AP_CXC(cen)
+ if (IS_INDEFR(AP_YERR(cen)))
+ AP_YCENTER(cen) = AP_CYC(cen)
+ }
+ AP_XCENTER(cen) = AP_XCENTER(cen) + AP_CXCUR(cen) - AP_CXC(cen)
+ AP_YCENTER(cen) = AP_YCENTER(cen) + AP_CYCUR(cen) - AP_CYC(cen)
+ AP_XSHIFT(cen) = AP_XCENTER(cen) - AP_CXCUR(cen)
+ AP_YSHIFT(cen) = AP_YCENTER(cen) - AP_CYCUR(cen)
+
+ case AP_GAUSS1D:
+
+ # Compute the 1D Gaussian x and y centers.
+ fier = ap_gctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen), AP_CNY(cen),
+ CONVERT * AP_FWHMPSF(ap) * AP_SCALE(ap), AP_CMAXITER(cen),
+ AP_XCENTER(cen), AP_YCENTER(cen), AP_XERR(cen), AP_YERR(cen))
+ AP_XCENTER(cen) = AP_XCENTER(cen) + AP_CXCUR(cen) - AP_CXC(cen)
+ AP_YCENTER(cen) = AP_YCENTER(cen) + AP_CYCUR(cen) - AP_CYC(cen)
+ AP_XSHIFT(cen) = AP_XCENTER(cen) - AP_CXCUR(cen)
+ AP_YSHIFT(cen) = AP_YCENTER(cen) - AP_CYCUR(cen)
+
+ case AP_OFILT1D:
+
+ # Compute the Goad 1D x and y centers.
+ fier = ap_lgctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen), AP_CNY(cen),
+ AP_CXC(cen), AP_CYC(cen), CONVERT * AP_FWHMPSF(ap) *
+ AP_SCALE(ap), AP_CMAXITER(cen), AP_EPADU(nse),
+ AP_SKYSIGMA(nse), AP_XCENTER(cen), AP_YCENTER(cen),
+ AP_XERR(cen), AP_YERR(cen))
+ AP_XCENTER(cen) = AP_XCENTER(cen) + AP_CXCUR(cen) - AP_CXC(cen)
+ AP_YCENTER(cen) = AP_YCENTER(cen) + AP_CYCUR(cen) - AP_CYC(cen)
+ AP_XSHIFT(cen) = AP_XCENTER(cen) - AP_CXCUR(cen)
+ AP_YSHIFT(cen) = AP_YCENTER(cen) - AP_CYCUR(cen)
+
+ default:
+
+ # do nothing gracefully
+ }
+
+ # Convert the coordinates.
+ switch (AP_WCSOUT(ap)) {
+ case WCS_WORLD, WCS_PHYSICAL:
+ call ap_ltoo (ap, AP_XCENTER(cen), AP_YCENTER(cen),
+ AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
+ call ap_ltoo (ap, AP_XCENTER(cen) - AP_XSHIFT(cen),
+ AP_YCENTER(cen) - AP_YSHIFT(cen), AP_OXINIT(cen),
+ AP_OYINIT(cen), 1)
+ AP_OXSHIFT(cen) = AP_OXCENTER(cen) - AP_OXINIT(cen)
+ AP_OYSHIFT(cen) = AP_OYCENTER(cen) - AP_OYINIT(cen)
+ case WCS_TV:
+ call ap_ltov (im, AP_XCENTER(cen), AP_YCENTER(cen),
+ AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
+ call ap_ltov (im, AP_XCENTER(cen) - AP_XSHIFT(cen),
+ AP_YCENTER(cen) - AP_YSHIFT(cen), AP_OXINIT(cen),
+ AP_OYINIT(cen), 1)
+ AP_OXSHIFT(cen) = AP_OXCENTER(cen) - AP_OXINIT(cen)
+ AP_OYSHIFT(cen) = AP_OYCENTER(cen) - AP_OYINIT(cen)
+ default:
+ AP_OXCENTER(cen) = AP_XCENTER(cen)
+ AP_OYCENTER(cen) = AP_YCENTER(cen)
+ AP_OXINIT(cen) = AP_XCENTER(cen) - AP_XSHIFT(cen)
+ AP_OYINIT(cen) = AP_YCENTER(cen) - AP_YSHIFT(cen)
+ AP_OXSHIFT(cen) = AP_XSHIFT(cen)
+ AP_OYSHIFT(cen) = AP_YSHIFT(cen)
+ }
+
+ # Return appropriate error code.
+ if (fier != AP_OK) {
+ AP_XCENTER(cen) = AP_CXCUR(cen)
+ AP_YCENTER(cen) = AP_CYCUR(cen)
+ AP_XSHIFT(cen) = 0.0
+ AP_YSHIFT(cen) = 0.0
+ AP_XERR(cen) = INDEFR
+ AP_YERR(cen) = INDEFR
+ switch (AP_WCSOUT(ap)) {
+ case WCS_WORLD, WCS_PHYSICAL:
+ call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
+ AP_OYINIT(cen), 1)
+ call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen),
+ AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
+ case WCS_TV:
+ call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
+ AP_OYINIT(cen), 1)
+ call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen),
+ AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
+ default:
+ AP_OXCENTER(cen) = AP_CXCUR(cen)
+ AP_OYCENTER(cen) = AP_CYCUR(cen)
+ AP_OXINIT(cen) = AP_CXCUR(cen)
+ AP_OYINIT(cen) = AP_CYCUR(cen)
+ }
+ AP_OXSHIFT(cen) = 0.0
+ AP_OYSHIFT(cen) = 0.0
+ return (fier)
+ } else if (ier == AP_CTR_BADDATA) {
+ return (AP_CTR_BADDATA)
+ } else if (ier == AP_CTR_LOWSNRATIO) {
+ return (AP_CTR_LOWSNRATIO)
+ } else if (abs (AP_XSHIFT(cen)) > (AP_MAXSHIFT(cen) * AP_SCALE(ap))) {
+ return (AP_CTR_BADSHIFT)
+ } else if (abs (AP_YSHIFT(cen)) > (AP_MAXSHIFT(cen) * AP_SCALE(ap))) {
+ return (AP_CTR_BADSHIFT)
+ } else if (ier == AP_CTR_OUTOFBOUNDS) {
+ return (AP_CTR_OUTOFBOUNDS)
+ } else
+ return (AP_OK)
+end
generated by cgit (git 2.34.1) at 2025-09-20 11:27:38 -0400