Jiang, J. F.
Jiang, J. F.
@greglucas The former answer may not work,but the latter one does work in utm projection.However, it does nor work for epsg=4326. Next, I will show you in details. The former...
``` def visualize_single_band_data(DataArray,label='Elevation(m)',map_crs=ccrs.PlateCarree(),cmap='terrain'): epsg = cn_dem.rio.crs.to_epsg() data_crs = CRS.from_epsg(epsg) # data_crs_str = DataArray.rio.crs.to_proj4() # data_crs = ccrs.Projection(data_crs_str) bounds = DataArray.rio.bounds() extent = bounds[0],bounds[2],bounds[1],bounds[3] nodata = DataArray._FillValue DataArray = DataArray.where(DataArray!= nodata)...
@rcomer Thanks.It would be nice if Gridliner had improvements in the future in terms of the ticker being displayed .