Question

我正在使用PyProj，并且试图找到一种方法来获取州平面坐标系的区域中心的中央子午线和纬度。任何帮助将不胜感激！

Answer 1

我相信您的回答是使用CRS类为https://pyproj4.github.io/pyproj/v2.2.0rel/api/crs.html的pyproj 2.2.0

>>> import pyproj
>>> pyproj.__version__
'2.2.0'
>>> crs = pyproj.CRS("ESRI:102719")
>>> crs
<Projected CRS: ESRI:102719>
Name: NAD_1983_StatePlane_North_Carolina_FIPS_3200_Feet
Axis Info [cartesian]:
- E[east]: Easting (US survey foot)
- N[north]: Northing (US survey foot)
Area of Use:
- name: USA - North Carolina
- bounds: (-84.33, 33.83, -75.38, 36.59)
Coordinate Operation:
- name: NAD_1983_StatePlane_North_Carolina_FIPS_3200_Feet
- method: Lambert Conic Conformal (2SP)
Datum: North American Datum 1983
- Ellipsoid: GRS 1980
- Prime Meridian: Greenwich
>>> crs.coordinate_operation
CONVERSION["NAD_1983_StatePlane_North_Carolina_FIPS_3200_Feet",
    METHOD["Lambert Conic Conformal (2SP)",
        ID["EPSG",9802]],
    PARAMETER["Latitude of false origin",33.75,
        ANGLEUNIT["degree",0.0174532925199433],
        ID["EPSG",8821]],
    PARAMETER["Longitude of false origin",-79,
        ANGLEUNIT["degree",0.0174532925199433],
        ID["EPSG",8822]],
    PARAMETER["Latitude of 1st standard parallel",34.3333333333333,
        ANGLEUNIT["degree",0.0174532925199433],
        ID["EPSG",8823]],
    PARAMETER["Latitude of 2nd standard parallel",36.1666666666667,
        ANGLEUNIT["degree",0.0174532925199433],
        ID["EPSG",8824]],
    PARAMETER["Easting at false origin",2000000.00261667,
        LENGTHUNIT["US survey foot",0.304800609601219],
        ID["EPSG",8826]],
    PARAMETER["Northing at false origin",0,
        LENGTHUNIT["US survey foot",0.304800609601219],
        ID["EPSG",8827]],
    ID["ESRI",102719]]
>>> crs.coordinate_operation.params
[Param(name=Latitude of false origin, auth_name=EPSG, code=8821, value=33.75, unit_name=degree, unit_auth_name=EPSG, unit_code=9102, unit_category=angular), Param(name=Longitude of false origin, auth_name=EPSG, code=8822, value=-79.0, unit_name=degree, unit_auth_name=EPSG, unit_code=9102, unit_category=angular), Param(name=Latitude of 1st standard parallel, auth_name=EPSG, code=8823, value=34.3333333333, unit_name=degree, unit_auth_name=EPSG, unit_code=9102, unit_category=angular), Param(name=Latitude of 2nd standard parallel, auth_name=EPSG, code=8824, value=36.1666666667, unit_name=degree, unit_auth_name=EPSG, unit_code=9102, unit_category=angular), Param(name=Easting at false origin, auth_name=EPSG, code=8826, value=2000000.00262, unit_name=US survey foot, unit_auth_name=EPSG, unit_code=9003, unit_category=linear), Param(name=Northing at false origin, auth_name=EPSG, code=8827, value=0.0, unit_name=US survey foot, unit_auth_name=EPSG, unit_code=9003, unit_category=linear)]
>>> crs.coordinate_operation.params[0]
Param(name=Latitude of false origin, auth_name=EPSG, code=8821, value=33.75, unit_name=degree, unit_auth_name=EPSG, unit_code=9102, unit_category=angular)
>>> crs.coordinate_operation.params[0].name
'Latitude of false origin'
>>> crs.coordinate_operation.params[0].value
33.75

获取中央子午线，PyProj？

1 个答案: