0

我有许多使用朗伯变换的图表。其中一些使用 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。这让我相信我错过了一些东西。

我从像素到地理坐标转换中遗漏了什么?

以下是其他一些 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
4

1 回答 1

0

对于像素到纬度/经度,输入向量经过 [x, y, 0, 1] 变换。将变换矩阵乘以向量后得到朗伯变换的 (x, y) 坐标。

反之,进行朗伯逆变换得到 (x, y) 坐标。然后,您必须求解 2x2 线性方程组才能获得像素位置。

于 2014-03-17T20:54:57.143 回答