Python 中的 pyproj 为 NAD_1927_StatePlane_Wisconsin_South_FIPS_4803

如何解决Python 中的 pyproj 为 NAD_1927_StatePlane_Wisconsin_South_FIPS_4803

我刚开始使用地理空间数据。我在使用 pyproj 从具有不同坐标系统的文件中获取正确的经度和纬度时遇到问题。我正在处理一个包含密尔沃基城市界限的文件，该文件位于州平面坐标 NAD_1927_StatePlane_Wisconsin_South_FIPS_4803 中，我想要 wgs84 坐标，它应该只是普通的经度和纬度。我想转换它，因为使用我正在使用的 folium 等工具会更容易工作。

https://data.milwaukee.gov/dataset/corporate-boundary shp文件中的第一个坐标是[2516440.160073474,440481.9301029742]

它有prj文件：

PROJCS["NAD_1927_StatePlane_Wisconsin_South_FIPS_4803",GEOGCS["GCS_north_American_1927",DATUM["D_north_American_1927",SPHEROID["Clarke_1866",6378206.4,294.9786982]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.017453292519943295]],PROJECTION["LAmbert_Conformal_Conic"],ParaMETER["False_Easting",2000000.0],ParaMETER["False_northing",ParaMETER["Central_Meridian",-90.0],ParaMETER["Standard_Parallel_1",42.73333333333333],ParaMETER["Standard_Parallel_2",44.06666666666667],ParaMETER["Latitude_Of_Origin",42.0],UNIT["Foot_US",0.30480060960121924]]

根据 https://spatialreference.org 这将具有代码 epsg6948，但 pyproj 无法识别，因此我构建了一个自定义字符串。查看文档和示例后，我想出了代码：

import pyproj
from pyproj import Transformer
cust = "+proj=lcc +ellps=clrk66 +datum=NAD27 +lat_1=42.73333333333333 +lat_2=44.06666666666667 +lat_0=42.0 +lon_0=-90.0 +x_0=2000000.0 +y_0=0.0 +units=us-ft +no_defs"
t = Transformer.from_crs(cust,"wgs84")
x1,y1 = [2516440.160073474,440481.9301029742]
lat,lon = t.transform(x1,y1)
print(lat,",lon)

根据站点 https://www.earthpoint.us/StatePlane.aspx，此代码应返回坐标： 纬度经度 43.1900328,-087.9451397

但是代码打印出来了： 42.206953018562416,-105.00939853775037 这是在怀俄明州的某个地方，而不是威斯康星州的密尔沃基。我已经检查了几次，但看不出我做错了什么。如果有人可以提供帮助，我将不胜感激。

解决方法

您可以使用 WKT 字符串来创建转换器：

>>> from pyproj import Transformer
>>> wkt = """PROJCS["NAD_1927_StatePlane_Wisconsin_South_FIPS_4803",...        GEOGCS["GCS_North_American_1927",...               DATUM["D_North_American_1927",...                     SPHEROID["Clarke_1866",6378206.4,294.9786982]],...               PRIMEM["Greenwich",0.0],...               UNIT["Degree",0.017453292519943295]],...        PROJECTION["Lambert_Conformal_Conic"],...        PARAMETER["False_Easting",2000000.0],...        PARAMETER["False_Northing",...        PARAMETER["Central_Meridian",-90.0],...        PARAMETER["Standard_Parallel_1",42.73333333333333],...        PARAMETER["Standard_Parallel_2",44.06666666666667],...        PARAMETER["Latitude_Of_Origin",42.0],...        UNIT["Foot_US",0.30480060960121924]]"""
>>> transformer = Transformer.from_crs(wkt,"WGS84")
>>> transformer.transform(2516440.160073474,440481.9301029742)
(43.19212933447268,-88.06334977811528)

旁注，我认为这是您想要的 EPSG 代码：

<Projected CRS: EPSG:32054>
Name: NAD27 / Wisconsin South
Axis Info [cartesian]:
- X[east]: Easting (US survey foot)
- Y[north]: Northing (US survey foot)
Area of Use:
- name: United States (USA) - Wisconsin - counties of Adams; Calumet; Columbia; Crawford; Dane; Dodge; Fond Du Lac; Grant; Green; Green Lake; Iowa; Jefferson; Juneau; Kenosha; La Crosse; Lafayette; Manitowoc; Marquette; Milwaukee; Monroe; Ozaukee; Racine; Richland; Rock; Sauk; Sheboygan; Vernon; Walworth; Washington; Waukesha; Waushara; Winnebago.
- bounds: (-91.43,42.48,-86.95,44.33)
Coordinate Operation:
- name: Wisconsin CS27 South zone
- method: Lambert Conic Conformal (2SP)
Datum: North American Datum 1927
- Ellipsoid: Clarke 1866
- Prime Meridian: Greenwich

>>> transformer = Transformer.from_crs("EPSG:32054",-88.06334977811528)