Changeset 5607 for trunk/psphot/src/pmObjects_EAM.c
- Timestamp:
- Nov 25, 2005, 5:00:41 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/pmObjects_EAM.c (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/pmObjects_EAM.c
r5593 r5607 6 6 * @author EAM, IfA: significant modifications. 7 7 * 8 * @version $Revision: 1.1 1$ $Name: not supported by cvs2svn $9 * @date $Date: 2005-11-2 5 23:13:43$8 * @version $Revision: 1.12 $ $Name: not supported by cvs2svn $ 9 * @date $Date: 2005-11-26 03:00:41 $ 10 10 * 11 11 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 723 723 for (psS32 row = 0; row < source->pixels->numRows ; row++) { 724 724 for (psS32 col = 0; col < source->pixels->numCols ; col++) { 725 if ((source->mask != NULL) && (source->mask->data.U8[row][col])) continue; 725 if ((source->mask != NULL) && (source->mask->data.U8[row][col])) { 726 continue; 727 } 726 728 727 729 psF32 xDiff = col + source->pixels->col0 - xPeak; … … 730 732 // XXX EAM : calculate xDiff, yDiff up front; 731 733 // radius is just a function of (xDiff, yDiff) 732 if (!VALID_RADIUS(xDiff, yDiff, R2)) continue; 734 if (!VALID_RADIUS(xDiff, yDiff, R2)) { 735 continue; 736 } 733 737 734 738 psF32 pDiff = source->pixels->data.F32[row][col] - sky; … … 736 740 // XXX EAM : check for valid S/N in pixel 737 741 // XXX EAM : should this limit be user-defined? 738 if (pDiff / sqrt(source->weight->data.F32[row][col]) < 1) continue; 742 if (pDiff / sqrt(source->weight->data.F32[row][col]) < 1) { 743 continue; 744 } 739 745 740 746 Sum += pDiff; … … 1031 1037 1032 1038 // XXX EAM : can we use the value of SATURATE if mask is NULL? 1033 // allArray = psRegionSet (0, 0, 0, 0);1034 1039 int Nsatpix = psImageCountPixelMask (tmpSrc->mask, allArray, PSPHOT_MASK_SATURATED); 1035 1040 … … 1670 1675 } 1671 1676 } 1672 # if (0)1673 else {1674 // try keeping the moments sky -- unsuccessful1675 for (int i = 1; i < paramMask->n; i++) {1676 paramMask->data.U8[i] = 0;1677 }1678 paramMask->data.U8[0] = 1;1679 }1680 # endif1681 1677 1682 1678 … … 1814 1810 return(p_pmSourceAddOrSubModel(image, mask, model, center, 1)); 1815 1811 } 1816 1817 1818 // XXX: Put this is psConstants.h1819 #define PS_VECTOR_CHECK_SIZE(VEC1, N, RVAL) \1820 if (VEC1->n != N) { \1821 psError(PS_ERR_BAD_PARAMETER_SIZE, true, \1822 "psVector %s has size %d, should be %d.", \1823 #VEC1, VEC1->n, N); \1824 return(RVAL); \1825 }
Note:
See TracChangeset
for help on using the changeset viewer.
