Closed JesJehle closed 1 month ago
Looks like the query coords are off. Hmmm
The issue is that there is no reprojection being done. Some orthophotos are in a CRS that is not EPSG:4326. The coordinates are then out of range for EPSG:4326.
As far as I remember, the consensus was that we store the incoming orthophotos in the CRS that we get it in for maximum accuracy and then reproject it to 4326 for the visualization.
For the bounds, which is where the OSM query gets its data from, we need to change this line...
...to bounds = rasterio.warp.transform_bounds(src_crs=src.crs, dst_crs="EPSG:4326", *src.bounds)
https://rasterio.readthedocs.io/en/latest/api/rasterio.warp.html#rasterio.warp.transform_bounds
For the COGs, there does not seem to be a reprojection hapenning here:
I think the correct flag is -t_crs=EPSG:4326
Perfect! We can rerun the process after the tiffs are imported.
A lot of these:
from id: 1143 upwards.