Not able to process certain GeoTIFFs

[kasra-keshavarz]

Problem

The problem is that for this specific raster, I cannot extract zonal statistics. I am not sure where the problem is, but included an MCVE. Let me know if you need more information.

MCVE

Here is a reproducible MCVE (using the data.zip file attached here):

wget https://github.com/user-attachments/files/18482547/data.zip
unzip data.zip
exactextract --polygons ./data/polygons/bcalgary_subbasins.shp --raster ./data/raster/nalcms_2020_landcover.tif --stat frac --fid COMID --output=fraction_landcover.csv --progress # results in an empty output
exactextract --polygons ./data/polygons/bcalgary_subbasins.shp --raster ./data/raster/nalcms_2020_landcover.tif --stat mean --fid COMID --output=mean_landcover.csv --progress # just a test - results in `nan` values

exactextrac version

$ ./exactextract 
Zonal statistics using exactextract: version 0.3.0-dev (d33bb54)
... 

Attachements:

data.zip

[kasra-keshavarz]

Forgot to mention that equivalent commands work well in exactextractr.

[dbaston]

The command-line interface does not check the CRS of its inputs. At a minimum, it should check the CRS of the inputs and produce a warning if they don't match. Some related discussion in #108

[dbaston]

Warning added in f312a11