New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Stac load with non-native CRS returns empty image in some cases #133

Open

Kirill888 opened this issue Nov 29, 2023 · 1 comment

Member

Kirill888 commented Nov 29, 2023

Source data:

low resolution tiff in epsg:4326 with tiny region of interest (4x4 pixel area kinda thing)

Desired output:

UTM projection, much higher resolution, 4x4px -> 100x100px
configured via geobox paramater

Loading same source file with projection change alone and without cropping works fine

Member Author

Kirill888 commented Nov 29, 2023

Culprit appears to be here:

odc-stac/odc/stac/_load.py

Lines 859 to 862 in 138846f

    
           geom = item.safe_geometry(gbt.base.crs) 
        
           if geom is None: 
        
               raise ValueError("Can not process items without defined footprint") 
        
           yield from gbt.tiles(geom)

Almost global footprint of the source dataset doesn't project well into UTM.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment