如何绘制两个具有不同扩展名的栅格

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 ，0 + units = m + no_defs

     

数据源：C：\ Users \ Usuario \ AppData \ Local \ Temp \ RtmpSYhw1w \ raster \ r_tmp_2018-11-29_153551_4484_16660.grd

     

名称：图层

     

值：1、3（最小，最大）

第二个栅格TS：

  

class：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 ，0 + units = m + no_defs

     

数据源：C：\ Users \ Usuario \ AppData \ Local \ Temp \ RtmpSYhw1w \ raster \ r_tmp_2018-11-29_154600_4484_99649.grd

     

名称：图层

     

值：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)

maps superimposed

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)