Карта Gadm со слоем точек в R

Я уже спрашивал о построении точек в зависимости от региона, но теперь мой вопрос касается наложения точки, заданные координатами по карте gadm. Я хочу показать разные метеостанции, используя данные на основе этого site или вот данные, используемые в коде Я пробовал использовать такой код:

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))

, следуя совету из уже упомянутого вопроса, но он заканчивается точками построения без слоя gadm введите описание изображения здесь , я думаю, это из-за неупоминания регионов

Извините, если вопрос глупый, но буду признателен за любую помощь.


r map cartography
person Ruvin Rafailov    schedule 25.06.2013    source источник
comment
См. этот перекрестный пост на gis.stackexchange.com для возможных решений: gis.stackexchange.com/questions/64448/   -  person fdetsch    schedule 26.06.2013

Ответы (1)


arrow_upward
0
arrow_downward

Вот ответ. Благодаря flowla

Немного отредактировал код, я мог построить другой график в зависимости от высоты.

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")

Сюжет

person Ruvin Rafailov    schedule 26.06.2013