this is really a WTF bug

This commit is contained in:
tth 2020-07-24 16:21:24 +02:00
parent 24976855d7
commit d8bb6f4ec6
1 changed files with 8 additions and 4 deletions

View File

@ -36,15 +36,19 @@ if (fits_create_file(&fptr, outname, &status)) {
} }
naxes[0] = src->width; naxes[1] = src->height; naxes[0] = src->width; naxes[1] = src->height;
array = calloc(src->height, sizeof(float *)); array = calloc(src->height, sizeof(float *));
if (verbosity) fimg_describe(src, "to be saved as FITS");
if (verbosity) fimg_describe(src, "to be saved as fits"); #define REVERSE 1
/* initialize pointers to the start of each row of the image */ /* initialize pointers to the start of each row of the image */
array[0] = src->R;
for( idx=0; idx<naxes[1]; idx++ ) { for( idx=0; idx<naxes[1]; idx++ ) {
/**** MAGIC CODE MUST COME HERE ****/ /**** MAGIC CODE MUST COME HERE ****/
#if REVERSE
k = naxes[1] - idx - 1; k = naxes[1] - idx - 1;
#else
k = idx;
#endif
array[idx] = src->R + (k*naxes[0]); array[idx] = src->R + (k*naxes[0]);
fprintf(stderr, " %6d %6d %p\n", idx, k, array[idx]); fprintf(stderr, " %6d %6d %p\n", idx, k, array[idx]);
} }