library(dplyr)
cdaspect = read.table("cdaspect.txt", header = T,sep=",")
cdaspect %>%
mutate(ASPECT = cut(RASTERVALU,
breaks = c(0, seq(1, 15, 2), 16) * 22.5,
labels = c(1L, 2L, 4L, 6L, 8L, 7L, 5L, 3L, 9L)
))->cdaspect
cdaspect$ASPECT[which(cdaspect$ASPECT==9)]=1
试了下,这样也是可以的,为cheetahfly兄弟提供的思路点赞!


雷达卡




京公网安备 11010802022788号







