0

我想创建扇区，然后为它们着色，下面的代码确实可以，但我还想添加一个图例来表示从提供的色调中选择的颜色光谱。它使用 colour.heat 中的颜色和为颜色提供的颜色始终介于 0 和 1 之间，作为其比率。

circs <- function(radii, sectors=4) {
 radii <- sort(radii)
 rads <- seq(0, 2*pi, length=2*length(radii)*sectors)# sample at these radians
do.call(rbind, lapply(radii, function(r)                   # points for drawing circles
data.frame(X=r*cos(rads), Y=r*sin(rads), 
           sector=rep(1:sectors, each=length(rads)/sectors),
           theta=rads, radius=r)))
}

drawCirc <- function(radii, sectors, hues=NULL, densities=NULL, ...) {
polys <- circs(radii, sectors)
 if (missing(hues)) {
 colors <- colorRampPalette(c("green","yellow","red","blue"))(sectors*length(radii))
 } 
else 
colors <- heat.colors(n=sectors*length(radii))
ind=0
plot(polys[,1:2], type="n" ,...)     # blank plot
for (i in seq_along(radii))  {  # add polygons
 for (j in 1:sectors) {
   ind <- ind+1
   color <- colors[ind]
  with(polys[polys$sector==j,],
       if (i == 1) {
         polygon(x=c(0, X[radius==radii[i]], 0), y=c(0, Y[radius==radii[i]], 0), 
                 col=color, density=densities[ind])
       } else
         polygon(x=c(X[radius==radii[i-1]], rev(X[radius==radii[i]])),
                 y=c(Y[radius==radii[i-1]], rev(Y[radius==radii[i]])), 
                 col=color, density=densities[ind]))
  }
 }
}

因此，当我运行以下代码时，它会绘制组成一个圆圈的扇区。

 drawCirc(radii=1:50, sectors=24, hues=ration, main="Ratio by Colors")
 ration = c( x_1,x_2... x_3600), where each x is between 0 and 1.

如何创建一个图例来表示颜色并将其添加到图表中？

我已经尝试过，但不幸的是无法添加图例。

4

1 回答 1

0

您可以在代码中函数的最末尾添加以下行drawCirc() ，就在最后一个大括号关闭之前：

 nColors <-sectors*length(radii)
 nLegend <- 5
 legend("topleft",legend=c(formatC(rev(seq(0,1,1/(nLegend-1))),format="f",digits=2)), fill=c(colors[seq(0,nColors,nColors/nLegend)]))

鉴于您发布的代码的复杂性，我假设您不会有困难根据需要调整和更改图例。

这是我得到的结果：

在此处输入图像描述

希望这可以帮助。

于 2015-06-28T14:57:55.227 回答