我的任务是读取fits
文件并在一个系统中显示图像(如果存在)。该系统接受图像数据作为像素阵列,该像素阵列既可以是8位的单个块(对于灰度),也可以是3个8位的块(对于RGB)。我正在尝试从CCfits
库读取测试中生成的文件,但对我来说,数据对我已检查的所有fits
查看器中的显示方式没有任何意义。
因此,这是在我正在使用的示例文件中生成数据数组的代码:
std::valarray<int> array(nelements);
for (int i = 0; i < numberOfRows; ++i)
{
array[std::slice(vectorLength*static_cast<int>(i),vectorLength,1)] = row + i;
}
这是标题:
SIMPLE = T / file does conform to FITS standard
BITPIX = 16 / number of bits per data pixel
NAXIS = 2 / number of data axes
NAXIS1 = 300 / length of data axis 1
NAXIS2 = 200 / length of data axis 2
EXTEND = T / FITS dataset may contain extensions
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
BZERO = 32768 / offset data range to that of unsigned short
BSCALE = 1 / default scaling factor
EXPOSURE= 1500 / Total Exposure Time
OMEGA = (-0.5, 0.8660254) / Complex cube root of 1
NEWVALUE= 42 / Test of adding keyword to different extension
STRING = ' Rope ' / trailing blank test 1
STRING2 = 'Cord ' / trailing blank test 2
END
这就是图像在gimp中的显示方式,例如:
因此,通过标题我可以理解,每个数组项将具有16位整数。具有该值的
BZERO
表示它是16位无符号整数。到目前为止没有问题。从图像生成开始,我假设每个数组项将代表一个像素,并且每个值都将代表灰度中的颜色。ccfits库中的以下代码生成的数字从0到最大498,这实际上与16位整数的边距很远。话虽这么说,这怎么从黑色变成白色?如何按比例缩放或变换这些数字以表示像素和颜色?
最佳答案
FITS专家在这里,但我不太了解您的问题。首先,您在此处发布的代码不完整-变量row
定义不明确。根本不清楚目的,但我想它是用一些渐变填充300x200的数组。 FITS阵列就是阵列。某些数据类型的某些值的N维数组。
尽管有些观看者是专为天文图像(如古老的ds9或更现代的Glue)设计的,但这些值的解释方式完全取决于特定的应用程序。
您写道,您只是在GIMP中打开了文件(实际上,我印象深刻的是它可以读取FITS文件)。由于缺少更好的选择,GIMP似乎正在做的所有事情就是在黑色(可能从0开始)和白色之间线性缩放数组中的值。
关于c++ - 如何在FITS文件中的数组数据中表示像素?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/49816222/