GIS中包括矢量与栅格两种数据,利用tmap可实现矢量数据的快速专题制图而raster可实现栅格数据的可视化及相关数据处理操作。如,波段运算。
数据下载地址:
需要注册:https://data.cdrc.ac.uk/dataset/introduction-spatial-data-analysis-and-visualisation-r
栅格数据:https://github.com/ourcodingclub/CC-spatial
(一)tmap制图
#读入数据
setwd("F:/test/camden/camdenoa11/Camden_oa11")
Census.Data
#----了解数据基本信息-----
#获取列数
ncol(Census.Data)
#获取行数
nrow(Census.Data)
#获取列名
names(Census.Data)
#mean, median, 25th and 75th quartiles,min, max
summary(Census.Data)
#----双变量-----
symbols(Census.Data$Unemployed,Census.Data$Qualification,
circles = Census.Data$White_British,fg="white",
bg ="purple", inches = 0.2)
# adds a regression line, sets the colourto red
abline(lm(Census.Data$Qualification~Census.Data$Unemployed),
col="red")
# a bubble plot with a dotted regressionline
symbols(Census.Data$Unemployed,Census.Data$Qualification,
circles = Census.Data$White_British,fg="white",
bg ="purple", inches = 0.2,
xlab="% in full time employmented",
ylab="% With a Qualification")+
abline(lm(Census.Data$Qualification~Census.Data$Unemployed),
col="red",lwd=2, lty=2)
#R中GIS制图
#加载包,rgdal- Bindingsfor the Geospatial Data Abstraction Library
#•rgeos- Interface to Geometry Engine -Open Source
library(rgdal)
library(rgeos)
#加载Camden_oa11.shp文件
Output.Areas
plot(Output.Areas)
#将census.data的数据连接到矢量属性表中
OA.Census
#设置坐标系
proj4string(OA.Census)
#快速地图可视化
library(tmap)
library(leaflet)
#分位图
qtm(OA.Census, fill ="Qualification")