非サンプリングポイントの値の推定

沿岸干潟のバイオマスを測定したい。ポリゴン内のポイントにしかアクセスできません。ポリゴン内部のポイントの値に基づいて、ポリゴン外部のポイントの値を推定できる方法はありますか？

set.seed(5)
x <- rnorm(50, -1.841, 0.01)
set.seed(50)
y <- rnorm(50, 55.663, 0.01)
xy <- data.frame(x,y, values=rnorm(50))
coordinates(xy) <- c("x", "y")
proj4string(xy) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
plot(xy)

makePolygons <- function(coordsx, coordsy){

  coords <- matrix(c(c(coordsx, coordsy)), ncol=2)
  p <- Polygon(coords)
  p <- Polygons(list(p), ID = "p")
  myPoly <- SpatialPolygons(list(p))
  spdf = SpatialPolygonsDataFrame(myPoly, data.frame(variable1 = c(2),
                                                     variable2 = c(3), row.names = c("p")))
  proj4string(spdf) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")

  print("polygon is in longlat!!!")

  spdf

}

myPoly <- makePolygons(coordsx=c(-1.841960, -1.843464, -1.888623, -1.841960), 
                      coordsy=c(55.633696, 55.68178, 55.63841, 55.633696))

plot(myPoly, add=T)

ここに画像の説明を入力してください

spatial-statistics interpolation sampling

— ルチアーノ
ソース

植生は潮間帯の干潟にあります-いくつかのエリアは到達するのに危険すぎます。植生バイオマスは予測変数

— ルチアーノ2013年

アクセスできないこととバイオマスが関連している可能性が非常に高いため、重大な問題があります。そのため、アクセス可能な場所から取得したデータをすべての場所に外挿することは無効です。有効なアプローチでは、アクセスできない領域の代表的な部分を測定する方法を（代理方法のみであっても）見つける必要があります。クリギング（および他のほとんどの輪郭作成手順）は問題を美しく描き、ソフトウェアは非常に詳細で信じられないほど間違った結果を喜んで与えます。

— whuber

私のアプローチは、おそらくランドサットデータに基づいて、地上ベースのバイオマス推定値をNDVI値に関連付けることです。危険ゾーン内のNDVIからバイオマスを予測するには、回帰を使用します。

— アーロン

@whuber赤いポリゴン内の領域はアクセス可能ですが、人々が使用することはありません。

— Luciano 2013年

どうしたの？それはあなたの研究の性質やサンプリング手順をどのように変えますか？

— whuber

私の推測では、ポリゴンをポイントグリッドに変換し、データを含むポイントと重なっていない各ポイントの値を推定します。これに関するかなりきちんとしたチュートリアルがあります。

— エリックリノ
ソース