Changeset 4223
- Timestamp:
- Jun 13, 2005, 10:09:53 AM (21 years ago)
- File:
-
- 1 edited
-
trunk/psModules/src/pmImageSubtract.c (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/pmImageSubtract.c
r4185 r4223 3 3 * This file will ... 4 4 * 5 * @author Paul Price, IfA (original prototype) 5 6 * @author GLG, MHPCC 6 * @author Paul Price, IfA7 7 * 8 * @version $Revision: 1. 8$ $Name: not supported by cvs2svn $9 * @date $Date: 2005-06- 09 06:18:19$8 * @version $Revision: 1.9 $ $Name: not supported by cvs2svn $ 9 * @date $Date: 2005-06-13 20:09:53 $ 10 10 * 11 11 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 37 37 #include "pmImageSubtract.h" 38 38 39 // XXX: Put this is psConstants.h 39 40 #define PS_ASSERT_INT_LESS_THAN(VAR1, VAR2, RVAL) \ 40 41 if (!(VAR1 < VAR2)) { \ … … 171 172 172 173 /******************************************************************************* 173 XXX: Get the types correct .174 XXX: Get the types correct (u, v, xOrder, yOrder). 174 175 175 176 XXX: For the various spatial orders, verify that the preCalc data does not … … 194 195 // 195 196 psS32 numSigmas = sigmas->n; 196 // psS32 nBF = 0;197 // for (psS32 s = 0 ; s < numSigmas ; s++) {198 // nBF+= ((orders->n + 1) * (orders->n + 2) / 2);199 // }200 // nBF*= ((spatialOrder + 1) * (spatialOrder + 2) / 2);201 202 // psS32 nBF = 1;203 // for (psS32 s = 0 ; s < numSigmas ; s++) {204 // nBF*= ((orders->n + 1) * (orders->n + 2) / 2);205 // }206 // nBF*= ((spatialOrder + 1) * (spatialOrder + 2) / 2);207 197 208 198 // XXX: Get rid of the sigma loop? We merely multiple nBF by numSigmas? … … 461 451 to ensure no seg faults. 462 452 ******************************************************************************/ 463 psF32 GeneralKernelConvolve(const psImage *input,464 const psSubtractionKernels *kernels,465 psS32 kernelID,466 psS32 col,467 psS32 row)453 static psF32 GeneralKernelConvolve(const psImage *input, 454 const psSubtractionKernels *kernels, 455 psS32 kernelID, 456 psS32 col, 457 psS32 row) 468 458 { 469 459 psS32 spatialOrder = kernels->p_spatialOrder; … … 500 490 501 491 /******************************************************************************* 502 XXX: How is the spatial order factor calculated? Is it simply a x^iy^j power503 evaluated at the (x, y) center for the stamp? Why bother with the 2-D504 polynomial in evaluating it?505 506 492 XXX: Should we assert that the footprint is equal to the kernel size, or does 507 493 they have nothing in common. … … 686 672 } 687 673 } else { 688 p rintf("XXX: put some warning message here (bad kernel->type).\n");674 psLogMsg(__func__, PS_LOG_WARN, "WARNING: unknown kernel->type.\n"); 689 675 return(false); 690 676 } // if-else on kernel type. … … 708 694 709 695 // 710 // XXX: understand this.711 696 // XXX: Why aren't they using numSolveParams instead of numKernels? 712 697 // XXX: is this POIS specific? 713 // XXX: What is config->penalty?714 698 // 715 699 #define XXX_CONFIG_PENALTY 1.0 … … 900 884 arbitrary image with either an ISIS or POIS set of kernel basis functions. 901 885 ******************************************************************************/ 902 psImage *ConvolveImage(const psImage *input,903 const psImage *mask,904 psU32 badStampMaskVal,905 const psVector *solution,906 const psSubtractionKernels *kernels)886 static psImage *ConvolveImage(const psImage *input, 887 const psImage *mask, 888 psU32 badStampMaskVal, 889 const psVector *solution, 890 const psSubtractionKernels *kernels) 907 891 { 908 892 PS_ASSERT_IMAGE_NON_NULL(input, NULL); … … 980 964 /******************************************************************************* 981 965 ******************************************************************************/ 982 psVector *CalculateDeviations(psVector *deviations,983 psArray *stamps,984 psS32 footprint,985 const psImage *refImage,986 const psImage *inImage,987 const psImage *mask,988 psU32 badStampMaskVal,989 const psSubtractionKernels *kernels,990 const psVector *solution)966 static psVector *CalculateDeviations(psVector *deviations, 967 psArray *stamps, 968 psS32 footprint, 969 const psImage *refImage, 970 const psImage *inImage, 971 const psImage *mask, 972 psU32 badStampMaskVal, 973 const psSubtractionKernels *kernels, 974 const psVector *solution) 991 975 { 992 976 PS_ASSERT_PTR_NON_NULL(stamps, NULL);
Note:
See TracChangeset
for help on using the changeset viewer.
