Как построить два растра с разными расширениями

Question

Я пытаюсь построить два растра с разными расширениями (из двух разных областей), и они накладываются друг на друга.

Первый растр TN:

класс: RasterLayer

размеры: 1785, 2363, 4217955 (nrow, ncol, ncell)

разрешение: 0.11766, 0.11766 (x, y)

экстент: 474953.5, 475231.5, 6539165, 6539375(xmin, xmax, ymin, ymax)

координ.ссылка: + proj = lcc + lat_1 = 58 + lat_2 = 59.33333333333335 + lat_0 = 57.51755393055556 + lon_0 = 24 + x_0 = 500000 + y_0 = 6375000 + ellps = GRS80 + towgs84 = 0,0,0, -0, -0, -0, 0 + единиц = m + no_defs

источник данных: C: \ Users \ Usuario \ AppData \ Local \ Temp \ RtmpSYhw1w \ raster \ r_tmp_2018-11-29_153551_4484_16660.grd

имена: слой

значения: 1, 3 (мин, макс)

Второй растр TS:

класс: RasterLayer

размеры: 7266, 5237, 38052042 (nrow, ncol, ncell)

разрешение: 0,1141, 0,1141 (x, y)

экстент: 474817, 475414,5, 6537431, 6538260 (xmin, xmax, ymin, ymax)

координ.ссылка: + proj = lcc + lat_1 = 58 + lat_2 = 59.33333333333335 + lat_0 = 57.51755393055556 + lon_0 = 24 + x_0 = 500000 + y_0 = 6375000 + ellps = GRS80 + towgs84 = 0,0,0, -0, -0, -0, 0 + единиц = m + no_defs

источник данных: C: \ Users \ Usuario \ AppData \ Local \ Temp \ RtmpSYhw1w \ raster \ r_tmp_2018-11-29_154600_4484_99649.grd

names: layer

значения: 1, 2 (мин, макс)

library(raster)
ext1 <- extent(6600000, 7000000, 66000000, 70000000) #total area extent

extent(TN) <- ext1
extent(TS) <- ext1 #rasters with the same extent (total area)

em = merge(extent(TN),extent(TS))
plot(em, type="n")
plot(TN,add=TRUE, legend=FALSE)
plot(TS, add=TRUE, legend=FALSE)

карты с наложением

Answer 1

Графики перекрываются только потому, что у вас одинаковые экстенты для каждого растра.Дайте им разные степени.Вот воспроизводимый пример:

library(raster)
logo <- raster(system.file("external/rlogo.grd", package="raster")) 
logo1 <- logo
logo2 <- logo

Назначьте правильные экстенты и объедините:

ext1 <- extent(0, 45, 0, 45)
ext2 <- extent(55, 100, 55, 100)
extent(logo1) <- ext1
extent(logo2) <- ext2
em <- merge(ext1, ext2)

Добавьте данные в график:

plot(em, type = "n")
plot(logo1, add = T, legend = F)
plot(logo2, add = T, legend = F)