0

我正在尝试使用cfitsio. 我不确定我在哪里做错了，但所有数据都没有写入 fit 文件。我用来binary vi检查数据和图像数据没有正确创建。我使用了来自NASA 网站的第二张图片 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;
}
4

1 回答 1

0

我发现代码有什么问题。我刚刚发布了答案，因为似乎没有很多社区在工作CFITSIO。代码在带有 variables 的双 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); 
    }
}
于 2017-11-17T23:25:18.837 回答