を使用して、画像をピクセルデータで新しい適合ファイルに書き込んで埋めようとしていますcfitsio
。どこが間違っているのかわかりませんが、すべてのデータが Fits ファイルに書き込まれていません。データbinary vi
をチェックアウトしたところ、画像データが正しく作成されませんでした。NASA の Web サイトの 2番目の画像 HST WFPC 2 適合ファイルを使用しました。これは私が書いたコードです:
#include <string.h>
#include <stdio.h>
#include "fitsio.h"
int main() {
fitsfile *fptr;
int status = 0, i;
fits_open_file(&fptr, "WFPC2ASSNu5780205bx.fits", READONLY, &status);
long naxes[2];
fits_get_img_size(fptr, 3, naxes, &status);
fitsfile *ofptr;
fits_create_file(&ofptr, "o_nasa.fits", &status);
fits_copy_header(fptr, ofptr, &status);
long fp[2] = {1, 1};
long nelements = naxes[0];
float arr[nelements];
for (i = 0; i < nelements; i++) {
arr[i] = 100;
}
int ii, jj, kk;
for (ii = 1; ii <= naxes[0]; ii++) {
for (jj = 1; jj <= naxes[1]; jj++) {
fits_write_pix(ofptr, TFLOAT, fp, nelements, arr, &status); // this is not working
}
}
fits_close_file(ofptr, &status);
fits_close_file(fptr, &status);
return status;
}