Open mdsumner opened 2 months ago
this almost works
library(vapour) library(ximage) dsn <- "NetCDF:../../rasterwise/extdata/dt_ref_global_merged_madt_uv_19921014_19921014_20060315.nc:Grid_0001" info <- vapour_raster_info(dsn) vrt <- vapour::vapour_vrt(sprintf("vrt://%s?a_gt=0,1,%i,0,%i,1&a_srs=EPSG:3857", dsn, info$dimension[1], -info$dimension[1])) dsn <- gdal_raster_dsn(vrt) #ximage(gdal_raster_data(dsn[[1]])) #ximage(gdal_raster_data(sprintf("vrt://%s?a_gt=40075017,-37106.5,0,-16957661,0,37065.93", dsn[[1]]))) dsn2 <- gdal_raster_dsn(sprintf("vrt://%s?a_gt=40075017,-37106.5,0,-16957661,0,37065.93", dsn[[1]]))
this almost works