o Dfc@s dZddlZddlZddlZddlZddlmZmZmZm Z ddl m Z m Z m Z mZmZmZmZmZmZmZmZmZddlmZmZmZmZmZmZddlmZddl m!Z!m"Z"m#Z#ddl$m%Z%m&Z&dd l'm(Z(dd l)m*Z*e+d ej,Z-Gd d d ej.Z/d3de0de1de2fddZ3de2de2fddZ4de2de2fddZ5de2de e2e6ffddZ7de e2e6ffddZ8dede1fdd Z9Gd!d"d"Z:Gd#d$d$e:Z;Gd%d&d&e;ZGd+d,d,e;Z?Gd-d.d.e;Z@Gd/d0d0e;ZAGd1d2d2e;ZBdS)4zx This module interfaces with PROJ to produce a pythonic interface to the coordinate reference system (CRS) information. N)AnyCallableOptionalUnion) _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 cs eZdZdZfddZZS)CRSLocalz Threading local instance for cython CRS class. For more details, see: https://github.com/pyproj4/pyproj/issues/782 csd|_tdSN)crssuper__init__self __class__G/home/ubuntu/webapp/venv/lib/python3.10/site-packages/pyproj/crs/crs.pyr$<szCRSLocal.__init__)__name__ __module__ __qualname____doc__r$ __classcell__r)r)r'r*r 4sr T projparams allow_jsonreturncCst|ts tdd|vrd|vr|rt|Sg}|D]:\}}t|ttfr3ddd|D}|dus=t |dkrF| d|qt |d krMq| d|d |qt d |S) NzCRS input is not a dictprojinit,cSg|]}t|qSr))str).0valr)r)r* Lz&_prepare_from_dict..True+False= ) isinstancedictrjsondumpsitemslisttuplejoinr7append_prepare_from_string)r0r1pjargskeyvaluer)r)r*_prepare_from_dictAs   rN in_crs_stringcCstdd|}d}||s5g}d}|D]}|s*||r*|d|d}q||qd|}d|vrFd |vrB|d 7}n|d 7}|d d  }|dr[t j dt dd|S)Nz [\s+]?=[\s+]?r?)+initz+projr4r3FrTr@ztype=crsr=z +type=crsz type=crsz +init=EPSGz +init=epsg)rPr4z'+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 startswithsplitinsertrIrHreplacestripwarningswarn FutureWarning)rOstarting_paramskvpairsfirst_item_insertedkvpairr)r)r*_prepare_from_proj_stringWs.        rbc Cst|ts td|std|d|vr=z tj|dd}Wnty0}ztd|d}~ww|s7tdt|}|St|rEt|}|S)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) rAr7rrCloads ValueErrorrNrrb)rOcrs_dicterrr)r)r*rJ}s$  rJ auth_name auth_codecCs|d|S)N:r))rirjr)r)r*_prepare_from_authoritysrlcCs td|S)NEPSG)rlrjr)r)r*_prepare_from_epsg rocCs@t|trdSt|tr|rdSt|dr|jdkrdSdS)NTshaper)F)rAintr7 isnumerichasattrrqrnr)r)r* _is_epsg_codes ruc @seZdZdZddeeddfddZeddZe d e d e e e fddfd d Z e d e e e fddfd dZe de ddfddZe de ddfddZe de ddfddZde fddZe deddfddZdeefddZe deddfd d!Ze d"e ddfd#d$Ze d%eddfd&d'Zdefd(d)Zejd*fd+e ee fd,edefd-d.Ze   dd/ed0eed1eed2eeddf d3d4Z de!efd5d6Z"d7edefd8d9Z#dd7ed:edefd;d<Z$ededfd=d>Z%ededfd?d@Z&ededfdAdBZ'ede!dfdCdDZ(edee fdEdFZ)ede fdGdHZ*ede fdIdJZ+ede!e,fdKdLZ-edee.fdMdNZ/edee0fdOdPZ1edee2fdQdRZ3edee4fdSdTZ5edee6fdUdVZ7edee8fdWdXZ9ede fdYdZZ:ede fd[d\Z;ejd*dfd]e ee fd^ed_eede fd`daZe?j@fd]e e?e fde fdhdiZAddke dee fdldmZBdd ee dke fdndoZC jdd ee dke de!eDfdpdqZEddree ddfdsdtZFddree ddfdudvZGedefdwdxZHedefdydzZIedefd{d|ZJedefd}d~ZKedefddZLedefddZMedefddZNeddZOd7edefddZPdee e ffddZQdee effddZRde fddZSde fddZTde fddZUdS)CRSaK A pythonic Coordinate Reference System manager. .. versionadded:: 2.0.0 See: :c:func:`proj_create` The functionality is based on other fantastic projects: * `rasterio `_ # noqa: E501 * `opendatacube `_ # noqa: E501 Attributes ---------- srs: str The string form of the user input used to create the CRS. Nr0r2cKsd}|rMt|tr |j}n@t|rt|}n7t|tr t|}n-t|tr*t|}n#t|t t fr>> 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 rQto_wktzInvalid CRS input: F)r1r@N)rArsrsruror7rJrBrNrFrGlenrlrtrxrrHr _localr")r&r0kwargs projstring projkwargsr)r)r*r$s0             z CRS.__init__cCs"|jjdur t|j|j_|jjS)zH Retrieve the Cython based _CRS object for this thread. N)r{r"rryr%r)r)r*_crs^s zCRS._crsricodecCs|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_inputrl)clsrirr)r)r*from_authoritygszCRS.from_authoritycC|t|S)zMake a CRS from an EPSG code Parameters ---------- code : int or str An EPSG code. Returns ------- CRS )rro)rrr)r)r* from_epsg{s z CRS.from_epsgin_proj_stringcC$t|s td||t|S)z .. versionadded:: 2.2.0 Make a CRS from a PROJ string Parameters ---------- in_proj_string : str A PROJ string. Returns ------- CRS zInvalid PROJ string: )rrrrb)rrr)r)r* from_proj4zCRS.from_proj4 in_wkt_stringcCr)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: )rrrrJ)rrr)r)r*from_wktrz CRS.from_wktrOcCr)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 )rrJ)rrOr)r)r* from_stringszCRS.from_stringcCs |jdd}|r d|S|jS)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 dmin_confidencerk) to_authorityrHry)r& auth_infor)r)r* to_strings  z CRS.to_stringrMcKst||r|S||fi|Sa 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 )rA)rrMr|r)r)r*rs zCRS.from_user_inputcCs(|jdurdSt|jj|jj|jjdS)zs Returns ------- pyproj.geod.Geod: Geod object based on the ellipsoid. N)arfb) ellipsoidrsemi_major_metreinverse_flatteningsemi_minor_metrer%r)r)r*get_geods z CRS.get_geod proj_dictcCr)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 )rrN)rrr)r)r* from_dictz CRS.from_dictcrs_jsoncCr)z .. versionadded:: 2.4.0 Create CRS from a CRS JSON string. Parameters ---------- crs_json: str CRS JSON string. Returns ------- CRS )rr)rrr)r)r* from_jsonrz CRS.from_jsonrgcCs|t|S)z .. versionadded:: 2.4.0 Create CRS from a JSON dictionary. Parameters ---------- crs_dict: dict CRS dictionary. Returns ------- CRS )rrCrD)rrgr)r)r*from_json_dict*szCRS.from_json_dictcCs`|}|dur iSdd}i}t|D]}|\}}|dur%||}|dur-|||<q|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 # noqa: E501 Returns ------- dict: PROJ params in dict format. NcSsb|dkrdS|dkrdSzt|WStyYnwzt|WSty0Yt|Sw)NtrueTfalseF)lowerrrrffloatr)r9r)r)r*_parseQs      zCRS.to_dict.._parseF)to_proj4_RE_PROJ_PARAMfinditergroups)r& proj_stringrrparamrLrMr)r)r*to_dict<s z CRS.to_dictF wkt_versionerrcheckc Cs@d||i}|jr.|jr.|jjr.|jr.|jj||d}|d|||jj|d<|S|jrJ|j D]}|j||d}|d||q4|S|j rh| }d|vr]|dd|d<|j rf|j j |d<|S|jr|j|jj|jj|jjd|jj |d <|jr|jj|d <|jj |d <|jr|jj |d <|jj r|jj j |d <|jr|jr|jjtvr|rtd|jjd|diS|t|jj|j|Sd|d<|Sd}|js|jr|j}|j |d<|sdn|jdd}|tvr|r|r td|jntdd|diS|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)rrtowgs84 geoid_modelname geoid_namegeopotential_datum_name)semi_major_axissemi_minor_axisrreference_ellipsoid_namelongitude_of_prime_meridianprime_meridian_namegeographic_crs_namehorizontal_datum_namez"Unsupported coordinate operation: latitude_longitudegrid_mapping_nameNprojected_crs_namer@_zCoordinate operation not found.)rxis_boundcoordinate_operationr source_crsto_cfpopupdate is_compound sub_crs_list is_vertical to_json_dictdatumrrrrrprime_meridian longitude geodetic_crs is_geographic method_namerrr[r\ is_projectedrYr) r&rrcf_dictsub_cfsub_crs vert_jsonrcoordinate_operation_namer)r)r*rjs                   z CRS.to_cfin_cfellipsoidal_cs cartesian_cs vertical_cscCsd}d|vr t|dSd|vrt|dS|d}|dur#tdt|}zt|}Wn ty8d}Ynw|d}|rJt|pDd||d } n%|rj||vrjt|} |duri| } t || d <t| } nt|d } |d kru| S|durt | |||d Szt |} Wntytd|dwt |dd| || |d} d} d|vrt | dt| jgt|dRd} d|vr| p| Std|d|d|d}td| p| |gdS)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_refrNz4CF projection parameters missing 'grid_mapping_name'rr)rrrcoordinate_system)rr)base_crs conversionrzUnsupported grid mapping name: r)rrrrrzWGS 84)r target_crstransformationrr)rrrr)r components)rvgetrrrKeyError GeographicCRSrr rDerivedGeographicCRSr ProjectedCRSBoundCRSrrr VerticalCRS CompoundCRS)rrrr unknown_namesrrgeographic_conversion_methodrgeographic_crsgeographic_crs_jsonconversion_method projected_crs bound_crs vertical_crsr)r)r*from_cfs(             z CRS.from_cfc Cs(g}dd}|jdkrZ|jrZ|j}|ddd}t|dkr)|d}|jd }d d |d d  d ddd}|j  d d}|rU|d||d<| ||j rk| |j j||d|S|jr|jr|jj r| |jj j||jd|S|jD] }| |q|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. cSs.z |jo |jjtvWStyYdSw)NF)rrrrr)r"r)r)r* rotated_pole|s  z"CRS.cs_to_cf..rotated_polez Temporal CRS time_originz 1875-05-20zZz-01-01rtimecalendarproleptic_gregorianr@rT) standard_name long_nameraxisrwz since units)r) type_namerrrrZrstriprz axis_inforrY unit_namerIrextendrrrrcs_to_cf) r& cf_axis_listr datum_jsonoriginrcf_temporal_axisrrr)r)r*rlsD         z CRS.cs_to_cfothercCs2zt|}Wn tyYdSw|j|jS)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)rvrrr is_exact_samer&r r)r)r*r s  zCRS.is_exact_sameignore_axis_ordercCs6zt|}Wn tyYdSw|jj|j|dS)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)rvrrrequals)r&r rr)r)r*rs  z CRS.equalscC|jjdurdS||jjSz .. versionadded:: 2.2.0 Returns ------- CRS: The geodeticCRS / geographicCRS from the CRS. N)rrr(r%r)r)r*r  zCRS.geodetic_crscCrz The base CRS of a BoundCRS or a DerivedCRS/ProjectedCRS, or the source CRS of a CoordinateOperation. Returns ------- CRS N)rrr(r%r)r)r*rs  zCRS.source_crscCrz .. versionadded:: 2.2.0 Returns ------- CRS: The hub CRS of a BoundCRS or the target CRS of a CoordinateOperation. N)rrr(r%r)r)r*rrzCRS.target_crscsfddjjDS) If the CRS is a compound CRS, it will return a list of sub CRS objects. Returns ------- list[CRS] csg|]}|qSr)r'r8rr%r)r*r:sz$CRS.sub_crs_list..rrr%r)r%r*r s zCRS.sub_crs_listcCsj|jr |jr |jjS|jr|jD] }|jr|jSqdS|jr3d|jjvr3|jjddSdS)z .. versionadded:: 2.6.0 Finds the UTM zone in a Projected CRS, Bound CRS, or Compound CRS Returns ------- Optional[str]: The UTM zone number and letter if applicable. zUTM ZONEz UTM ZONE N)rrutm_zonerrrupperrW)r&rr)r)r*rs   z CRS.utm_zonecC|jjS)zx Returns ------- str: The name of the CRS (from :cpp:func:`proj_get_name`). )rrr%r)r)r*r/zCRS.namecCr)zj Returns ------- str: The name of the type of the CRS object. )rrr%r)r)r*r9rz CRS.type_namecCr)a; Retrieves all relevant axis information in the CRS. If it is a Bound CRS, it gets the axis list from the Source CRS. If it is a Compound CRS, it gets the axis list from the Sub CRS list. Returns ------- list[Axis]: The list of axis information. )rrr%r)r)r*rCs z CRS.axis_infocCr)z{ Returns ------- AreaOfUse: The area of use object with associated attributes. )r area_of_user%r)r)r*rQrzCRS.area_of_usecCr)z .. versionadded:: 2.2.0 Returns ------- Ellipsoid: The ellipsoid object with associated attributes. )rrr%r)r)r*r[ z CRS.ellipsoidcCr)z .. versionadded:: 2.2.0 Returns ------- PrimeMeridian: The prime meridian object with associated attributes. )rrr%r)r)r*rgrzCRS.prime_meridiancCr)zX .. versionadded:: 2.2.0 Returns ------- Datum )rrr%r)r)r*rs z CRS.datumcCr)zc .. versionadded:: 2.2.0 Returns ------- CoordinateSystem )rrr%r)r)r*r~rzCRS.coordinate_systemcCr)zf .. versionadded:: 2.2.0 Returns ------- CoordinateOperation )rrr%r)r)r*rrzCRS.coordinate_operationcCr)zy .. versionadded:: 2.4.0 Returns ------- str: Remarks about object. )rremarksr%r)r)r*r rz CRS.remarkscCr)zt .. versionadded:: 2.4.0 Returns ------- str: Scope of object. )rscoper%r)r)r*r!rz CRS.scopeversionprettyoutput_axis_rulecCs.|jj|||d}|durtd|d|S)aA Convert the projection to a WKT string. Version options: - WKT2_2015 - WKT2_2015_SIMPLIFIED - WKT2_2019 - WKT2_2019_SIMPLIFIED - WKT1_GDAL - WKT1_ESRI .. versionadded:: 3.6.0 output_axis_rule Parameters ---------- version: pyproj.enums.WktVersion, optional The version of the WKT output. Default is :attr:`pyproj.enums.WktVersion.WKT2_2019`. pretty: bool, default=False If True, it will set the output to be a multiline string. output_axis_rule: bool, optional, default=None If True, it will set the axis rule on any case. If false, never. None for AUTO, that depends on the CRS and version. Returns ------- str )r"r#r$Nz.CRS cannot be converted to a WKT string of a 'zG' version. Select a different version of a WKT string or edit your CRS.)rrxr)r&r"r#r$wktr)r)r*rxs" z CRS.to_wktrQ indentationcCs$|jj||d}|durtd|S)a .. versionadded:: 2.4.0 Convert the object to a JSON string. Parameters ---------- pretty: bool, default=False If True, it will set the output to be a multiline string. indentation: int, default=2 If pretty is True, it will set the width of the indentation. Returns ------- str )r#r&Nz.CRS cannot be converted to a PROJ JSON string.)rto_jsonr)r&r#r& proj_jsonr)r)r*r'sz CRS.to_jsoncCs |jS)z .. versionadded:: 2.4.0 Convert the object to a JSON dictionary. Returns ------- dict )rrr%r)r)r*rs zCRS.to_json_dictcCs"|jj|d}|durtd|S)a< Convert the projection to a PROJ string. .. 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 # noqa: E501 Parameters ---------- version: pyproj.enums.ProjVersion The version of the PROJ string output. Default is :attr:`pyproj.enums.ProjVersion.PROJ_4`. Returns ------- str )r"Nz)CRS cannot be converted to a PROJ string.)rrr)r&r"r3r)r)r*rsz CRS.to_proj4FrcCs|jj|dS)ac Return the EPSG 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_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 ------- Optional[int]: The best matching EPSG code matching the confidence level. r)rto_epsg)r&rr)r)r*r*s#z CRS.to_epsgcC|jj||dS)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. rir)rrr&rirr)r)r*r7s&zCRS.to_authoritycCr+)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. r,)rlist_authorityr-r)r)r*r.as(zCRS.list_authorityrcC||jj|dSa .. 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(rto_3dr&rr)r)r*r2sz CRS.to_3dcCr/)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 r1)r(rto_2dr3r)r)r*r4sz CRS.to_2dcCr)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. )rrr%r)r)r*r zCRS.is_geographiccCr)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. )rrr%r)r)r*rr5zCRS.is_projectedcCr)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. )rrr%r)r)r*rszCRS.is_verticalcCr)zY Returns ------- bool: True if CRS is bound. )rrr%r)r)r*rrz CRS.is_boundcCr)z} .. versionadded:: 3.1.0 Returns ------- bool: True if CRS is compound. )rrr%r)r)r*rrzCRS.is_compoundcCr)z .. versionadded:: 2.2.0 Returns ------- bool: True if CRS is local/engineering. )ris_engineeringr%r)r)r*r6rzCRS.is_engineeringcCr)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%r)r)r*r7s zCRS.is_geocentriccCr)z .. versionadded:: 3.2.0 Returns ------- bool: True if CRS is a Derived CRS. )r is_derivedr%r)r)r*r8rzCRS.is_derivedcCs ||Sr!)rr r)r)r*__eq__rpz CRS.__eq__cCs d|jiS)Nryryr%r)r)r* __getstate__!rpzCRS.__getstate__statecCs|j|t|_dSr!)__dict__rr r{)r&r<r)r)r* __setstate__$s  zCRS.__setstate__cCs t|Sr!)hashrxr%r)r)r*__hash__(s z CRS.__hash__cCs|jSr!r:r%r)r)r*__str__+sz CRS.__str__c Csg}|jD] }|dt|dgqd|}d}d}|jr(|jjr(t|j}n;|jr=|jr=t|jj}d|jjd}n&g}dg}|j D]} | t| j|d| jdgqEd|}d|}d} |j rxdddt|j d |j j dg} | } t| d kr| n d | dd d g} |pd }d|jd| d|jd|pdd|d|jpdd| d|jd|jpdd|jpdd||S)Nz-  rwz Source CRS: z Sub CRS: |zCoordinate Operation: z- name: z - method: 2r@z...z - undefined  Name: z Axis Info [rz]: z Area of Use: z - undefinedzDatum: z - Ellipsoid: z - Prime Meridian: )rrr7rHr axis_listrrrrrIrrrrzrrrrr) r&axis_info_listr axis_info_strsource_crs_repr sub_crs_reprcoordinate_system_namecoordinate_system_namessub_crs_repr_listrrsrs_reprr)r)r*__repr__.sl         &   z CRS.__repr__r!)NNN)F)FrQ)r))Nr))Vr+r,r-r.rrr$propertyr classmethodr7rrrrrrrrrrrrrBrrrrr WKT2_2019boolr staticmethodrrFrr rrrrrrrrr rrrr rrrr rr rr rr r!rxr'rrPROJ_5rr*rrr.r2r4rrrrrr6r7r8r9r;r>r@rArOr)r)r)r*rvs(  0   ~>              , %+ ,     rvcseZdZdZedeedffddZddZe de dd ffd d Z ede d fd d Z ede d fddZede d fddZeded fddZdde edd fddZZS)CustomConstructorCRSz This class is a base class for CRS classes that use a different constructor than the main CRS class. .. versionadded:: 3.2.0 See: https://github.com/pyproj4/pyproj/issues/847 r2.cCst)zt These are the type names of the CRS class that are expected when using the from_* methods. )NotImplementedErrorr%r)r)r*_expected_typesrsz$CustomConstructorCRS._expected_typescCs*|j|jvrtd|jd|jddS)zl This validates that the type of the CRS is expected when using the from_* methods. z Invalid type z . Expected .N)rrXrr%r)r)r* _check_typezs z CustomConstructorCRS._check_typerMrvc s<t||r|S||}tt|j|fi|||Sr)rA__new__r#rVr$rZ)rrMr|r"r'r)r*rs  z$CustomConstructorCRS.from_user_inputcC|jjdurdSt|jjSr)rrrvr%r)r)r*r z!CustomConstructorCRS.geodetic_crscCr\r)rrrvr%r)r)r*rs zCustomConstructorCRS.source_crscCr\r)rrrvr%r)r)r*rr]zCustomConstructorCRS.target_crscCsdd|jjDS)rcSr6r))rvrr)r)r*r:r;z5CustomConstructorCRS.sub_crs_list..rr%r)r)r*rs z!CustomConstructorCRS.sub_crs_listNrcCst|jj|dSr0)rvrr2r3r)r)r*r2szCustomConstructorCRS.to_3dr!)r+r,r-r.rPrGr7rXrZrQrrrrrrrFrr2r/r)r)r'r*rVhs       rVc sBeZdZdZdZ   d dededeed dffd d ZZ S) rzR .. versionadded:: 2.5.0 This class is for building a Geographic CRS )zGeographic CRSzGeographic 2D CRSzGeographic 3D CRSrurn:ogc:def:ensemble:EPSG::6326Nrrrr2csZt|}dd|t|ptd}|ddkr!||d<n||d<t|dS) 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)$schematyperrra DatumEnsembledatum_ensemblerN)r rrr rr#r$)r&rrrrr'r)r*r$s  zGeographicCRS.__init__)rr^N r+r,r-r.rXr7rrr$r/r)r)r'r*rsrc sDeZdZdZdZ  d dededeeded df fd d ZZ S) rzZ .. versionadded:: 2.5.0 This class is for building a Derived Geographic CRS )zDerived Geographic CRSzDerived Geographic 2D CRSzDerived Geographic 3D CRSNrrrrrr2csFdd|t|t|t|ptd}t|dS)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_rr`rarrrrN)rvrrr r rr#r$)r&rrrrderived_geographic_crs_jsonr'r)r*r$$s  zDerivedGeographicCRS.__init__)Nr) r+r,r-r.rXrrr7r$r/r)r)r'r*rs  rcs8eZdZdZdZ  d dededdffd d ZZS) GeocentricCRSzR .. versionadded:: 3.2.0 This class is for building a Geocentric CRS )zGeocentric CRSrurn:ogc:def:datum:EPSG::6326rrr2Nc sTdd|t|ddddddd d d ddd d dddgdd}t|dS)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)subtyperr`rarrrN)r rrr#r$)r&rrgeocentric_crs_jsonr'r)r*r$Ss0 zGeocentricCRS.__init__)rrh) r+r,r-r.rXr7rr$r/r)r)r'r*rgJsrgc sJeZdZdZdZ   d dededeedeed df fd d ZZ S) rzR .. versionadded:: 2.5.0 This class is for building a Projected CRS. )z Projected CRSzDerived Projected CRSrNrrrrr2csLdd|t|p tt|t|ptd}t|dS)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_rreN) rvrrrr r rr#r$)r&rrrr proj_crs_jsonr'r)r*r$s"zProjectedCRS.__init__)rNN) r+r,r-r.rXrr7rr$r/r)r)r'r*rs"rc sHeZdZdZdZ  d dededeedeeddf fd d ZZ S) rz .. versionadded:: 2.5.0 This class is for building a Vetical CRS. .. warning:: geoid_model support only exists in PROJ >= 6.3.0 )z Vertical CRSNrrrrr2csNdd|t|t|ptd}|durd|i|d<t|dS)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_rrvNrr)r rrr rr#r$)r&rrrr vert_crs_jsonr'r)r*r$s  zVerticalCRS.__init__)NNrdr)r)r'r*rs  rcs6eZdZdZdZdedeeddffdd ZZ S) rzQ .. versionadded:: 2.5.0 This class is for building a Compound CRS. )z Compound CRSrrr2Ncs(dd|dd|Dd}t|dS)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_rcSsg|] }t|qSr))rvrr)r8 componentr)r)r*r:s z(CompoundCRS.__init__..)r`rarrN)r#r$)r&rrcompound_crs_jsonr'r)r*r$s  zCompoundCRS.__init__) r+r,r-r.rXr7rFrr$r/r)r)r'r*r&rcs6eZdZdZdZdedededdffdd ZZS) rzN .. versionadded:: 2.5.0 This class is for building a Bound CRS. )z Bound CRSrrrr2Ncs>ddt|t|t|d}t|dS)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`rarrrN)rvrrr r#r$)r&rrrbound_crs_jsonr'r)r*r$ s    zBoundCRS.__init__)r+r,r-r.rXrr$r/r)r)r'r*rr|r)T)Cr.rCrS threadingr[typingrrrr pyproj._crsrrrr r r r r rrrrpyproj.crs._cf1x8rrrrrrpyproj.crs.coordinate_operationrpyproj.crs.coordinate_systemrrr pyproj.enumsrrpyproj.exceptionsr pyproj.geodrcompilerkrlocalr rBrSr7rNrbrJrrrlrorurvrVrrrgrrrrr)r)r)r*sZ8      & P,3710