IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Dec 16, 2009, 12:22:44 PM (16 years ago)
Author:
bills
Message:

switch to using auto tools

Location:
trunk/magic/remove
Files:
2 edited
1 moved

Legend:

Unmodified
Added
Removed
  • trunk/magic/remove

    • Property svn:ignore set to
      configure
      Makefile.in
      Doxyfile
      config.log
      depcomp
      config.status
      config.guess
      ltmain.sh
      config.sub
      autom4te.cache
      libtool
      missing
      Makefile
      aclocal.m4
      install-sh

  • trunk/magic/remove/src

    • Property svn:ignore
      •  

        old new  
        44streakscompare
        55streaksrelease
        6 makefile
         6Makefile
         7Makefile.in
         8config.h
         9.deps
         10streaksVersionDefinitions.h
         11config.h.in
  • trunk/magic/remove/src/diffedpixels.c

    • Property svn:mergeinfo set to (toggle deleted branches)
      /branches/eam_branches/20091113/magic/remove/src/warpedpixels.c26119-26255
    r26406 r26438  
    99static int xRight(psPlane *pt, int y);
    1010
     11// Examine the set of diff skycells and compute the pixels that were
     12// included in difference processing
    1113bool
    12 computeWarpedPixels(streakFiles *sf)
     14computeDiffedPixels(streakFiles *sf)
    1315{
    1416    bool status;
     
    2931    psFree(bounds);
    3032
    31     sf->warpedPixels = psImageAlloc(width, height, PS_TYPE_U8);
    32     psImageInit(sf->warpedPixels, 0);
     33    sf->diffedPixels = psImageAlloc(width, height, PS_TYPE_U8);
     34    psImageInit(sf->diffedPixels, 0);
    3335
    3436    for (int i=0; i<n; i++) {
     
    4244    bool writeTouchedPixels = false; // XXX: make this an argument to the program
    4345    if (writeTouchedPixels) {
    44         // write out the warped pixels
     46        // write out the diffed pixels
    4547        psMetadata * header = psMetadataAlloc();
    4648        if (sf->inAstrom->fpa->toSky->type != PS_PROJ_DIS) {
    4749            pmAstromWriteWCS(header, sf->inAstrom->fpa, sf->chip, 0.001);
    4850        }
    49         psFits *fits = psFitsOpen("warpedpixels.fits", "w");
    50 
    51         psFitsWriteImage(fits, header, sf->warpedPixels, 0, NULL);
     51        psFits *fits = psFitsOpen("diffedpixels.fits", "w");
     52
     53        psFitsWriteImage(fits, header, sf->diffedPixels, 0, NULL);
    5254        psFitsClose(fits);
    5355        psFree(header);
     
    8486    /* now set up our wrapper to the chip astrometry to apply to the whole chip */
    8587    sf->astrom = streakSetAstrometry(sf->astrom, sf->stage, sf->inAstrom->fpa, sf->chip, false, NULL,
    86         sf->warpedPixels->numCols, sf->warpedPixels->numRows);
     88        sf->diffedPixels->numCols, sf->diffedPixels->numRows);
    8789
    8890    // convert corners of skycell to sky coordinates
     
    138140    // Now set the touched pixels
    139141    int ymin = fmax(0, pt[1].y );
    140     int ymax = fmin(pt[3].y + 0.5, sf->warpedPixels->numRows - 1);
     142    int ymax = fmin(pt[3].y + 0.5, sf->diffedPixels->numRows - 1);
    141143#if (DEBUG_PRINT > 1)
    142144    printf("\nymin: %d ymax: %d\n", ymin, ymax);
     
    149151            xleft = 0;
    150152        }
    151         if (xleft > sf->warpedPixels->numCols) {
     153        if (xleft > sf->diffedPixels->numCols) {
    152154            continue;
    153155        }
     
    155157            continue;
    156158        }
    157         if (xright >= sf->warpedPixels->numCols) {
    158             xright = sf->warpedPixels->numCols - 1;
     159        if (xright >= sf->diffedPixels->numCols) {
     160            xright = sf->diffedPixels->numCols - 1;
    159161        }
    160162#if (DEBUG_PRINT > 1)
     
    162164#endif
    163165
    164         psU8 *rowPixels = sf->warpedPixels->data.U8[y];
     166        psU8 *rowPixels = sf->diffedPixels->data.U8[y];
    165167
    166168        int n = xright - xleft + 1;
Note: See TracChangeset for help on using the changeset viewer.