diff options
| -rw-r--r-- | .gitignore | 1 | ||||
| -rw-r--r-- | python-pycrs.spec | 225 | ||||
| -rw-r--r-- | sources | 1 |
3 files changed, 227 insertions, 0 deletions
@@ -0,0 +1 @@ +/PyCRS-1.0.2.tar.gz diff --git a/python-pycrs.spec b/python-pycrs.spec new file mode 100644 index 0000000..d663aa5 --- /dev/null +++ b/python-pycrs.spec @@ -0,0 +1,225 @@ +%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: + +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): + +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]]' + +### 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: + +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): + +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]]' + +### 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: + +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): + +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]]' + +### 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 <Python_Bot@openeuler.org> - 1.0.2-1 +- Package Spec generated @@ -0,0 +1 @@ +8ad01d22bf1b114e885b98567ea2e528 PyCRS-1.0.2.tar.gz |