0

我有一些从奥地利的雷达中提取的图像和一个投影公式(练习说投影公式将像素位置映射到地理坐标),它将像素转换成(phi, lambda)练习状态为球角的像素。

PRODUCT-PH-FORMS   =     "(1.570796326795-2*atan(0.404026225835*pow((1/sin(acotan(5899179.2/\
((5676531.8*pow((tan((1.570796326795-(v1*0.017453292520))/2)/\
0.373884679485),0.737361597616)*sin(0.012869387656*(v0-13.333333333))-\
v3*v4)+COL*v4)-((6380000.0*(0.924636243305-0.889738536848*pow((tan((\
1.570796326795-(v1*0.017453292520))/2)/0.373884679485),0.737361597616)*\
cos(0.012869387656*(v0-13.333333333)))+v2*v5)-ROW*v5)/((5676531.8*\
pow((tan((1.570796326795-(v1*0.017453292520))/2)/0.373884679485),\
0.737361597616)*sin(0.012869387656*(v0-13.333333333))-v3*v4)+COL*v4))))*\
((5676531.8*pow((tan((1.570796326795-(v1*0.017453292520))/2)/\
0.373884679485),0.737361597616)*sin(0.012869387656*(v0-13.333333333))-\
v3*v4)+COL*v4)/6380000.0*0.737361597616/sin(0.767944870878),1/\
0.737361597616)))"

PRODUCT-LA-FORMS   = f(pixel row, pixel col, v0, v1, v2, v3, v4, v5)

但这给出了非常小的值,这是不可能(latitude, longitude)的,所以我正在寻找球面角之间的(phi,lambda)转换(latitude,longitude)

对于我需要提供的结果(latitude, longitude),如何进行这种转换?该练习还为我提供了一些我不理解的数据,这可能是必要的:

v0=8.194 (reference lon)
v1=50.437 (reference lat)
v2 = 0
v3 = 0
v4=1018.18 (meter / pixel)
v5=1018.18 (meter / pixel)

此变量用于转换公式,但也可用于获得(latitude, longitude)


编辑

我试图直接将我从转换中获得的值视为(纬度,经度),*180/pi但我得到的坐标不在奥地利

4

0 回答 0