Changeset 4424
- Timestamp:
- Jun 28, 2005, 3:17:23 PM (21 years ago)
- File:
-
- 1 edited
-
trunk/psModules/src/pmImageSubtract.c (modified) (35 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/pmImageSubtract.c
r4417 r4424 7 7 * @author GLG, MHPCC 8 8 * 9 * @version $Revision: 1.1 1$ $Name: not supported by cvs2svn $10 * @date $Date: 2005-06-2 8 23:23:55$9 * @version $Revision: 1.12 $ $Name: not supported by cvs2svn $ 10 * @date $Date: 2005-06-29 01:17:23 $ 11 11 * 12 12 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 98 98 int spatialOrder) 99 99 { 100 psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 1,100 psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 3, 101 101 "Calling pmSubtractionKernelsAllocPOIS(%d, %d)\n", size, spatialOrder); 102 102 PS_ASSERT_INT_POSITIVE(size, NULL); … … 170 170 } 171 171 172 psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 1,172 psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 3, 173 173 "Exiting pmSubtractionKernelsAllocPOIS(%d, %d)\n", size, spatialOrder); 174 174 return(tmp); … … 187 187 PS_ASSERT_VECTOR_NON_NULL(sigmas, NULL); 188 188 PS_ASSERT_VECTOR_NON_NULL(orders, NULL); 189 psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 1,189 psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 3, 190 190 "Calling pmSubtractionKernelsAllocISIS(%d, %d, %d, %d)\n", 191 191 sigmas->n, orders->n, size, spatialOrder); … … 295 295 } 296 296 297 psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 1,297 psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 3, 298 298 "Exiting pmSubtractionKernelsAllocISIS(%d, %d, %d, %d)\n", 299 299 sigmas->n, orders->n, size, spatialOrder); … … 322 322 ) 323 323 { 324 psTrace("ImageSubtract.pmSubtractionFindStamps", 1,324 psTrace("ImageSubtract.pmSubtractionFindStamps", 3, 325 325 "Calling pmSubtractionFindStamps(%d, %f, %d, %d, %d)\n", 326 326 maskVal, threshold, xNum, yNum, border); … … 444 444 } 445 445 } 446 psTrace("ImageSubtract.pmSubtractionFindStamps", 1,446 psTrace("ImageSubtract.pmSubtractionFindStamps", 3, 447 447 "Exiting pmSubtractionFindStamps(%d, %f, %d, %d, %d)\n", 448 448 maskVal, threshold, xNum, yNum, border); … … 462 462 psF32 y) 463 463 { 464 psTrace("ImageSubtract.GenSpatialOrder", 5,464 psTrace("ImageSubtract.GenSpatialOrder", 4, 465 465 "Calling GenSpatialOrder(%d, %f, %f)\n", spatialOrder, x, y); 466 466 … … 478 478 } 479 479 480 psTrace("ImageSubtract.GenSpatialOrder", 5,480 psTrace("ImageSubtract.GenSpatialOrder", 4, 481 481 "Exiting GenSpatialOrder(%d, %f, %f)\n", spatialOrder, x, y); 482 482 … … 497 497 498 498 499 psTrace("ImageSubtract.IsisKernelConvolve", 5,499 psTrace("ImageSubtract.IsisKernelConvolve", 4, 500 500 "Calling IsisKernelConvolve(%d, %d, %d)\n", kernelID, col, row); 501 501 psS32 spatialOrder = kernels->p_spatialOrder; … … 529 529 psFree(polyValues); 530 530 531 psTrace("ImageSubtract.IsisKernelConvolve", 5,531 psTrace("ImageSubtract.IsisKernelConvolve", 4, 532 532 "Exiting IsisKernelConvolve(%d, %d, %d)\n", kernelID, col, row); 533 533 return(conv); … … 551 551 psS32 row) 552 552 { 553 psTrace("ImageSubtract.ConvolvePixelPois", 5,553 psTrace("ImageSubtract.ConvolvePixelPois", 4, 554 554 "Calling ConvolvePixelPois(%d, %d)\n", col, row); 555 555 psS32 nBF = kernels->u->n; … … 592 592 } 593 593 594 psTrace("ImageSubtract.ConvolvePixelPois", 5,594 psTrace("ImageSubtract.ConvolvePixelPois", 4, 595 595 "Exiting ConvolvePixelPois(%d, %d)\n", col, row); 596 596 return(conv); … … 616 616 psS32 row) 617 617 { 618 psTrace("ImageSubtract.ConvolvePixelIsis", 5,618 psTrace("ImageSubtract.ConvolvePixelIsis", 4, 619 619 "Calling ConvolvePixelIsis(%d, %d)\n", col, row); 620 620 psF32 background = solution->data.F64[solution->n-1]; … … 628 628 } 629 629 630 psTrace("ImageSubtract.ConvolvePixelIsis", 5,630 psTrace("ImageSubtract.ConvolvePixelIsis", 4, 631 631 "Exiting ConvolvePixelIsis(%d, %d)\n", col, row); 632 632 return(conv); … … 643 643 const psSubtractionKernels *kernels) 644 644 { 645 psTrace("ImageSubtract.ConvolveImage", 5, "Calling ConvolveImage()\n");645 psTrace("ImageSubtract.ConvolveImage", 4, "Calling ConvolveImage()\n"); 646 646 PS_ASSERT_IMAGE_NON_NULL(input, NULL); 647 647 PS_ASSERT_IMAGE_NON_EMPTY(input, NULL); … … 712 712 } 713 713 714 psTrace("ImageSubtract.ConvolveImage", 5, "Exiting ConvolveImage()\n");714 psTrace("ImageSubtract.ConvolveImage", 4, "Exiting ConvolveImage()\n"); 715 715 return convolved; 716 716 } … … 731 731 ) 732 732 { 733 psTrace("ImageSubtract.pmSubtractionCalculateEquation", 1,733 psTrace("ImageSubtract.pmSubtractionCalculateEquation", 3, 734 734 "Calling pmSubtractionCalculateEquation()\n"); 735 735 PS_ASSERT_PTR_NON_NULL(stamps, false); … … 793 793 for (psS32 s = 0; s < stamps->n; s++) { 794 794 pmStamp *stamp = (pmStamp *) stamps->data[s]; 795 psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d\n", s);795 psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d\n", s); 796 796 if (stamp->status == PM_STAMP_RECALC) { 797 psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d: status is PM_STAMP_RECALC.\n", s);797 psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: status is PM_STAMP_RECALC.\n", s); 798 798 psImage *stampMatrix = stamp->matrix; 799 799 psVector *stampVector = stamp->vector; … … 816 816 } 817 817 PS_VECTOR_SET_F64(stampVector, 0.0); 818 psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d: allocate matrix and vector.\n", s);818 psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: allocate matrix and vector.\n", s); 819 819 820 820 // … … 827 827 ((psF64) (stamp->y - numHalfRows)) / ((psF64) numHalfRows)); 828 828 829 psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d: generated spatial order terms.\n", s);829 psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: generated spatial order terms.\n", s); 830 830 831 831 if (kernels->type == PM_SUBTRACTION_KERNEL_POIS) { … … 835 835 for (psS32 y = stamp->y - footprint; y < stamp->y + footprint; y++) { 836 836 for (psS32 x = stamp->x - footprint; x < stamp->x + footprint; x++) { 837 psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): pixel (%d, %d).\n", y, x);837 psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x); 838 838 839 839 // The inverse of the noise, squared. … … 919 919 for (psS32 y = stamp->y - footprint; y < stamp->y + footprint; y++) { 920 920 for (psS32 x = stamp->x - footprint; x < stamp->x + footprint; x++) { 921 psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): pixel (%d, %d).\n", y, x);921 psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x); 922 922 psF32 invNoise2 = 1.0/reference->data.F32[y][x]; // The inverse of the noise, squared. 923 923 … … 985 985 } 986 986 } 987 psTrace("ImageSubtract.pmSubtractionCalculateEquation", 1,987 psTrace("ImageSubtract.pmSubtractionCalculateEquation", 3, 988 988 "Exiting pmSubtractionCalculateEquation()\n"); 989 989 return(true); … … 999 999 ) 1000 1000 { 1001 psTrace("ImageSubtract.pmSubtractionSolveEquation", 1,1001 psTrace("ImageSubtract.pmSubtractionSolveEquation", 3, 1002 1002 "Calling pmSubtractionSolveEquation()\n"); 1003 1003 PS_ASSERT_PTR_NON_NULL(stamps, NULL); … … 1067 1067 1068 1068 // XXX: psTrace() 1069 if ( 1) {1069 if (0) { 1070 1070 PS_IMAGE_PRINT_F64(sumMatrix); 1071 1071 } … … 1081 1081 } 1082 1082 // XXX: psTrace() 1083 if ( 1) {1083 if (0) { 1084 1084 PS_IMAGE_PRINT_F64(luMatrix); 1085 1085 } … … 1087 1087 solution = psMatrixLUSolve(solution, luMatrix, sumVector, permutation); 1088 1088 // XXX: psTrace() 1089 if (1) { 1089 // XXX: should we be checking for NAN's in the solution vector? 1090 if (0) { 1090 1091 PS_VECTOR_PRINT_F64(solution); 1091 1092 } … … 1104 1105 psFree(permutation); 1105 1106 1106 psTrace("ImageSubtract.pmSubtractionSolveEquation", 1,1107 psTrace("ImageSubtract.pmSubtractionSolveEquation", 3, 1107 1108 "Exiting pmSubtractionSolveEquation()\n"); 1108 1109 return(solution); … … 1122 1123 const psVector *solution) 1123 1124 { 1124 psTrace("ImageSubtract.CalculateDeviations", 5,1125 psTrace("ImageSubtract.CalculateDeviations", 4, 1125 1126 "Calling CalculateDeviations()\n"); 1126 1127 PS_ASSERT_PTR_NON_NULL(stamps, NULL); … … 1203 1204 psFree(subStamp); 1204 1205 1205 psTrace("ImageSubtract.CalculateDeviations", 5,1206 psTrace("ImageSubtract.CalculateDeviations", 4, 1206 1207 "Exiting CalculateDeviations()\n"); 1207 1208 return deviations; … … 1221 1222 ) 1222 1223 { 1223 psTrace("ImageSubtract.pmSubtractionRejectStamps", 5,1224 psTrace("ImageSubtract.pmSubtractionRejectStamps", 3, 1224 1225 "Calling pmSubtractionRejectStamps()\n"); 1225 1226 PS_ASSERT_PTR_NON_NULL(stamps, false); … … 1296 1297 1297 1298 psFree(deviations); 1298 psTrace("ImageSubtract.pmSubtractionRejectStamps", 5,1299 psTrace("ImageSubtract.pmSubtractionRejectStamps", 3, 1299 1300 "Exiting pmSubtractionRejectStamps()\n"); 1300 1301 return(true); … … 1310 1311 ) 1311 1312 { 1312 psTrace("ImageSubtract.pmSubtractionKernelImage", 5,1313 psTrace("ImageSubtract.pmSubtractionKernelImage", 3, 1313 1314 "Calling pmSubtractionKernelImage()\n"); 1314 1315 PS_ASSERT_VECTOR_NON_NULL(solution, NULL); … … 1384 1385 psFree(polyValues); 1385 1386 1386 psTrace("ImageSubtract.pmSubtractionKernelImage", 5,1387 psTrace("ImageSubtract.pmSubtractionKernelImage", 3, 1387 1388 "Exiting pmSubtractionKernelImage()\n"); 1388 1389 return(out);
Note:
See TracChangeset
for help on using the changeset viewer.
