Записать данные пикселей в файл FITS Cfitsio

Question

Я пытаюсь записать и заполнить изображение данными о пикселях в файл нового соответствия, используя cfitsio, Я не уверен, где я делаю неправильно, но все данные не записываются в файл соответствий. я использовал binary vi проверить данные и данные изображения не созданы должным образом. Я использовал второе изображение 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;
} 

Answer 1

Я узнал, что не так с кодом. Я только что опубликовал ответ, так как, кажется, не так много сообщества работает над CFITSIO, Код неверен в строках в двойном цикле for с переменными ii, jj, я думал fpпервый пиксельный массив, указывающий, где массив изображений начинает записываться, используется один раз, когда начинается вся итерация, но вместо этого write_pix использует его на протяжении итераций, чтобы знать, куда писать. Мой код также пытается написать / обновить 1, 1 pixel и кажется NULL, Правильный код:

for (fp[1] = 1; fp[1] <= naxes[0]; fp[1]++) { 
    for (fp[0] = 1; fp[0] <= naxes[1]; fp[0]++) { 
        fits_write_pix(ofptr, TFLOAT, fp, nelements, arr, &status); 
    }
}