R中带有点图层的Gadm地图
我的问题是绘制依赖于区域的点,但现在我的问题是由坐标在gadm地图上给出的叠加点。我想使用基于此的数据显示不同的计量,或者这里是我尝试使用的代码:R中带有点图层的Gadm地图,r,map,cartography,R,Map,Cartography,我的问题是绘制依赖于区域的点,但现在我的问题是由坐标在gadm地图上给出的叠加点。我想使用基于此的数据显示不同的计量,或者这里是我尝试使用的代码: require(ggplot2) library(maptools) library(rgdal) library(RColorBrewer) library(gpclib) library(rgeos) library(PBSmapping) gpclibPermit() rus<-load("C://RUS_adm1.RData") pro
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:3413 +lon_0=105")
gadm.prj <- spTransform(gadm, proj4.str)
rus<-gadm.prj
met <- read.csv2('C:\\meteo.txt')
cds <- data.frame(
Longitude=met$Long,
Latitude=met$Lat)
k<-as.matrix(cds)
popSP <- SpatialPointsDataFrame(coords=k,met["Elevation"], proj4string=proj4.str)
spplot(popSP, sp.layout=list('sp.polygons', gadm.prj))
require(ggplot2)
图书馆(地图工具)
图书馆(rgdal)
图书馆(RColorBrewer)
图书馆(gpclib)
图书馆(rgeos)
图书馆(PBSmapping)
gpclibPermit()
rus。多亏了
编辑了一点代码,我可以根据海拔绘制不同的图
library(rgdal)
library(raster)
rus <- load("C://RUS_adm1.RData")
popSP <- met <- read.csv2("C://meteo.txt")
for (i in c(3, 4))
popSP[, i] <- popSP[, i] / 1000
coordinates(popSP) <- ~ Long + Lat
projection(popSP) <- projection(gadm)
# Reprojection to EPSG:3413 (see http://www.spatialreference.org/ref/epsg/3413/)
proj4.str <- CRS("+proj=stere +lat_0=90 +lat_ts=70 +lon_0=105 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs")
gadm.prj <- spTransform(gadm, proj4.str)
popSP.prj <- spTransform(popSP, proj4.str)
col_no <- as.factor(as.numeric(cut(popSP.prj$Elevation,
c(0,2500,5000,10000,15000,20000,100000))))
levels(col_no) <- c("<2500", "2500-5000","5000-10000", "10000-15000","15000-20000",">20000" )
color = rainbow(6)
popSP.prj$col_no<-col_no
spplot(popSP.prj, "col_no", sp.layout=list('sp.polygons',cols="ID_0", gadm.prj), col.regions=color, main="Meteostations in Russia", key.space="right")
库(rgdal)
图书馆(光栅)
rus可在gis.stackexchange.com上查看此交叉柱,了解可能的解决方案: