Próbuję wyrównać dwie siatki rastrowe w R. Po wyrównaniu chciałbym móc je dodać.
Próbowałem sprawdzić, czy wykonanie stack
działa:
grid_snap <- stack(habi_sdw, Pop_sdw)
I pojawia się następujący błąd:
Błąd w porównaniu Raster (x): inny zakres
Siatki rastrowe mają następujące właściwości:
show(habi_sdw)
# class : RasterLayer
# dimensions : 9187, 9717, 89270079 (nrow, ncol, ncell)
# resolution : 0.00892857, 0.00892857 (x, y)
# extent : -28.83706, 57.92186, -36.02464, 46.00214 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs
# data source : C:\Users\di39\AppData\Local\Temp\R_raster_di39\raster_tmp_2015-08-12_172902_12860_17067.grd
# names : layer
# values : 0, 333707.6 (min, max)
show(Pop_sdw)
# class : RasterLayer
# dimensions : 10143, 8858, 89846694 (nrow, ncol, ncell)
# resolution : 0.008333333, 0.008333333 (x, y)
# extent : -17.53524, 56.28143, -46.97893, 37.54607 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
# data source : C:\Users\di39\AppData\Local\Temp\R_raster_di39\raster_tmp_2015-08-12_170421_12860_12760.grd
# names : pop2010ppp
# values : 0, 128925.9 (min, max)
Korzystanie alignExtent()
z pakietu rastrowego wydaje się niewłaściwe.
Czy muszę ponownie próbkować, ponieważ rozdzielczości są nieco inne?
(0,00892857 x 0,00892857) vs (0,008333333 vs 0,008333333)