Changeset 3671 in Sophya for trunk/AddOn/TAcq/minifits.cc
- Timestamp:
- Nov 13, 2009, 11:54:04 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/AddOn/TAcq/minifits.cc
r3658 r3671 97 97 if (fip != NULL) throw MiniFITSException("MiniFITSFile::Open() - fip != NULL"); 98 98 if (rwm == MF_Write) { 99 for(int i=0; i<MFITSHLEN; i++) header[i]=' '; 100 nkeya_ = 0; 99 101 FillHeader(); 100 102 fip = fopen(nom, "w"); … … 119 121 { 120 122 if (fip) { 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 123 if (rwmode == MF_Write) { 124 // on remplit avec des zeros pour avoir une longueur multiple de 2880 125 size_t padsz = MFITSHLEN-(totwsz%MFITSHLEN); 126 char zeros[160]; 127 for(size_t k=0; k<160; k++) zeros[k]=0; 128 while(padsz>160) { 129 fwrite(zeros, 1, 160, fip); 130 padsz-=160; 131 } 132 if (padsz>0) fwrite(zeros, 1, padsz, fip); 133 // On reecrit l'entete 134 FillHeader(); 135 fseek(fip, 0, SEEK_SET); 136 fwrite(header, 1, MFITSHLEN, fip); 137 } 138 fclose(fip); 137 139 } 138 140 fip = NULL; … … 187 189 else if (dtype == MF_Float32) bpix = -32; 188 190 char * buff = header+80; 189 sprintf(buff, "BITPIX = %20d / 191 sprintf(buff, "BITPIX = %20d / number of bits per data pixel", bpix); 190 192 buff[strlen(buff)] = ' '; 191 193 buff = header+160; … … 193 195 buff[strlen(buff)] = ' '; 194 196 buff = header+240; 195 sprintf(buff, "NAXIS1 = %20ld / number of bits per data pixel", (long)nax1);197 sprintf(buff, "NAXIS1 = %20ld / nb of pixels along X = PaquetSize", (long)nax1); 196 198 buff[strlen(buff)] = ' '; 197 199 buff = header+320; 198 sprintf(buff, "NAXIS2 = %20ld / number of bits per data pixel", (long)nax2);200 sprintf(buff, "NAXIS2 = %20ld / nb of rows = NumberOfPaquets", (long)nax2); 199 201 buff[strlen(buff)] = ' '; 200 202 buff = header+400+nkeya_*80; … … 248 250 sprintf(buff,"%s %20lg / %s", cle, val, tcom); 249 251 } 250 else sprintf(buff,"%s %20l d/ ", cle, val);252 else sprintf(buff,"%s %20lg / ", cle, val); 251 253 buff[strlen(buff)] = ' '; 252 254 nkeya_++;
Note:
See TracChangeset
for help on using the changeset viewer.