यह एक ऐसी प्रक्रिया है जो जीआईएस सॉफ्टवेयर में बस कुछ सेकंड का समय लेती है। R में इसे करने का मेरा प्रयास स्मृति की एक बड़ी मात्रा का उपयोग करता है, फिर विफल हो जाता है। क्या मेरे कोड में कुछ गड़बड़ है, या यह ऐसा कुछ है जो R नहीं कर सकता है? मैंने पढ़ा है आर ग्रास के अंदर काम कर सकता है, क्या मैं आर के अंदर ग्रास फ़ंक्शन का उपयोग कर सकता हूं?
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)
बस जिज्ञासु, जो पैकेज आप उपयोग कर रहे हैं?
—
djq
@सेलीनियस: इस पैकेज को कहा जाता है
—
जे विन।
raster