如何在Matlab中确定.fits格式天文图像的配色方案？

问题背景 #

• 用Unistellar eVscope2拍摄的夜空星场图像保存为FITS格式
• IrfanView打开显示为正常灰度，但通过Matlab执行data=fitsread('stars.fits');读取为二维数组后，调用imagesc(data)显示时呈现蓝色调
• 图像数据值范围为0到32767（2^15-1），多数数值集中在4000到12000区间
• 推测Unistellar采用16位RGB配色映射，IrfanView可正确识别，但Matlab默认处理逻辑不匹配

解决方法 #

1. 强制切换为灰度配色映射 #

Matlab的imagesc默认使用parula配色（偏蓝），直接指定灰度色表即可解决：

data = fitsread('stars.fits');
imagesc(data);
colormap(gray); % 切换到灰度配色
colorbar; % 可选，添加色阶参考

2. 优化动态范围提升显示细节 #

由于多数数据集中在4000-12000，手动限定显示范围可以避免细节被压缩：

data = fitsread('stars.fits');
imagesc(data, [4000 12000]); % 仅展示4000到12000区间的数值
colormap(gray);
axis image; % 保持图像原始比例

3. 检查并读取多通道RGB数据 #

如果FITS文件实际是16位RGB格式，可先查看文件结构确认通道数：

info = fitsinfo('stars.fits');
disp(info.Image); % 查看图像维度、通道数等信息

若为三通道RGB，需分别读取并组合：

red = fitsread('stars.fits', 'Image', 1);
green = fitsread('stars.fits', 'Image', 2);
blue = fitsread('stars.fits', 'Image', 3);
rgb_image = cat(3, red, green, blue);
imagesc(rgb_image);
axis image;

内容的提问来源于stack exchange，提问作者Tom Mozdzen