R に 2 つのラスター (Landsat slc-off 画像) があります。どちらも一部のデータが欠落していますが、ギャップの位置は完全にオフセットされています。例として、以下の 2 つのラスターを作成 r1しr2ます。
r1 <- raster(system.file("external/test.grd", package="raster"))
r1_mat <- as.matrix(r1)
r1[which(!is.na(as.matrix(r1)))[1:600]] <- NA
par(mfrow=c(3,1))
plot(r1)
r2 <- raster(system.file("external/test.grd", package="raster"))
r2[which(!is.na(as.matrix(r2)))[900:1400]] <- NA
plot(r2)
ただし、2 番目の画像は、異なる大気条件、たとえば、より良い太陽照度で撮影されています。この効果をシミュレートするには (非常に簡単に):
r2 <- r2 + 200
ここで、r1 の欠落データを埋めたいと思います。r2 ピクセルは、自明なギャップに重なっています。
r3 <- cover(r1, r2)
plot(r3)
の結果を見ると、このページの図 2 にも示されplot(r3)ているいわゆる「ストリッピング効果」が明らかです。推奨される解決策の 1 つは、ギャップを埋める前にヒストグラムを照合して 2 つの画像を正規化することです。この手法は、Grass 7 で行われるように、候補画像の累積分布関数に基づいています。
Rでこれをどのように正確に達成できますか? ありがとう！