%global _empty_manifest_terminate_build 0 Name: python-PyCRS Version: 1.0.2 Release: 1 Summary: GIS package for reading, writing, and converting between CRS formats. License: MIT URL: http://github.com/karimbahgat/PyCRS Source0: https://mirrors.nju.edu.cn/pypi/web/packages/50/0b/33c6ab39701d982eabfdc732d920862d79d7598893c7291c5f6ec1e66d3c/PyCRS-1.0.2.tar.gz BuildArch: noarch %description ## Recipes ### Modifying the CS Instance In most case you will only ever need to load a CRS and convert it to some format. Sometimes, however, you may want to tweak the parameters of your CS type instance. Knowing the composition of your CS type instance, this is as easy as setting/replacing the desired attributes. Let's demonstrate some examples using the World Robinson projection: >>> crs = pycrs.parse.from_esri_code(54030) # Robinson projection from esri code >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["WGS_1984", SPHEROID["WGS_1984", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' Here is a map of the default Robinson projection: ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_orig.png "Defualt Robinson") Let's say we wanted to switch its datum from WGS84 to NAD83, we could do it like so: >>> crs.geogcs.datum = pycrs.elements.datums.NAD83() >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 0], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' Or let's say we wanted to switch its prime meridian, so that the longitude axis is centered closer to the Pacific instead of over Greenwhich: >>> crs.geogcs.prime_mer.value = 160.0 >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 160], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' And here is what that map would look like (the odd-looking lines is just a rendering issue due to polygons that cross the meridian): ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_tweak2.png "Modified Robinson") Or if we just switch the projection type alltogether: >>> crs.proj = pycrs.elements.projections.Sinusoidal() >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 160], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Sinusoidal"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_tweak3.png "Modified Sinusoidal") ### Coordinate Transformations A common reason for wanting to convert between CRS formats, is if you want to transform coordinates from one coordinate system to another. In Python this is typically done with the PyProj module, which only takes proj4 format. Using PyCRS we can easily define the original coordinate system that we want to convert and get its proj4 representation: >>> fromcrs = pycrs.parse.from_epsg_code(4326) # WGS84 projection from epsg code >>> fromcrs_proj4 = fromcrs.to_proj4() We can then use PyCRS to define our target projection from the format of your choice, before converting it to the proj4 format that PyProj expects: >>> tocrs = pycrs.parse.from_esri_code(54030) # Robinson projection from esri code >>> tocrs_proj4 = tocrs.to_proj4() With the source and target projections defined in the proj4 crs format, we are ready to transform our data coordinates with PyProj: >>> import pyproj >>> fromproj = pyproj.Proj(fromcrs_proj4) >>> toproj = pyproj.Proj(tocrs_proj4) >>> lng,lat = -76.7075, 37.2707 # Williamsburg, Virginia :) >>> pyproj.transform(fromproj, toproj, lng, lat) (-6766170.001635834, 3985755.032695593) ### Writing a Shapefile .prj file After you transform your data coordinates you may also wish to save the data back to file along with the new crs. With PyCRS you can do this in a variety of crs format. For instance, to write a shapefile .prj file: %package -n python3-PyCRS Summary: GIS package for reading, writing, and converting between CRS formats. Provides: python-PyCRS BuildRequires: python3-devel BuildRequires: python3-setuptools BuildRequires: python3-pip %description -n python3-PyCRS ## Recipes ### Modifying the CS Instance In most case you will only ever need to load a CRS and convert it to some format. Sometimes, however, you may want to tweak the parameters of your CS type instance. Knowing the composition of your CS type instance, this is as easy as setting/replacing the desired attributes. Let's demonstrate some examples using the World Robinson projection: >>> crs = pycrs.parse.from_esri_code(54030) # Robinson projection from esri code >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["WGS_1984", SPHEROID["WGS_1984", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' Here is a map of the default Robinson projection: ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_orig.png "Defualt Robinson") Let's say we wanted to switch its datum from WGS84 to NAD83, we could do it like so: >>> crs.geogcs.datum = pycrs.elements.datums.NAD83() >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 0], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' Or let's say we wanted to switch its prime meridian, so that the longitude axis is centered closer to the Pacific instead of over Greenwhich: >>> crs.geogcs.prime_mer.value = 160.0 >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 160], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' And here is what that map would look like (the odd-looking lines is just a rendering issue due to polygons that cross the meridian): ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_tweak2.png "Modified Robinson") Or if we just switch the projection type alltogether: >>> crs.proj = pycrs.elements.projections.Sinusoidal() >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 160], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Sinusoidal"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_tweak3.png "Modified Sinusoidal") ### Coordinate Transformations A common reason for wanting to convert between CRS formats, is if you want to transform coordinates from one coordinate system to another. In Python this is typically done with the PyProj module, which only takes proj4 format. Using PyCRS we can easily define the original coordinate system that we want to convert and get its proj4 representation: >>> fromcrs = pycrs.parse.from_epsg_code(4326) # WGS84 projection from epsg code >>> fromcrs_proj4 = fromcrs.to_proj4() We can then use PyCRS to define our target projection from the format of your choice, before converting it to the proj4 format that PyProj expects: >>> tocrs = pycrs.parse.from_esri_code(54030) # Robinson projection from esri code >>> tocrs_proj4 = tocrs.to_proj4() With the source and target projections defined in the proj4 crs format, we are ready to transform our data coordinates with PyProj: >>> import pyproj >>> fromproj = pyproj.Proj(fromcrs_proj4) >>> toproj = pyproj.Proj(tocrs_proj4) >>> lng,lat = -76.7075, 37.2707 # Williamsburg, Virginia :) >>> pyproj.transform(fromproj, toproj, lng, lat) (-6766170.001635834, 3985755.032695593) ### Writing a Shapefile .prj file After you transform your data coordinates you may also wish to save the data back to file along with the new crs. With PyCRS you can do this in a variety of crs format. For instance, to write a shapefile .prj file: %package help Summary: Development documents and examples for PyCRS Provides: python3-PyCRS-doc %description help ## Recipes ### Modifying the CS Instance In most case you will only ever need to load a CRS and convert it to some format. Sometimes, however, you may want to tweak the parameters of your CS type instance. Knowing the composition of your CS type instance, this is as easy as setting/replacing the desired attributes. Let's demonstrate some examples using the World Robinson projection: >>> crs = pycrs.parse.from_esri_code(54030) # Robinson projection from esri code >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["WGS_1984", SPHEROID["WGS_1984", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' Here is a map of the default Robinson projection: ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_orig.png "Defualt Robinson") Let's say we wanted to switch its datum from WGS84 to NAD83, we could do it like so: >>> crs.geogcs.datum = pycrs.elements.datums.NAD83() >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 0], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' Or let's say we wanted to switch its prime meridian, so that the longitude axis is centered closer to the Pacific instead of over Greenwhich: >>> crs.geogcs.prime_mer.value = 160.0 >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 160], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Robinson"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' And here is what that map would look like (the odd-looking lines is just a rendering issue due to polygons that cross the meridian): ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_tweak2.png "Modified Robinson") Or if we just switch the projection type alltogether: >>> crs.proj = pycrs.elements.projections.Sinusoidal() >>> crs.to_ogc_wkt() 'PROJCS["Unknown", GEOGCS["Unknown", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980", 6378137.0, 298.257222101]], PRIMEM["Greenwich", 160], UNIT["degree", 0.017453292519943295], AXIS["Lon", EAST], AXIS["Lat", NORTH]], PROJECTION["Sinusoidal"], PARAMETER["Central_Meridian", 0], PARAMETER["false_easting", 0], PARAMETER["false_northing", 0], UNIT["Meters", 1.0], AXIS["X", EAST], AXIS["Y", NORTH]]' ![](https://github.com/karimbahgat/pycrs/raw/master/testrenders/docs_tweak3.png "Modified Sinusoidal") ### Coordinate Transformations A common reason for wanting to convert between CRS formats, is if you want to transform coordinates from one coordinate system to another. In Python this is typically done with the PyProj module, which only takes proj4 format. Using PyCRS we can easily define the original coordinate system that we want to convert and get its proj4 representation: >>> fromcrs = pycrs.parse.from_epsg_code(4326) # WGS84 projection from epsg code >>> fromcrs_proj4 = fromcrs.to_proj4() We can then use PyCRS to define our target projection from the format of your choice, before converting it to the proj4 format that PyProj expects: >>> tocrs = pycrs.parse.from_esri_code(54030) # Robinson projection from esri code >>> tocrs_proj4 = tocrs.to_proj4() With the source and target projections defined in the proj4 crs format, we are ready to transform our data coordinates with PyProj: >>> import pyproj >>> fromproj = pyproj.Proj(fromcrs_proj4) >>> toproj = pyproj.Proj(tocrs_proj4) >>> lng,lat = -76.7075, 37.2707 # Williamsburg, Virginia :) >>> pyproj.transform(fromproj, toproj, lng, lat) (-6766170.001635834, 3985755.032695593) ### Writing a Shapefile .prj file After you transform your data coordinates you may also wish to save the data back to file along with the new crs. With PyCRS you can do this in a variety of crs format. For instance, to write a shapefile .prj file: %prep %autosetup -n PyCRS-1.0.2 %build %py3_build %install %py3_install install -d -m755 %{buildroot}/%{_pkgdocdir} if [ -d doc ]; then cp -arf doc %{buildroot}/%{_pkgdocdir}; fi if [ -d docs ]; then cp -arf docs %{buildroot}/%{_pkgdocdir}; fi if [ -d example ]; then cp -arf example %{buildroot}/%{_pkgdocdir}; fi if [ -d examples ]; then cp -arf examples %{buildroot}/%{_pkgdocdir}; fi pushd %{buildroot} if [ -d usr/lib ]; then find usr/lib -type f -printf "/%h/%f\n" >> filelist.lst fi if [ -d usr/lib64 ]; then find usr/lib64 -type f -printf "/%h/%f\n" >> filelist.lst fi if [ -d usr/bin ]; then find usr/bin -type f -printf "/%h/%f\n" >> filelist.lst fi if [ -d usr/sbin ]; then find usr/sbin -type f -printf "/%h/%f\n" >> filelist.lst fi touch doclist.lst if [ -d usr/share/man ]; then find usr/share/man -type f -printf "/%h/%f.gz\n" >> doclist.lst fi popd mv %{buildroot}/filelist.lst . mv %{buildroot}/doclist.lst . %files -n python3-PyCRS -f filelist.lst %dir %{python3_sitelib}/* %files help -f doclist.lst %{_docdir}/* %changelog * Tue Apr 11 2023 Python_Bot - 1.0.2-1 - Package Spec generated