r - 使用 colorRampPalette 合并调色板并使用传单进行绘图

Question

我正在尝试合并两个colorRampPalette要使用的方案，leaflet并且一直在关注这个很好的 例子。该示例运行良好，但我似乎无法让它适用于我的工作，下面是可重现的示例。我正在使用RdYlGn调色板，我希望低于阈值的数字为深绿色，高于阈值的数字为更红色（跳过一些内部颜色）。

对于我的示例，我的截止值是nc$PERIMETER< 1.3，所以我希望低于此值的数字为绿色，高于红色的所有数字（颜色#FDAE61向前）。

library(sf)  
library(leaflet)
library(RColorBrewer)

#palette im using
palette <- rev(brewer.pal(11, "RdYlGn"))
# [1] "#006837" "#1A9850" "#66BD63" "#A6D96A" "#D9EF8B" "#FFFFBF" "#FEE08B" "#FDAE61" "#F46D43" "#D73027" "#A50026"
previewColors(colorNumeric(palette = palette, domain = 0:10), values = 0:10)


# preparing the shapefile
nc <- st_read(system.file("gpkg/nc.gpkg", package="sf"), quiet = TRUE) %>% 
  st_transform(st_crs(4326)) %>% 
  st_cast('POLYGON')
nc

x <- sum(nc$PERIMETER < 1.3)  
x # number of values below threshold = 21


### Create an asymmetric color range
## Make vector of colors for values smaller than 1.3 (21 colors)
rc1 <- colorRampPalette(colors = c("#006837", "#1A9850"), space = "Lab")(x)    #21 

## Make vector of colors for values larger than 1.3 
rc2 <- colorRampPalette(colors = c("#FDAE61", "#A50026"), space = "Lab")(length(nc$PERIMETER) - x)

## Combine the two color palettes
rampcols <- c(rc1, rc2)

mypal <- colorNumeric(palette = rampcols, domain = nc$PERIMETER)
previewColors(colorNumeric(palette = rampcols, domain = NULL), values = 1:length(nc$PERIMETER))

查看预览，它似乎有效（1.3 以下的 21 个值应该是绿色的）：

绘制它：

leaflet() %>%
  addTiles() %>%
  addPolygons(data = nc,
              fillOpacity = 0.7,
              fillColor = ~mypal(PERIMETER),
              popup = paste("PERIMETER: ", nc$PERIMETER) )

绘图正常但没有给出正确的颜色，突出显示的颜色高于阈值（1.3），因此不应该是绿色，但它是：

我认为我创建调色板的方式是错误的，但预览似乎表明我做得对？

有人有想法么？谢谢

标签： rcolorsleafletmapscolor-palette

Answer 1

自从我写了那个答案以来，我有点对这个问题负责。我不知道传单如何为多边形分配颜色。但我认为我们见证了您的方法行不通。根据我之前的想法，我为您做了以下工作。我创建了一个新的连续变量（即ranking）。此信息是 中的值的顺序PERIMETER。PERIMETER这样， （即0.999）的最小值肯定会得到第一个颜色。在我之前的回答中，我建议使用colorFactor()，但这让您很难创建传奇。所以这里有额外的信息。当我创建一个图例时，我使用ranking并colorNumeric()创建了一个调色板，它是mypal2. 我们使用相同的信息来填充多边形并添加图例，但我们使用不同的函数（colorFactor 或 colorNumeric）。一旦我们有了图例，我们就必须改变标签格式。因此我们使用labelFormat(). 我ranking用作索引并在PERIMETER.

library(sf)  
library(leaflet)
library(RColorBrewer)

#palette im using
palette <- rev(brewer.pal(11, "RdYlGn"))
# [1] "#006837" "#1A9850" "#66BD63" "#A6D96A" "#D9EF8B" "#FFFFBF" "#FEE08B" "#FDAE61" "#F46D43" "#D73027" "#A50026"
previewColors(colorNumeric(palette = palette, domain = 0:10), values = 0:10)


# preparing the shapefile
nc2 <- st_read(system.file("gpkg/nc.gpkg", package="sf"), quiet = TRUE) %>% 
       st_transform(st_crs(4326))


# Add sequence information in order to create 108 categories for
# colorFactor(). I sorted the data and added the sequence information.

arrange(nc2, PERIMETER) %>% 
mutate(ranking = 1:n()) -> nc2

x <- sum(nc2$PERIMETER < 1.3)   
x # number of values below threshold = 21


### Create an asymmetric color range
## Make vector of colors for values smaller than 1.3 (21 colors)
rc1 <- colorRampPalette(colors = c("#006837", "#1A9850"), space = "Lab")(x)    #21 

## Make vector of colors for values larger than 1.3 
rc2 <- colorRampPalette(colors = c("#FDAE61", "#A50026"), space = "Lab")(length(nc2$PERIMETER) - x)

## Combine the two color palettes
rampcols <- c(rc1, rc2)

# Create a palette to fill in the polygons
mypal <- colorFactor(palette = rampcols, domain = factor(nc2$ranking))
previewColors(colorNumeric(palette = rampcols, domain = NULL), values = 1:length(nc$PERIMETER))


# Create a palette for a legend with ranking again. But this time with
# colorNumeric()

mypal2 <- colorNumeric(palette = rampcols, domain = nc2$ranking)

leaflet() %>%
addTiles() %>%
addPolygons(data = nc2,
            fillOpacity = 0.7,
            fillColor = ~mypal(nc2$ranking),
            popup = paste("PERIMETER: ", nc2$PERIMETER)) %>% 
addLegend(position = "bottomright", pal = mypal2, values = nc2$ranking,
          title = "PERIMETER",
          opacity = 0.7,
          labFormat = labelFormat(transform = function(x) nc2$PERIMETER[x]))

如果我将阈值级别设置为 2.3（小于 2.3），我会得到这个。