使用 healpix 读取 FITS 文件时出现分段错误

Question

我正在尝试使用 C 代码通过 healpix 打开一个 FITS 文件：

#define _GNU_SOURCE
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include "chealpix.h"
#include "fitsio.h"
#include <string.h>
#include "fitsio2.h"

int main(int argc, char **argv){

        char *coordsys="G";
        char *order="RING";
        long nside;
        long *p_nside=&nside;
        nside=512.;
        float *map;
        map=read_healpix_map("file.fits",p_nside,coordsys,order);
         
    return(0);
}

这个代码在使用gdb时返回一个segmentation fault (core dumped)，我明白了

Program received signal SIGSEGV, Segmentation fault.
ffgkys (status=0x7fffffffdaf4, comm=<optimized out>, value=0x55555555a004 "G", keyname=<optimized out>, 
    fptr=<optimized out>) at getkey.c:808
808     value[0] = '\0';

这是我第一次使用 healpix，我不知道是什么导致了问题。有谁知道可能出了什么问题？

Answer 1

根据read_healpix_map()的文档：

ordering 和 coordsys 是 output 参数 - 它们是指向函数将写入结果的空间的指针。您已将指针传递给内存中的字符串常量，这些常量通常会被标记为只读，并且在任何情况下都可能不足以接收结果。

nside 也是一个输出参数，将其初始化为 512 没有任何意义。你不需要p_nside；你可以简单地传递&nside

它不是最清晰的接口；具体而言，尚不清楚coordsys 是字符串还是单个char。该文档使用单引号，因此您可能会假设 char 但它也会为 ordering 这样做，这显然是一个字符串。假设coordsys 是一个字符串是最安全的——它不会造成任何伤害。

    char  coordsys[] = "G";    // This will be overwritten
    char  order[] = "NESTED" ; // This will be overwritten, 
                                  initialise to largest expected string
                                  to ensure sufficient space is allocated
    long nside = 0 ;
    float* map = read_healpix_map( "file.fits", &nside, coordsys, order ) ;

    printf( "nside:%ld, coordsys:%s, order:%s", nside, coordsys, order ) ;

---

后置脚本：查看source code，coordsys 和ordering 都传递给fits_read_key()，因此coordsys 被视为字符串。