在C#中基于GeoTransform计算GMapControl叠加GeoTIFF的变换参数

我有一份GeoTIFF文件，需要叠加到使用GMap.NET的GMapControl组件展示的地图上。目前通过自定义GMapMarker，以BitmapImage作为源实现了图像叠加，但现在用的是固定变换参数，无法实现精准对齐。希望利用GeoTIFF的GeoTransform元数据，计算出RotateTransform、ScaleTransform、TranslateTransform这三个变换的正确参数，实现精准叠加。

当前使用固定参数的示例代码：

_marker = new GMapMarker(new PointLatLng(topLeftY, topLeftX))
{
    Shape = new System.Windows.Controls.Image
    {
        Source = new BitmapImage(new Uri("teste.tif")),
        Width = 100,
        Height = 100,
        
        RenderTransform = new TransformGroup
        {
            Children = new TransformCollection
            {
                new RotateTransform(rotationAngle), // 固定旋转角度
                new ScaleTransform(1, -1), // 固定缩放
                new TranslateTransform(0.5, -1) // 固定平移
            }
        }
    }
};

此前尝试过GdalWarp和OpenCV，但C#环境下GdalWarp调用卡顿，原生OpenCV不支持GeoTIFF且读取卡顿，同时不想将文件转为PNG格式，因此希望通过调整GMapMarker的变换参数解决问题。

1. 读取GeoTIFF的GeoTransform元数据 #

GeoTIFF的GeoTransform是一个6元素数组（[a, b, c, d, e, f]），对应仿射变换公式：

X = a + b*col + c*row
Y = d + e*col + f*row

各参数含义：

• a: 左上角像素的X坐标（地理坐标系）
• b: X方向的像素分辨率（每个像素代表的地理单位长度，正方向为东）
• c: X方向旋转分量（非0表示图像有旋转）
• d: 左上角像素的Y坐标（地理坐标系）
• e: Y方向旋转分量（非0表示图像有旋转）
• f: Y方向的像素分辨率（每个像素代表的地理单位长度，正方向为北，通常为负，因图像行号向下增加时Y坐标减小）

可以用GDAL的C#绑定轻量读取元数据（仅获取GeoTransform，不处理图像内容）：

using OSGeo.GDAL;

// 初始化GDAL（全局只做一次）
Gdal.AllRegister();

// 只读打开GeoTIFF文件
Dataset ds = Gdal.Open("teste.tif", Access.GA_ReadOnly);
double[] geoTransform = new double[6];
ds.GetGeoTransform(geoTransform);
// 提取图像宽高用于设置Marker尺寸
int imgWidth = ds.RasterXSize;
int imgHeight = ds.RasterYSize;
ds.Dispose();

// 拆分参数
double a = geoTransform[0];
double b = geoTransform[1];
double c = geoTransform[2];
double d = geoTransform[3];
double e = geoTransform[4];
double f = geoTransform[5];

2. 计算ScaleTransform参数 #

缩放比例需将地理单位转换为GMapControl的屏幕像素单位，先获取当前地图分辨率（每像素代表的地理单位长度）：

// 获取当前缩放级别下的地面分辨率（米/像素，假设地图用WGS84投影）
double mapResolution = gMapControl.MapProvider.Projection.GetGroundResolution(gMapControl.Zoom, gMapControl.CurrentPosition.Lat);

再计算X、Y方向的缩放因子：

• X方向：scaleX = b / mapResolution（保留b的符号，匹配地理坐标系X方向）
• Y方向：scaleY = f / mapResolution（保留f的负号，实现图像Y轴与地理坐标系Y轴的翻转）

3. 计算RotateTransform参数 #

利用GeoTransform的旋转分量计算旋转角度，公式为：

// 计算旋转角度（弧度转角度）
double rotationAngle = Math.Atan2(c, b) * (180 / Math.PI);

• 正值为逆时针旋转，负值为顺时针旋转，与GMapControl的旋转方向匹配。

4. 计算TranslateTransform参数 #

平移参数需将地理坐标偏差转换为屏幕像素偏移：

// Marker的地理位置（原代码中的左上角坐标）
PointLatLng markerGeoPos = new PointLatLng(topLeftY, topLeftX);
// Marker对应的屏幕坐标
Point markerScreenPos = gMapControl.FromLatLngToLocal(markerGeoPos);

// GeoTransform定义的图像左上角实际地理坐标
PointLatLng actualGeoPos = new PointLatLng(d, a);
// 实际左上角对应的屏幕坐标
Point actualScreenPos = gMapControl.FromLatLngToLocal(actualGeoPos);

// 计算平移偏移量
double translateX = actualScreenPos.X - markerScreenPos.X;
double translateY = actualScreenPos.Y - markerScreenPos.Y;

5. 整合变换到GMapMarker #

注意：WPF的TransformGroup变换应用顺序为从后往前，需按「平移→缩放→旋转」的顺序添加（或根据实际逻辑调整），同时设置旋转原点为图像中心避免偏移：

_marker = new GMapMarker(markerGeoPos)
{
    Shape = new System.Windows.Controls.Image
    {
        Source = new BitmapImage(new Uri("teste.tif")),
        Width = imgWidth,
        Height = imgHeight,
        // 设置旋转原点为图像中心
        RenderTransformOrigin = new Point(0.5, 0.5),
        RenderTransform = new TransformGroup
        {
            Children = new TransformCollection
            {
                new RotateTransform(rotationAngle),
                new ScaleTransform(scaleX, scaleY),
                new TranslateTransform(translateX, translateY)
            }
        }
    }
};

// 添加Marker到地图
gMapControl.Markers.Add(_marker);

关键注意事项 #

• 确保GeoTIFF与GMapControl的坐标系一致（如均为WGS84经纬度），不一致需先做坐标转换。
• 当地图缩放级别变化时，需在gMapControl.ZoomChanged事件中重新计算缩放和平移参数。
• GDAL初始化仅需执行一次，避免重复初始化影响性能。

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