Mais

Comparar e encontrar imprecisões em duas camadas raster

Comparar e encontrar imprecisões em duas camadas raster


Eu tenho duas camadas raster com células urbanas e células rurais. Usei uma função de sobreposição para comparar essas duas camadas raster e ver as diferenças. Antes da sobreposição, dei à célula urbana o valor 1 e à célula rural o valor 0. Após a sobreposição, obtenho o resultado: 0,1,2, NA: Claramente, 2 e 0 significam que as células urbanas e rurais correspondem e 1 significa que não. t. no entanto, apenas olhando para 1 não consigo ver se o erro ocorreu devido a célula urbana ou célula rural. Alguém sabe como posso ver onde está a imprecisão nesta camada raster se eu tomar uma das camadas raster como referência?

classe b: dimensões RasterLayer: 190, 333, 63270 (nrow, ncol, ncell) resolução: 0,008344257, 0,008344257 (x, y) extensão: 12,10615, 14,88479, 51,68836, 53,27377 (xmin, xmax, ymin, ymax) coord. ref. : + proj = longlat + datum = WGS84 + no_defs + ellps = WGS84 + towgs84 = 0,0,0 c classe: RasterLayer dimensões: 190, 333, 63270 (nrow, ncol, ncell) resolução: 0,008344257, 0,008344257 (x, y ) extensão: 12,10615, 14,88479, 51,68836, 53,27377 (xmin, xmax, ymin, ymax) coord. ref. : + proj = longlat + datum = WGS84 + no_defs + ellps = WGS84 + towgs84 = 0,0,0 uc <- overlay (c, b, fun = function (x, y) {return (x + y)}) urx <- overlay (c, b, fun = function (x, y) {return (x * y)}) cellStats (uc, stat = "sum", na.rm = TRUE) exclusivo (getValues ​​(uc)) [1 ] NA 0 1 2

Portanto, você tem dois rasters booleanos contendo as mesmas classes, urbano e rural, e gostaria de saber as classes combinadas, retendo informações tanto para a classe (rural e urbana) quanto para o raster de origem (1 ou 2). Eu sugeriria o uso de números de base 2 (binários) para seus valores de classe iniciais. Reclassifique seus rasters iniciais de forma que:

Raster1 Rural = 1 (00000001)

Raster1 Urbano = 2 (00000010)

Raster2 Rural = 4 (00000100)

Raster2 Urban = 8 (00001000)

Agora, quando você soma as duas imagens, os pixels conterão apenas os seguintes valores possíveis:

Acordo Rural = 5 (00000101)

Acordo Urbano = 10 (00001010)

Raster Rural1 / Raster Urbano2 = 9 (00001001)

Raster Urbano1 / Raster Rural2 = 6 (00000110)

Você pode então reclassificar esses valores para uma sequência mais lógica (por exemplo, 1, 2, 3, 4) conforme necessário. Esse é o truque para representar mais de uma informação no valor numérico único de um pixel. A chave é que os próprios valores numéricos nada mais são do que sinalizadores, conforme indicam sua representação binária.


Assista o vídeo: Mosaic or Combine Multiple Raster Datasets - DEM or Satellite Imagery Using ArcGIS