如何使用spplot函数和RColorBrewer填充区域
spplot函数是R语言中用于绘制空间数据的函数，而RColorBrewer是一个R包，用于生成颜色调色板。通过结合这两个工具，可以实现对区域进行填充。
使用spplot函数和RColorBrewer填充区域的步骤如下：
- 安装和加载所需的R包：install.packages("sp") install.packages("RColorBrewer") library(sp) library(RColorBrewer)
- 准备空间数据：
假设我们有一个空间数据对象
spatial_data，其中包含了要绘制的区域的几何信息和属性数据。 - 选择颜色调色板：
使用RColorBrewer包中的
brewer.pal函数选择一个合适的颜色调色板。该函数的第一个参数是调色板的名称，第二个参数是调色板中颜色的数量。例如，选择一个包含5种颜色的调色板：palette <- brewer.pal(5, "Set3") - 绘制填充区域：
使用spplot函数绘制填充区域。该函数的第一个参数是空间数据对象，第二个参数是要绘制的属性字段，第三个参数是调色板。例如，绘制属性字段
value的填充区域：spplot(spatial_data, "value", col.regions = palette)
完整的代码示例：
# 安装和加载所需的R包
install.packages("sp")
install.packages("RColorBrewer")
library(sp)
library(RColorBrewer)
# 准备空间数据
spatial_data <- readShapePoly("path/to/shapefile.shp")
# 选择颜色调色板
palette <- brewer.pal(5, "Set3")
# 绘制填充区域
spplot(spatial_data, "value", col.regions = palette)以上代码中的"path/to/shapefile.shp"需要替换为实际的shapefile文件路径。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云地理信息服务：https://cloud.tencent.com/product/gis
- 腾讯云数据万象：https://cloud.tencent.com/product/ci
- 腾讯云云服务器：https://cloud.tencent.com/product/cvm
- 腾讯云云数据库MySQL版：https://cloud.tencent.com/product/cdb_mysql
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云物联网平台：https://cloud.tencent.com/product/iotexplorer
- 腾讯云移动开发：https://cloud.tencent.com/product/mobdev
- 腾讯云对象存储：https://cloud.tencent.com/product/cos
- 腾讯云区块链服务：https://cloud.tencent.com/product/baas
- 腾讯云游戏多媒体引擎：https://cloud.tencent.com/product/gme
- 腾讯云视频直播：https://cloud.tencent.com/product/lvb
- 腾讯云音视频处理：https://cloud.tencent.com/product/mps
- 腾讯云元宇宙：https://cloud.tencent.com/product/mu
相关·内容
我有一个数据文件，有两列，Column1-状态；Column2-饱和度(%)。使用包sp和RColorBrewer；制作地图，我下载gadm - level1印度。但是我不能把百分比放在印度地区，我该怎么做呢，请检查我的R代码代码如下所示
ind1 = readRDS("IND_adm1.rds")
spplot(ind1, "NAME_1", scales=list(draw=T), colorkey=F, main="India")
ind1$NAME_1 = as.factor(ind1$NAME_1)
ind1$adhaar = runi
浏览 17提问于2018-01-05得票数 0
我想合并gadm数据中的一些区域，然后绘制地图。到目前为止，我有以下几点：
#install.packages("sp",dependencies=TRUE)
#install.packages("RColorBrewer",dependencies=TRUE)
#install.packages("maptools",dependencies=TRUE)
library(sp)
library(maptools)
#library(RColorBrewer)
# get spatial data
con <- url("http
浏览 7提问于2013-03-04得票数 0
回答已采纳
2回答
我试图用GADM数据绘制印度各州/地区的详细地图，并使用R。我使用了以下代码
# Load required libraries
library(sp)
library(RColorBrewer)
# ---------------------------------------------------------------------------
# load level 2 india data downloaded from http://gadm.org/country
load("IND_adm2.RData")
ind2 = gadm
# plotting dis
浏览 3提问于2015-05-04得票数 4
回答已采纳
1回答
我有一张俄罗斯地图，上面有区域分区
library(raster)
data <- getData('GADM', country='RUS', level=1)
这个链接是一个Word.doc，上面有俄罗斯地区犯罪率的数据(表)。我可以提取这些数据并在R中使用它，我想用2015年，并根据犯罪率在地图上着色区域(还添加了一个传说)。我该怎么做？问题是，区域名称有时在形状文件(NL_NAME_1)和来自www.gks.ru的数据中是不同的。我也有我需要的图形代码，除了这里我们有无意义的颜色：
library(sp)
library(RColorBrewer
浏览 1提问于2018-04-08得票数 1
回答已采纳
2回答
我在R中绘制栅格数据，当我调整到我正在处理的区域时，R显示不想被切断的ylim。
我试过了：
# set lon lat
ylim=c(-4,2)
xlim=c(118,126)
plot(pm10_mean,xlim=xlim, ylim=ylim)
plot(shp, add=TRUE)
plot(shp2, add=TRUE)
但我得到的照片是这样的
如何删除大于2和低于-4的空闲空间？我只想在xlim和ylim中绘制区域
浏览 0提问于2020-04-06得票数 4
我开始使用R的map()函数绘制地图。我注意到，当我调整绘图窗口的大小时，图像不会缩放来填充窗口。根据拖动窗口的大小，如何使地图图像自动调整大小？
我在MacOS上使用RVersion3.0.2。
例如，这里有一张地图，我在这里拖动了一个越来越大的绘图窗口。注意，地图图像的大小没有变化。
library(maps)
map("state")
另一方面，通常的plot()命令做调整图形大小以适应窗口。
plot(1:100, 201:300)
浏览 5提问于2014-02-18得票数 4
回答已采纳
当我'spplot‘一个用readOGR导入的SHP多边形时，补丁的东西总是用颜色填充的。颜色栏显示它使用的是具有单个值的shapefile ID值。
有没有一种方法可以去掉ID值并填充？
谢谢!
浏览 0提问于2013-06-11得票数 0
回答已采纳
1回答
人口点阵图
、、
我想要像在这个中一样绘制一种isarhytmoc点，但是要使用人口数据。然而，不幸的是，我没有足够的技能，只带了这样的代码：
require(sp)
require(rgdal)
require(RColorBrewer)
require(ggplot2)
library(plyr)
library(maptools)
library(rgeos)
require(gpclib)
gpclibPermitStatus()
gpclibPermit()
rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_a
浏览 0提问于2013-06-24得票数 1
回答已采纳
1回答
在R中使用spplot {sp}时，一切工作正常。然而，当我在一个函数中使用它时，spplot似乎什么都不做。 下面是我打算做的一个例子： A)这是可行的： #load packages
require(gplots)
require(sp)
require(profvis)
#create spatial polygon with categories-attribute
grd <- GridTopology(c(1,1), c(1,1), c(10,10))
polys <- as.SpatialPolygons.GridTopology(grd)
myAttribut
浏览 12提问于2019-05-02得票数 0
回答已采纳
1回答
我有一个包含多个序列(T1、T2、T3 )的sf多边形数据帧，它们都在相同的尺度上:它们是不同时间点的观测值。我可以用来描述T1
ggplot(map)+geom_sf(aes(fill=T1))
我想要做的是将所有这三个(T1，T2和T3)作为面(单独的地图)绘制在同一个图形上。我确信有办法做到这一点，但我找不到。有人能告诉我是怎么做到的吗？谢谢!
添加:关于此问题的两个附加说明。首先，上面描述的数据结构可以使用spplot绘制，其中T是spplot的zcol参数的参数。因此，在这方面，我的问题相当于询问如何将spplot结构转换为geom_sf可用的结构。
其次，假设我使用sf读入一个sh
浏览 0提问于2021-02-06得票数 0
1回答
我正在创建一座房子的平面图，以下是我现在所拥有的：
是一种jpg文件，是3 shp文件的扫描版本，该文件显示了房屋、房间和床。
我想了一下如何使用sp.layout的3个shp文件创建一个spplot，现在我还需要在spplot的背面绘制jpg。我该怎么做呢？谢谢。
浏览 1提问于2012-04-04得票数 1
2回答
在使用简单函数图时，我知道如何绘制两个图：
old.par <- par(mfrow=c(1, 2))
plot(faithful, main="Faithful eruptions")
plot(large.islands, main="Islands", ylab="Area")
par(old.par)
这会把某物归还：
我也需要对一个相当复杂的spplot函数做同样的操作。我想要的是3x3平方。
我要绘制9次的函数是：
labelat = fivenum(gwr.res$SDF$Unempl)
labeltext
浏览 2提问于2015-09-03得票数 4
回答已采纳
我是R的新手，我正在尝试使用adm1 (或州)分区来绘制阿根廷蚂蚁物种的分布图。
我已经从GADM网站上下载了数据，我创建了一个csv文件，其中包含了每个adm1中物种是否存在的信息。即使我没有梯度，我还能做一个合唱吗？如果没有，我还可以使用哪些其他类型的地图？
我看过几个网站，包括，和Choropleth Challenge，它们都很有帮助，但它们都有数字数据，我使用的是现在(1)或缺席(0)列。我尝试过的不同包是sp(带有RColorBrewer)、ggplot2、rgeos和maptools。
以下是我到目前为止拥有的代码：
library(sp)
library(RColorBrewe
浏览 0提问于2013-06-04得票数 2
回答已采纳
我用一些数据在R中做了一个空间插值，然后被困在了创建地图上。为此，我必须使用。当我从那里测试代码示例时：
data(meuse)
coordinates(meuse) <- ~x+y
l2 = list("SpatialPolygonsRescale", layout.north.arrow(), offset =
c(181300,329800), scale = 400)
l3 = list("SpatialPolygonsRescale", layout.scale.bar(), offset =
c(180500,329800), sc
浏览 2提问于2015-09-14得票数 1
回答已采纳
1回答
我试图在地块空间数据中使用R中的spplot函数。虽然以我的标准来看，我的情节看上去很好，但我想弄清楚如何改变我的可视化范围，并显示它们的X和Y轴。
我目前正在使用spplot函数，如下所示：
spplot(no2_surf.idw['var1.pred'],main = 'interpolated no2 surface concentrations',ylim<- c(latv[200],latv[40]), xlim <- c(lonv[60],lonv[180]))
使用xlim和ylim来指定输出的范围，但它根本没有改变我绘制的输出。
我
浏览 1提问于2019-05-15得票数 1
回答已采纳
1回答
在sp图的不同面板中绘制等高线
、、、、
我有3个不同的栅格在一个堆栈。我需要绘制一个面板图，并在每个面板上添加不同的shapefile。到目前为止，我设法做到了以下几点；
## read the libraries
library(raster)
library(rgdal)
library(sp)
library(rworldmap)
library(OceanView)
##random raster object
r <- raster(ncol=40, nrow=20)
r[] <- rnorm(n=ncell(r))
# Crea
浏览 4提问于2018-07-28得票数 0
1回答
我想用美国各州的冠状病毒感染者做一张地图。因此，我们的想法是用美国的所有状态可视化一张地图，并在所有地图中看到一系列受感染的人(即500-2000，等等)。这应该用一种颜色的不同色调来表示。深色将是有更多冠状病毒(州)病例的州。
这是我的密码：
install.packages("sp")
library(sp)
install.packages("sf")
library(sf)
install.packages("maptools")
library(maptools)
install.packages("spdep")
浏览 2提问于2020-03-12得票数 0
回答已采纳
1回答
向多个映射图中添加加权点
、、、、
如何使用spplot或ggplot生成以下绘图。这两种方法我都试过，但不幸的是，它们中的任何一种都没有得逞。我面临的问题是与批准的连接，在地图上添加点。重现问题所需的所有必要数据(.shp,.prj,.dbf,.shx)都是可用的
我的代码如下：
library(sp)
library(maptools)
library(RColorBrewer)
germK <- readShapePoly("C:/Users/XYZ/Dropbox/R Skripts/PolygonG/vg250_krs.shp")
germK <- g
浏览 5提问于2015-02-21得票数 0
回答已采纳
1回答
R在地图上绘制示意图
、、、
大家好，任何人都可以帮助我升级我的代码，有可能插入额外的数据到我的地图。这是为我绘制迁移强度地图的代码，我试图添加每个地区的动态信息(许多小饼图)。
绘制地图的步骤
con <- url("http://biogeo.ucdavis.edu/data/gadm2/R/UKR_adm1.RData")
print(load(con))
close(con)
name<-gadm$VARNAME_1
value<-c(4,2,5,2,1,2,4,2,2,4,1,1,1,4,3,3,1,1,3,1,2,4,5,3,4,2,1)
gadm$VARNAME_1<-
浏览 0提问于2015-03-19得票数 0
2回答
我想用spplot + sp.lines (格)代替plot + segments。你知道一个简单的方法来实现这一点吗，比如
library(dismo)
require(rgdal)
require(FNN)
laurus <- gbif("Laurus", "nobilis")
locs <- subset(laurus, !is.na(lat) & !is.na(lon),
select = c("country", "lat", "lon"))
lo
浏览 1提问于2015-04-03得票数 1
回答已采纳
我正在试图绘制美国的热图，并分别给出它的国名和价值。我得到热图的代码是：
library(maps)
library(maptools)
library(sp)
txt <- "AB AK AL AN AR AZ CA CO CT DC DE EN FL GA HI IA ID IL IN KS
1 21 31 1 12 56 316 53 31 16 7 1 335 63 11 42 29 73 40 2"
dat <- stack(read.table(text = txt,
浏览 5提问于2015-01-28得票数 1
回答已采纳
2回答
我想为spplot()上的区域添加名称标签。
示例：
load(url('http://gadm.org/data/rda/FRA_adm0.RData'))
FR <- gadm
FR <- spChFIDs(FR, paste("FR", rownames(FR), sep = "_"))
load(url('http://gadm.org/data/rda/CHE_adm0.RData'))
SW <- gadm
SW <- spChFIDs(SW, paste("SW", rowna
浏览 0提问于2011-02-28得票数 4
回答已采纳
3回答
我已经绘制了一个基本的地图，使用sp()(包sp)。
小数点栏图例在0到1之间的小数位，但我实际上希望它能显示百分比。
如果我对一个简单的图使用Plot()函数，我可以通过手动添加它们来改变比例。是否有一种类似的方法可以在spplot中手动定义单元类型？
示例：
library(sp)
library(lattice)
trellis.par.set(sp.theme()) # sets bpy.colors() ramp
data(meuse)
coordinates(meuse) <- ~x+y
spplot(meuse, "ffreq")
故事情节以1,2,3作为传
浏览 3提问于2015-01-14得票数 2
回答已采纳
1回答
带有面板图的地图
、、、
我正在学习这个很好的例子。一切都正常直到网格绘制，
library("lattice")
spplot(sp_grd, "id", colorkey=FALSE,
panel = function(...) {
panel.gridplot(..., border="black", col.regions="white")
sp.polygons(shp)
sp.points(poi, cex=1.5)
panel.text(..., col=&
浏览 2提问于2015-12-16得票数 4
回答已采纳
我试图使用sp包中的函数spplot绘制多边形数据，但我的数据帧中有一些缺失值(NA)。当我绘制此数据帧时，缺少的值具有透明色。我想把它们画成黑色。我该怎么做呢？
library(sp)
spplot(TestData,12)
这是我的
浏览 0提问于2015-06-17得票数 5
我有一堆数据是我在县一级绘制的，。我想加上一些州的边界。我有一个状态形状文件(多边形)，但是spplot似乎没有任何方法可以添加到之前的地图之上。有没有办法做到这一点，而不用重写面板函数来接受两个SPDF(对于其他人可能遇到的问题，这似乎是非常专门的)？
下面是一个可重复使用的示例：
library(sp)
Srs1 = Polygons(list(Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))), "s1")
Srs2 = Polygons(list(Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))), "
浏览 0提问于2012-10-13得票数 8
回答已采纳
我一直在尝试为spplot创建的地图添加一个标题(通过main)，该地图由两行文本组成，但第二行是斜体的。每行中的文本都是从变量名中提取的。我尝试过不同的引用和表达式，但都没有成功。
下面是我的起始点的一个基本示例，尽管在我的实际使用中，a和b的文本是作为循环的一部分迭代地从向量中提取的。任何关于如何根据示例保持两行分开的建议，但与第二行文本斜体将非常感谢。我希望这个例子是清楚的。提前谢谢。
require(sp)
data(meuse)
coordinates(meuse)<-~x+y
a<-"line 1"
b<-"line 2"
tit
浏览 1提问于2013-05-04得票数 1
回答已采纳
我需要ggplot的帮助。我有一个很大的数据框架，包括不同年份收集的不同物种。对于每个物种和年份，有两个样本，每个样本代表一个长度频率分布。我需要创建一个ggplot函数来为不同的物种、年份和样本绘制相同的图。下面是我多次使用的一个图的代码：
ggplot(df, aes(x=lenClass, y=Ntot, fill=Sex)) +
geom_bar(stat="identity") +
facet_grid(Year~DayTime) + theme_bw() +
xlab('CL mm') + ylab('No.individuals'
浏览 0提问于2016-07-07得票数 0
1回答
我有一个从RasterStack的八个地图的面板，并希望添加一个100公里比例尺的地块。我在spplot函数中使用了sp.layout将SpatialLines对象添加到栅格之上，并假设使用比例尺(或指北针)也可以做到这一点，但我什么也做不到。
我使用的代码是：
library(sp)
library(raster)
library(grid)
spplot(rasterstack, layout=c(2,4),
names.attr=c("a", "b", "c", "d", "e", "
浏览 4提问于2016-08-13得票数 3
我想根据我提供的六角码颜色来更改颜色刻度。
例如，color_scale参数下面的值基于来自RColorBrewer库的可用调色板
library(bupaR)
patients %>%
process_map(type = frequency("relative_case", color_scale = "Purples"))
我是否可以创建一个自定义调色板，它不是RColorBrewer调色板(RColorBrewer::brewer.pal.info)的一部分，以便它可以作为process_map color_scale参数的输入？
浏览 3提问于2021-03-18得票数 1
回答已采纳
1回答
为多个空间图形绘制相同的图例
、、、、
我有一个for循环，它通过sp图进行插值和绘制栅格图。它给出了一个单一的颜色栏，但标签在每个情节中都会发生变化。你可以看到下面的区别。
我想使用相同的颜色栏，相同的标签，但我不能按照相同的地图图例绘制图形。
下面是代码的最后一部分
为(.){
..。WElev.IDW =idw(公式=变量~1，位置= spdf2，newdata = r.pts)
mypath <- file.path("C:“、”.“、粘贴(”WElevMonth“、冒号(变量)、".png”、sep = ""))
png(文件名= mypath)
打印(spplot(WEle
浏览 1提问于2018-09-18得票数 2
1回答
我正在使用spplot命令在晶格中并排绘制2张地图，使用的是简单的默认颜色填充。我使用的语法与此示例相同：
library(sp)
library(rgdal)
library(lattice)
spplot(nc, c("SID74", "SID79"), names.attr = c("1974","1979"),
colorkey=list(space="bottom"), scales = list(draw = TRUE),
main = "SIDS (sudden infant
浏览 0提问于2015-06-11得票数 0
我正在努力理解scale_fill_brewer与scale_fill_manual of package ggplot2之间的联系。
首先，生成一个填充了颜色的ggplot：
library(ggplot2)
p <- ggplot(data = mtcars, aes(x = mpg, y = wt,
group = cyl, fill = factor(cyl))) +
geom_area(position = 'stack')
# apply ready-made palette with scale_fill_brewer from ggp
浏览 7提问于2015-01-14得票数 2
回答已采纳
2回答
我正在使用seg软件包来获得西雅图州人口普查区的种族隔离值。函数spseg()计算四种偏析度量，spplot函数映射每一束在其中心点的偏析结果。
如何获得在每个点上创建的单个域值的向量，比如在spplot输出中？
(下载有关的 )
library(seg)
library(rgdal)
# load data
sea.shp <- readOGR(".","seattle")
# run segregation measure
sea.seg <- spseg(x=sea.shp, data=sea.shp@data[,2:5])
print(s
浏览 0提问于2015-09-21得票数 0
2回答
我想画一张有两个变量的地图。一个是多边形，另一个是点。
多边形-省边界(name1，name2，name3，name4 )，
点类(a，b，c，d)
因为我必须把它们画在一起，所以我需要两个图例。下面是我获取多边形特征图例的代码。请建议如何为点要素添加图例。
#importing shp files
states = readOGR ( dsn = "C:/.../shp", layer = "states")
river = readOGR ( dsn = "C:/.../shp", layer = "Province"
浏览 1提问于2015-03-05得票数 1
1回答
我关于绘制依赖于区域的点，但现在我的问题是关于在gadm地图上覆盖由坐标给出的点。我想展示不同的气象站，使用基于这个的数据，或者这里是我尝试使用这样的代码的：
require(ggplot2)
library(maptools)
library(rgdal)
library(RColorBrewer)
library(gpclib)
library(rgeos)
library(PBSmapping)
gpclibPermit()
rus<-load("C://RUS_adm1.RData")
proj4.str <- CRS("+init=epsg:341
浏览 3提问于2013-06-26得票数 0
1回答
我已经使用Raster软件包上传西班牙地图的各省(2级)，我想填补他们的颜色，根据收入的人均。这是一个文件，上面有每股收入和分配给每个省的号码。请注意，ID_2是包Raster分配的号码，变量PROV是西班牙政府指定的官方编号。
library(raster)
esp<-getData('GADM', country="ESP", level=2)
espPols <- unionSpatialPolygons(esp, esp$ID_2)
renta <- read.table("renta.csv",sep = ";
浏览 1提问于2016-03-13得票数 3
回答已采纳
1回答
我有一个地图，它总结了一个指标的饱和百分比的房地产按邻里在巴黎(房地产的观察价格/法律规定的最高价格)。我想在我的数据集中包含少于5个观测值的社区添加孵化。
我找遍了，但我找不到办法。任何正确方向的建议都会受到欢迎。谢谢。
这是我的代码：
library(sp)
library(sf)
library(rgdal)
library(RColorBrewer)
library(raster)
library(classInt)
library(cartography)
#Importation
setwd("path")
shp <- readOGR(dsn="p
浏览 2提问于2021-08-25得票数 0
我有一个使用以下命令绘制的SpatialPolygonsDataFrame：
spplot(milanoNILNew, "percPolVita",at = c(0,0.0696,0.08979,0.0974,0.1116,0.181),col.regions =pal,main="Quota Clienti Unipol con Polizze Vita ")
我已经使用以下命令创建了一个spatialPointsDataFrame
agenziePoints<-SpatialPointsDataFrame(coords=datiAgenzie[,c(&
浏览 0提问于2014-05-09得票数 6
回答已采纳
1回答
如何设置spplot的默认主题？
、、、、
比方说，我有一些空间数据，我想使用sp包中的函数绘制它：
library('sp')
library('lattice')
demo(meuse, ask = FALSE, echo = FALSE)
spplot(meuse.grid[,'dist'])
很容易更改主题，例如col.regions，只需一个情节：
spplot(meuse.grid[,'dist'], col.regions=rainbow(100))
但是，如果我想为我的所有情节更改col.regions怎么办？我尝试过一些功能，例如：
t
浏览 0提问于2016-01-11得票数 3
回答已采纳
1回答
首先是样本数据：
library(raster)
# download a sample shape file
shape.file <- getData('GADM', country='FRA', level=2)
ID_2 <- rep(shape.file@data$ID_2,times=2)
group.id <- rep(c(100,200),each=96)
set.seed(1)
cat1<- runif(192,min=0,max=100)
set.seed(2)
cat2<- runif(192,min=0,
浏览 4提问于2017-12-15得票数 0
回答已采纳
1回答
编辑
x <- raster::getData('worldclim', var='tmin', res = 10)
s <- stack(x$tmin1,x$tmin12)
sp <- as(s, 'SpatialGridDataFrame')
fra <- raster::getData('GADM',country = 'FRA', level = 2)
my.layer <- list("sp.polygons", fra, col = "black&
浏览 0提问于2018-04-16得票数 0
回答已采纳
我正试图使用R中的spplot函数创建一个映射。
我加载了执行此功能所必需的不同库：
> library(mosaic)
> library(latticeExtra)
我读过两个数据帧。第一个是人口普查，它包含了一个城市33个病房的8个不同的特征。其中之一是汤森贫困指数。第二个区是包含这些病房的多边形数据的SpatialPolygonsDataFrame。
> DISTRICTS=readShapePoly("DISTRICTS.shp")
> CENSUS=read.csv("CENSUS.csv")
我想为不同的病房绘制汤森指数。
浏览 1提问于2018-05-02得票数 1
我有一个这种类型的SPDF：
class : SpatialPolygonsDataFrame
features : 11723
extent : -2294145, -1102155, 2191605, 3506145 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0
vari
浏览 0提问于2017-11-29得票数 0
2回答
我有一个类SpatialPointsDataFrame的对象，它看起来像：
coordinates div score
1 (-86.2905, 32.7131) 1 0.73
2 (-113.17, 34.4462) 2 3.00
3 (-112.769, 35.1567) 4 0.94
4 (-113.193, 33.9167) 5 4.09
5 (-113.265, 34.407) 1 1.50
6 (-113.585, 34.8241) 2 5.98
7 (-113.793, 34.7597)
浏览 2提问于2014-01-13得票数 3
回答已采纳
我使用Venn来可视化values_one和values_two之间的重叠。
默认结果是黑白的，但我想用颜色填充它.基于，我使用RColorBrewer；但是，我得到以下错误消息：
length(x[])，中的错误(area1= area2 = length(x[2]) )：“fill”的意外参数长度
我做错了什么，我怎么解决它？
谢谢你的帮助!
# ------------------
# reproducible values
# ------------------
values_one <- c("Physical Review B", "IFAC
浏览 2提问于2021-01-20得票数 0
回答已采纳
1回答
我试图在R中的位图中添加轴，但它没有显示这些轴。我使用OSX10.10和Rv3.3.0
set.seed(0)
#Sample matrix
bitmap<-matrix(rnorm(150000,mean=1:500),nrow = 300, ncol = 500)
image(bitmap,col = RColorBrewer::brewer.pal(9,"Greys"), axes=FALSE,
useRaster = TRUE)
axis(1,at = seq(from=1000,to = 10000,length.out = 19),
labe
浏览 2提问于2017-04-09得票数 0
回答已采纳
我使用org.rosuda.jri库运行一个R脚本，该脚本通过java实现普通的kirging算法。我使用了ubuntu13.04和版本1.7-3的REngine，虽然所有的结果都完美地生成了，但是我不能创建这些情节并存储它们。
当以下行在R控制台中完美地执行时
png('/home/panisis/Desktop/plots/tralala.png');
spplot(df);
dev.off();
这些都被忽略了
re.eval("png('/home/panisis/Desktop/plots/tralala.png');");
re.e
浏览 4提问于2015-01-09得票数 1
回答已采纳
2回答
我在用R编写程序时遇到了一些问题，通常我会在另一种语言的for循环中这样做……
我现在试图使用sp包(这是一个叶绿体映射)中的spplot绘制一个图，但是我想对一堆变量这样做。在for循环中，我会这样编程：
for (TRAIT in c("height","bmi","whr","body_fat","income")){
pdf("plot_of_TRAIT.pdf")
spplot(MyData, "TRAIT", col.regions = my.palette, cuts
浏览 3提问于2017-11-07得票数 2
回答已采纳
使用了两个数据集：
3列空间数据(x，y，data)
2列(x，y)中的网格数据
automap包autoKrige可以计算克里格，并可以用no、x和y勾标和标签绘制：
plot(kriging_result)
automapPlot(kriging_result$krige_output, "var1.pred", sp.layout = list("sp.points", shimadata), main="OK without grids", xlab="x", ylab="y")
当我使用g
浏览 0提问于2014-08-15得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
云直播活动推荐
腾讯云开发者
