AssertionError when trying to use Netherlands DEMs about demquery

[kylebarron]

Can you create an minimal example of using demquery directly? E.g. try using query_points(points) directly and find the point where it fails?

[craigxchen]

Hi Kyle,

Thanks for responding.

I tried it using a minimal query and it's failing on every single point of the graph.

It seems like the assertion error doesn't have anything to do with the specific points right? It's coming from the _check_bounds method. My best guess is that somehow the TIFs aren't getting tiled together correctly?

[kylebarron]

I think it's likely your DEM is not in the same projection as your input points. This package won't do any automatic reprojection of the raster file, and if the raster data is in some local projection, it won't work and _check_bounds will fail.

[craigxchen]

Do you have any suggestions on how I can fix the projection issue? Do you think there's a quick fix or is it something more involved?

[kylebarron]

You can look into GDAL/rasterio for reprojecting your data. https://gdal.org/programs/gdalwarp.html