Hi,
I suspect you are running into
https://gdal.org/en/stable/tutorials/osr_api_tut.html#crs-and-axis-order
Concretely try adding:
source_srs.SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER)
at line 34.
Even
Le 10/01/2025 à 17:49, ky via gdal-dev a écrit :
Hello all, hope you are well
Got some unexpected results computing area using GDAL==3.4.1 and
official Brazilian projection, where results are about half of what is
reported as correct. Does anyone knows what is causing it ?
Just to compare shapely/geopandas yields correct results using same
projections.
code can be found on this gist:
https://gist.github.com/rupestre-campos/a82a1403cb47a8153f07c775d7902f50
the linked scrip outputs
Area ha farm: 6.0
Area ha calculated (gdal): 3.23
Area ha calculated (geopandas): 6.0
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
http://www.spatialys.com
My software is free, but my time generally not.
Butcher of all kinds of standards, open or closed formats. At the end, this is
just about bytes.
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
