এটি এমন একটি প্রক্রিয়া যা জিআইএস সফ্টওয়্যারটিতে কয়েক সেকেন্ড সময় নেয়। আর-তে এটি করার চেষ্টা আমার প্রচুর পরিমাণে স্মৃতি ব্যবহার করে ব্যর্থ হয়। আমার কোডে কিছু ভুল আছে, নাকি এটি আর কিছু করতে পারে না? আমি পড়েছি আর গ্রাসের ভিতরে আর কাজ করতে পারে, আমি কি আর এর ভিতর থেকে গ্রাস ফাংশন ব্যবহার করতে পারি?
library(raster)
# I have many environmental rasters in this format
new_r <- raster(ncol=615, nrow=626, xmn=-156.2, xmx=-154.8, ymn=18.89, ymx=20.30)
res(new_r) <- 0.00225
projection(new_r) <- "+proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs +towgs84=0,0,0"
R> new_r ### not too big with a few hundred cells per side
class : RasterLayer
dimensions : 627, 622, 1 (nrow, ncol, nlayers)
ncell : 389994
resolution : 0.00225, 0.00225 (x, y)
projection : +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs +towgs84=0,0,0
extent : -156.2, -154.8, 18.89, 20.3 (xmin, xmax, ymin, ymax)
values : none
# I get the DEM at much higher resolution (zipfile is 182Mb)
zipurl <- "ftp://soest.hawaii.edu/coastal/webftp/Hawaii/dem/Hawaii_DEM.zip"
DEMzip <- download.file(zipurl, destfile = "DEMzip")
unzip("DEMzip", exdir = "HIDEM")
HIDEM <- raster("HIDEM/hawaii_dem")
R> HIDEM ### 10m resolution, file is way too big
class : RasterLayer
dimensions : 15067, 13136, 1 (nrow, ncol, nlayers)
ncell : 197920112
resolution : 10, 10 (x, y)
projection : +proj=utm +zone=5 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0
extent : 179066, 310426, 2093087, 2243757 (xmin, xmax, ymin, ymax)
values : HIDEM/hawaii_dem
min value : 0
max value : 4200
# the following line fails (after a long time)
new_HIDEM <- projectRaster(HIDEM, new_r)
শুধু কৌতূহল, আপনি যে প্যাকেজটি ব্যবহার করছেন তা কী?
—
ডিজেকিউ
@ এসেনিয়াস: এই প্যাকেজটি বলা হয়
—
জে উইন।
raster