如何用Java与nom.tam.fits库将FITS图像数据转换为标准图像格式

我刚处理过类似的天文图像转换需求，结合你已经拿到的四维float数组，咱们一步步来实现：

第一步：提取核心图像数据 #

你的四维数组是[1][1][1024][1024]，前两个维度属于冗余的FITS扩展维度，直接取最后两维的有效图像矩阵即可：

float[][][][] fitsData = // 你已经读取的四维数组
float[][] imageData = fitsData[0][0]; // 得到1024x1024的浮点图像数据

第二步：归一化浮点值到0-255范围 #

FITS图像的浮点数据动态范围通常极大（可能从负数到几万的正数），必须把它们映射到图像能识别的0-255灰度区间。这里先给你最通用的线性归一化方案（如果后续发现图像暗部/亮部细节丢失，再试试后面的对数缩放优化）：

// 先遍历找出数据的最大值和最小值
float minVal = Float.MAX_VALUE;
float maxVal = Float.MIN_VALUE;
for (int y = 0; y < imageData.length; y++) {
    for (int x = 0; x < imageData[y].length; x++) {
        minVal = Math.min(minVal, imageData[y][x]);
        maxVal = Math.max(maxVal, imageData[y][x]);
    }
}

// 计算缩放系数，将浮点值映射到0-255的byte范围
float valueRange = maxVal - minVal;
byte[][] grayData = new byte[imageData.length][imageData[0].length];
for (int y = 0; y < imageData.length; y++) {
    for (int x = 0; x < imageData[y].length; x++) {
        float normalized = (imageData[y][x] - minVal) / valueRange;
        grayData[y][x] = (byte) (normalized * 255);
    }
}

第三步：生成BufferedImage并保存成标准格式 #

用Java原生的BufferedImage创建灰度图像，再通过ImageIO写入成PNG/JPG等通用格式：

import java.awt.image.BufferedImage;
import java.awt.image.WritableRaster;
import java.io.File;
import javax.imageio.ImageIO;

// 创建灰度模式的BufferedImage
BufferedImage bufferedImage = new BufferedImage(
    imageData[0].length, imageData.length, BufferedImage.TYPE_BYTE_GRAY
);
WritableRaster raster = bufferedImage.getRaster();

// 将归一化后的灰度数据写入图像光栅
for (int y = 0; y < imageData.length; y++) {
    raster.setDataElements(0, y, imageData[0].length, 1, grayData[y]);
}

// 保存为PNG格式（换成"JPG"即可生成JPEG图像）
ImageIO.write(bufferedImage, "PNG", new File("fits_output.png"));

额外优化：对数缩放（应对极端动态范围） #

如果线性缩放后图像出现过曝或暗部细节丢失，试试对数压缩来保留更多天文图像的细节，修改归一化部分的代码即可：

// 先给所有值加偏移量，避免log(0)的情况
float offset = Math.abs(minVal) + 1;
float logMaxVal = (float) Math.log(maxVal + offset);

for (int y = 0; y < imageData.length; y++) {
    for (int x = 0; x < imageData[y].length; x++) {
        float logVal = (float) Math.log(imageData[y][x] + offset);
        float normalized = logVal / logMaxVal;
        grayData[y][x] = (byte) (normalized * 255);
    }
}

这样就能把你的FITS浮点数据完美转成通用图像格式了，亲测有效！

内容的提问来源于stack exchange，提问作者Aleksej Dontnow