xqj z d Z ddlZddlZddlZddlZddlmZ ddlmZm Z ddl
m Z m Z m
Z
mZmZmZmZmZmZmZmZmZ ddlmZmZmZmZm Z m Z ddl m Z ddl m!Z!m"Z"m#Z# dd l$m%Z%m&Z& dd
l'm(Z( dd l)m*Z* ejV d ejX Z- G d
dej\ Z/d-de0de1de2fdZ3de2de2fdZ4de2de2fdZ5de2de2e6z fdZ7de2e6z fdZ8dede1fdZ9 G d d Z: G d d e: Z; G d d e; Z< G d! d"e; Z= G d# d$e; Z> G d% d&e; Z? G d' d(e; Z@ G d) d*e; ZA G d+ d,e; ZBy).zx
This module interfaces with PROJ to produce a pythonic interface
to the coordinate reference system (CRS) information.
N)Callable)AnyOptional) _CRS AreaOfUseAuthorityMatchInfoAxisCoordinateOperationCoordinateSystemDatum Ellipsoid
PrimeMeridian_load_proj_jsonis_projis_wkt)!_GEOGRAPHIC_GRID_MAPPING_NAME_MAP_GRID_MAPPING_NAME_MAP)_INVERSE_GEOGRAPHIC_GRID_MAPPING_NAME_MAP _INVERSE_GRID_MAPPING_NAME_MAP _horizontal_datum_from_params_try_list_if_string)ToWGS84Transformation)
Cartesian2DCSEllipsoidal2DCS
VerticalCS) ProjVersion
WktVersion)CRSError)Geoda9
\+ # parameter starts with '+' character
(?P\w+) # capture parameter name
\=? # match both key only and key-value parameters
(?P\S+)? # capture all characters up to next space (None if no value)
\s*? # consume remaining whitespace, if any
c " e Zd ZdZ fdZ xZS )CRSLocalz
Threading local instance for cython CRS class.
For more details, see:
https://github.com/pyproj4/pyproj/issues/782
c 0 d | _ t | y N)crssuper__init__)self __class__s J/var/www/html/truck-me/venv/lib/python3.12/site-packages/pyproj/crs/crs.pyr& zCRSLocal.__init__> s
)__name__
__module__ __qualname____doc__r&
__classcell__r( s @r) r! r! 6 s r* r!
projparams
allow_jsonreturnc t | t s t d d| vrd| vr|rt j | S g }| j D ] \ }}t |t t f r(dj |D cg c]
}t | c} }|t | dk( r|j d| it | dk( rx|j d| d| t d j | S c c}w )
NzCRS input is not a dictprojinit,True+False= )
isinstancedictr jsondumpsitemslisttuplejoinstrappend_prepare_from_string)r1 r2 pjargskeyvaluevals r) _prepare_from_dictrL C s j$ '011
Z F*$<zz*%%
F &&(
-
U edE] +HH% 83c#h 89E =CJ& 0 MMAcU) $
Z7
" MMAcU!E7+ ,
- 0 11 9s 8C8
in_crs_stringc t j dd| j } d}| j | scg }d}| j D ]; }|s&|j | r|j d| d}+|j
| = dj | } d| vrd | v r| d
z
} n| d z
} | j d d
j } | j d r t j dt d | S )Nz
[\s+]?=[\s+]?r; )+initz+projr6 r5 Fr Tr< ztype=crsr9 z
+type=crsz type=crsz
+init=EPSGz
+init=epsg)rO r6 z'+init=:' syntax is deprecated. ':' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6 )
stacklevel)
resublstrip
startswithsplitinsertrF rD replacestripwarningswarn
FutureWarning)rM starting_paramskvpairsfirst_item_insertedkvpairs r) _prepare_from_proj_stringra Y s FF+S-2F2F2HIM8O # #O 4 ##))+ 'F&F,=,=o,Nq&)&*#v& ' )
& - \ )M [ (M ")), EKKMM 12
,
r* c . t | t s t d | st d| d| v r2 t j | d }|s t d t
| } | S t | r t | } | S # t
$ r}t d |d }~ww xY w)NzCRS input is not a stringz CRS string is empty or invalid: {F)strictz'CRS appears to be JSON but is not validzCRS is empty JSON) r= rE r r? loads
ValueErrorrL r ra )rM crs_dicterrs r) rG rG s mS )233 9-9JKLL
m Ozz->H ./ /*84
1-@
ODE3 N Os A: : B BB auth_name auth_codec | d| S )N: )ri rj s r) _prepare_from_authorityrn s [) %%r* c t d| S )NEPSG)rn rj s r) _prepare_from_epsgrr s "69 55r* c t | t ryt | t r| j ryt | d r| j
dk( ryy)NTshaperm F)r= intrE isnumerichasattrrt rq s r)
_is_epsg_coderx sB )S!)S!i&9&9&;y'"y"'< r* c e Zd ZdZdZdedz ddfdZed Zede de e
z dd fd Z ede e
z dd fd
Z ed e dd fd Z
ed
e dd fd Zede dd fd Zde fdZededd fd Zdedz fdZededd fd Zede dd fd Zededd fd ZdefdZej4 d fd ee z d edefd Z e d[d ed!edz d"edz d#edz dd f
d$ Z de e fd%Z d&edefd'Z!d\d&ed(edefd)Z"ede#d fd* Z$ede#d fd+ Z%ede#d fd, Z&ede d fd- Z'ede dz fd. Z(ede fd/ Z)ede fd0 Z*ede e+ fd1 Z,ede-dz fd2 Z.ede/dz fd3 Z0ede1dz fd4 Z2ede3dz fd5 Z4ede5dz fd6 Z6ede7dz fd7 Z8ede fd8 Z9ede fd9 Z:ej4 d dfd:ee z d;ede
de fd?Zj~ fd:e>e
z de fdAZ@d^dBe
de
dz fdCZAd_de dz dBe
fdDZB d_de dz dBe
de eC fdEZDdZdFe dz dd fdGZEdZdFe dz dd fdHZFedefdI ZGedefdJ ZHedefdK ZIedefdL ZJedefdM ZKedefdN ZLedefdO ZMedP ZNedefdQ ZOde d fdRZPd&eQdefdSZRdee e f fdTZSdUee ef fdVZTde
fdWZUde fdXZVde fdYZWy)`CRSa/
A pythonic Coordinate Reference System manager.
.. versionadded:: 2.0.0
See: :c:func:`proj_create`
The functionality is based on other fantastic projects:
* `rasterio `_
* `opendatacube `_
Attributes
----------
srs: str
The string form of the user input used to create the CRS.
Nr1 r3 c d}|rt |t r
|j }nt | r t | }nt |t
r t
| }ntt |t r t | }nXt |t t f rt | dk( r t | }n+t |d r|j }nt d| |r)t |d }t
dj! ||f }|| _ t# | _ t |t r|| j$ _ yt | j | j$ _ y) u
Initialize a CRS class instance with:
- PROJ string
- Dictionary of PROJ parameters
- PROJ keyword arguments for parameters
- JSON string with PROJ parameters
- CRS WKT string
- An authority string [i.e. 'epsg:4326']
- An EPSG integer code [i.e. 4326]
- A tuple of ("auth_name": "auth_code") [i.e ('epsg', '4326')]
- An object with a `to_wkt` method.
- A :class:`pyproj.crs.CRS` class
Example usage:
>>> from pyproj import CRS
>>> crs_utm = CRS.from_user_input(26915)
>>> crs_utm
Name: NAD83 / UTM zone 15N
Axis Info [cartesian]:
- E[east]: Easting (metre)
- N[north]: Northing (metre)
Area of Use:
- name: North America - 96°W to 90°W and NAD83 by country
- bounds: (-96.0, 25.61, -90.0, 84.0)
Coordinate Operation:
- name: UTM zone 15N
- method: Transverse Mercator
Datum: North American Datum 1983
- Ellipsoid: GRS 1980
- Prime Meridian: Greenwich
>>> crs_utm.area_of_use.bounds
(-96.0, 25.61, -90.0, 84.0)
>>> crs_utm.ellipsoid
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1],
ID["EPSG",7019]]
>>> crs_utm.ellipsoid.inverse_flattening
298.257222101
>>> crs_utm.ellipsoid.semi_major_metre
6378137.0
>>> crs_utm.ellipsoid.semi_minor_metre
6356752.314140356
>>> crs_utm.prime_meridian
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8901]]
>>> crs_utm.prime_meridian.unit_name
'degree'
>>> crs_utm.prime_meridian.unit_conversion_factor
0.017453292519943295
>>> crs_utm.prime_meridian.longitude
0.0
>>> crs_utm.datum
DATUM["North American Datum 1983",
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]],
ID["EPSG",6269]]
>>> crs_utm.coordinate_system
CS[Cartesian,2],
AXIS["(E)",east,
ORDER[1],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]],
AXIS["(N)",north,
ORDER[2],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]
>>> print(crs_utm.coordinate_operation.to_wkt(pretty=True))
CONVERSION["UTM zone 15N",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",-93,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",0.9996,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",500000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",0,
LENGTHUNIT["metre",1],
ID["EPSG",8807]],
ID["EPSG",16015]]
>>> crs = CRS(proj='utm', zone=10, ellps='WGS84')
>>> print(crs.to_wkt(pretty=True))
PROJCRS["unknown",
BASEGEOGCRS["unknown",
DATUM["Unknown based on WGS84 ellipsoid",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1],
ID["EPSG",7030]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8901]]],
CONVERSION["UTM zone 10N",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",-123,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",0.9996,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",500000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",0,
LENGTHUNIT["metre",1],
ID["EPSG",8807]],
ID["EPSG",16010]],
CS[Cartesian,2],
AXIS["(E)",east,
ORDER[1],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]],
AXIS["(N)",north,
ORDER[2],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]]
>>> geod = crs.get_geod()
>>> f"+a={geod.a:.0f} +f={geod.f:.8f}"
'+a=6378137 +f=0.00335281'
>>> crs.is_projected
True
>>> crs.is_geographic
False
rP to_wktzInvalid CRS input: F)r2 r< N)r= r srsrx rr rE rG r> rL rB rC lenrn rw r} r rD r! _localr$ )r' r1 kwargs
projstring
projkwargss r) r& z CRS.__init__ s T
*d+ '^^
z* /
;
J, 1* =
J- /
;
Ju
63z?a;O 4j A
X. ' . . 0
!4ZN CDD +FuEJ-chh
J7O.PQJ j j$ ' (DKKO "488nDKKOr* c | j j $t | j | j _ | j j S )zH
Retrieve the Cython based _CRS object for this thread.
)r r$ r r~ r' s r) _crszCRS._crs` s3
;;?? " "488nDKKO{{ r* ri codec 8 | j t || S )a?
.. versionadded:: 2.2.0
Make a CRS from an authority name and authority code
Parameters
----------
auth_name: str
The name of the authority.
code : int or str
The code used by the authority.
Returns
-------
CRS
)from_user_inputrn )clsri r s r) from_authorityzCRS.from_authorityi s $ ""#:9d#KLLr* c 6 | j t | S )zMake a CRS from an EPSG code
Parameters
----------
code : int or str
An EPSG code.
Returns
-------
CRS
)r rr )r r s r) from_epsgz
CRS.from_epsg} s ""#5d#;<2BCD D""#<^#LMMr*
in_wkt_stringc h t | st d| | j t | S )z
.. versionadded:: 2.2.0
Make a CRS from a WKT string
Parameters
----------
in_wkt_string : str
A WKT string.
Returns
-------
CRS
zInvalid WKT string: )r r r rG )r r s r) from_wktz CRS.from_wkt s6 m$ 1- AB B""#7
#FGGr* rM c 6 | j t | S )az
Make a CRS from:
Initialize a CRS class instance with:
- PROJ string
- JSON string with PROJ parameters
- CRS WKT string
- An authority string [i.e. 'epsg:4326']
Parameters
----------
in_crs_string : str
An EPSG, PROJ, or WKT string.
Returns
-------
CRS
)r rG )r rM s r) from_stringzCRS.from_string s ( ""#7
#FGGr* c d | j d }|rdj | S | j S )a
.. versionadded:: 2.2.0
Convert the CRS to a string.
It attempts to convert it to the authority string.
Otherwise, it uses the string format of the user
input to create the CRS.
Returns
-------
str
d min_confidencerl ) to_authorityrD r~ )r' auth_infos r) to_stringz
CRS.to_string s2 %%S%9 88I& &xxr* rJ c 0 t || r|S | |fi |S a
Initialize a CRS class instance with:
- PROJ string
- Dictionary of PROJ parameters
- PROJ keyword arguments for parameters
- JSON string with PROJ parameters
- CRS WKT string
- An authority string [i.e. 'epsg:4326']
- An EPSG integer code [i.e. 4326]
- A tuple of ("auth_name": "auth_code") [i.e ('epsg', '4326')]
- An object with a `to_wkt` method.
- A :class:`pyproj.crs.CRS` class
Parameters
----------
value : obj
A Python int, dict, or str.
Returns
-------
CRS
)r= )r rJ r s r) r zCRS.from_user_input s# 0 eS !L5#F##r* c | j yt | j j | j j | j j S )zs
Returns
-------
pyproj.geod.Geod:
Geod object based on the ellipsoid.
N)arfb) ellipsoidr semi_major_metreinverse_flatteningsemi_minor_metrer s r) get_geodz CRS.get_geod sE >> !nn--~~00nn--
r* proj_dictc 6 | j t | S )z
.. versionadded:: 2.2.0
Make a CRS from a dictionary of PROJ parameters.
Parameters
----------
proj_dict : str
PROJ params in dict format.
Returns
-------
CRS
)r rL )r r s r) from_dictz
CRS.from_dict s ""#5i#@AAr* crs_jsonc 6 | j t | S )z
.. versionadded:: 2.4.0
Create CRS from a CRS JSON string.
Parameters
----------
crs_json: str
CRS JSON string.
Returns
-------
CRS
)r r )r r s r) from_jsonz
CRS.from_json s ""?8#<==r* rg c J | j t j | S )z
.. versionadded:: 2.4.0
Create CRS from a JSON dictionary.
Parameters
----------
crs_dict: dict
CRS dictionary.
Returns
-------
CRS
)r r? r@ )r rg s r) from_json_dictzCRS.from_json_dict, s ""4::h#788r* c | j }|i S d }i }t j | D ]) }|j \ }}| || }|dus%|||< + |S )a
.. versionadded:: 2.2.0
Converts the CRS to dictionary of PROJ parameters.
.. warning:: You will likely lose important projection
information when converting to a PROJ string from
another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
Returns
-------
dict:
PROJ params in dict format.
c | j dk( ry| j dk( ry t | S # t $ r Y nw xY w t | S # t $ r
Y t | S w xY w)NtrueTfalseF)lowerru rf floatr )rK s r) _parsezCRS.to_dict.._parseS sp yy{f$yy{g%
3x
Sz!
&s+ +
s!
5 A A
A A&%A&F)to_proj4_RE_PROJ_PARAMfinditergroups)r' proj_stringr r paramrI rJ s r) to_dictz CRS.to_dict> sw " mmo I
, #,,[9 'E JC u
E!!& # ' r* F wkt_versionerrcheckc d| j | i}| j r| j r|| j j rf| j rZ| j j || }|j
d |j | | j j |d< |S | j rH| j D ]7 }|j || }|j
d |j | 9 |S | j rF| j }d|v r |d d |d< | j r| j j |d< |S | j ri|j | j j | j j | j j" | j j |d < | j$ r2| j$ j& |d
< | j$ j |d < | j( rR| j( j |d < | j( j r#| j( j j |d
< | j* r| j r| j j, j/ t0 vr5|r,t3 j4 d| j j, d|d iS |j t1 | j j, j/ | j |S d|d< |S d}| j s'| j6 r| j }| j |d< |sdn)|j, j/ j9 dd }|t: vrC|r:|r#t3 j4 d|j, nt3 j4 d d|d iS |j t; | | |S )aS
.. versionadded:: 2.2.0
This converts a :obj:`pyproj.crs.CRS` object
to a Climate and Forecast (CF) Grid Mapping Version 1.8 dict.
:ref:`build_crs_cf`
Parameters
----------
wkt_version: str or pyproj.enums.WktVersion
Version of WKT supported by CRS.to_wkt.
Default is :attr:`pyproj.enums.WktVersion.WKT2_2019`.
errcheck: bool, default=False
If True, will warn when parameters are ignored.
Returns
-------
dict:
CF-1.8 version of the projection.
crs_wkt)r r towgs84 geoid_modelname
geoid_namegeopotential_datum_name)semi_major_axissemi_minor_axisr reference_ellipsoid_namelongitude_of_prime_meridianprime_meridian_namegeographic_crs_namehorizontal_datum_namez"Unsupported coordinate operation: latitude_longitudegrid_mapping_nameNprojected_crs_namer< _z Coordinate operation not found.) r} is_boundcoordinate_operationr
source_crsto_cfpopupdate is_compound sub_crs_list is_vertical to_json_dictdatumr r r r r prime_meridian longitude geodetic_crs
is_geographic method_namer r rZ r[ is_projectedrX r ) r' r r cf_dictsub_cfsub_crs vert_jsonr coordinate_operation_names r) r z CRS.to_cfl s 8 $-dkk+.F"G
MM))))11%)__%:%: '! &; &F
JJy ! NN6 "!%!:!:!B!BGI N ,,
' ; R
9%v&
' N ))+I )(1-(@(H %zz59ZZ__12N >> NN $ ? ? $ ? ?#'>>#D#D
37..2E2EG. / 595H5H5R5RG1 2-1-@-@-E-EG) * -1->->-C-CG) * &&373D3D3J3J3O3O/0 ((--99??ADE
@ # 8 8 D D E G &wy'9::= 11==CCE//1 N 0D+,N $ }}!2!2#'#<#< ,0IIG( ) (
%11779AA#sK "
%,J J' MM</;;<>
MM"CD wy 12 2 *+D E$
r* in_cfellipsoidal_cs cartesian_cs vertical_csc d}d| v rt | d S d| v rt | d S | j d }| t d t | } t | }| j d }|rt
|xs d|| } na|rS||vrOt | } |N| j }
t j | j |
d
< t |
} n t
| } |d k( r| S |t | || |
S t | } t | j dd | | | | } d}
d| v r.t | dt | j gt! | d }
d| vr|
xs | S t# d| d | j d | }t% d|
xs | |g S # t
$ r d}Y Sw xY w# t
$ r t d| dw xY w)a
.. versionadded:: 2.2.0
.. versionadded:: 3.0.0 ellipsoidal_cs, cartesian_cs, vertical_cs
This converts a Climate and Forecast (CF) Grid Mapping Version 1.8
dict to a :obj:`pyproj.crs.CRS` object.
:ref:`build_crs_cf`
Parameters
----------
in_cf: dict
CF version of the projection.
ellipsoidal_cs: Any, optional
Input to create an Ellipsoidal Coordinate System.
Anything accepted by :meth:`pyproj.crs.CoordinateSystem.from_user_input`
or an Ellipsoidal Coordinate System created from :ref:`coordinate_system`.
cartesian_cs: Any, optional
Input to create a Cartesian Coordinate System.
Anything accepted by :meth:`pyproj.crs.CoordinateSystem.from_user_input`
or :class:`pyproj.crs.coordinate_system.Cartesian2DCS`.
vertical_cs: Any, optional
Input to create a Vertical Coordinate System accepted by
:meth:`pyproj.crs.CoordinateSystem.from_user_input`
or :class:`pyproj.crs.coordinate_system.VerticalCS`
Returns
-------
CRS
)unknown undefinedr spatial_refr Nz4CF projection parameters missing 'grid_mapping_name'r r )r r r coordinate_system)r r )base_crs
conversionr z Unsupported grid mapping name: r )r r r r r zWGS 84)r
target_crstransformationr r )r r r r )r
components)rz getr r r KeyError
GeographicCRSr r r DerivedGeographicCRSr ProjectedCRSBoundCRSr r r VerticalCRS CompoundCRS)r r r r
unknown_namesr r geographic_conversion_methodr geographic_crsgeographic_crs_jsonconversion_method
projected_crs bound_crs vertical_crss r) from_cfz CRS.from_cf sY P 1
uY'( ( E !u]+, , !II&9 : $QR R .e4 012CD
) $ii(= > "/(7K -#N
!%8
%M !4 5N )&4&A&A&C#$44^DQQS $$78 "%%8!9 *. IN 4 4! ! ' 3''7> -
67H I
%/=(/'%
(# 4!.. 1DU9EU1V I %E 1 -
- # 12 ,/#
)*D}l)S
C 0+/ ( 0@ 12C1DE
s F ? F1 F.-F.1G
c g }d }| j dk( r| j r| j j }|j dd j j d }t
| dk( r| d}| j d }d d |j d
d j j d d
dd}|j j j d
d j }|r
| d| |d< |j | | j r3|j | j j || |S | j ri| j r]| j j rG|j | j j j || j |S | j" D ]! }|j |j% # |S )a1
.. versionadded:: 3.0.0
This converts all coordinate systems (cs) in the CRS
to a list of Climate and Forecast (CF) Version 1.8 dicts.
:ref:`build_crs_cf`
Returns
-------
list[dict]:
CF-1.8 version of the coordinate systems.
c | j xr* | j j j t v S # t $ r Y yw xY w)NF)r r r r r )r$ s r) rotated_polez"CRS.cs_to_cf..rotated_pole~ sP
,, A00<<BBD@A
s 7: AAz Temporal CRS time_originz
1875-05-20zZ z-01-01r timecalendarproleptic_gregorianr< r T)
standard_name long_namer axisr| z since units)r ) type_namer r r rY rstripr axis_infor rX unit_namerF r extendr r r r cs_to_cf) r' cf_axis_listr
datum_jsonoriginr cf_temporal_axisr r s r) r z CRS.cs_to_cfn s >>^ +
002J ^^M<@FFHOOPTUF6{a "86*>>!$D!' # NN:/DEUW WS#& ,,.66z2FLLNI/8k,I ) 0 1 ! ! &&,,,t:L,M
]]t4??3T3T 1177!-doo!> 8
,,
8 ##G$4$4$67
8r* otherc t j | }| j j |j S # t $ r Y yw xY w)ab
Check if the CRS objects are the exact same.
Parameters
----------
other: Any
Check if the other CRS is the exact same to this object.
If the other object is not a CRS, it will try to create one.
On Failure, it will return False.
Returns
-------
bool
F)rz r r r
is_exact_samer' r s r) r! zCRS.is_exact_same sF ''.E yy&&uzz22 s < AAignore_axis_orderc t j | }| j j |j | S # t $ r Y yw xY w)a
.. versionadded:: 2.5.0
Check if the CRS objects are equivalent.
Parameters
----------
other: Any
Check if the other object is equivalent to this object.
If the other object is not a CRS, it will try to create one.
On Failure, it will return False.
ignore_axis_order: bool, default=False
If True, it will compare the CRS class and ignore the axis order.
Returns
-------
bool
F)r# )rz r r r equals)r' r r# s r) r% z
CRS.equals sL ( ''.E yy
>O PP s > A
A
c | | j j dS | j | j j S z
.. versionadded:: 2.2.0
Returns
-------
CRS:
The geodeticCRS / geographicCRS from the CRS.
N)r r r( r s r) r zCRS.geodetic_crs s; yy%%-
6 67
r* c | | j j dS | j | j j S z
The base CRS of a BoundCRS or a DerivedCRS/ProjectedCRS,
or the source CRS of a CoordinateOperation.
Returns
-------
CRS
N)r r r( r s r) r zCRS.source_crs s; yy##+
4 45
r* c | | j j dS | j | j j S z
.. versionadded:: 2.2.0
Returns
-------
CRS:
The hub CRS of a BoundCRS or the target CRS of a CoordinateOperation.
N)r r r( r s r) r zCRS.target_crs s; yy##+
4 45
r* c r | j j D cg c] }| j | c}S c c}w z
If the CRS is a compound CRS, it will return a list of sub CRS objects.
Returns
-------
list[CRS]
)r r r( r' r s r) r zCRS.sub_crs_list
s+ 8>> from pyproj import CRS
>>> ccs = CRS("EPSG:4328")
>>> ccs.to_epsg()
4328
If the CRS is bound, you can attempt to get an epsg code from
the source CRS:
>>> from pyproj import CRS
>>> ccs = CRS("+proj=geocent +datum=WGS84 +towgs84=0,0,0")
>>> ccs.to_epsg()
>>> ccs.source_crs.to_epsg()
4978
>>> ccs == CRS.from_epsg(4978)
False
Parameters
----------
min_confidence: int, default=70
A value between 0-100 where 100 is the most confident.
:ref:`min_confidence`
Returns
-------
int | None:
The best matching EPSG code matching the confidence level.
r )r to_epsg)r' r s r) rN z CRS.to_epsg s F yy ??r* c < | j j || S )a!
.. versionadded:: 2.2.0
Return the authority name and code best matching the CRS
or None if it a match is not found.
Example:
>>> from pyproj import CRS
>>> ccs = CRS("EPSG:4328")
>>> ccs.to_authority()
('EPSG', '4328')
If the CRS is bound, you can get an authority from
the source CRS:
>>> from pyproj import CRS
>>> ccs = CRS("+proj=geocent +datum=WGS84 +towgs84=0,0,0")
>>> ccs.to_authority()
>>> ccs.source_crs.to_authority()
('EPSG', '4978')
>>> ccs == CRS.from_authorty('EPSG', '4978')
False
Parameters
----------
auth_name: str, optional
The name of the authority to filter by.
min_confidence: int, default=70
A value between 0-100 where 100 is the most confident.
:ref:`min_confidence`
Returns
-------
tuple(str, str) or None:
The best matching (, ) for the confidence level.
ri r )r r r' ri r s r) r zCRS.to_authority9 s& L yy%% &
r* c < | j j || S )aT
.. versionadded:: 3.2.0
Return the authority names and codes best matching the CRS.
Example:
>>> from pyproj import CRS
>>> ccs = CRS("EPSG:4328")
>>> ccs.list_authority()
[AuthorityMatchInfo(auth_name='EPSG', code='4326', confidence=100)]
If the CRS is bound, you can get an authority from
the source CRS:
>>> from pyproj import CRS
>>> ccs = CRS("+proj=geocent +datum=WGS84 +towgs84=0,0,0")
>>> ccs.list_authority()
[]
>>> ccs.source_crs.list_authority()
[AuthorityMatchInfo(auth_name='EPSG', code='4978', confidence=70)]
>>> ccs == CRS.from_authorty('EPSG', '4978')
False
Parameters
----------
auth_name: str, optional
The name of the authority to filter by.
min_confidence: int, default=70
A value between 0-100 where 100 is the most confident.
:ref:`min_confidence`
Returns
-------
list[AuthorityMatchInfo]:
List of authority matches for the CRS.
rP )r list_authorityrQ s r) rS zCRS.list_authorityc s& P yy'' (
r* r c X | j | j j | S a
.. versionadded:: 3.1.0
Convert the current CRS to the 3D version if it makes sense.
New vertical axis attributes:
- ellipsoidal height
- oriented upwards
- metre units
Parameters
----------
name: str, optional
CRS name. Defaults to use the name of the original CRS.
Returns
-------
CRS
r )r( r to_3dr' r s r) rW z CRS.to_3d s" ( ~~diioo4o 899r* c X | j | j j | S )a#
.. versionadded:: 3.6.0
Convert the current CRS to the 2D version if it makes sense.
Parameters
----------
name: str, optional
CRS name. Defaults to use the name of the original CRS.
Returns
-------
CRS
rV )r( r to_2drX s r) rZ z CRS.to_2d s" ~~diioo4o 899r* c . | j j S )aW
This checks if the CRS is geographic.
It will check if it has a geographic CRS
in the sub CRS if it is a compound CRS and will check if
the source CRS is geographic if it is a bound CRS.
Returns
-------
bool:
True if the CRS is in geographic (lon/lat) coordinates.
)r r r s r) r zCRS.is_geographic s yy&&&r* c . | j j S )a6
This checks if the CRS is projected.
It will check if it has a projected CRS
in the sub CRS if it is a compound CRS and will check if
the source CRS is projected if it is a bound CRS.
Returns
-------
bool:
True if CRS is projected.
)r r r s r) r zCRS.is_projected s yy%%%r* c . | j j S )aS
.. versionadded:: 2.2.0
This checks if the CRS is vertical.
It will check if it has a vertical CRS
in the sub CRS if it is a compound CRS and will check if
the source CRS is vertical if it is a bound CRS.
Returns
-------
bool:
True if CRS is vertical.
)r r r s r) r zCRS.is_vertical s yy$$$r* c . | j j S )zY
Returns
-------
bool:
True if CRS is bound.
)r r r s r) r z CRS.is_bound s yy!!!r* c . | j j S )z}
.. versionadded:: 3.1.0
Returns
-------
bool:
True if CRS is compound.
)r r r s r) r zCRS.is_compound s yy$$$r* c . | j j S )z
.. versionadded:: 2.2.0
Returns
-------
bool:
True if CRS is local/engineering.
)r is_engineeringr s r) ra zCRS.is_engineering r: r* c . | j j S )z
This checks if the CRS is geocentric and
takes into account if the CRS is bound.
Returns
-------
bool:
True if CRS is in geocentric (x/y) coordinates.
)r
is_geocentricr s r) rc zCRS.is_geocentric yy&&&r* c . | j j S )z
.. versionadded:: 3.2.0
Returns
-------
bool:
True if CRS is a Derived CRS.
)r
is_derivedr s r) rf zCRS.is_derived s yy###r* c . | j j S )z
.. versionadded:: 3.7.0
Check if the CRS is deprecated
Returns
-------
bool
)r
is_deprecatedr s r) rh zCRS.is_deprecated rd r* c 6 | j j S )z
.. versionadded:: 3.7.0
Return a list of non-deprecated objects related to this.
Returns
-------
list[CRS]
)r get_non_deprecatedr s r) rj zCRS.get_non_deprecated- s yy++--r* c $ | j | S r# )r% r" s r) __eq__z
CRS.__eq__9 s {{5!!r* c d| j iS )Nr~ r~ r s r) __getstate__zCRS.__getstate__< s txx r* statec X | j j | t | _ y r# )__dict__r r! r )r' rp s r) __setstate__zCRS.__setstate__? s
U# j r* c 4 t | j S r# )hashr} r s r) __hash__z CRS.__hash__C s DKKM""r* c | j S r# rn r s r) __str__z CRS.__str__F s xxr* c g }| j D ] }|j dt | dg ! dj | }d}d}| j r,| j j
rt | j }n| j rF| j r:t | j j }d| j j d}nzg }dg}| j D ]D } |j t | j |j d| j dg F dj | }dj | }d}
| j r>dj ddt | j d | j j dg }
| j } t | d
k r| nd j | d d
d g } |xs d
}d| j d| d| j d|xs d d| d| j xs d d|
d| j" d| j$ xs d d| j& xs d d| | S )Nz-
r| z Source CRS: z Sub CRS:
|zCoordinate Operation:
z- name: z
- method: 2 r< z...z - undefined
Name: z
Axis Info [r z]:
z
Area of Use:
z - undefinedzDatum: z
- Ellipsoid: z
- Prime Meridian: )r r rE rD r axis_listr r r r rF r r r r r r7 r r r ) r' axis_info_listr
axis_info_strsource_crs_repr sub_crs_reprcoordinate_system_namecoordinate_system_namessub_crs_repr_listr r srs_reprs r) __repr__z CRS.__repr__I sT $&NN ;D ! !4TD"9 : ;/
! !d&<&<&F&F%()?)?%@ "
]]t%( )J)J%K " ,T__-A-A,B" EO&( #!- ,,
E'..s73L3L/MN!(($ d)CD
E &)XX.E%F " 77#45L " $ $#%77- 112"--99
$ >># "8} 28(3B-QVAW8X%8
r( ,YYK 0?K@ o 0=1#$djj\ " NN 9k: ;!!%!4!4!C DB n
r* r# )NNN)F)FrP )F )Nr )Xr+ r, r- r. r r& propertyr classmethodrE ru r r r r r r r r r r> r r r r r WKT2_2019boolr staticmethodr rB r r! r% r r r r r r1 r r r r r r7 r
r r r r r r r r
r r? rA r} rI r r PROJ_5r rN r r rS rW rZ r r r r r ra rc rf rh rj objectrl ro rs rv rx r rm r* r) rz rz s &e-3: e-4 e-N Ms M#) M M M& =S3Y =5 = = N N N N& HS HU H H& H H H H*3 & $C $e $ $6
$+
B$ B5 B B" > > > >" 9d 9u 9 9", ,` )3(<(< A #%A A
AF &*#'"& |
|
d
|
Dj|
4Z |
|
|
|< $t* < |33 34 3*QC QD QT Q4
huo
HUO
HUO
Od5k O O #* 0 c #3 # # #4: # # %Y - % % #9t+ # # (
4 ( ( ut| +#3d#: + + .&9D&@ . . ! ! ! s %/$8$8 (, *c!* * + *
*Xd S ,
(d
( 5@4F4F c 1 3 0#@c #@3: #@J(
cDj (
(
V CE*
t*
>!5!5 5 /{4;O;O:PPQR
6r* rJ rz c t || r|S | j | }t t | |fi | |j |S r )r= __new__r% r r& r )r rJ r r$ r( s r) r z$CustomConstructorCRS.from_user_input sF 0 eS !Lkk#
"C1%B6B
r* c p | j j dS t | j j S r' )r r rz r s r) r z!CustomConstructorCRS.geodetic_crs s- yy--5tV3tyy?U?U;VVr* c p | j j dS t | j j S r) )r r rz r s r) r z CustomConstructorCRS.source_crs s- yy++3tRTYY=Q=Q9RRr* c p | j j dS t | j j S r+ )r r rz r s r) r z CustomConstructorCRS.target_crs s- yy++3tRTYY=Q=Q9RRr* c f | j j D cg c]
}t | c}S c c}w r- )r r rz r. s r) r z!CustomConstructorCRS.sub_crs_list s& -1II,B,BCG CCCs .Nr c L t | j j | S rU )rz r rW rX s r) rW zCustomConstructorCRS.to_3d s ( 499???-..r* r# )r+ r, r- r. r rC rE r r r r r r r r r rB r rW r/ r0 s @r) r r s "sCx " " C e <
Whuo
W
W SHUO S S
SHUO
S
S Dd5k D D/#* / /r* r c F e Zd ZdZdZ d dedededz ddf fd
Z xZS )
r zR
.. versionadded:: 2.5.0
This class is for building a Geographic CRS
)zGeographic CRSzGeographic 2D CRSzGeographic 3D CRSNr r r r3 c t j | j }dd|t j |xs
t j d}|d dk( r||d< n||d< t
| | y) aC
Parameters
----------
name: str, default="undefined"
Name of the CRS.
datum: Any, default="urn:ogc:def:ensemble:EPSG::6326"
Anything accepted by :meth:`pyproj.crs.Datum.from_user_input` or
a :class:`pyproj.crs.datum.CustomDatum`.
ellipsoidal_cs: Any, optional
Input to create an Ellipsoidal Coordinate System.
Anything accepted by :meth:`pyproj.crs.CoordinateSystem.from_user_input`
or an Ellipsoidal Coordinate System created from :ref:`coordinate_system`.
2https://proj.org/schemas/v0.2/projjson.schema.jsonr )$schematyper r r
DatumEnsembledatum_ensembler N)r r r r r r% r& )r' r r r r r( s r) r& zGeographicCRS.__init__ s & %%e,99;K#!1!A!A 3/"3"ln
=O +49 0 1+0 (
,-r* )r z urn:ogc:def:ensemble:EPSG::6326N r+ r, r- r. r rE r r& r/ r0 s @r) r r sO SO 6%) . . . d
.
. .r* r c H e Zd ZdZdZ d
dedededz deddf
fd
Z xZS ) r zZ
.. versionadded:: 2.5.0
This class is for building a Derived Geographic CRS
)zDerived Geographic CRSzDerived Geographic 2D CRSzDerived Geographic 3D CRSNr r r r r3 c dd|t j | j t j | j t j |xs
t j d}t | | y)a
Parameters
----------
base_crs: Any
Input to create the Geodetic CRS, a :class:`GeographicCRS` or
anything accepted by :meth:`pyproj.crs.CRS.from_user_input`.
conversion: Any
Anything accepted by :meth:`pyproj.crs.CoordinateSystem.from_user_input`
or a conversion from :ref:`coordinate_operation`.
ellipsoidal_cs: Any, optional
Input to create an Ellipsoidal Coordinate System.
Anything accepted by :meth:`pyproj.crs.CoordinateSystem.from_user_input`
or an Ellipsoidal Coordinate System created from :ref:`coordinate_system`.
name: str, default="undefined"
Name of the CRS.
r r r r r r r r N)rz r r r
r r r% r& )r' r r r r derived_geographic_crs_jsonr( s r) r& z DerivedGeographicCRS.__init__? sv 0 L*++H5BBD-==ln!1!A!A 3/"3"ln '
# 45r* )Nr r+ r, r- r. r r rE r& r/ r0 s @r) r r 2 sU O &* #6#6 #6 d
#6
#6
#6 #6r* r c : e Zd ZdZdZ ddededdf fd
Z xZS )
GeocentricCRSzR
.. versionadded:: 3.2.0
This class is for building a Geocentric CRS
)zGeocentric CRSr r r3 Nc dd|t j | j ddddddd d
d ddd d
dddgdd}t | | y)a0
Parameters
----------
name: str, default="undefined"
Name of the CRS.
datum: Any, default="urn:ogc:def:datum:EPSG::6326"
Anything accepted by :meth:`pyproj.crs.Datum.from_user_input` or
a :class:`pyproj.crs.datum.CustomDatum`.
r GeodeticCRS Cartesianz Geocentric XX geocentricXmetre)r abbreviation directionunitz Geocentric YY geocentricYz Geocentric ZZ geocentricZ)subtyper r r r r r N)r r r r% r& )r' r r geocentric_crs_jsonr( s r) r& zGeocentricCRS.__init__n s M!**51>>@& !/(+%2 ' !/(+%2 ' !/(+%2 ' "
: ,-r* )r z urn:ogc:def:datum:EPSG::6326r r0 s @r) r r e s= *O 3+.+. +.
+. +.r* r c P e Zd ZdZdZ d
dedededz dedz ddf
fd
Z xZS ) r zR
.. versionadded:: 2.5.0
This class is for building a Projected CRS.
)z
Projected CRSzDerived Projected CRSNr r r r r3 c , dd|t j |xs
t j t j | j t j |xs
t
j d}t | ! | y)a
Parameters
----------
conversion: Any
Anything accepted by :meth:`pyproj.crs.CoordinateSystem.from_user_input`
or a conversion from :ref:`coordinate_operation`.
name: str, optional
The name of the Projected CRS. Default is undefined.
cartesian_cs: Any, optional
Input to create a Cartesian Coordinate System.
Anything accepted by :meth:`pyproj.crs.CoordinateSystem.from_user_input`
or :class:`pyproj.crs.coordinate_system.Cartesian2DCS`.
geodetic_crs: Any, optional
Input to create the Geodetic CRS, a :class:`GeographicCRS` or
anything accepted by :meth:`pyproj.crs.CRS.from_user_input`.
r r r N) rz r r r r
r r r% r& )r' r r r r
proj_crs_jsonr( s r) r& zProjectedCRS.__init__ s{ 0 L"++ /
ln-==ln!1!A!A /
"ln
'r* )r NNr r0 s @r) r r s] AO
#'#' %(%( %( Dj %(
Dj %(
%( %(r* r c N e Zd ZdZdZ d
dedededz dedz ddf
fd
Z xZS ) r z
.. versionadded:: 2.5.0
This class is for building a Vertical CRS.
.. warning:: geoid_model support only exists in PROJ >= 6.3.0
)z Vertical CRSNr r r r r3 c dd|t j | j t j |xs
t j d}|d|i|d< t
| | y)a!
Parameters
----------
name: str
The name of the Vertical CRS (e.g. NAVD88 height).
datum: Any
Anything accepted by :meth:`pyproj.crs.Datum.from_user_input`
vertical_cs: Any, optional
Input to create a Vertical Coordinate System accepted by
:meth:`pyproj.crs.CoordinateSystem.from_user_input`
or :class:`pyproj.crs.coordinate_system.VerticalCS`
geoid_model: str, optional
The name of the GEOID Model (e.g. GEOID12B).
r r r Nr r )r r r r r r% r& )r' r r r r
vert_crs_jsonr( s r) r& zVerticalCRS.__init__ sp , L!**51>>@!1!A!A+z|"ln
",2K+@M- (
'r* )NNr r0 s @r) r r sY (O #'"& !(!( !( 4Z !(
4Z !(
!( !(r* r c : e Zd ZdZdZdedee ddf fd Z xZ S )r zQ
.. versionadded:: 2.5.0
This class is for building a Compound CRS.
)z Compound CRSr r r3 Nc dd||D cg c]% }t j | j ' c}d}t | | yc c}w )a#
Parameters
----------
name: str
The name of the Compound CRS.
components: list[Any], optional
List of CRS to create a Compound Coordinate System.
List of anything accepted by :meth:`pyproj.crs.CRS.from_user_input`
r r )r r r r N)rz r r r% r& )r' r r componentcompound_crs_jsonr( s r) r& zCompoundCRS.__init__ sX L! ", ##I.;;=
*+
s *A
)
r+ r, r- r. r rE rB r r& r/ r0 s @r) r r s2 (O,S ,d3i ,D , ,r* r c 8 e Zd ZdZdZdedededdf fd Z xZS ) r zN
.. versionadded:: 2.5.0
This class is for building a Bound CRS.
)z Bound CRSr r r r3 Nc ddt j | j t j | j t j | j d}t | | y)a
Parameters
----------
source_crs: Any
Input to create a source CRS.
target_crs: Any
Input to create the target CRS.
transformation: Any
Input to create the transformation.
r r )r r r r r N)rz r r r
r% r& )r' r r r bound_crs_jsonr( s r) r& zBoundCRS.__init__& sh L --j9FFH --j9FFH 1 A A ln
(r* )r+ r, r- r. r r r& r/ r0 s @r) r r s6 %O)3 )C ) )QU ) )r* r )T)Cr. r? rR threadingrZ collections.abcr typingr r pyproj._crsr r r r r
r r r
r r r r pyproj.crs._cf1x8r r r r r r pyproj.crs.coordinate_operationr pyproj.crs.coordinate_systemr r r pyproj.enumsr r pyproj.exceptionsr pyproj.geodr compiler r localr! r> r rE rL ra rG ru rn rr rx rz r r r r r r r r rm r* r) r sw $
B S S 0 & DD
y
24 2T 2S 2,#S #S #L (&s &sSy &6#) 6S T [
[
|.@/3 @/F).( ).X06/ 06f4.( 4.n.(' .(b-(& -(` ,& ,@ )# )r*