通過在R中使用SpatialDataPoint提取值形式的堆棧柵格

Question

我對Landsat的柵格數據感興趣，我想通過在R中使用SpatialPointData來提取它們的值，然后在SpatialPointData中將提取的值與相關變量繪制在一起，最后，我想將提取的數據與變量一起導出。空間點對象的屬性。 我已經使用了提取函數來執行此操作，但是問題是，提取后每次都會遇到一些錯誤，有時它會起作用，但是它只給我一個數據幀作為矩陣，無法將其與觀察點匹配。

我的劇本

      #raster
      lsat <- stack(b1,b2,b3,b4,b5,b6_1,b6_2,b7)
      #SpatialPoint
      soil_sp=SpatialPoints(cbind(soil.clean2$x,soil.clean2$y))

從堆棧層提取值

     soil_sp$ref<- extract(lsat2, soil_sp)

在每個波段的點數據中繪制拉伸值和觀測值

     plot( soil_sp$ref ~., data=soil_sp)

最后，我想將提取的值與點數據中的變量一起導出到單個數據幀或SpatialPointsDataFrame中。

Answer 1

解決方案是使用SpatialPointDataFrame提取堆棧柵格的值，然后將結果寫入CSV文件或任何您想要的文件。

實現：從SpatialPointData中創建SpatialPointDataFram

  soil_spdf=SpatialPointsDataFrame(coords = soil_sp,data = soil.clean2, 
  proj4string = soil_crs)

提取值。

soil_spdf$ref <- extract(lsat2 , soil_spdf)

將數據以CSV文件格式寫入PC上所需的目錄中。

write.csv(x =soil_spdf,file ="C:/lsat2.csv")