如何将数字映射到颜色？

Question

我正在试图将此文件（无符号字符，像素= 1440和行= 720）视为一张地图。我试了一下给出的代码段。首先，我下载的问题可能是这个代码使用了连续的配色方案，即使我有离散数据（这是一个分类方案）。如何将数字映射到颜色？请通缉规模的 例子如下所示：

所有，没有代码可以使用“连续”的颜色方案

conne <- file("C:\\landcover.bin", "rb") 
dfr<- readBin(conne, integer(), size=1, n=720*1440, signed=F) 

y<-matrix((data=dfr), ncol=1440, nrow=720) 
image(y) 

Answer 1

raster包提供了一种使用分类 数据的方法。详情请阅读ratify的帮助页面。

首先，让我们创建与你的数据RasterLayer：

library(raster) 

dfr <- readBin('biome1440s.bin', integer(), size=1, n=720*1440, signed=F) 
r <- raster(nrow=720, ncol=1440) 
r[] <- dfr 

现在我们定义这个RasterLayer与ratify一个因素。您应该 使用您的信息，而不是LETTERS改变其等级：

r <- ratify(r) 
rat <- levels(r)[[1]] 
rat$soil <- LETTERS[1:15] 
levels(r) <- rat 

最后，这绝对RasterLayer可以 的rasterVispackage的levelplot方法来显示。

library(rasterVis) 
myPal <- c('lightblue', terrain.colors(14)) 
## using par.settings 
levelplot(r, par.settings=rasterTheme(region=myPal)) 
## or with col.regions 
levelplot(r, col.regions=myPal) 

编辑：levelplot使用lattice图形，而plot使用 base图形。它们不能一起使用（除非您使用 gridBase包）。但是，可以使用 latticeExtra程序包中的+.trellis和layer函数轻松地覆盖其他 信息。由于wrld_simpl是 SpatialPolygonsDataFrame您可以使用sp包中的sp.polygons函数 对其进行绘制。

library(maptools) ## needed for wrld_simpl 
data(wrld_simpl) ## a SpatialPolygonsDataFrame 
levelplot(r, col.regions=myPal) + 
    layer(sp.polygons(wrld_simpl, lwd=0.5)) 

Answer 2

首先，因为所有的数据是数字和量化。
现在，R中有许多内置颜色模式。有关用于地图等的一些常用调色板的讨论见?image。如果您的地图数据对于不同的土壤类型具有离散的数值范围，则使用小比例数学可以为每个类别应用不同的颜色。

Answer 3

或者你可以使用优秀的ggplot2软件包。这有geom_raster几何，它应该适用于您的数据。请注意，ggplot2要求其数据的格式与矩阵不同，可以使用reshape2包中的melt执行此转换。另请参阅我对an earlier question的回答，它与您拥有很多相似之处。