Question

我有两张图片（一张RGB和一张NDVI）。我用R下的光栅加载它们。我想将NDVI的特定部分（NDVI值<0.2）添加到RGB图像中。但是，两幅图像的范围是不同的。

RGB图像：

class       : RasterBrick 
dimensions  : 25223, 19462, 490890026, 4  (nrow, ncol, ncell, nlayers)
resolution  : 0.02973, 0.02973  (x, y)
extent      : 179774.1, 180352.7, 2583827, 2584577  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
data source : /Users/shu-yunchen/Desktop/UVA_related_articals/Images_for_testing/Image/1005_rgb_a7_transparent_mosaic_group1.tif 
names       : X1005_rgb_a7_transparent_mosaic_group1.1, X1005_rgb_a7_transparent_mosaic_group1.2, X1005_rgb_a7_transparent_mosaic_group1.3, X1005_rgb_a7_transparent_mosaic_group1.4 
min values  : 0, 0, 0, 0 
max values  : 255, 255, 255, 255

NDVI图像（NDVI值<0.2提取）：

class       : RasterLayer 
dimensions  : 7311, 6124, 44772564  (nrow, ncol, ncell)
resolution  : 0.08541, 0.08541  (x, y)
extent      : 225357.8, 225880.8, 2585532, 2586156  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=51 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : layer 
values      : 1, 1  (min, max)

以下是我用于将NDVI添加到RGB图像的代码：

plotRGB(RGB, r = 1, g = 2, b = 3, axes = TRUE, stretch = "lin", 
        main = "Landsat False Color Composite")
plot(NDVI,add = TRUE,col = rev(terrain.colors(4)))

它没有将NDVI图像添加到RGB文件中，我认为原因可能是x轴是如此不同，有什么建议吗？如果我应该提供进一步的信息，请告诉我。非常感谢您的帮助。

最佳，陈淑云

匹配两个图像之间的x轴

0 个答案: