vapour icon indicating copy to clipboard operation
vapour copied to clipboard

vapour_warp_raster() is clamping ??

Open mdsumner opened this issue 2 years ago • 0 comments

src <- "/vsicurl/https://gebco2023.s3.valeria.science/gebco_2023_land_cog.tif"
library(terra)
r <- rast(src)
ext <- c(-2165880,  1768374, -1014018,  2262496)

dm <- c(512, 512)
proj <- "+proj=stere +lon_0=80.5196726901466 +lat_0=-67.7777429448851 +lat_ts=71 +datum=WGS84"
t <- rast(ext(ext), ncols = dm[1], nrows= dm[1], crs = proj)
project(r, t, by_util = T)

range(vapour::vapour_warp_raster(src, dimension = dm, extent = ext, projection = proj)[[1]], na.rm = TRUE)
#[1]    1 3730

range(vapour::gdal_raster_data(src, target_dim = dm, target_ext = ext, target_crs = proj)[[1]], na.rm = TRUE)
#[1] -5614  3730

mdsumner avatar Nov 22 '23 00:11 mdsumner