2015-12-04 14 views
0

計算輪廓我試圖生成,使用R錯誤作爲R

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

但是一個LST(地表溫度)文件的輪廓,我得到的錯誤

錯誤contourLines( LST,20):沒有適當的 'Z' 矩陣指定

據我所知,命令需要xy(座標)和z(在這種情況下是LST值)。但是,我無法使用這種格式輸入數據。我嘗試使用as.matrix將我的文件轉換爲矩陣,但問題仍然保持不變。

回答

0

雖然contourLines沒有爲我工作,但我用rasterToContour,它的工作完美。

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