如何使用模型变换标记矩阵对Lambert Geotiff进行纬度/经度？

Question

我有很多使用朗伯变换的图表。其中一些使用Tie Point标签。为此，我找到了Lat/Long<=>X,Y转换。

我也有一些使用模型转换标签的图表。到目前为止，我无法弄清楚如何进行点转换。

Geotiff文档说Model-Coordinates = Matrix * Image Coordinates

一个图表有这个矩阵。

798.320542  0.009593987 0   -6084843.93
0.020321671 -798.350061 0   3716484.45
0                     0 0   0
0                     0 0   1 

通过Photoshop我找到了

9572    2107 (X/Y) corresponds to 50N on the prime meridian and that
4957    410 corresponds to 60N 60W

当我乘以矩阵（矢量转置）时，我得到了

M * [9572   2107 0 1] = [1556700.513 2034555.391 0 1]
M * [4957   410 0 1] = [-2127565.07 3389261.66 0 1]

我没有得到像我期待的度数（弧度）。也没有扭转X / Y.这让我相信我错过了一些东西。

我从Pixel to Geo坐标转换中遗漏了什么？

以下是其他一些Geotiff数据

Image Width                     : 24000
Image Height                    : 7998
Proj Linear Units               : Linear Meter
Proj Std Parallel 1             : 45
Proj Std Parallel 2             : 33
Proj Nat Origin Long            : -95
Proj False Origin Long          : -95
Proj False Origin Lat           : 39
Proj False Origin Easting       : 0
Proj False Origin Northing      : 0

Answer 1

对于像素到纬度/经度，输入矢量为[x，y，0,1]变换。将变换矩阵乘以向量后，给出朗伯变换的（x，y）坐标。

反过来，执行逆lambert变换以获得（x，y）坐标。然后你必须解决一个2x2线性方程组来获得像素位置。