Open-EO / Open-EO/openeo-python-driver
geometry_to_crs reprojection is incorrect
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 12
- Forks
- 8
- Avg merge
- 6d 3h
- Merged PRs (30d)
- 4
Description
The projection from geometry_to_crs differs from the reprojection we get in qgis. This is especially visible when you create polygons near the north pole, such as FeatureCollection09.json. The qgis reprojection appears to be correct because there the area remains identical, where our reprojected polygons have a larger area compared to the original (area is calculated using WGS84 ellipsoid).
For nearly all use-cases the difference is negligible but it is still worth investigating if we are using the pyproj library wrong in some way.
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start with geometry_to_crs in openeo_driver/util/utm.py around line 59 and reproduce the issue using tests/data/geojson/FeatureCollection09.json. Compare the pyproj reprojection with the QGIS result, focusing on the preserved area of polygons near the north pole. Done means identifying whether pyproj is used incorrectly and correcting the reprojection so the result matches the expected area.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- backend
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Quiet
- Clarity
- Mostly clear
- Newbie friendliness
- 48/100