Reading GeoTIFFs with rioxarray/rasterio

chpolste commented 4 weeks ago

A basic reader for GeoTIFF files using rioxarray under the hood. To address #120.

Implements GeoTIFFReader, GeoTIFFFieldList, GeoTIFFField, GeoTIFFMetadata and GeoTIFFGeography.
The unit test comes with a small multi-band test file.
New geotiff optional dependencies in pyproject.toml.

There is not much in terms of metadata extraction other than the projection information yet, but many GeoTIFFs don't provide much metadata in the first place, unfortunately.

I have configured the reader to apply masking and scaling to the values by default, even though rioxarray doesn't.

Small things adjusted along the way:

Turned shape of the core.geography.Geography abstract class into a method (was a property), since that is how it seems to be used in practice.
Added a TransverseMercator projection, which seems to be common for GeoTIFF data.
Changed utils.projections.cf.CF_PARAMS_TO_GLOBE_KWARGS["inverse_flattening"] to "inverse_flattening" (was "flattening").

FussyDuck commented 4 weeks ago

All committers have signed the CLA.

sandorkertesz commented 3 weeks ago

@chpolste , thank you very much for this development. I tested it and works nicely. I made a few comments, once these and the branch conflict are sorted out the PR can be merged.

As for the documentation, what notebook example do you think should be added? What should it highlight?

codecov-commenter commented 1 week ago

Codecov Report

Attention: Patch coverage is 85.71429% with 3 lines in your changes missing coverage. Please review.

Project coverage is 91.08%. Comparing base (acda540) to head (2598847). Report is 7 commits behind head on develop.

Files with missing lines	Patch %	Lines
tests/readers/test_geotiff_reader.py	85.71%	2 Missing and 1 partial :warning:

Additional details and impacted files

```diff @@ Coverage Diff @@ ## develop #503 +/- ## =========================================== - Coverage 91.09% 91.08% -0.02% =========================================== Files 147 148 +1 Lines 10744 10765 +21 Branches 510 511 +1 =========================================== + Hits 9787 9805 +18 - Misses 794 796 +2 - Partials 163 164 +1 ```

:umbrella: View full report in Codecov by Sentry.
:loudspeaker: Have feedback on the report? Share it here.

🚨 Try these New Features:

Flaky Tests Detection - Detect and resolve failed and flaky tests

ecmwf / earthkit-data

Reading GeoTIFFs with rioxarray/rasterio #503

Codecov Report