Changeset 8246 for trunk/psModules/src/imcombine/pmImageSubtract.c
- Timestamp:
- Aug 8, 2006, 4:37:08 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psModules/src/imcombine/pmImageSubtract.c (modified) (32 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/imcombine/pmImageSubtract.c
r7769 r8246 7 7 * @author GLG, MHPCC 8 8 * 9 * @version $Revision: 1. 3$ $Name: not supported by cvs2svn $10 * @date $Date: 2006-0 6-30 23:59:49$9 * @version $Revision: 1.4 $ $Name: not supported by cvs2svn $ 10 * @date $Date: 2006-08-09 02:37:07 $ 11 11 * 12 12 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 97 97 int spatialOrder) 98 98 { 99 psTrace(" ImageSubtract.pmSubtractionKernelsAllocPOIS", 3,99 psTrace("psModules.imcombine", 3, 100 100 "Calling pmSubtractionKernelsAllocPOIS(%d, %d)\n", size, spatialOrder); 101 101 PS_ASSERT_INT_POSITIVE(size, NULL); … … 173 173 } 174 174 175 psTrace(" ImageSubtract.pmSubtractionKernelsAllocPOIS", 3,175 psTrace("psModules.imcombine", 3, 176 176 "Exiting pmSubtractionKernelsAllocPOIS(%d, %d)\n", size, spatialOrder); 177 177 return(tmp); … … 190 190 PS_ASSERT_VECTOR_NON_NULL(sigmas, NULL); 191 191 PS_ASSERT_VECTOR_NON_NULL(orders, NULL); 192 psTrace(" ImageSubtract.pmSubtractionKernelsAllocISIS", 3,192 psTrace("psModules.imcombine", 3, 193 193 "Calling pmSubtractionKernelsAllocISIS(%d, %d, %d, %d)\n", 194 194 sigmas->n, orders->n, size, spatialOrder); … … 306 306 } 307 307 308 psTrace(" ImageSubtract.pmSubtractionKernelsAllocISIS", 3,308 psTrace("psModules.imcombine", 3, 309 309 "Exiting pmSubtractionKernelsAllocISIS(%d, %d, %d, %d)\n", 310 310 sigmas->n, orders->n, size, spatialOrder); … … 333 333 ) 334 334 { 335 psTrace(" ImageSubtract.pmSubtractionFindStamps", 3,335 psTrace("psModules.imcombine", 3, 336 336 "Calling pmSubtractionFindStamps(%d, %f, %d, %d, %d)\n", 337 337 maskVal, threshold, xNum, yNum, border); … … 456 456 } 457 457 } 458 psTrace(" ImageSubtract.pmSubtractionFindStamps", 3,458 psTrace("psModules.imcombine", 3, 459 459 "Exiting pmSubtractionFindStamps(%d, %f, %d, %d, %d)\n", 460 460 maskVal, threshold, xNum, yNum, border); … … 474 474 psF32 y) 475 475 { 476 psTrace(" ImageSubtract.GenSpatialOrder", 4,476 psTrace("psModules.imcombine", 4, 477 477 "Calling GenSpatialOrder(%d, %f, %f)\n", spatialOrder, x, y); 478 478 … … 490 490 } 491 491 492 psTrace(" ImageSubtract.GenSpatialOrder", 4,492 psTrace("psModules.imcombine", 4, 493 493 "Exiting GenSpatialOrder(%d, %f, %f)\n", spatialOrder, x, y); 494 494 … … 509 509 510 510 511 psTrace(" ImageSubtract.IsisKernelConvolve", 4,511 psTrace("psModules.imcombine", 4, 512 512 "Calling IsisKernelConvolve(%d, %d, %d)\n", kernelID, col, row); 513 513 psS32 spatialOrder = kernels->spatialOrder; … … 541 541 psFree(polyValues); 542 542 543 psTrace(" ImageSubtract.IsisKernelConvolve", 4,543 psTrace("psModules.imcombine", 4, 544 544 "Exiting IsisKernelConvolve(%d, %d, %d)\n", kernelID, col, row); 545 545 return(conv); … … 563 563 psS32 row) 564 564 { 565 psTrace(" ImageSubtract.ConvolvePixelPois", 4,565 psTrace("psModules.imcombine", 4, 566 566 "Calling ConvolvePixelPois(%d, %d)\n", col, row); 567 567 psS32 nBF = kernels->u->n; … … 604 604 } 605 605 606 psTrace(" ImageSubtract.ConvolvePixelPois", 4,606 psTrace("psModules.imcombine", 4, 607 607 "Exiting ConvolvePixelPois(%d, %d)\n", col, row); 608 608 return(conv); … … 628 628 psS32 row) 629 629 { 630 psTrace(" ImageSubtract.ConvolvePixelIsis", 4,630 psTrace("psModules.imcombine", 4, 631 631 "Calling ConvolvePixelIsis(%d, %d)\n", col, row); 632 632 psF32 background = solution->data.F64[solution->n-1]; … … 640 640 } 641 641 642 psTrace(" ImageSubtract.ConvolvePixelIsis", 4,642 psTrace("psModules.imcombine", 4, 643 643 "Exiting ConvolvePixelIsis(%d, %d)\n", col, row); 644 644 return(conv); … … 655 655 const psSubtractionKernels *kernels) 656 656 { 657 psTrace(" ImageSubtract.ConvolveImage", 4, "Calling ConvolveImage()\n");657 psTrace("psModules.imcombine", 4, "Calling ConvolveImage()\n"); 658 658 PS_ASSERT_IMAGE_NON_NULL(input, NULL); 659 659 PS_ASSERT_IMAGE_NON_EMPTY(input, NULL); … … 724 724 } 725 725 726 psTrace(" ImageSubtract.ConvolveImage", 4, "Exiting ConvolveImage()\n");726 psTrace("psModules.imcombine", 4, "Exiting ConvolveImage()\n"); 727 727 return convolved; 728 728 } … … 743 743 ) 744 744 { 745 psTrace(" ImageSubtract.pmSubtractionCalculateEquation", 3,745 psTrace("psModules.imcombine", 3, 746 746 "Calling pmSubtractionCalculateEquation()\n"); 747 747 PS_ASSERT_PTR_NON_NULL(stamps, false); … … 805 805 for (psS32 s = 0; s < stamps->n; s++) { 806 806 pmStamp *stamp = (pmStamp *) stamps->data[s]; 807 psTrace("p mSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d\n", s);807 psTrace("psModules.imcombine", 6, "subCalcEqn(): stamp %d\n", s); 808 808 if (stamp->status == PM_STAMP_RECALC) { 809 psTrace("p mSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: status is PM_STAMP_RECALC.\n", s);809 psTrace("psModules.imcombine", 6, "subCalcEqn(): stamp %d: status is PM_STAMP_RECALC.\n", s); 810 810 psImage *stampMatrix = stamp->matrix; 811 811 psVector *stampVector = stamp->vector; … … 829 829 } 830 830 psVectorInit(stampVector, 0.0); 831 psTrace("p mSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: allocate matrix and vector.\n", s);831 psTrace("psModules.imcombine", 6, "subCalcEqn(): stamp %d: allocate matrix and vector.\n", s); 832 832 833 833 // … … 840 840 ((psF64) (stamp->y - numHalfRows)) / ((psF64) numHalfRows)); 841 841 842 psTrace("p mSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: generated spatial order terms.\n", s);842 psTrace("psModules.imcombine", 6, "subCalcEqn(): stamp %d: generated spatial order terms.\n", s); 843 843 844 844 if (kernels->type == PM_SUBTRACTION_KERNEL_POIS) { … … 848 848 for (psS32 y = stamp->y - footprint; y < stamp->y + footprint; y++) { 849 849 for (psS32 x = stamp->x - footprint; x < stamp->x + footprint; x++) { 850 psTrace("p mSubtractionCalculateEquation", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x);850 psTrace("psModules.imcombine", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x); 851 851 852 852 // The inverse of the noise, squared. … … 932 932 for (psS32 y = stamp->y - footprint; y < stamp->y + footprint; y++) { 933 933 for (psS32 x = stamp->x - footprint; x < stamp->x + footprint; x++) { 934 psTrace("p mSubtractionCalculateEquation", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x);934 psTrace("psModules.imcombine", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x); 935 935 psF32 invNoise2 = 1.0/reference->data.F32[y][x]; // The inverse of the noise, squared. 936 936 … … 998 998 } 999 999 } 1000 psTrace(" ImageSubtract.pmSubtractionCalculateEquation", 3,1000 psTrace("psModules.imcombine", 3, 1001 1001 "Exiting pmSubtractionCalculateEquation()\n"); 1002 1002 return(true); … … 1012 1012 ) 1013 1013 { 1014 psTrace(" ImageSubtract.pmSubtractionSolveEquation", 3,1014 psTrace("psModules.imcombine", 3, 1015 1015 "Calling pmSubtractionSolveEquation()\n"); 1016 1016 PS_ASSERT_PTR_NON_NULL(stamps, NULL); … … 1120 1120 psFree(permutation); 1121 1121 1122 psTrace(" ImageSubtract.pmSubtractionSolveEquation", 3,1122 psTrace("psModules.imcombine", 3, 1123 1123 "Exiting pmSubtractionSolveEquation()\n"); 1124 1124 return(solution); … … 1138 1138 const psVector *solution) 1139 1139 { 1140 psTrace(" ImageSubtract.CalculateDeviations", 4,1140 psTrace("psModules.imcombine", 4, 1141 1141 "Calling CalculateDeviations()\n"); 1142 1142 PS_ASSERT_PTR_NON_NULL(stamps, NULL); … … 1220 1220 psFree(subStamp); 1221 1221 1222 psTrace(" ImageSubtract.CalculateDeviations", 4,1222 psTrace("psModules.imcombine", 4, 1223 1223 "Exiting CalculateDeviations()\n"); 1224 1224 return deviations; … … 1238 1238 ) 1239 1239 { 1240 psTrace(" ImageSubtract.pmSubtractionRejectStamps", 3,1240 psTrace("psModules.imcombine", 3, 1241 1241 "Calling pmSubtractionRejectStamps()\n"); 1242 1242 PS_ASSERT_PTR_NON_NULL(stamps, false); … … 1313 1313 1314 1314 psFree(deviations); 1315 psTrace(" ImageSubtract.pmSubtractionRejectStamps", 3,1315 psTrace("psModules.imcombine", 3, 1316 1316 "Exiting pmSubtractionRejectStamps()\n"); 1317 1317 return(true); … … 1327 1327 ) 1328 1328 { 1329 psTrace(" ImageSubtract.pmSubtractionKernelImage", 3,1329 psTrace("psModules.imcombine", 3, 1330 1330 "Calling pmSubtractionKernelImage()\n"); 1331 1331 PS_ASSERT_VECTOR_NON_NULL(solution, NULL); … … 1401 1401 psFree(polyValues); 1402 1402 1403 psTrace(" ImageSubtract.pmSubtractionKernelImage", 3,1403 psTrace("psModules.imcombine", 3, 1404 1404 "Exiting pmSubtractionKernelImage()\n"); 1405 1405 return(out);
Note:
See TracChangeset
for help on using the changeset viewer.
