将纬度/经度转换为Alberta 10 TM Projection

我需要在Alberta 10 TM Projection中将纬度/经度坐标转换为Easting/Northing坐标.

10 TM投影类似于UTM,但它是加拿大阿尔伯塔省的定制投影.我认为(付出一些努力)我可以自己编码,但如果已经完成,我宁愿不重新发明轮子.

抓住PROJ.4制图投影库 - 开源库.

10TM的建议参数:

+ proj = tmerc + lon_0 = -115 + k_0 = 0.9992 + x_0 = 500000 + datum = NAD27

根据这篇文章,您可能需要:

如果您的10TM数据被引用到NAD83数据(而不是NAD27/clrk66),则将椭圆更改为GRS80.如果您的Alberta的10TM坐标小于5,000,000(AltaLIS"标准"),您可能还需要将false northing(y_0)更改为-5000000.

我应该指出,PROJ.4是图书馆可以为任何类型的地理坐标系统改造.几乎没有任何转型它无法做到.

如果你喜欢这些东西,我还建议阅读John Snyder的Map Projections-A Working Manual(平装本) ......这是经典之作.:)(修复了链接)