如何用VB/C#修改TIFF元数据，实现3D球面地形数据转换与缩放

核心逻辑 #

无需修改TIFF像素数据，仅通过修改地理空间元数据标签就能实现经纬度范围、高度范围的缩放，以及球面半径的转换。重点操作以下几类TIFF标签：

• ModelPixelScaleTag（标签号33550）：控制像素对应的地理坐标比例
• ModelTiepointTag（标签号33922）：控制地理坐标的基准点
• GeoKeyDirectoryTag（标签号34735）：包含椭球半径等地理参考参数

C# 实现代码 #

依赖LibTiff.Net库（可通过NuGet直接安装），代码示例如下：

using LibTiff.Net;

// 定义转换参数
double originalSphereRadius = 6371000; // 原南极球面半径（根据实际数据调整）
double targetEquatorRadius = 6378137; // 目标赤道球面半径
double heightScaleFactor = 0.5; // 高度缩放系数（按需调整）

using (Tiff tiff = Tiff.Open("antarctic_terrain.tif", "r+"))
{
    // 修改像素地理比例
    float[] pixelScale = new float[3];
    if (tiff.GetField(TiffTag.ModelPixelScale, out pixelScale))
    {
        double radiusRatio = targetEquatorRadius / originalSphereRadius;
        pixelScale[0] *= (float)radiusRatio; // 经度方向缩放
        pixelScale[1] *= (float)radiusRatio; // 纬度方向缩放
        pixelScale[2] *= (float)heightScaleFactor; // 高度缩放
        tiff.SetField(TiffTag.ModelPixelScale, pixelScale);
    }

    // 修改地理基准点
    double[] tiePoints = new double[6];
    if (tiff.GetField(TiffTag.ModelTiepoint, out tiePoints))
    {
        double radiusRatio = targetEquatorRadius / originalSphereRadius;
        tiePoints[3] *= radiusRatio; // 基准经度调整
        tiePoints[4] *= radiusRatio; // 基准纬度调整
        tiePoints[5] *= heightScaleFactor; // 基准高度调整
        tiff.SetField(TiffTag.ModelTiepoint, tiePoints);
    }

    // 修改椭球半径参数
    ushort[] geoKeys = null;
    if (tiff.GetField(TiffTag.GeoKeyDirectory, out geoKeys))
    {
        // 查找并替换椭球长半轴（GeoKey 3008）
        for (int i = 0; i < geoKeys.Length; i += 4)
        {
            if (geoKeys[i] == 3008)
            {
                geoKeys[i + 3] = (ushort)targetEquatorRadius;
                break;
            }
        }
        tiff.SetField(TiffTag.GeoKeyDirectory, geoKeys);
    }

    // 保存修改
    tiff.Checkpoint();
}

VB.NET 实现代码 #

Imports LibTiff.Net

' 定义转换参数
Dim originalSphereRadius As Double = 6371000 ' 原南极球面半径（按需调整）
Dim targetEquatorRadius As Double = 6378137 ' 目标赤道球面半径
Dim heightScaleFactor As Double = 0.5 ' 高度缩放系数

Using tiff As Tiff = Tiff.Open("antarctic_terrain.tif", "r+")
    ' 修改像素地理比例
    Dim pixelScale(2) As Single
    If tiff.GetField(TiffTag.ModelPixelScale, pixelScale) Then
        Dim radiusRatio As Double = targetEquatorRadius / originalSphereRadius
        pixelScale(0) *= CSng(radiusRatio)
        pixelScale(1) *= CSng(radiusRatio)
        pixelScale(2) *= CSng(heightScaleFactor)
        tiff.SetField(TiffTag.ModelPixelScale, pixelScale)
    End If

    ' 修改地理基准点
    Dim tiePoints(5) As Double
    If tiff.GetField(TiffTag.ModelTiepoint, tiePoints) Then
        Dim radiusRatio As Double = targetEquatorRadius / originalSphereRadius
        tiePoints(3) *= radiusRatio
        tiePoints(4) *= radiusRatio
        tiePoints(5) *= heightScaleFactor
        tiff.SetField(TiffTag.ModelTiepoint, tiePoints)
    End If

    ' 修改椭球半径参数
    Dim geoKeys() As UShort = Nothing
    If tiff.GetField(TiffTag.GeoKeyDirectory, geoKeys) Then
        For i As Integer = 0 To geoKeys.Length - 1 Step 4
            If geoKeys(i) = 3008 Then
                geoKeys(i + 3) = CUShort(targetEquatorRadius)
                Exit For
            End If
        Next
        tiff.SetField(TiffTag.GeoKeyDirectory, geoKeys)
    End If

    ' 保存修改
    tiff.Checkpoint()
End Using

注意事项 #

• 先通过TiffView等工具查看目标TIFF的元数据标签，确认存在上述标签后再执行修改
• 若数据使用自定义地理参考标签，需调整代码中对应的标签编号和处理逻辑
• 操作前务必备份原TIFF文件，避免数据损坏

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