This would be quite useful to shorten tests, mirroring `raster_equal()` with a tolerance, using `np.allclose` on `.data` and `np.array_equal()` on .mask