使用 Xarray.where 进行屏蔽会引发错误“维度 'y' 不相等”

Question

I am trying to mask a raster using another raster. Both rasters have same dimensions, but when i apply the mask the xarraydata.where throws error

f[0].where(mask,other=np.nan)

ValueError: indexes along dimension 'y' are not equal

Answer 1

I struggled with the same problem for a few hours due to a lack of helpful answers online.

This line:

self.file[var] = self.file[var].where(self.not_rectangle, other=val)

was giving an error similar to yours:

indexes along dimension 'lsmlat' are not equal

I solved it by adding the extension .data in my definition of self.rectangle:

self.rectangle = self._landmask_file.landmask_diff.data

I hope this helps!