Changeset 27478 for trunk/Ohana/src/relastro
- Timestamp:
- Mar 26, 2010, 1:53:22 PM (16 years ago)
- Location:
- trunk/Ohana/src/relastro/src
- Files:
-
- 4 edited
-
ImageOps.c (modified) (6 diffs)
-
MosaicOps.c (modified) (2 diffs)
-
relastro.c (modified) (1 diff)
-
select_images.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/relastro/src/ImageOps.c
r27435 r27478 1 1 # include "relastro.h" 2 # define TESTING 2 3 3 # define USE_IMAGE_ID 1 4 4 … … 7 7 static off_t **mlist; // measure reference for measurement on image 8 8 static off_t *Nlist; // number of measurements on image 9 static off_t *NLIST; // allocated number of measurements on image 9 static off_t *NLIST; // allocated number of measurements on image 10 10 11 11 static Image *image; // list of available images 12 12 static off_t Nimage; // number of available images 13 13 14 // if we search by image ID, we sort (imageIDs, imageIdx) by imageIDs to get a sorted 15 // index 16 14 17 # if USE_IMAGE_ID 15 // we sort (imageIDs, imageIdx) by imageIDs to get a sorted index16 18 static off_t *imageIDs; // list of all image IDs 17 19 static off_t *imageIdx; // list of index for image IDs … … 55 57 } 56 58 # endif 57 58 59 } 59 60 … … 132 133 for (i = 0; i < Ncatalog; i++) { 133 134 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; 134 138 matchImage (catalog, j, i); 135 139 } … … 144 148 145 149 # if USE_IMAGE_ID 146 / * modify this function to use the measure->imageID field */150 // this is the imageID-based match 147 151 void matchImage (Catalog *catalog, off_t meas, int cat) { 148 152 … … 178 182 179 183 # else 180 / * modify this function to use the measure->imageID field */184 // this is the time-based match 181 185 void matchImage (Catalog *catalog, off_t meas, int cat) { 182 186 -
trunk/Ohana/src/relastro/src/MosaicOps.c
r27435 r27478 2 2 3 3 // array of mosaic definition structures 4 static off_t Nmosaic;4 static off_t Nmosaic; 5 5 static Mosaic *mosaic; 6 6 … … 161 161 } 162 162 163 free (startMos); 164 free (stopMos); 165 free (indexMos); 163 166 return; 164 167 } -
trunk/Ohana/src/relastro/src/relastro.c
r27435 r27478 13 13 FITS_DB db; 14 14 struct timeval start, stop; 15 16 15 SkyList *skylist = NULL; 17 16 -
trunk/Ohana/src/relastro/src/select_images.c
r27435 r27478 20 20 Image *image; 21 21 off_t i, j, k, m, nStart, iSky, nimage, NIMAGE; 22 off_t *line_number; 22 23 int InRange, found; 23 24 double Ri[5], Di[5], Xi[5], Yi[5], dx, dy; 24 off_t *line_number;25 25 Coords tcoords; 26 26 SkyRegionCoords *skycoords; … … 127 127 found = FALSE; 128 128 129 /* transform to ra,dec */129 /* transform corners to ra,dec */ 130 130 double RminImage = 360.0; 131 // double RmaxImage = 0.0;132 131 for (j = 0; j < 5; j++) { 133 132 XY_to_RD (&Ri[j], &Di[j], Xi[j], Yi[j], &timage[i].coords); 134 133 RminImage = MIN(RminImage, Ri[j]); 135 // RmaxImage = MIN(RmaxImage, Ri[j]);136 134 } 137 135 … … 143 141 144 142 m = index[iSky]; 145 // if (RmaxImage < skylist[0].regions[m][0].Rmin) {146 // break;147 // }148 143 149 144 /* we make positional comparisons in the projection of catalog */ … … 180 175 } 181 176 if (RESET) { 182 // XXX do we need this in relastro?177 // XXX do we need / want to do this in relastro? 183 178 assignMcal (&image[nimage], (double *) NULL, -1); 184 179 image[nimage].dMcal = NAN;
Note:
See TracChangeset
for help on using the changeset viewer.
