PyCRS
PyCRS copied to clipboard
Projection creation, and conversion to misc formats
Is there any maintenance activity for this package? I found some mistakes. Thx!
It seems one could install `pycrs` on [Pyodide](https://pyodide.org) and hence [JupyterLite](https://github.com/jtpio/jupyterlite) using [micropip](https://pyodide.org/en/stable/usage/api/micropip-api.html) if it would have a source wheel package on PyPI. Maybe worth a thought?
the specified projection name 'Gauss_Kruger' could not be found in the list of supported projections
NotImplementedError: Unsupported projection: The specified projection name 'Gauss_Kruger' could not be found in the list of supported projections Do we have a plan to implement Gauss_Kruger projection?
Hello. I'm very excited about this project. It's what the World needs now. I'm just getting started but have hit an error pretty early. I'm trying to load this projection:...
One more projection format we could support is the CF metadata standard, which appears to be used by NetCDF files and possibly elsewhere. Links below could be a good starting...
As mentioned in > https://www.reddit.com/r/gis/comments/bdnir7/prj2epsgorgsearch_not_working_anymoreanyone_have/ It seems the `prj2epsg/search` will give 500 internal error, which makes `utils.wkt_to_epgs()` function doesn't work anymore, is it possible to change to another searching web...
I'm converting an Albers projection from WKT to Proj4. My WKT is: `PROJCS["NAD_1983_Albers",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.2572221010042,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["standard_parallel_1",29.5],PARAMETER["standard_parallel_2",45.5],PARAMETER["latitude_of_center",23],PARAMETER["longitude_of_center",-96],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]` Then this line: `pycrs.parse.from_ogc_wkt(wkt).to_proj4()` Gives: `'+proj=aea +datum=NAD83 +ellps=GRS80 +a=6378137.0 +rf=298.2572221010042 +pm=0 +lat_1=29.5 +lat_2=45.5 +lonc=-96.0 +x_0=0.0 +y_0=0.0...
Can PyCRS pass this test? ```python import rasterio import pycrs def test_pycrs_parses_ogc_wkt(): # data from https://spatialreference.org/ref/epsg/4326/ epsg_4326_to_ogc_wkt = 'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]' rasterio_ogc_wkt = rasterio.crs.CRS.from_string('EPSG:4326').to_wkt(morph_to_esri_dialect=False) epsg_4326_ogc_crs = pycrs.parse.from_ogc_wkt(epsg_4326_to_ogc_wkt) rasterio_ogc_crs = pycrs.parse.from_ogc_wkt(rasterio_ogc_wkt)...
http://geotiff.maptools.org/proj_list/lambert_conic_conformal_2sp.html An error related to this arose while trying to parse WKT, e.g. http://spatialreference.org/ref/esri/102740/ogcwkt/ ``` >>> esri_102740_ogc_wkt = 'PROJCS["NAD_1983_StatePlane_Texas_South_Central_FIPS_4204_Feet",GEOGCS["GCS_North_American_1983",DATUM["North_American_Datum_1983",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["False_Easting",1968500],PARAMETER["False_Northing",13123333.33333333],PARAMETER["Central_Meridian",-99],PARAMETER["Standard_Parallel_1",28.38333333333333],PARAMETER["Standard_Parallel_2",30.28333333333334],PARAMETER["Latitude_Of_Origin",27.83333333333333],UNIT["Foot_US",0.30480060960121924],AUTHORITY["EPSG","102740"]]' >>> pycrs.parse.from_ogc_wkt(esri_102740_ogc_wkt) ... File "/home/joe/miniconda3/envs/gis-dataprocessing/lib/python3.7/site-packages/pycrs/parse.py", line 303, in _parse_top raise NotImplementedError("Unsupported...
I'm getting an unexpected exception. It's easy to reproduce: ``` >>> import pycrs >>> pycrs.__version__ '1.0.0' >>> crs = pycrs.parse.from_epsg_code(2230) Traceback (most recent call last): File "", line 1, in...