r - 如何在 R 中将 SpatialPointsDataFrame 更改为 SpatialPolygonsDataFrame 以在 ggplot2 中使用它？-6ren

r - 如何在 R 中将 SpatialPointsDataFrame 更改为 SpatialPolygonsDataFrame 以在 ggplot2 中使用它？

转载作者：行者123 更新时间：2023-12-04 11:34:31

我使用 R 使用 ggplot2 形式的大圆圈在世界地图上投影一些数据流。
我还想在我的 map 上投影城市区域:
http://www.naturalearthdata.com/downloads/

然而，这些在 SpatialPointsDataFrame 中。也许我的问题是微不足道的，但我不知道如何将文件更改为 SpatialPolygons。

我的代码如下:

urbanareasin  <- readShapePoly("//....//ne_10m_populated_places//ne_10m_populated_places.shp") 
simp<-gSimplify(urbanareasin, tol=0.1)

urbanareas<-ggplot2:::fortify(simp)

我也试过:

urbanareas<-fortify.SpatialPolygonsDataFrame(simp)

和:

urbanareas<-ggplot2:::fortify.SpatialPolygonsDataFrame(simp)

但它们都不起作用。我必须错过一些东西......
我是 R 的初学者，如果有任何建议，我将不胜感激。

先感谢您!

PS。查找以下数据信息:

str(urbanareasin) # to get info about the object

Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots

..@ data :'data.frame': 7322 obs. of 92 variables:

.. ..$ SCALERANK : int [1:7322] 10 10 10 10 10 10 10 10 10 10 ...

.. ..$ NATSCALE : int [1:7322] 1 1 1 1 1 1 1 1 1 1 ...

.. ..$ LABELRANK : int [1:7322] 8 8 8 8 8 8 8 7 7 7 ...

.. ..$ FEATURECLA: Factor w/ 8 levels "Admin-0 capital",..: 4 4 4 4 4 4 4 4 4 4 ...

.. ..$ NAME : Factor w/ 7069 levels "'Ataq","'s-Hertogenbosch",..: 1453 6358 2017 1135 1973 612 5894 3924 2991 6136 ...

.. ..$ NAMEPAR : Factor w/ 81 levels "Adi Ugri","Alleppey",..: NA NA NA NA NA NA NA NA NA NA ...

.. ..$ NAMEALT : Factor w/ 454 levels "\xdcr\xfcmqi|Wulumqi",..: NA NA NA NA NA NA NA NA NA NA ...

.. ..$ DIFFASCII : int [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ NAMEASCII : Factor w/ 7063 levels "'Ataq","'s-Hertogenbosch",..: 1441 6355 2008 1125 1964 605 5887 3912 2981 6127 ...

.. ..$ ADM0CAP : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ CAPALT : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ CAPIN : Factor w/ 23 levels "Administrative",..: NA NA NA NA NA NA NA NA NA NA ...

.. ..$ WORLDCITY : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MEGACITY : int [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ SOV0NAME : Factor w/ 200 levels "Afghanistan",..: 191 191 191 191 191 177 177 180 180 180 ...

.. ..$ SOV_A3 : Factor w/ 201 levels "AFG","AGO","ALB",..: 189 189 189 189 189 175 175 182 182 182 ...

.. ..$ ADM0NAME : Factor w/ 223 levels "Afghanistan",..: 214 214 214 214 214 200 200 203 203 203 ...

.. ..$ ADM0_A3 : Factor w/ 223 levels "ABW","AFG","AGO",..: 210 210 210 210 210 196 196 203 203 203 ...

.. ..$ ADM1NAME : Factor w/ 2477 levels "?li Bayramli",..: 535 718 1856 431 719 1065 460 1318 1100 2202 ...

.. ..$ ISO_A2 : Factor w/ 225 levels "-99","AD","AE",..: 212 212 212 212 212 197 197 202 202 202 ...

.. ..$ NOTE : Factor w/ 6 levels "1","From 1996 as a summer only station",..: NA NA NA NA NA NA NA NA NA NA ...

.. ..$ LATITUDE : num [1:7322] -34.5 -33.5 -33.1 -34.5 -34.1 ...

.. ..$ LONGITUDE : num [1:7322] -57.8 -56.9 -58.3 -56.3 -56.2 ...

.. ..$ CHANGED : num [1:7322] 4 4 4 4 4 4 4 4 4 4 ...

.. ..$ NAMEDIFF : int [1:7322] 1 1 1 1 1 1 1 1 1 1 ...

.. ..$ DIFFNOTE : Factor w/ 52 levels "Added","Added from GeoNames for UN mega cities.",..: 9 9 9 9 9 9 9 6 9 9 ...

.. ..$ POP_MAX : int [1:7322] 21714 21093 23279 19698 32234 61845 21054 61705 19875 62577 ...

.. ..$ POP_MIN : int [1:7322] 21714 21093 23279 19698 32234 61845 21054 61705 19875 62577 ...

.. ..$ POP_OTHER : int [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ RANK_MAX : int [1:7322] 7 7 7 6 7 8 7 8 6 8 ...

.. ..$ RANK_MIN : int [1:7322] 7 7 7 6 7 8 7 8 6 8 ...

.. ..$ GEONAMEID : num [1:7322] 3443013 3439749 3442568 3443413 3442585 ...

.. ..$ MEGANAME : Factor w/ 588 levels "\xdcr\xfcmqi (Wulumqi)",..: NA NA NA NA NA NA NA NA NA NA ...

.. ..$ LS_NAME : Factor w/ 6559 levels "25 de Mayo","28 de Noviembre",..: NA NA NA NA NA NA NA NA NA NA ...

.. ..$ LS_MATCH : int [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ CHECKME : int [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_POP10 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_POP20 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_POP50 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_POP300: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_POP310: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_NATSCA: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_AREAKM: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_AREAKM: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_AREAMI: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_AREAMI: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_PERKM : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_PERKM : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_PERMI : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_PERMI : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_BBXMIN: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_BBXMIN: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_BBXMAX: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_BBXMAX: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_BBYMIN: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_BBYMIN: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MIN_BBYMAX: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MAX_BBYMAX: num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MEAN_BBXC : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ MEAN_BBYC : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ COMPARE : int [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ GN_ASCII : Factor w/ 5960 levels "'Ayoun el 'Atrous",..: 1236 5382 1692 979 1656 529 5000 NA 2536 5197 ...

.. ..$ FEATURE_CL: Factor w/ 1 level "P": 1 1 1 1 1 1 1 NA 1 1 ...

.. ..$ FEATURE_CO: Factor w/ 10 levels "PPL","PPLA","PPLA2",..: 1 1 1 1 1 1 1 NA 1 1 ...

.. ..$ ADMIN1_COD: num [1:7322] 4 6 12 2 7 5 22 0 31 34 ...

.. ..$ GN_POP : num [1:7322] 21714 21093 23279 19698 32234 ...

.. ..$ ELEVATION : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ GTOPO30 : num [1:7322] 28 134 43 29 74 463 374 0 49 247 ...

.. ..$ TIMEZONE : Factor w/ 319 levels "Africa/Abidjan",..: 121 121 121 121 121 32 32 NA 49 49 ...

.. ..$ GEONAMESNO: Factor w/ 16 levels "Added from GeoNames.",..: 2 2 2 2 2 5 2 13 2 5 ...

.. ..$ UN_FID : int [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ UN_ADM0 : Factor w/ 119 levels "Afghanistan",..: NA NA NA NA NA NA NA NA NA NA ...

.. ..$ UN_LAT : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ UN_LONG : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1950 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1955 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1960 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1965 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1970 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1975 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1980 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1985 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1990 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP1995 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP2000 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP2005 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP2010 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP2015 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP2020 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP2025 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ POP2050 : num [1:7322] 0 0 0 0 0 0 0 0 0 0 ...

.. ..$ CITYALT : Factor w/ 140 levels "Ad Damman","Al Hudaydah",..: NA NA NA NA NA NA NA NA NA NA ...

..@ coords.nrs : num(0)

..@ coords : num [1:7322, 1:2] -57.8 -56.9 -58.3 -56.3 -56.2 ...

.. ..- attr(*, "dimnames")=List of 2

.. .. ..$ : NULL

.. .. ..$ : chr [1:2] "coords.x1" "coords.x2"

..@ bbox : num [1:2, 1:2] -179.6 -90 179.4 82.5

.. ..- attr(*, "dimnames")=List of 2

.. .. ..$ : chr [1:2] "coords.x1" "coords.x2"

.. .. ..$ : chr [1:2] "min" "max"

..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots

.. .. ..@ projargs: chr "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

最佳答案

ggplot不适用于形状文件或光栅图像等空间对象。您需要将它们转换为 data.frames。

如果您正在使用 SpatialPointsDataFrame ，这很简单:

mapdata <- data.frame(yourshapefile)

# now create the map
  ggplot() +  geom_point( data= mapdata, aes(x=long, y=lat), color="red")

如果您正在使用 SpatialPolygonsDataFrame , 你需要 fortify对象，像这样:

yourshapefile_df <- fortify(yourshapefile, region ="id")

# now create the map
  ggplot() + geom_point(data= yourshapefile_df, aes(x=long, y=lat, group=group), color="red")

这个答案来自 ZevRoss webpage ，使用 R 进行空间分析的非常有用的来源.

关于r - 如何在 R 中将 SpatialPointsDataFrame 更改为 SpatialPolygonsDataFrame 以在 ggplot2 中使用它？，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/20320377/

文章推荐： r - 根据组内另一列中的唯一值排列列中的值

文章推荐： ssh - 缺少开始标记的错误是什么意思

文章推荐： javascript - 读取 cookie 数据并重定向

Python ggplot 和 ggplotly
前 R 用户，我曾经通过 ggplotly() 函数广泛地结合 ggplot 和 plot_ly 库来显示数据。刚到 Python 时，我看到 ggplot 库可用，但在与 plotly 的简单组合
r - ggplotly 从 ggplot 中删除图例
ggplotly 使用 ggplot 删除 geom_line 图的图例。见例如以下: library(plotly) g % ggplotly() 关于r - ggplotly 从 gg
r - 设置带有端点的 ggplot 网格线/ggplot 的中断计算
我有一个 ggplot我试图以非常简约的外观制作线图的问题。我已经摆脱了图例，转而使用每行右侧的文本标签。如果标签不是那么长，它可能不会那么明显，但如果网格线停在最大 x 值(在这种情况下，在 201
r - 在一个 ggplot() 中生成多个 ggplot 图形
我想使用相同的 ggplot 代码以我的数据框中的数字为条件生成 8 个不同的数字。通常我会使用 facet_grid，但在这种情况下，我希望最终得到每个单独数字的 pdf。例如，我想要这里的每一行一
r - ggplot : conflict between geom_text and ggplot(fill)
当我在 ggplot 上使用 geom_text 时，与 ggplot 的“填充”选项发生冲突。这是问题的一个明显例子: library(ggplot2) a=ChickWeight str(a)
r - 将 ggplotly 和 ggplot 与拼凑而成？
是否可以结合使用 ggplot ly 和拼凑而成的ggplot？例子这将并排显示两个图 library(ggplot2) library(plotly) library(patchwork) a
r - ggplot、ggplotly、scale_y_连续、ylim 和百分比
我想绘制一个图表，其中 y 轴以百分比表示: p = ggplot(test, aes(x=creation_date, y=value, color=type)) + geom_line(aes
R ggplot，删除 ggsave/ggplot 中的白边
如何去除ggsave中的白边距？我的问题和Remove white space (i.e., margins) ggplot2 in R一模一样。然而，那里的答案对我来说并不理想。我不想对固定但未知
r - 文本层在 ggplot 中工作，但用 ggplotly 删除
我有一个带有一些文本层的条形图，在 ggplot 库中一切正常，但现在我想添加一些与 ggplotly 的交互性，但它无法显示文本层我更新了所有软件包但问题仍然存在 df = read.table(
r - ggplot 到 ggplotly 不适用于自定义的 geom_boxplot 宽度
当我尝试在 ggplot 中为我的箱线图设置自定义宽度时，它工作正常: p=ggplot(iris, aes(x = Species,y=Sepal.Length )) + geom_boxplot(
r - 如何通过从 ggplot 中的不同数据帧映射 aes_string 在 ggplot 中生成图例？
我正在尝试为 ggplot 密度创建一个图例，将一个组与所有组进行比较。使用此示例 - R: Custom Legend for Multiple Layer ggplot - 我可以使用下面的代码成
r - ggplot 在多面图上有一些错误。尝试使用多面 ggplot 协调 y 值
所以我试图在一个多面的 ggplot 上编辑 y 值，因为我在编织时在情节上有几个不准确之处。我对 R 和 R Markdown 很陌生，所以我不太明白为什么，例如，美国的 GDP PPP 在美元金额
python-ggplot - 如何在 Python Ggplot 上格式化 x 轴？
我需要在 python 条形图的 x 轴 ggplot 上格式化日期。我该怎么做？最佳答案使用 scale_x_date() 格式化 x 轴上的日期。 p = ggplot(aes(x='dat
r - 为什么 ggplotly 在 rmarkdown 中不能像 ggplot 一样工作
我想使用 ggplotly因为它的副作用相同ggplot甚至graphics做。我的意思是当我 knitr::knit或 rmarkdown::render我期望的 Rmd 文档 print(obj)
r - 在 Shiny 的应用程序中显示 ggplot 时，如何捕获控制台中出现的 ggplot 警告并显示在应用程序中？
我在下面有一个简单的应用程序，它显示了一个 ggplot。 ggplot 在控制台中生成警告(见底部图片)。我想捕获警告，并将其显示在应用程序的情节下方。这是我的代码: library(shiny)
r - 在 Shiny 的应用程序中缓存基本 ggplot 并允许动态修改图层(与 ggplot 等效的leafletProxy)
如果显示的基本数据集很大(下面的示例工作代码)，则在 Shiny 的应用程序中向/从 ggplot 添加/删除图层可能需要一段时间。问题是: 有没有办法缓存 ggplot(基本图)并添加/删除/修改
r - ggplot 和网格 : Find the relative x and y positions of a point in a ggplot grob
我正在组合 ggplot 的多个绘图，使用网格视口(viewport)，这是必要的(我相信)，因为我想旋转绘图，这在标准 ggplot 中是不可能的，甚至可能是 gridExtra 包。我想在两个图
R中的相对频率直方图，ggplot
我可以使用 lattice 在 R 中绘制相对频率直方图包裹: a <- runif(100) library(lattice) histogram(a) 我想在 ggplot 中获得相同的图形.我试
ggplot geom_area的R堆叠区域顺序
我需要重新安装 R，但我现在遇到了 ggplot 的一个小问题。我确信有一个简单的解决方案，我感谢所有提示! 我经常使用堆叠面积图，通常我通过定义因子水平并以相反的顺序绘制来获得所需的堆叠和图例顺序。
ggplot 中的数据重新排序
新的并且坚持使用ggplot: 我有以下数据: tribe rho preference_watermass 1 Luna2 -1.000 hypolimnic 2 OP10I-A1

行者123

个人简介

我是一名优秀的程序员,十分优秀！

作者热门文章

滴滴打车优惠券免费领取

全站热门文章

首页

博学

6Ren·AI

商城

r - 如何在 R 中将 SpatialPointsDataFrame 更改为 SpatialPolygonsDataFrame 以在 ggplot2 中使用它？