Deltares / hydromt_sfincs

https://deltares.github.io/hydromt_sfincs
GNU General Public License v3.0
20 stars 15 forks source link

Plotting extent #123

Open lgrimley opened 1 year ago

lgrimley commented 1 year ago

I have a model that covers more than one UTM zone. When I use the plot_basemap() function it cuts off and only shows the area of the model that is in the UTM zone that the model as been assigned to. Is there a way for it to clip to the extent of the region instead? region

DirkEilander commented 1 year ago

Thanks for the suggestion @lgrimley! I'll need to figure out why this happens, but agree this should be changed. If it is urgent, as a work around, you can copy the plot_basemaps method and adapt/delete the line where the plot extent is set https://github.com/Deltares/hydromt_sfincs/blob/main/hydromt_sfincs/plots.py#L189

DirkEilander commented 1 year ago

Currently the approach is to use the projection of the data. This could lead to negative coordinates if a model spans multiple UTM zones with cartopy doesn't seem to like (this assumption needs to be checked). Instead of a UTM zone it would perhaps be better if the model uses another CRS that can describe the entire model in that case.

@lgrimley can you check if this error persists if you reproject the model grid to e.g. WGS84 and then plot it?

sf = SfincsModel(<root>, mode='r')
sf._grid = sf.grid.raster.reproject('epsg:4326')
sf.plot_basemap()