【lidR】"ttops"の保存方法

ttopsの保存方法がわからなかったので、その備忘録。
これがベストかわからないが。そもそも保存後に元の形式に戻す必要があるのかどうかもわからん。

コード

ttops <- locate_trees(nlas, lmf(ws = 5))　#nlasはnormalization済み

#ttops保存
ttops_new <- sf::st_zm(ttops) #Zの重複を無くす
st_write(ttops_new, "ttops_zdropped.shp") #保存

#ttops復活
ttops_again <- st_read("ttops_zdropped.shp") #上で書き出したファイルを読み込み
ttops_again <- st_zm(ttops_again, drop = FALSE, what = "Z") #もともとのttopsと同じ形式

参考サイト

R lidR write tree points to shapefile

st_zm: Drop or add Z and/or M dimensions from feature geometries in sf: Simple Features for R