sf | 使用plot函数绘制地图

2024-03-01 15:40

本文主要是介绍sf | 使用plot函数绘制地图,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

前面已经介绍过绘图基础包graphics中的plot函数是一个比较全能的函数,本篇就介绍如何使用plot函数绘制地图。在这里,plot是一个经过sf改装过的函数,但除了新增几个与地图相关的参数外,与原本用法差别并不大。

首先加载sf包,将plot(x, y, ...)中的x参数设置为sf对象,y缺省就可以绘制出一张地图。sf对象主要包括下图几种类型:

如果sf对象是点要素,那么就可以调用plotpar函数中有关点要素的美化参数,如果是线要素亦然。下面主要介绍面要素POLYGON的绘制方法。

先加载工具包和示例数据:

library(sf)
nc = st_read(system.file("gpkg/nc.gpkg", package="sf"), quiet = TRUE)
summary(nc)
# 部分输出结果
summary(nc)
##       AREA          PERIMETER         CNTY_         CNTY_ID            NAME   
##  Min.   :0.0420   Min.   :0.999   Min.   :1825   Min.   :1825   Alamance : 1  
##  1st Qu.:0.0910   1st Qu.:1.324   1st Qu.:1902   1st Qu.:1902   Alexander: 1  
##  Median :0.1205   Median :1.609   Median :1982   Median :1982   Alleghany: 1  
##  Mean   :0.1263   Mean   :1.673   Mean   :1986   Mean   :1986   Anson    : 1  
##  3rd Qu.:0.1542   3rd Qu.:1.859   3rd Qu.:2067   3rd Qu.:2067   Ashe     : 1  
##  Max.   :0.2410   Max.   :3.640   Max.   :2241   Max.   :2241   Avery    : 1  
##                                                                 (Other)  :94  
##       FIPS        FIPSNO         CRESS_ID          BIR74           SID74      
##  37001  : 1   Min.   :37001   Min.   :  1.00   Min.   :  248   Min.   : 0.00  
##  37003  : 1   1st Qu.:37051   1st Qu.: 25.75   1st Qu.: 1077   1st Qu.: 2.00  
##  37005  : 1   Median :37100   Median : 50.50   Median : 2180   Median : 4.00  
##  37007  : 1   Mean   :37100   Mean   : 50.50   Mean   : 3300   Mean   : 6.67  
##  37009  : 1   3rd Qu.:37150   3rd Qu.: 75.25   3rd Qu.: 3936   3rd Qu.: 8.25  
##  37011  : 1   Max.   :37199   Max.   :100.00   Max.   :21588   Max.   :44.00  
##  (Other):94                                                                   
##     NWBIR74           BIR79           SID79          NWBIR79       
##  Min.   :   1.0   Min.   :  319   Min.   : 0.00   Min.   :    3.0  
##  1st Qu.: 190.0   1st Qu.: 1336   1st Qu.: 2.00   1st Qu.:  250.5  
##  Median : 697.5   Median : 2636   Median : 5.00   Median :  874.5  
##  Mean   :1050.8   Mean   : 4224   Mean   : 8.36   Mean   : 1352.8  
##  3rd Qu.:1168.5   3rd Qu.: 4889   3rd Qu.:10.25   3rd Qu.: 1406.8  
##  Max.   :8027.0   Max.   :30757   Max.   :57.00   Max.   :11631.0  
##                                                                    
##             geom    
##  MULTIPOLYGON :100  
##  epsg:4267    :  0  
##  +proj=long...:  0  
##                     
##                     
##                     
## 

nc中有多个属性变量,直接使用plot(nc)可以绘制至多10个变量的分布情况:

plot(nc)

通过调整max.plot参数绘制全部变量的分布情况:

plot(nc, max.plot = 14)

前面也介绍过,使用st_geometry函数可以提取对象的几何信息(sfc格式),再用plot就可以只绘制对象的几何形状,而不进行属性填充:

plot(st_geometry(nc))

绘制单个属性变量的分布情况:

plot(nc["BIR79"])# 或者
library(tidyverse)
plot(select(nc, BIR79))

通过以下参数调整图例:

  • key.pos:位置,1-4分别表示下、左、上、右;

  • key.length:长度;

  • key.width:宽度。

plot(nc["BIR79"], key.pos = 4, key.length = 1, key.width = 0.1)

通过以下参数可以绘制出常见的重分类地图:

  • nbreaks:分类个数;

  • breaks:分类方法,继承自classInt::classIntervals函数中的style参数,包括fixed、sd、pretty、equal、quantile、jenks、kmeans、hclust、bclust、fisher、dpih、headtails等方法,也可通过向量形式手动指定分割点;

plot(nc["BIR79"], nbreaks = 5, breaks = "quantile")

如果绘图变量是因子变量,nbreaksbreaks将失效:

nc$factor <- factor(round(runif(100, 1, 6)))
plot(nc["factor"], nbreaks = 5, breaks = "quantile")
  • 尽管nbreaks参数设定的是5分类,但是由于factor变量是包含6个水平的因子变量,因此绘制出的地图仍然是6分类。

如果不满意默认配色,可以使用pal参数指定颜色组:

library(RColorBrewer)
plot(nc["BIR79"], nbreaks = 5, breaks = "hclust",pal = brewer.pal(5, "OrRd"))

col参数与pal参数的区别:

  • palsf新增的参数,用在重分类时指定色卡;

  • colplot函数自有的参数,是按对象所包含要素的行序依次给要素指定颜色,与重分类无关;

  • 两参数不能同时存在。

nc$n <- 1:100 # 行序号
plot(st_geometry(nc), col = brewer.pal(9, "Blues")[1:2],nbreaks = 5, breaks = "hclust")
text(st_coordinates(st_centroid(nc)), labels = nc$n)
  • 从图中可以看出,行序号为奇数和偶数的分别同色,与重分类结果无关。

通过以下参数实现图层的叠加:

  • reset:后续图层是否重设坐标系、尺寸等绘图属性,默认为TRUE;叠加图层时需要在第一个plot函数中将其设置为FALSE;

  • add:同boxplot等函数的add参数。

nc2 <- dplyr::filter(nc, AREA > 0.2)
pt2 <- st_centroid(nc2)plot(nc["BIR79"], nbreaks = 5, breaks = "jenks",pal = brewer.pal(5, "OrRd"), reset = F)
plot(pt2, add = T, col = "black", pch = 16) # 添加面积超过0.2的县质心
  • 从图中可以看出,基础绘图系统中用于设置点要素的pch参数仍然可以用于sf的点要素,其他参数亦然。

其他新增的参数还有:

  • border:地图边界颜色,设置为NA时隐藏边界;

  • graticule:经纬网;如果要显示经纬度数值,需要将plot的坐标轴显示参数axes设置为TRUE;

  • col_graticule:经纬网颜色;

  • extend:绘图范围;

  • logz:对属性变量对数化。

plot函数原有的参数仍然适用,如使用main参数修改标题:

plot(nc["BIR79"], nbreaks = 5, breaks = "hclust",pal = brewer.pal(5, "OrRd"), border = "blue", graticule = T, axes = T,main = "北卡罗莱纳州")

本篇只介绍plot函数绘制地图的基本用法,后续推文会介绍更精细的调整方法。

plot函数绘制地图的缺点在于不能添加比例尺和指北针。更专业的绘制地图工具包是tmap,后续将会重点介绍。


这篇关于sf | 使用plot函数绘制地图的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/762962

相关文章

鸿蒙中@State的原理使用详解(HarmonyOS 5)

《鸿蒙中@State的原理使用详解(HarmonyOS5)》@State是HarmonyOSArkTS框架中用于管理组件状态的核心装饰器,其核心作用是实现数据驱动UI的响应式编程模式,本文给大家介绍... 目录一、@State在鸿蒙中是做什么的?二、@Spythontate的基本原理1. 依赖关系的收集2.

Python基础语法中defaultdict的使用小结

《Python基础语法中defaultdict的使用小结》Python的defaultdict是collections模块中提供的一种特殊的字典类型,它与普通的字典(dict)有着相似的功能,本文主要... 目录示例1示例2python的defaultdict是collections模块中提供的一种特殊的字

C++ Sort函数使用场景分析

《C++Sort函数使用场景分析》sort函数是algorithm库下的一个函数,sort函数是不稳定的,即大小相同的元素在排序后相对顺序可能发生改变,如果某些场景需要保持相同元素间的相对顺序,可使... 目录C++ Sort函数详解一、sort函数调用的两种方式二、sort函数使用场景三、sort函数排序

C语言函数递归实际应用举例详解

《C语言函数递归实际应用举例详解》程序调用自身的编程技巧称为递归,递归做为一种算法在程序设计语言中广泛应用,:本文主要介绍C语言函数递归实际应用举例的相关资料,文中通过代码介绍的非常详细,需要的朋... 目录前言一、递归的概念与思想二、递归的限制条件 三、递归的实际应用举例(一)求 n 的阶乘(二)顺序打印

Java String字符串的常用使用方法

《JavaString字符串的常用使用方法》String是JDK提供的一个类,是引用类型,并不是基本的数据类型,String用于字符串操作,在之前学习c语言的时候,对于一些字符串,会初始化字符数组表... 目录一、什么是String二、如何定义一个String1. 用双引号定义2. 通过构造函数定义三、St

Pydantic中Optional 和Union类型的使用

《Pydantic中Optional和Union类型的使用》本文主要介绍了Pydantic中Optional和Union类型的使用,这两者在处理可选字段和多类型字段时尤为重要,文中通过示例代码介绍的... 目录简介Optional 类型Union 类型Optional 和 Union 的组合总结简介Pyd

Vue3使用router,params传参为空问题

《Vue3使用router,params传参为空问题》:本文主要介绍Vue3使用router,params传参为空问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐... 目录vue3使用China编程router,params传参为空1.使用query方式传参2.使用 Histo

使用Python自建轻量级的HTTP调试工具

《使用Python自建轻量级的HTTP调试工具》这篇文章主要为大家详细介绍了如何使用Python自建一个轻量级的HTTP调试工具,文中的示例代码讲解详细,感兴趣的小伙伴可以参考一下... 目录一、为什么需要自建工具二、核心功能设计三、技术选型四、分步实现五、进阶优化技巧六、使用示例七、性能对比八、扩展方向建

使用Python实现一键隐藏屏幕并锁定输入

《使用Python实现一键隐藏屏幕并锁定输入》本文主要介绍了使用Python编写一个一键隐藏屏幕并锁定输入的黑科技程序,能够在指定热键触发后立即遮挡屏幕,并禁止一切键盘鼠标输入,这样就再也不用担心自己... 目录1. 概述2. 功能亮点3.代码实现4.使用方法5. 展示效果6. 代码优化与拓展7. 总结1.

使用Python开发一个简单的本地图片服务器

《使用Python开发一个简单的本地图片服务器》本文介绍了如何结合wxPython构建的图形用户界面GUI和Python内建的Web服务器功能,在本地网络中搭建一个私人的,即开即用的网页相册,文中的示... 目录项目目标核心技术栈代码深度解析完整代码工作流程主要功能与优势潜在改进与思考运行结果总结你是否曾经