IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Mar 26, 2010, 1:53:22 PM (16 years ago)
Author:
eugene
Message:

minor changes (comments, format)

Location:
trunk/Ohana/src/relastro/src
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana/src/relastro/src/ImageOps.c

    r27435 r27478  
    11# include "relastro.h"
    2 # define TESTING
     2
    33# define USE_IMAGE_ID 1
    44
     
    77static off_t       **mlist;   // measure reference for measurement on image
    88static off_t        *Nlist;   // number of measurements on image
    9 static off_t        *NLIST;   // allocated number of measurements on image
     9static off_t        *NLIST;   // allocated number of measurements on image   
    1010
    1111static Image        *image;   // list of available images
    1212static off_t        Nimage;   // number of available images
    1313
     14// if we search by image ID, we sort (imageIDs, imageIdx) by imageIDs to get a sorted
     15// index
     16
    1417# if USE_IMAGE_ID
    15 // we sort (imageIDs, imageIdx) by imageIDs to get a sorted index
    1618static off_t        *imageIDs; // list of all image IDs
    1719static off_t        *imageIdx; // list of index for image IDs
     
    5557  }
    5658# endif
    57 
    5859}
    5960
     
    132133  for (i = 0; i < Ncatalog; i++) {
    133134    for (j = 0; j < catalog[i].Nmeasure; j++) {
     135      // XXX : update to optionally restrict by photcode equiv a la relphot
     136      // ecode = GetPhotcodeEquivCodebyCode (catalog[i].measure[j].photcode);
     137      // if (photcode[0].code != ecode) continue;
    134138      matchImage (catalog, j, i);
    135139    }
     
    144148
    145149# if USE_IMAGE_ID
    146 /* modify this function to use the measure->imageID field */
     150// this is the imageID-based match
    147151void matchImage (Catalog *catalog, off_t meas, int cat) {
    148152
     
    178182
    179183# else
    180 /* modify this function to use the measure->imageID field */
     184// this is the time-based match
    181185void matchImage (Catalog *catalog, off_t meas, int cat) {
    182186
  • trunk/Ohana/src/relastro/src/MosaicOps.c

    r27435 r27478  
    22
    33// array of mosaic definition structures
    4 static off_t    Nmosaic;
     4static off_t   Nmosaic;
    55static Mosaic *mosaic;
    66
     
    161161  }
    162162
     163  free (startMos);
     164  free (stopMos);
     165  free (indexMos);
    163166  return;
    164167}
  • trunk/Ohana/src/relastro/src/relastro.c

    r27435 r27478  
    1313  FITS_DB db;
    1414  struct timeval start, stop;
    15 
    1615  SkyList *skylist = NULL;
    1716
  • trunk/Ohana/src/relastro/src/select_images.c

    r27435 r27478  
    2020  Image *image;
    2121  off_t i, j, k, m, nStart, iSky, nimage, NIMAGE;
     22  off_t *line_number;
    2223  int InRange, found;
    2324  double Ri[5], Di[5], Xi[5], Yi[5], dx, dy;
    24   off_t *line_number;
    2525  Coords tcoords;
    2626  SkyRegionCoords *skycoords;
     
    127127    found = FALSE;
    128128
    129     /* transform to ra,dec */
     129    /* transform corners to ra,dec */
    130130    double RminImage = 360.0;
    131     // double RmaxImage =   0.0;
    132131    for (j = 0; j < 5; j++) {
    133132      XY_to_RD (&Ri[j], &Di[j], Xi[j], Yi[j], &timage[i].coords);
    134133      RminImage = MIN(RminImage, Ri[j]);
    135       // RmaxImage = MIN(RmaxImage, Ri[j]);
    136134    }
    137135
     
    143141
    144142      m = index[iSky];
    145       // if (RmaxImage < skylist[0].regions[m][0].Rmin) {
    146       // break;
    147       // }
    148143
    149144      /* we make positional comparisons in the projection of catalog */
     
    180175      }
    181176      if (RESET) {
    182         // XXX do we need this in relastro?
     177        // XXX do we need / want to do this in relastro?
    183178        assignMcal (&image[nimage], (double *) NULL, -1);
    184179        image[nimage].dMcal = NAN;
Note: See TracChangeset for help on using the changeset viewer.