我正在编写一个python脚本,该脚本读取包含x和y坐标的多个XML文件,并将它们全部合并为一个csv文件。纬度和经度是csv中的必填字段,但我很难将俄亥俄州北州飞机usFt中的x,y坐标转换为WGS84。
>>> p = Proj(r'+proj=lcc +lat_1=41.7 +lat_2=40.43333333333333 +lat_0=39.66666666666666 +lon_0=-82.5 +x_0=600000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs') #Nad83 State Plane Ohio North US Feet Proj object using parameters
>>> p(739400.91,2339327.3,inverse=True)
(-80.138057868777224, 60.278230707978487)
>>> p1 = Proj(init="epsg:3734") #Nad83 State Plane Ohio North US Feet Proj object using EPSG code
>>> p1(739400.91,2339327.3,inverse=True)
(-80.138057868777224, 60.278230707978487)
上面的两种方法都返回相同的结果,但是这很长一段时间都位于哈德逊湾中。当我在ArcMap中绘制坐标时,正确的纬度为:-81.142311,41.688205。
*请注意,所有纬度均以经纬度提供,因为这是Proj使用的顺序
有谁知道为什么我会从Proj.4和pyproj中得到错误的坐标?