R中计算轮廓时出错

Question

我正在尝试使用R

为LST（地表温度）文件生成轮廓

lst<-raster("test_lst.tif")
cl <- contourLines(lst,20)
shp <- ContourLines2SLDF(cl)

但是，我收到错误

  

contourLines（lst，20）出错：没有指定正确的'z'矩阵

我知道该命令需要x，y（坐标）和z（在这种情况下是LST值）。但是，我无法以这种格式提供数据。我尝试使用as.matrix将文件转换为矩阵，但问题仍然存在。

Answer 1

虽然contourLines对我不起作用，但我使用rasterToContour并且效果很好。

cnt <- rasterToContour(lst) writeOGR(cnt, "lst_contour", layer="cnt", driver="ESRI Shapefile")