保持r/ sf中点顺序的多边形边界st_intersection线性化
保持r/sf中点顺序的多边形边界st_intersection线性化是指将多边形的边界线性化，使得多边形的边界上的点按照一定的顺序排列，以便进行后续的计算和处理。
多边形边界线性化的目的是为了方便进行多边形的空间分析、几何计算和可视化等操作。通过线性化，可以将多边形的边界表示为一系列有序的线段或曲线，从而更容易进行各种几何运算，如求交、求并、求差等。
在保持r/sf中点顺序的多边形边界st_intersection线性化过程中，可以采用以下步骤：
- 确定多边形的边界点：遍历多边形的所有顶点，将顶点按照一定的顺序连接起来，形成多边形的边界线段。
- 线段化边界：将多边形的边界线段进行线段化处理，将曲线段拆分为一系列线段，使得每个线段都是直线段或简单曲线段。
- 线性化边界：对线段化后的边界进行线性化处理，将曲线段转化为一系列直线段，使得多边形的边界可以用一系列有序的线段来表示。
- 保持r/sf中点顺序：在线性化边界的过程中，需要保持多边形边界上的点的顺序不变，即按照原始多边形的边界点顺序进行线性化处理。
保持r/sf中点顺序的多边形边界st_intersection线性化可以应用于许多领域，如地理信息系统（GIS）、计算机图形学、计算机辅助设计（CAD）等。它可以用于多边形的空间分析、路径规划、碰撞检测、可视化等应用场景。
腾讯云提供了一系列与云计算相关的产品，其中包括云服务器、云数据库、云存储、人工智能服务等。具体推荐的产品和产品介绍链接地址可以根据具体需求和场景进行选择。
相关·内容
2回答
我正在尝试使用sf包中的sf_intersection()函数来识别两个不同多边形的公共边界。 我在我的数据中尝试了这个simple approach，它来自一个shapefile，但它并不完全像我预期的那样工作。 我的数据是来自this repository的shapefile名称"zones“，这是我尝试过的： library(sf)
library(ggplot2)
zones <- st_read('./Data/zones.shp')
zones$id <- seq(nrow(zones))
borders <- st_intersect
浏览 29提问于2019-12-19得票数 1
回答已采纳
1回答
我画了一张地图，其中点被表示为多边形。这些点是在靠近南部国家边界的地方发现的。守则是：
library(leaflet)
library(sf)
long <- c( 4.676119175, 4.53172103 , 4.939782877, 5.074127987, 5.072757119)
lat <- c(51.477299959, 51.589766239, 51.624436295, 51.520707997, 51.631483055)
labs <- c("A", "B", "C", &#
浏览 11提问于2022-10-11得票数 0
回答已采纳
如何使用R中的sf包检查点与多边形是否与布尔结果相交？
library(sf)
dub<- st_as_sf(data.frame(lon=72.917,lat=19.112), coords=c("lon","lat"),crs=4326)
a<- st_as_sf(data.frame(lon=72.930,lat=19.112), coords=c("lon","lat"), crs=4326)
buf<- st_buffer(a,0.01)
我试过交叉口的长度，但不起作用
length(st_inters
浏览 1提问于2018-10-23得票数 0
回答已采纳
1回答
我有两个空间数据集。一个数据集包含许多多边形(总计超过150 K)，指定不同的特征，如河流、植被。其他数据集包含的指定不同区域的多边形(500)要少得多。我需要交叉这两个数据集，以获得不同区域的特性。我可以根据不同的特性对第一个数据集进行子集。如果我使用一个小特征的子集(2,500个多边形)，与区域的交集非常快(5分钟)。但是如果我想让一个更大的特征子集(20,000个多边形)感兴趣，计算就会运行很长时间(两个小时后我就结束了)。这甚至不是我需要交叉的最大功能(5万多边形)。
这是我运行的代码片段：
clean_intersect_save = function(geo_feature
浏览 2提问于2019-11-12得票数 3
回答已采纳
1回答
我试图通过使用sf图书馆和传单来找出哪个市政厅区代表了纽约的哪个邮政编码。我有两个区的空间多边形和邮政编码。我的问题是，每当我试图找到多边形之间的交点时，我就会得到额外的邮政编码，这些邮政编码与地区边界重叠。该地区没有一部分属于边界以外的邮政编码，但使用空间交叉口包括那些与邮政编码相邻的区域。
是否有一个空间操作函数只返回多边形的交集，而不是多边形只共享边界的地方？
以下是我的观点(见邮政编码(灰色)和议会区(黑色)重叠的边框)：
我的密码是：
library(sf)
#NYC Zip Codes
nyc_zips <- read_sf('ZIP_CODE_04011
浏览 1提问于2022-03-30得票数 0
1回答
我试着用sf::st_union在三个多边形之间联合。在下图中显示了ArcGIS "Overlay，Union，All“的结果，我希望通过使用R中的sf包获得与输出中的五个不同多边形相似的结果。
library(sf)
a1 <- st_polygon(list(rbind(c(0, 10), c(45, 10), c(45, 90), c(0, 90), c(0, 10))))
a2 <- st_polygon(list(rbind(c(45, 10), c(90,10), c(90, 90), c(45, 90), c(45, 10))))
b <- s
浏览 0提问于2019-02-15得票数 6
回答已采纳
1回答
我正在使用sf库来绘制聚合数据，并且我想将我的多边形限制为另一组多边形的边界。 下面的图片是我所拥有的 ? 我是sf的新手，所以如果我没有使用正确的术语，我很抱歉，但我想将图像“限制”为灰色边框，同时保留其所有内容。有没有一种干净利落的方法呢？换句话说，丢弃边界之外的一切。 在df_map中定义新多边形还是在ggplot调用中格式化图像本身更容易？我更倾向于后者，但如果是前者，如何做到这一点。 关于如何生成上面的图形，我可以分享的代码如下： library(tidyverse)
library(sf)
library(ggplot2)
ggplot() +
geom_sf(
浏览 16提问于2021-01-11得票数 0
回答已采纳
2回答
我是PostgreSQL和PostGIS (一般情况下的数据库)的绝对初学者，但在R方面有相当好的工作经验。我有两个来自两个不同来源的印度脆弱地区的多多边形数据集--一个是12 2gb左右，它是.gdb格式的(我们称之为mygdb)，另一个是2GB左右的shapefile (让我们称之为myshp)。我希望比较这两组脆弱性映射，并利用映射之间的交(I)、差(D)和并(U)生成一些状态匹配度量。我想利用PostGIS功能(通过R)既不使用R (崩溃！)qgis (太慢)对此也是有效的。首先，我已经将这两个数据集上载到我的PostGIS数据库中。我用ogr2ogr在R上传mygdb。但我在这点上有
浏览 10提问于2022-02-01得票数 0
我正在使用包装传单和sf制作一个voronoi地图，如下所示：
library(leaflet)
library(sf)
library(rnaturalearth)
library(rnaturalearthdata)
long <- c(4.35556 , 5.83745, 4.63683 , 6.06389, 6.41111, 5.639722)
lat <- c(52.00667, 53.09456, 52.38084 , 52.475 , 52.15917, 53.440278)
labs <- c("Delft" , &
浏览 10提问于2022-09-15得票数 0
回答已采纳
1回答
我已经构造了一些缓冲多边形，我想计算如下：
与缓冲区相交的多边形缓冲面积负面积
为了更好地理解我的问题，我附上了一张图片。在图片中，您可以看到带有缓冲区的红色多边形。在这些缓冲区中有其他多边形或其他多边形的一部分。我现在需要计算红色多边形缓冲区内其他多边形的面积。是否也可以只计算部分相交多边形(因为一些多边形只是部分与缓冲区相交)？
我使用的软件包sf，sp和rgdal，我考虑使用intersect命令？但我完全不知道如何开始这个问题，因为我是一个新的蜜蜂R一般(否则我可以提供我的一些失败的代码)。以下是有关我的数据的一些信息：
class(polygons)
[1] "Spatia
浏览 3提问于2021-08-04得票数 0
回答已采纳
1回答
我是{tmap}包的新手，我在创建适当的路线图方面遇到了困难。下面是我的代码(不能立即运行，因为我在磁盘上有一个导出的OpenStreetMap .osm文件)。 # Load OpenStreetMap (OSM) polygon data of the Holden area from disk
holden_bbox_shapes = sf::read_sf('./localdata/holden.osm', 'multipolygons')
# Load OSM line data from disk. These are almost all ro
浏览 9提问于2020-12-23得票数 1
回答已采纳
我正在用tmap包绘制几个多边形的边界。下面的代码是一个基本示例。
library(sf)
library(tmap)
nc <- st_read(system.file("shape/nc.shp", package = "sf"))
tm_shape(nc) +
tm_borders()
这个看起来不错。然而，如果我改变了边界线的样式，多边形之间的边界看起来与轮廓不同。下面是一个例子。我把线的类型改为虚线。一些线段看起来很坚固，或者有很多点。
tm_shape(nc) +
tm_borders(lwd = 1, lty =
浏览 1提问于2019-12-04得票数 3
回答已采纳
我有个棘手的问题。
我正试图将一些数据可视化为一种“漂亮”的元数据浏览器。它是以下格式的基本点数据：
> print(tempdata[1:5, ])
Station Lat_dec Long_dec Surface_T
1 247 50.33445 -2.240283 15.19
2 245 50.58483 -2.535217 14.11
3 239 50.16883 -2.509250 15.41
4 225 50.32848 -2.765967 15.34
5 229 50.63900 -2.9648
浏览 0提问于2019-03-18得票数 2
回答已采纳
我有两个sf对象，每个对象包含3个多边形(圆圈)。每个对象中的每个多边形都有一个唯一的一对一的标识符(“名称”)。我想交叉这两个对象，以便它返回一个相同长度的对象，其中每个几何学都是两个具有相同名称的多边形的交集。例如，产生的几何A将是从第一个sf对象到第二个对象的几何A的交集。这在QGIS中很容易实现，但我正在尝试自动化该工作流。
下面是在一个简单的xy空间中创建三个圆圈的代码。我如何才能像描述的那样将pts_a与pts_b相交？
我试过st_intersects和aggregate
(我刚进入R的几天，所以我真的不知道如何正确地编排这些问题--谢谢你的帮助)
library(sf)
#cr
浏览 1提问于2022-03-03得票数 2
举个例子，我试图给英国Wetwang村的地图涂上颜色，这是基于与该村两座教堂之一的距离。这是没有Voronoi图的基本情节的代码，它显示了村庄和两个教堂的多边形作为点。
library(osmdata)
library(sf)
library(tidyverse)
bb <- getbb("Wetwang", featuretype = "settlement", format_out = "polygon")
wetwang <- getbb("Wetwang", featuretype = "settl
浏览 9提问于2022-10-25得票数 2
回答已采纳
2回答
我在试着测量海岸线的长度。这将用作分析中沿海岸线位置的度量。例如，假设我有一个区域中所有公共海滩的位置数据，我想通过测量它们与海岸上参考点的距离来描述它们在空间上的分布情况。 我跟随this extremely helpful tutorial使用不同长度的尺子计算海岸线的长度。但是，只有当您想要测量整个多边形的长度，即您感兴趣的地理对象是一个岛屿时，它才是准确的。 要获取海岸线的shapefile (请注意，在ne_countries调用中，我故意使用粗略的比例，以使海岸线更平滑，并且只保留返回的第一个形状-- "scalerank“名称并不重要)： library(raster)
浏览 12提问于2019-04-23得票数 1
回答已采纳
我不确定我是否完全理解了创建voronoi多边形的帮助页面。
library(sf)
# function to get polygon from boundary box
bbox_polygon <- function(x) {
bb <- sf::st_bbox(x)
p <- matrix(
c(bb["xmin"], bb["ymin"],
bb["xmin"], bb["ymax"],
bb["xmax"], bb["ymax
浏览 15提问于2017-08-17得票数 12
回答已采纳
我需要将LINESTRING与边界POLYGON功能相交，同时保留生成的POINT功能的顺序。背景是，我需要找出汽车进入/离开特定国家的边境口岸，但边境口岸的顺序很重要。
我实现了以下方法：
# setup test data
poly <-
list(matrix(c(0,0,10,0,10,10,0,10,0,0),ncol=2, byrow=TRUE)) %>%
sf::st_polygon() %>%
sf::st_sfc() %>%
sf::st_sf()
line1 <- matrix(c(-1, 10
浏览 20提问于2019-07-25得票数 0
回答已采纳
2回答
我有一个多边形(zones)和一组坐标(points)。我想为整个多边形创建一个空间核密度栅格，并按区域提取密度之和。多边形以外的点应该被丢弃。
library(raster)
library(tidyverse)
library(sf)
library(spatstat)
library(maptools)
load(url("https://www.dropbox.com/s/iv1s5butsx2v01r/example.RData?dl=1"))
# alternatively, links to gists for each object
# https://gis
浏览 2提问于2020-01-22得票数 4
回答已采纳
我有多边形，我想要计算它们之间重叠的百分比。其思想是，当一个多边形与另一个多边形相交时，可以根据一个多边形或另一个多边形的透视图(即最大值)计算百分比。因此，我想要制作一个脚本，在多边形之间生成百分比覆盖率--从一个多边形到另一个多边形的重叠百分比，并将所有的结果放在一个数据帧中。
下面是我目前的代码：
set.seed(131)
library(sf)
library(mapview)
m = rbind(c(0,0), c(1,0), c(1,1), c(0,1), c(0,0))
p = st_polygon(list(m))
n = 5
l = vector("list"
浏览 11提问于2021-11-17得票数 2
回答已采纳
1回答
给定一组具有全球分布的坐标，我想在世界地图上绘制一个多边形(最好是一个矩形)，然后提取出这些图中的所有坐标。多边形的位置将根据点密度手动放置，但它们的大小必须保持不变。
在这里，我提供了一个在欧洲随机分布的点的最小可重现性示例。我还添加了一个定向图像，希望它有助于理解所需的输出。首先，我想在地图上添加多边形，然后提取该区域内的所有点。谢谢你提前提供任何可能的帮助。
#Load libraries
library(rgdal) #v1.5-28
library(rgeos) #v.0.5-9
library(ggplot2) # 3.3.5
library(rworldmap) #plot wo
浏览 6提问于2022-01-29得票数 1
回答已采纳
2回答
我对地理数据很陌生，所以这也许是个愚蠢的问题。我有两个多角形，我想要做的是测试“选区”中的哪个多边形，每个“选区”的多边形都包含在其中。基本上，看看每个病房的情况。
但是，我对来自sf的函数sf有问题，因为只有当第一个几何图形完全在第二个几何图形内时，它才返回TRUE。当较小的单位触及较大单位的边界时，这是一个问题。下面是一个例子
我有一个多多边形，叫做“called”和一个叫做"wards“。你可以在这里找到他们：
选区里有两个物体：
plot(constituencies$geometry)
然后，我试着和他们一起探讨病房是否在选区之内。
wards <- st_tra
浏览 3提问于2021-06-04得票数 2
我正在研究一些地理数据，看看波士顿街区的边界，并试图找出哪些社区获得了建筑许可。
到目前为止我已经：
在shapefile中阅读并使用readshapePoly将其转换为纬度和经度的数据。
把一个名字和每一个邻里界限联系起来--布莱顿，唐人街等等。
long lat order hole piece id group name 1 -71.12593 42.27201 1 FALSE 1 0 0.1 Roslindale 2 -71.12575 4
浏览 1提问于2019-02-24得票数 2
回答已采纳
1回答
我有一个名为Operation的表，多个点作为geom，我创建了一个100 m的缓冲区，并将其作为一个新列添加到同一个表中。我有另一个表名为住宅，它有多个多边形，目标是找到每个缓冲区内的多边形百分比，并将其作为一个列添加到操作表中。我不知道如何处理这件事。
SELECT AddGeometryColumn ('public','operations','buffer',4326,'POLYGON',2);
UPDATE operations SET buffer = ST_Buffer(geom::geography,100)::
浏览 1提问于2021-06-08得票数 1
回答已采纳
1回答
我有两个多边形层。我想对它们运行st_intersection，给出它们重叠的区域作为新图层的结果。新层应该包含来自两个输入层的属性。我找到了这张图片，它似乎说明了我想要的最终结果。 ? 我的两个输入层都是多边形： SELECT st_geometrytype(geom),
COUNT(*)
FROM a
GROUP BY st_geometrytype(geom)
-- Result is 1368 st_polygons
SELECT st_geometrytype(geom),
COUNT
浏览 40提问于2020-08-17得票数 1
回答已采纳
1回答
我想将共享颜色的颜色更改为不同的颜色，比如红色。到目前为止，我正在策划德国联邦州巴伐利亚和接触奥地利各州。我从那里得到数据-
德国2级-
德国1级-
奥地利2级-
奥地利1级-
library("sf")
library("raster")
library("dplyr")
library("spData")
library("spDataLarge")
library("ggplot2")
library("patchwork")
library(tmap) # fo
浏览 7提问于2021-12-14得票数 3
回答已采纳
提前感谢大家的帮助。
我有两个重叠的复杂多边形。
我需要得到重叠区域的边界。
我有这两个多边形的边界坐标，我想要取回黄色区域的边界坐标。
我试着在另一个中获取$geoWithing坐标，但是坐标没有排序，我不能正确地绘制多边形。
使用$geoIntersects，我只能得到边界交叉点的坐标，这不足以绘制所需的区域。
一定有人以前做过这件事，所以我希望有人能帮上忙。
//Get first zone coordinates within zone 2
var zoneOneArea = db.zone.find(
{
"name":"Zone 1",
浏览 0提问于2015-07-28得票数 1
1回答
我正试图将"LINESTRING“sf格式的2行拆分成一个dataframe格式，在另一个dataframe中用"MULTIPOLYGON”格式的2个圆圈分隔sf格式。
# make data frame with sf points
lndf <- data.frame(
x = c(40, 55, 60, 70),
y = c(5, 20, 30, 35),
attr_data = c(10,10,10,10),
var = c("abc", "abc", "bac", "bac")
)
浏览 1提问于2018-03-05得票数 2
回答已采纳
1回答
我可以读一个rastre和exrec多边形:如果我有这个代码来读取光栅和shapefile：
library(raster)
library(geojsonsf)
library(sf)
library(exactextractr)
r <- raster(matrix(rnorm(10*12), nrow=10), xmn = -180, xmx= 180, ymn = -90, ymx= 90)
myurl <- "http://eric.clst.org/assets/wiki/uploads/Stuff/gz_2010_us_050_00_5
浏览 5提问于2022-04-29得票数 0
1回答
我想限制一个基于geojson中定义的多边形的绘图，以便它只显示带蓝色阴影的区域。
也就是说，只需绘制内部特征并包括环路。
geojson是可用的。
在边缘周围添加一个缓冲区以包括环路也是很棒的。
我绘制所有特征的代码(不受geojson的限制如下所示)。
library(tidyverse)
library(osmdata)
bounding_box <- getbb("Birmingham", featuretype = "city")
streets <- bounding_box %>%
opq()%>%
add_os
浏览 20提问于2020-03-07得票数 1
回答已采纳
1回答
网格设计中的随机地理空间采样
、、、、
我试图创建一个空间调查设计，首先在一个研究地点设置一个网格，然后在每个网格单元中随机抽样。下面是一个虚构的例子，它在蒙大拿州创建了一个六边形网格。问题是st_sample在整个形状上提供随机点，而不是在网格中。我在st_sample中尝试了另一个st_sample，它们不起作用。我怀疑st_sample不认为我的网格是多个多边形？
有什么想法，或者使用这些包/函数或其他在R？
States <- raster::getData("GADM", country = "United States", level = 1)
Montana <- S
浏览 15提问于2022-07-23得票数 2
2回答
R中空间目标面积的密度曲线
、、、、
我试图在R中观察不同形状的纬向和纵向分布，我试图找出一个sf物体在纬度范围内的区域，基本上是横向平缓它的大部分质量分布在哪里。我该怎么做？
我想过使用st_cast()将多个多边形重新投影成多个点，然后在纬度上对这些点进行密度曲线，但这只会使点围绕形状的周长，而不是内部。也许我能找到同一纬度多个点之间的距离，然后把这些距离画成一个伪直方图？
我认为，如果我转换为栅格，并计算每个纬度行中的细胞数，也许会更容易一些？任何想法都值得赞赏。
下面的代码可以得到一个小的shapefile，开始于：
# try with the country of brazil
brazil <- rnatura
浏览 12提问于2022-09-08得票数 0
回答已采纳
我需要从不同形状的文件中提取信息，用于4兆网格--1公顷的单元。目前，我在所有单元格上的每一层都使用st_crop，但这将永远运行。我想通过使用“data.table”(DT)--按坐标裁剪形状文件--来加速这个过程。让我们考虑下面的示例，其中我正在寻找感兴趣区域中多边形边缘的范围：
require(sf)
require(data.table)
require(ggplot2)
require(tidyverse)
# load shapefile
nc = st_read(system.file("shape/nc.shp", package="sf"))
浏览 3提问于2019-11-29得票数 3
回答已采纳
我是新来PostgreSQL的，请耐心听我说。
我有两个表，一个叫"polys“，有几个多边形，另一个叫"box”，只包含一个多边形(我的边界框)。我的查询选择掉在边界框" box“内的"polys”中的所有多边形--如果你愿意的话，这是一个剪辑。两个表都有两列，一列包含它们的ID，另一列包含它们的名为"the_geom“的GeoJSON。
我想要的是一列包含边界框内多边形的ID，另一列包含这些多边形的GeoJSON，称为"the_geom_webmercator"，另一列称为"polygonarea“，包含每个多边形的面积，然后
浏览 4提问于2012-06-06得票数 0
回答已采纳
1回答
我有一个excel文件与加油站和另一个形状文件，有多边形(县)和相应的人口每个县。
我使用st_as_sf命令将excel文件转换为地理空间文件。县形文件包含市所有县的多边形，每个县有一列表示县内的人口。
我能够在加油站周围创建一个缓冲区，并使用地图视图将县、缓冲区和加油站划分为一起。
我的问题是，离这个加油站有多少人住在x米以内(我的缓冲区)。在图中，蓝色区域是加油站周围的缓冲区，下面的一层是用“人口”颜色表示的各县。搜索的重点是计算多边形中的点数，但我的查询有点不同。
我基本上是想从我的县形状文件中的人口列中提取，并且说400 ppl生活在加油站1等的缓冲区中，etc..for所有的加油站
浏览 6提问于2022-03-08得票数 0
1回答
我遇到一个计算多边形和多多边形之间的交集多边形的问题，如下所示：
st_intersection(poly1,poly2)
此命令获取错误：
Error in `[[<-.data.frame`(`*tmp*`, attr(x, "sf_column"), value = list( :
replacement has 2 rows, data has 1
poly1看起来是这样的：
Simple feature collection with 1 feature and 0 fields
Geometry type: MULTIPOLYGON
Dimension:
浏览 15提问于2022-06-29得票数 0
回答已采纳
我正在学习用sf来处理空间数据。
当我试图生成以下情节时，会出现一个错误：
library(sf)
library(ggplot2)
library(rnaturalearth)
library(rnaturalearthdata)
world <- ne_countries(scale = "medium", returnclass = "sf")
ggplot(data = world) +
geom_sf() +
coord_sf(xlim = c(-102.15, -74.12), ylim = c(0, 33.97), expan
浏览 1提问于2022-01-15得票数 3
我有两个多边形和一个等间距的网格。我的目标是让每个单元计算每个单元中多边形的总面积。
library(sf)
library(dplyr)
# Create two polygons of equal size area
pol_A = st_polygon(list(rbind(c(0,0), c(1,1), c(0,1), c(0,0))))
polys = st_sfc(pol_A, pol_A + c(-0.97, .05)) %>%
st_as_sf() %>%
mutate(poly_ID = c("A","B"))
# C
浏览 0提问于2020-05-18得票数 0
我试图‘剪辑’voronoi多边形(与sf-包创建)到一个边界框.但它给我带来了一个我无法定义的错误。我对R的空间世界不是很有经验。
所有的帮助都是感激的。
样本数据
stations <- structure(list(STN = c(209L, 210L, 215L, 225L, 235L, 240L, 242L,
248L, 249L, 251L, 257L, 258L, 260L, 267L, 269L, 270L, 273L, 275L,
277L, 278L, 279L, 280L, 283L, 285L, 286L, 290L, 308L, 310L, 311L,
3
浏览 0提问于2018-09-24得票数 3
回答已采纳
1回答
我想识别野火烟雾多边形(shapefile)下的诊所(csv文件与经纬度信息和诊所ID)，这是从2005年至2019年的日常水平。
我可以用一个带有R代码的野火烟雾多边形来识别诊所。
hms191030 <-readOGR("/Users/jin/Documents/HMS/hms_smoke20191030.shp")
clinics <-read_csv('/Users/jin/Documents/ESRD/Clinics_cont.csv')
clinicsW <-subset(clinics, PHYS_STATE_NM=="
浏览 15提问于2022-11-11得票数 0
回答已采纳
我有一个名为p_1，p_2，.，p_n的多边形(和多个多边形)的列表，我想得到它们都相交的区域。由于st_intersection()不接受列表作为参数，所以我尝试了以下三种方法。他们没有提供一个令人满意的解决方案，这就是为什么我正在寻找替代的，更有效的技术。
(i)我可以循环浏览列表
for(i in P) p_1 <- st_intersection(p_1, i)
其中P是一个包含多边形的列表，p_2到p_n，但这是相当慢的。
(ii) do.call()方法，即
p <- do.call(st_intersection, P)
其中P是从p_1到p_n的多边形列表，其中只计
浏览 2提问于2020-09-11得票数 2
回答已采纳
我试图在另一个多边形中找到多边形的百分比，并希望将其写入现有的表中。我计算了面积，但不能把它加到现有的表格中。
SELECT
ST_Intersection(ST_MakeValid(r.geom),o.buffer) AS intersection,
ST_Area(ST_Intersection(ST_MakeValid(r.geom),o.buffer))/ST_Area(r.geom)*100
FROM public."Residential" r, public."Operation" o
WHERE ST_Intersects(o
浏览 10提问于2021-06-21得票数 1
回答已采纳
1回答
我有一个数据have (addresses)，它有6个地点的纬度/经度信息。我还添加了(为了说明)这个位置实际上属于哪个国家(我在我的真实数据中没有这个)。然后，我使用addresses_sf函数将数据转换为sf对象( st_as_sf )。
然后通过rnaturalearth包获取加拿大边界的sf对象，并将其保存为can。
我想确定address_sf中哪些条目在加拿大，所以我使用sf包中的st_intersection函数。
不幸的是，当它运行时，它只识别出第二个条目(45.41126，-75.70591)位于加拿大边界，并且无法识别结果中的其余3个条目。
我在下面附上了一些可重复的代码，
浏览 2提问于2021-07-19得票数 2
回答已采纳
1回答
我在地图上有一组空间坐标，还有一个贯穿它们的多行字符串。我需要计算出线条在每种颜色的多边形中所花费的长度。
注意事项：
我的最小表达式不是很好，我的实际数据是带有颜色列(在实际数据中为组)的data.frame，以及一些标识列和一列多角形几何图形。我真的不知道怎么做，所以我的reprex只有3个单独的多边形。
有时两种颜色的多边形重叠，在这种情况下，我希望区分的长度，线相交的形状，与当这条线只相交一个形状。
最小折旧率：
poly1 <-
# create list of matrices and the first point same as last point
浏览 2提问于2021-09-17得票数 1
我想从数据框中的一个变量计算多边形的两两相交。 这就是我要找的： set.seed(131)
library(sf)
m = rbind(c(0,0), c(1,0), c(1,1), c(0,1), c(0,0))
p = st_polygon(list(m))
n = 5
l = vector("list", n)
for (i in 1:n)
l[[i]] = p + 2 * runif(2)
s = st_sfc(l)
s.f = st_sf(s)
s.f$id = c(1,1,2,2,3)
plot(s.f.2, col = s.f.2$id)
s.
浏览 14提问于2021-11-15得票数 2
回答已采纳
2回答
按组查找多个多边形的交点
、、、
因此，使用dplyr中的st_union和group_by()可以很容易地找到两个以上多边形的并集。由于st_intersection只适用于查找两个对象的交集，因此使用st_intersection查找交集并不那么简单。按组查找多个多边形的交点的最佳方法是什么？以下是一些示例数据 s1 <- rbind(c(1, 1), c(10, 1), c(10, 10), c(1, 10), c(1, 1))
s2 <- s1 + 4
s3 <- s1 - 4
group <- c(rep(1, 3), rep(2, 3))
df <- data.frame(
浏览 42提问于2021-10-19得票数 2
我有一个数据，显示在格拉斯哥的公共汽车站的位置和另一个数据，显示格拉斯哥的数据多边形。我正在使用的sf包，并已使两个数据空间。我想要做一个空间连接，以创建一个新的数据create (joined_ds)，以匹配每个公交站位置与数据交换多边形及其相关特征(剥夺分数)。我使用的是st_intersection，它为我提供了一个新的数据格式，它包含所有正确的列，但有0条观察。
joined_ds <- st_intersection(st_buffer(bus_stop_data,0), st_buffer(datazones,0))
这两个数据集都使用适当的CRS (英国国家电网的EPSG:
浏览 4提问于2021-08-17得票数 0
回答已采纳
1回答
我想知道在下面的rstats示例中，计算底层网格中多边形形状文件的(小数)区域覆盖率是否更简单。到目前为止，我还没有找到一个为我做这项工作的函数。以下是我的例子：
library(sf)
## create a multipolygon
p1 <- rbind(c(7, 48), c(8, 52), c(11, 49), c(9, 48), c(7, 48))
p2 <- rbind(c(8, 50), c(8, 51), c(9, 50), c(8, 50))
mpol <- st_sfc(st_multipolygon(list(list(p1, p2))), crs=
浏览 0提问于2018-01-25得票数 3
1回答
我正在尝试执行空间表中所有多边形之间的递归交集，并获得结果(多个)多边形以及每个多边形的每个交集的信息。
一个图像(不是真正的规模)来解释它：
假设表中有A, B, C正方形。我希望输出中有A, B, C, A+B, A+C, B+C, A+B+C多边形，并且我需要知道A+B是A和B的交集，依此类推。
到目前为止，我有一个执行交集的查询，但它不会“切断”原始多边形的相交部分。例如：
Polygon A should be A - (A+B) - (A+C) - (A+B+C)
Polygon A+C should be A+C - (A+B+C)
这是我现在得到的A和A+
浏览 0提问于2013-03-26得票数 12
回答已采纳
我是SF和stack的新手，希望我的问题足够清楚。我能够创建一组线，将一个点连接到全美各地的一组点。我可以把美国的各个县读成多个多边形。我的目标是找到并对我创建的跨越县边界的线的所有点进行地理定位。 到目前为止，我可以从这些点创建直线： points_to_lines <- dt %>%
st_as_sf(coords = c("lon", "lat"), crs = 4326) %>%
group_by(lineid) %>%
summarize(do_union = FALSE) %>% lineid
st
浏览 20提问于2019-03-12得票数 4
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
