Open jbpaduan opened 2 years ago
Looking at gdal utilities… Looks like dgal_translate can be used to convert a geotiff to cylindrical equidistant, e.g.:
gdal_translate -a_srs EPSG:4087 <src.tif> <dst.tif>
And then gdal_polygonize.py can be used to generate Polygons representing the area where ther's valid data in a grid.
The Geopandas package has a .area() method that we could use:
https://geopandas.org/en/stable/docs/reference/api/geopandas.GeoSeries.area.html#geopandas.GeoSeries.area
It would require generating a Polygon representing the boundary of the irregular area surveyed.
Here's another example implementation: https://stackoverflow.com/a/69220519/1281657