绘制每个段属性的颜色线

Question

我需要为每个分段组件绘制具有离散颜色的线，并给出以下数据框：

head(profilo1, 10)
   dx TPI      DEM
1   0   5 2159.219
2   2   5 2157.735
3   4   5 2156.269
4   6   5 2154.752
5   8   5 2153.282
6  10   5 2151.925
7  12   5 2150.660
8  14   3 2149.402
9  16   3 2148.277
10 18   2 2147.149

dx＆amp; DEM是x和y，而TPI是颜色属性（从1:10开始）我需要链接到颜色。我希望使用RColorBrewer中的一些颜色表（例如，配对）。 虽然我正在寻找解决方案，但我觉得这很有用

color_paired = brewer.pal(n = 10, "Paired")

plot(profilo1$dx, profilo1$DEM,t='p', 
     xlab = "Profile distance (m)", 
     ylab = "Elevation (m.s.l.m)",
     col = color_paired[profilo1$TPI], cex= .3)

但我更倾向于指出。

Answer 1

您可以先创建一个空白的图，然后使用segements函数。

library(RColorBrewer)
color_paired=brewer.pal(n = 10, "Paired")

plot(profilo1$dx,profilo1$DEM, xlab="Profile distance (m)",
     ylab="Elevation (m.s.l.m)", type="n")
segments(head(profilo1$dx, -1), head(profilo1$DEM, -1), 
    profilo1$dx[-1], profilo1$DEM[-1], 
    col=color_paired[head(profilo1$TPI, -1)])

数据

profilo1 = read.table(text="dx TPI      DEM
1   0   5 2159.219
2   2   5 2157.735
3   4   5 2156.269
4   6   5 2154.752
5   8   5 2153.282
6  10   5 2151.925
7  12   5 2150.660
8  14   3 2149.402
9  16   3 2148.277
10 18   2 2147.149",
header=TRUE)