将法定描述(部分,乡镇,范围)转换为纬度和经度

时间:2018-10-10 20:12:39

标签: python coordinates gdal data-conversion gis

问题在于如何将乡镇系统中的法律说明转换为纬度和经度。

这里是Township System描述的链接。下面是我想做的一个例子。

该图显示了乡镇系统中的包裹:Township Example。假设我要在红线圈出的小镇中找到的位置在第23节的东北四分之一处。“法律说明”可能看起来像NE23-5-3W。 “ W”表示在本初子午线的西侧。

我接下来要做的是从描述中提取信息。例如,我可以将其拆分为列表['NE',23、5、3,'W']。我可以计算出偏移量(距本初子午线和加拿大和美国边界的距离-这是乡镇系统的起点-原点)。在Google Map上,原点将位于红色圆圈Origin周围。

在此GitHub回购中,这是我执行此转换的脚本。我希望脚本的结果与我发现的here地块的位置相匹配。包裹是t05r03w。

该脚本可以基本定位乡镇,但准确性不如预期。如果您有任何想法可以解决精度问题,请给我一些建议吗?

非常感谢您。

0 个答案:

没有答案