BUG: interpolate_na function loses crs information 
### HydroMT version checks

- [X] I have checked that this issue has not already been reported.

- [X] I have checked that this bug exists on the latest version of HydroMT.   

- [X] I have checked that this bug exists on the main branch of HydroMT.            


### Reproducible Example

when running an update of setup_lulc in hydromt_wflow using a landuse map with a utm coordinate system, I noticed that the crs was lost in the step interpolate_na of landuse. 

I initially fixed it in hydromt_wflow 

```
    da_crs = da.raster.crs
    da = da.raster.interpolate_na(method="nearest")
    da.raster.set_crs(da_crs)
```

but @hboisgon correctly pointed out that this should be fixed in interpolate_na of hydromt-core. 
it is probably related by the creation of an empty array which then loses the crs.

### Issue Description

see above. not sure if the landuse map is needed to fix the error, but I can provide.

### Installed Versions

sorry not the last version of hydromt-core, but just made sure the code has not changed for interpolate_na and this is the case. 
the version I used to test is still:  0.6.1.dev0

### Additional context

_No response_
diff --git a/hydromt/raster.py b/hydromt/raster.py
index 3115bb6..4218dad 100644
--- a/hydromt/raster.py
+++ b/hydromt/raster.py
@@ -1850,6 +1850,8 @@ class RasterDataArray(XRasterBase):
             data=interp_data,
             attrs=self._obj.attrs,
         )
+        interp_array.raster.set_nodata(self.nodata)
+        interp_array.raster.set_crs(self.crs)
         return interp_array
 
     def to_xyz_tiles(
['hydromt/raster.py']