Changeset 7326
- Timestamp:
- Jun 2, 2006, 6:04:10 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotApResid.c (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotApResid.c
r6918 r7326 98 98 if (!strcasecmp (ApTrendOption, "CONSTANT")) { 99 99 stats->clipIter = 2; 100 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);101 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 102 103 stats->clipIter = 3; 104 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);100 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 101 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 102 103 stats->clipIter = 3; 104 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 105 105 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 106 106 } … … 110 110 // first clip out objects which are too far from the median 111 111 stats->clipIter = 2; 112 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);113 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 114 115 // apply the fit 116 stats->clipIter = 3; 117 pmPSF_MaskApTrend (psf , PM_PSF_SKYBIAS);112 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 113 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 114 115 // apply the fit 116 stats->clipIter = 3; 117 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKYBIAS); 118 118 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 119 119 } … … 122 122 // first clip out objects which are too far from the median 123 123 stats->clipIter = 2; 124 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);125 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 126 127 // apply the fit 128 stats->clipIter = 2; 129 pmPSF_MaskApTrend (psf , PM_PSF_SKYBIAS);130 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 131 132 // apply the fit 133 stats->clipIter = 3; 134 pmPSF_MaskApTrend (psf , PM_PSF_SKYSAT);124 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 125 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 126 127 // apply the fit 128 stats->clipIter = 2; 129 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKYBIAS); 130 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 131 132 // apply the fit 133 stats->clipIter = 3; 134 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKYSAT); 135 135 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 136 136 } … … 140 140 // first clip out objects which are too far from the median 141 141 stats->clipIter = 2; 142 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);143 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 144 145 // apply the fit 146 stats->clipIter = 3; 147 pmPSF_MaskApTrend (psf , PM_PSF_XY_LIN);142 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 143 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 144 145 // apply the fit 146 stats->clipIter = 3; 147 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_XY_LIN); 148 148 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 149 149 } … … 153 153 // first clip out objects which are too far from the median 154 154 stats->clipIter = 2; 155 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);156 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 157 158 // apply the fit 159 stats->clipIter = 3; 160 pmPSF_MaskApTrend (psf , PM_PSF_XY_QUAD);155 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 156 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 157 158 // apply the fit 159 stats->clipIter = 3; 160 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_XY_QUAD); 161 161 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 162 162 } … … 166 166 // first clip out objects which are too far from the median 167 167 stats->clipIter = 2; 168 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);169 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 170 171 // apply the fit 172 stats->clipIter = 3; 173 pmPSF_MaskApTrend (psf , PM_PSF_SKY_XY_LIN);168 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 169 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 170 171 // apply the fit 172 stats->clipIter = 3; 173 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKY_XY_LIN); 174 174 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 175 175 } … … 179 179 // first clip out objects which are too far from the median 180 180 stats->clipIter = 2; 181 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);182 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 183 184 // apply the fit 185 stats->clipIter = 3; 186 pmPSF_MaskApTrend (psf , PM_PSF_SKYBIAS);187 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 188 189 // apply the fit 190 stats->clipIter = 3; 191 pmPSF_MaskApTrend (psf , PM_PSF_SKYSAT_XY_LIN);181 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 182 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 183 184 // apply the fit 185 stats->clipIter = 3; 186 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKYBIAS); 187 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 188 189 // apply the fit 190 stats->clipIter = 3; 191 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKYSAT_XY_LIN); 192 192 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 193 193 } … … 196 196 // first clip out objects which are too far from the median 197 197 stats->clipIter = 2; 198 pmPSF_MaskApTrend (psf , PM_PSF_CONSTANT);198 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_CONSTANT); 199 199 psf->ApTrend = psVectorClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 200 200 201 201 // fit just SkyBias and clip out objects which are too far from the median 202 202 stats->clipIter = 2; 203 pmPSF_MaskApTrend (psf , PM_PSF_SKYBIAS);203 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_SKYBIAS); 204 204 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 205 205 206 206 // finally, fit x, y, SkyBias and clip out objects which are too far from the median 207 207 stats->clipIter = 3; 208 pmPSF_MaskApTrend (psf , PM_PSF_ALL);208 pmPSF_MaskApTrend (psf->ApTrend, PM_PSF_ALL); 209 209 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 210 210 }
Note:
See TracChangeset
for help on using the changeset viewer.
