From ab309b4c330b244966efffc885c5cbd441eb79e3 Mon Sep 17 00:00:00 2001 From: cnavacch Date: Thu, 2 Feb 2023 01:31:13 +0100 Subject: [PATCH] introduced rounding of spatial coordinates when accessing an xarray dataset --- CHANGELOG.rst | 5 +++++ src/veranda/raster/mosaic/netcdf.py | 4 +++- 2 files changed, 8 insertions(+), 1 deletion(-) diff --git a/CHANGELOG.rst b/CHANGELOG.rst index ca92f29..605b822 100644 --- a/CHANGELOG.rst +++ b/CHANGELOG.rst @@ -2,6 +2,11 @@ Changelog ========= +Version 1.0.1 +============= + +- introduced rounding of spatial coordinates when accessing an xarray dataset + Version 1.0.0 ============= - completely new package and module structure diff --git a/src/veranda/raster/mosaic/netcdf.py b/src/veranda/raster/mosaic/netcdf.py index 3039ecc..68e4bb2 100644 --- a/src/veranda/raster/mosaic/netcdf.py +++ b/src/veranda/raster/mosaic/netcdf.py @@ -10,6 +10,7 @@ from geospade.crs import SpatialRef from geospade.raster import Tile from geospade.raster import MosaicGeometry +from geospade import DECIMALS from veranda.utils import to_list from veranda.raster.native.netcdf import NetCdf4File @@ -777,7 +778,8 @@ def write(self, data, use_mosaic=False, data_variables=None, encoder=None, encod if not src_tile.intersects(data_geom): continue dst_tile = data_geom.slice_by_geom(src_tile, inplace=False) - data_write = data_filt.sel(**{space_dims[0]: dst_tile.y_coords, space_dims[1]: dst_tile.x_coords}) + data_write = data_filt.sel(**{space_dims[0]: np.around(dst_tile.y_coords, decimals=DECIMALS), + space_dims[1]: np.around(dst_tile.x_coords, decimals=DECIMALS)}) else: dst_tile = data_geom src_tile = data_geom