This PR fixes the bugs in the cases where the bursts are under the anti-meridian crossing.
When testing with some cases (i.e., MS_1_59), the DEM got zero values only where the input DEM is geographic coordinates (EPSG4326). This is because the original code expects the input DEM to have the UTM coordinate, but the Copernicus DEM is the EPSG 4326.
To be able to process even in the EPSG 4326 DEM, the PR replaced the old function (interpolate_gdal) with the modified one (_warp) .
This function came from DSWx-HLS after some modifications. : https://github.com/nasa/PROTEUS/blob/e5b411a91d5e78bc1479ba6d146ce62a29823e49/src/proteus/dswx_hls.py#L3185
This PR fixes the bugs in the cases where the bursts are under the anti-meridian crossing. When testing with some cases (i.e., MS_1_59), the DEM got zero values only where the input DEM is geographic coordinates (EPSG4326). This is because the original code expects the input DEM to have the UTM coordinate, but the Copernicus DEM is the EPSG 4326. To be able to process even in the EPSG 4326 DEM, the PR replaced the old function (interpolate_gdal) with the modified one (_warp) . This function came from DSWx-HLS after some modifications. : https://github.com/nasa/PROTEUS/blob/e5b411a91d5e78bc1479ba6d146ce62a29823e49/src/proteus/dswx_hls.py#L3185