R/sf:将多边形边界约束到重叠多边形层
R/sf是一个R语言的包，用于空间数据分析和地理信息系统 (GIS) 相关的操作。它是基于sf包的扩展，提供了更多的功能和操作。
多边形边界约束到重叠多边形层是指将一个多边形的边界限制在与其重叠的其他多边形层上。这个过程可以帮助我们解决一些空间分析和地理数据处理的问题。
具体来说，R/sf可以通过以下步骤实现多边形边界约束到重叠多边形层：
- 导入所需的包：使用library函数导入R/sf包，确保已经安装了sf包和其他所需的依赖包。
- 读取数据：使用st_read函数读取相关的空间数据。可以读取包括多边形层和边界多边形层的数据。
- 重叠检测：使用st_intersects函数检测多边形层与边界多边形层之间的重叠情况。该函数会返回一个逻辑值，表示两个多边形是否有重叠。
- 边界约束：使用st_intersection函数将多边形的边界限制在与其重叠的多边形上。该函数会返回一个新的多边形层，其中包含了被约束的多边形。
- 结果可视化：使用plot函数将结果可视化，以便查看多边形边界约束的效果。
R/sf包的优势在于它提供了强大的地理空间数据处理能力，可以进行空间查询、空间连接、缓冲区分析等常见的GIS操作。它还与其他R语言的数据分析和可视化包完美集成，可以方便地进行数据分析和可视化展示。
在实际应用中，多边形边界约束到重叠多边形层可以用于许多场景，如城市规划中的用地分析、自然资源管理中的生态分区划定等。通过将多边形边界限制在重叠的多边形上，我们可以更好地理解和分析地理空间数据，为决策提供更准确的支持。
对于使用腾讯云的用户，腾讯云提供了一系列与云计算和地理空间数据分析相关的产品和服务，可以帮助用户完成各种任务。具体推荐的产品包括腾讯云地理空间分析、腾讯云地图、腾讯云云服务器等。更多关于腾讯云相关产品的介绍和详情，请参考腾讯云官方网站（https://cloud.tencent.com/）。
相关·内容
我试图使用st_centroid函数在sf库中显示R中GIS多边形特征的标签。不幸的是，虽然head()函数似乎显示每个多边形具有与其相关的不同的x和y坐标，但在地图上的单个点(这显然是一个特定多边形的质心)，所有标签都会被呈现重叠。我在这里做错什么了？
当前代码设置：
library("ggplot2")
library("sf")
sf::sf_use_s2(FALSE) #makes centroids not break
world <- st_read("C:/prgrm/gis/source/10m_land_and_islands
浏览 2提问于2021-11-25得票数 0
回答已采纳
1回答
用sf查找R中无重叠区域的多边形
、、、、
我有两套多边形在形状文件，进口使用R的sf包。我会叫他们A和B。
我想得到B中的所有多边形，其中a) 0%与A或b中的任何多边形重叠，或b) 中的一些与A中的多边形不重叠。
例如，对于B中给定的多边形B，我做而不是想要的区域，其中b是100%包含在两个(或三个或四个.)内。A中的多边形。
但是，假设b包含在两个(或三个或四个.)中的50%。A中的多边形，但其中50%不包含在A中的任何其他多边形中。在这种情况下，我想选择b。或者如果b 100%在A中的所有多边形之外，我想选择b。
是否有可能使用st_intersects进行此操作？请让我知道一个图表或其他东西是否有助于澄清。
======
更明
浏览 3提问于2021-07-12得票数 0
回答已采纳
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回答
有办法将ggplot2映射分割成两个单独的映射吗？我有一张很大的地图，上面有难以辨认的标识。我想把地图垂直分割成两个不同的地图，最好有一个重叠的区域，这样每个多边形至少会在一张地图中显示整体。
下面是一个可重复的例子。我想把地图划分为北纬35度的北地图和北纬35.5度的南面地图(两者重叠在35到35.5之间)。(虽然我意识到用这个例子将另一种方式分开可能更有意义，但我的实际地图是垂直长的。)
library(sf)
library(ggplot2)
sf_nc <- sf::st_read(system.file("shape/nc.shp", package = &#
浏览 2提问于2020-09-11得票数 2
回答已采纳
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
回答已采纳
我正在用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
回答已采纳
2回答
我正在尝试来自：的多边形重叠(POST请求)
下面的内容返回一个空列表：
library(sf)
library(geojsonsf)
library(httr)
# Create a polygon
lon = c(-114.88984, -114.05445)
lat = c(48.99889,49.46062)
poly_df = data.frame(lon, lat)
poly <- poly_df %>%
st_as_sf(coords = c("lon", "lat"),
crs = 4326) %&g
浏览 12提问于2022-10-07得票数 0
回答已采纳
我在使用Geoserver和OpenLayers 2时遇到地图问题
我有一个基本层是带标签的多边形。我有一个河层重叠的基础层。
当多边形层与河层重叠时，如何使多边形层的标签不显示(或显示在顶层)？
这是map 的图像
浏览 11提问于2016-05-18得票数 0
1回答
我有一个典型的数据库，其中包含数百万个多边形作为地块，我想将这些多边形拆分成线，然后删除彼此重叠的线。这些线将纯粹用于在mapnik/和/或geoserver中渲染，因为此时每个地块边界都会渲染两次。
我建议将地块多边形拆分到一个新的表("boundary_lines")中，然后搜索并删除重叠的线。我该如何删除postgis中的这些重叠行？
浏览 0提问于2010-08-05得票数 3
回答已采纳
1回答
当我在没有添加多边形的情况下向我的闪亮的应用程序添加一个简单的地图时，缩放参数可以正常工作，并且地图按它应该的方式初始化。
另一方面，当我用add_polygon()添加多边形时，映射的初始化比它应该的要大。
为什么会发生这种情况？
这是我的密码：
library(mapdeck)
library(sf)
sf = st_as_sf(my_LargeSpatialPolygonDataframe)
output$my_map = renderMapdeck({
mapdeck(token = mytoken, location = c(a, b), zoom = 12,
浏览 1提问于2019-05-16得票数 3
回答已采纳
我有许多栅格，我想检查它们是否完全包含在空间多边形中，完全没有空间多边形，或者与空间多边形交叉(这可能意味着多边形完全在栅格内，或者多边形和栅格重叠)。我这样做是为了在可能的情况下把我从时间密集的伪装中拯救出来。
下面是一个示例：
# create 3 example rasters
r <- raster()
r[] <- rnorm(n = ncell(r))
e1 <- extent(c(45,55,45,50))
r1 <- crop(r,e1)
e2 <- extent(c(20
浏览 3提问于2016-01-08得票数 4
回答已采纳
1回答
这些(k-d树，r-树)中的哪一个适合于搜索和索引多边形。我的用例是，我得到了一些经度较长的点(有效多边形的最小长度为3)，我需要从这些点中找到最小的多边形。我说的最小的意思是，如果一个多边形在另一个多边形中，那么应该返回内部多边形。如果多边形重叠，则不应选择它们。我认为找到位置，然后该地区可能会尝试，但我不确定。
我还想了解一下哪种数据结构更有用。我认为postgis使用了R树索引。
浏览 11提问于2019-10-23得票数 0
1回答
计算无多边形重叠的面积
、、、、
我有个问题。我正在试着计算下面图层的面积(见图)
我使用交集工具找到了4个重叠缓冲区的层与另一个多边形(从栅格转换而来，因此由许多其他多边形组成)之间的交点。这些层现在由200多个多边形组成，并且大多数多边形相互重叠。我实际上想要计算这个层的2D面积，所以我实际上想把这个由多个多边形组成的层转换成一个多边形，这样你就可以计算出这个多边形的面积。因此，我的问题是，有没有可能将这个图层转换为彼此相邻的多边形，并且不再有重叠的多边形，以便我可以计算面积？也许还有另一种方法可以做到这一点？
浏览 19提问于2020-05-07得票数 1
回答已采纳
我在一个PIXI.Graphics中有一个PIXI.Container (还有其他一些东西，包括一个掩码和一个边框)。图形对象用于绘制各种多边形。容器的alpha属性设置为0.5。结果如下：
明亮的正方形是两个多边形之间的重叠。看起来，即使两个多边形都被绘制到相同的不透明图形对象上，但它们似乎是具有自己的alpha通道的独立对象。
是否有任何方法将所有多边形合并在一起，以使生成的图形具有统一的α，尽管有一些重叠的多边形？
Pixi版本为4.7.3。
浏览 1提问于2018-07-05得票数 1
回答已采纳
我在QGIS中有一个有趣的问题，我有一个矢量层的点(以紫色显示)和一个多边形层(绿色)，并附加了许多元数据属性(例如，日期，id号等)。
我想将多边形中的数据附加到它们下面的点数据上。问题是我在点数据被两个多边形覆盖的地方有一些重叠。
在这些情况下，我基本上想要复制点数据，因此在图层的属性表中，我为每个重叠的多边形都有一个点，并具有相应的元数据。这意味着我不能简单地使用“按位置连接属性”工具，因为它只会分配最近多边形的数据，或者两个重叠的总和/平均值等。
下面是数据的插图。任何关于如何解决这个问题的建议都将不胜感激。
浏览 6提问于2015-09-04得票数 4
回答已采纳
3回答
二维平面上的重叠多边形
、、、
我想要构建一个动态数据结构，它可以容纳一个多边形列表，并返回一个与指定矩形重叠的多边形列表。
我观察了bst树(和四叉树)，但是当多边形重迭时，这些树似乎不能很好地工作。
在我胡说八道之前有什么好主意吗？
编辑
假设所有多边形都是正常的、非旋转的矩形。我愿意在点测试期间(无论如何我也是这么做)接受打击(多边形测试中的点)，而在区域测试中，获得他们的边界框也同样好。其中只有一小部分实际上不会与有关区域重叠。
浏览 4提问于2010-08-10得票数 1
回答已采纳
1回答
我有一个带有几个重叠路径的shapefile GIS层。它们包括从多个起源到一个目的地的大地测量/最短路径计算。我想把这些路径想象成直线或者3D中的缓冲多边形，因为这在2D中是不可能的。理想情况下，我希望在R中这样做，以保持一致的配色方案和一个工作流程。ArcGIS Pro能够做到这一点-参见附件。
我试图用光线着色器实现这一点，但这似乎不能画堆叠线或多边形(附)。
plot=ggplot() +
geom_sf(data=Example,aes(colour = value), size=0.8) +
scale_color_viridis(option = "H"
浏览 7提问于2022-10-30得票数 1
回答已采纳
1回答
无缝匹配两个sf多边形
、、、
问题所在
假设我们有两个应该无缝边界的shapefiles。只是，他们没有。有什么办法可以迫使他们彼此坚持下去吗？
具体情况
我有两种形状:一种是欧洲地区的REG，另一种是邻国的NEI。这两个shapefiles都是从中提取的，应该很好地结合在一起，但是有很小的差距。此外，我需要简化多边形，然后差距变得非常显著。
我能想到的最好的
我尝试过几种方法，但都没有成功。实现我所看到的预期结果的唯一方法需要以下步骤：
创建一个线sf与只是边界之间的形状文件；
从这一行创建一个足够大的缓冲区多边形来覆盖所有的空白；
加入并将此缓冲区溶解到邻居的shapefile中- NEI；
浏览 0提问于2018-01-20得票数 23
回答已采纳
1回答
我试图通过使用sf图书馆和传单来找出哪个市政厅区代表了纽约的哪个邮政编码。我有两个区的空间多边形和邮政编码。我的问题是，每当我试图找到多边形之间的交点时，我就会得到额外的邮政编码，这些邮政编码与地区边界重叠。该地区没有一部分属于边界以外的邮政编码，但使用空间交叉口包括那些与邮政编码相邻的区域。
是否有一个空间操作函数只返回多边形的交集，而不是多边形只共享边界的地方？
以下是我的观点(见邮政编码(灰色)和议会区(黑色)重叠的边框)：
我的密码是：
library(sf)
#NYC Zip Codes
nyc_zips <- read_sf('ZIP_CODE_04011
浏览 1提问于2022-03-30得票数 0
1回答
我有一个包含多边形的集合和一个包含点位置的集合，我想看看每个多边形中有哪些位置，并用多边形的ID标记位置。
为了实现这一点，我想在多边形集合中实现一个存储过程，它检索所有多边形，然后对每个多边形调用位置集合中的存储过程。第二个过程将以多边形ID和多边形坐标作为输入，并使用ST_WITHIN方法检查多边形中的哪些位置。
这样可以解决吗?这是个好办法吗？
浏览 2提问于2017-07-10得票数 1
回答已采纳
2回答
从光栅中提取值
、、、
我试图从一层随机点的光栅中提取数据。
输入的数据是光栅，在那里我必须提取值和多边形的形状文件。通过这个形状文件，我得到了多边形内部点的随机抽样。这是通过SF包完成的，我得到了一个层sfc_POINTS。然后，我尝试用这些点提取光栅的值，使用光栅包。
我得到了这个错误: Error in (函数(classes，fdef，mtable) )：无法找到用于签名“RasterLayer”、“sfc_POINT”的函数‘提取’的继承方法
以下是代码：
# Clean environment
rm(list = ls())
#Import packages
library(sf)
library(r
浏览 2提问于2022-01-22得票数 -1
回答已采纳
我对处理多边形数据(在使用NetTopologySuite的C#中)有点陌生，我想获得几个多边形的边界框，这取决于每个多边形的边界框是否与另一个多边形重叠(多边形聚类)。 在这个演示中，我有3个多边形，它们的边界框与其他多边形重叠，并希望最终结果是红色边界框。 ? 基本上，我已经被多边形的边界框卡住了。我试过Geometry.Boundary，但那只会把外环还给我... 最后，我可以遍历坐标，但我想知道Geometry或Polygon类是否内置了此功能(或者库是否内置了此功能)。
浏览 59提问于2020-12-07得票数 4
回答已采纳
1回答
如何在没有用户交互的情况下在google地图上检索和对象。这意味着用户不单击对象。我有10个多边形，我希望用户导航到其中的每一个，操纵可能的颜色。而不点击多边形本身。
我确实有一个" next“按钮，当用户单击它时，它应该导航到下一个多边形，设置它是活动的，并且充满了渴望的颜色。
多边形以json格式存储在数组中。
我在使用数据层。我不知道如何实施。我可以点击多边形本身和改变颜色，但我需要遍历他们不点击。
谢谢
绘制地图：
函数loadJson(url，样式){ var jsonfile =$.getJSON(Url).done(函数(数据){map.data.addGeoJson
浏览 1提问于2014-06-26得票数 0
回答已采纳
2回答
我有代表下面两个多边形的图像精灵。这些多边形大致表示精灵区域。我想要做的是使用这些多边形来检测精灵的重叠(或碰撞)。但是，重叠应在绿色正方形内有效。(这是一个拼图游戏，我正在尝试实现的是当拼图碎片被移动到更近的地方时，它们会被抓取)
我尝试了Intersector.overlapConvexPolygons(adjacentPiece.polygon, currentPiece.polygon);，但是这个检测到了整个多边形的重叠。
我可以在这里做的任何聪明的事情来检测重叠。
浏览 0提问于2017-09-26得票数 2
我有一个带有西班牙干旱数据的sf对象，我想把这些数据投影到一张传单地图上。数据源是：
我正在使用函数esp_get_country()从图书馆mapSpain，以获得国家的边界。
问题是，在sf物体中，加那利群岛被代表在远离边界地图的地方。我如何从sf投射加那利群岛的多边形，因为它们在边界地图中？
我试过函数st_crs()从边界映射获得crs，函数st_transform()重新设计sf对象，但它没有做任何事情。
我使用的代码如下：
# get the boundary map from package 'mapSpain'
spain_coun_sf <- esp_g
浏览 2提问于2022-05-09得票数 0
回答已采纳
1回答
我有一个光栅和sf多边形，如下所示： library(raster)
libary(sf)
# Create raster r
r = raster(ncol=1000, nrow=1000, xmn=0, xmx=1000, ymn=0, ymx=1000)
values(r) = round(runif(ncell(r),1,10))
# Create sf polygon
poly_sf = st_sfc(st_polygon(list(cbind(c(0,10,50,100,0),c(0,70,300,500,0))))) 栅格包含值在1到10之间的单元格。我希望能够生成一个数据
浏览 28提问于2020-04-27得票数 2
回答已采纳
我正在尝试在罗宾逊投影中投影一张世界地图，其中中央子午线不同于0。根据this StackOverFlow thread的说法，这应该是一件容易的事情(尽管示例中使用了sp)。 下面是我的可重现代码： library(sf)
library(ggplot2)
library(rnaturalearth)
world <- ne_countries(scale = 'small', returnclass = 'sf')
# Notice +lon_0=180 instead of 0
world_robinson <- st_transform
浏览 41提问于2019-05-15得票数 0
回答已采纳
我很难把传说画成一幅格格图。我已经在网上搜索了好几个小时了，但还没有找到任何解决方案。
我正在尝试创建一个具有不同的线条和多边形形状的ggplot。我用函数geom_sf()绘制shapefile。这幅画不成问题。我的问题和传说的创造有关。
目标是创建一个有一个图例的情节，并有所有不同的形状文件类型，包括它们的颜色和或填充优雅。这意味着，线条形状应在图例中表示为简单直线，多边形形状应表示为简单多边形。我已经做到了这一切。
主要的困难在于，一个多边形(为了更好的可读性)不仅使用“填充”参数作为astetic，而且使用“颜色”参数。
如果我创建一个所有多边形都使用填充参数作为astetic的情节，
浏览 6提问于2022-02-25得票数 3
回答已采纳
3回答
我有一堆来自3D模型的重叠三角形投影到2D平面上。我需要将每个相互接触的三角形岛合并成一个闭合的非凸多边形。
生成的多边形不应该有任何洞(因为源数据没有)。
许多源三角形与源数据中的其他三角形共享(浮点相同)边。
做这件事最简单的方法是什么？性能并不是特别重要，因为这将在设计时完成。
浏览 1提问于2010-05-12得票数 8
回答已采纳
我们有一个应用程序，它有一个充满多边形(当前存储为点)的数据库，.net应用程序拉出这些多边形并检查它们是否重叠。
我突然想到，将这些点数组转换为数据库中的多边形/折线对象，并使用sql来了解它们是否重叠会更好。
我已经看到了不同的方法建议这样做，但没有一个例子是非常符合我的需要。
我将非常高兴收到那些好心的人的意见，以提供他们的经验。
其他：
回答问题:它确实是2D的。是的，这两者的任何交叉都被认为是正确的。多边形有n个点，可以是凹面的。多边形将保存为每行1个(在数据转换任务之后)作为多边形(即多边形类型。它可能被称为其他空间/ geom我的记忆现在不在我这边)
浏览 0提问于2012-01-10得票数 2
回答已采纳
我想知道两个矩形是否使用SFML相互重叠。这是我的代码：
if (ball.getLocalBounds().intersects(paddle.getLocalBounds()))
{
//perform action
}
如果两个矩形相互重叠，代码将执行该操作。但是，它总是在每一种情况下返回true，即使在两个矩形甚至不相交的情况下也是如此，如下所示：
从图片中可以看出，左边的桨和球甚至不是很近，但是控制台仍然说它们是重叠的。这里发生了什么，你怎么解决的？
编辑:我试着用足够的代码来创建一个测试项目来重现问题。从两个矩形的位置上看，它们绝对不是相交的。但是，它仍然说它们
浏览 0提问于2018-03-03得票数 0
回答已采纳
1回答
我在一个叫做data.frame的walktime中有一些多边形。
> class(walktime)
[1] "sf" "data.frame"
我想把多边形组合成一个。据我所知这是有可能的
walktime2 <- sf::st_combine(walktime)
> class(walktime2)
[1] "sfc_MULTIPOLYGON" "sfc"
但是如果我绘制walktime2，我的组合多边形上有洞。
我怎么才能避免这种情况？
谢谢你的帮助。
可复制代码在这里：
浏览 2提问于2022-12-01得票数 2
大家好，
我正在努力解决这个问题，希望有人能拿出一个简单的解决方案。
我的目标是在多边形的范围上创建一个规则多边形网格，但是由用户定义的角度旋转。
我知道我可以很容易地在sf中创建一个南北多边形网格，例如：
library(sf)
#> Linking to GEOS 3.6.2, GDAL 2.2.3, proj.4 4.9.3
inpoly <- st_read(system.file("shape/nc.shp", package="sf"))[1,] %>%
sf::st_transform(3857) %>%
sf:
浏览 0提问于2018-07-11得票数 14
回答已采纳
1回答
我使用以下代码将一个栅格读入我的R会话：
raster <- stack("raster.tif")
现在我想创建一个简单的特征(sf)对象来表示该栅格的轮廓。我不能使用边界框，因为光栅是多部分的，所以边界框会比光栅大得多。因此，轮廓线也需要是多部分特征(sf多多边形)。
如果能帮上忙我将不胜感激。谢谢!
标记
浏览 2提问于2020-05-29得票数 0
1回答
单叶多边形排序
、、、
现在有没有一种方法，如何使用Leaflet对多边形进行z索引？当地图启动时，我是可以的，但当向现有地图添加新多边形时，我需要根据面积将新多边形排序到现有多边形中-这样较大的多边形就不会与较小的多边形重叠。我找到了这个解决方案：
但是当我的地图包含更多的要素时，它是非常慢的。有什么想法吗？
浏览 0提问于2015-11-21得票数 0
我有一组多边形，它们可能共享公共边和节点。所有这些多边形都是严格不重叠的，尽管它们可以共享一个公共顶点或边。
我想批量三角剖分所有这些多边形，因此，我能想到的解决方案是约束delaunay三角剖分。但是约束Delaunay三角剖分的输出将生成不在原始多边形中的三角形。
有没有办法识别这些多边形外的三角形？
编辑：；我正在寻找这样一个独立于语言的算法。
浏览 0提问于2011-11-30得票数 2
1回答
我对R和sf包还是比较陌生的。
我有两组要分析的多面体数据。我的第一组多边形(火灾)包含数百个野火边界。第二个集合(城镇)包含数百个城市区域边界。
对于每个火灾，我想计算到最近城镇的距离(火灾多边形边缘到最近的城镇多边形边缘)，并将其作为字段添加到每个火灾中。
到目前为止，我主要使用sf包来处理空间数据。在我的搜索中，我只能找到多边形到点、点到点、线到点等的最小距离方法，但似乎找不到多边形到多边形的示例。如果能帮助我找到正确的方向，我们将不胜感激！谢谢。
浏览 2提问于2018-12-19得票数 6
我正在开发一个2D程序，它需要对多边形(由用户生成或导入，可能是凸的还是凹的)与AABB的交叉(重叠)进行大量的判断。我用kdtree实现的数据结构部分，现在需要一种方法来确定多边形是否与AABB盒相交。
我在网上见过很多讨论，但它们都是关于凸多边形的，我也见过用多边形的每个线段来确定它是否与AABB相交的方法，但由于它是AABB盒，恐怕这种方法太昂贵了。
1多边形由一系列按逆时针顺序存储的Vector2列表定义，AABB由Listminx、miny、maxx、maxy定义。
我只关心它们是否相交，我不需要返回重叠范围或交点。
3任何多边形的预处理都是可能的。
4所使用的语言是C#，欢迎任何可
浏览 4提问于2022-06-04得票数 1
我试图将一个1度1度的普通地理网格转换为多边形，同时使用R软件包维护与每个网格点相关的数据。
示例代码：
library(tidyverse)
library(sf)
library(sfheaders)
library(sfheaders)
# latitude
lldata <- tibble(lat=as.numeric(rep(seq(40,49.5),10)))
# Longitude + add data
lldata <- lldata %>%
group_by(lat) %>%
mutate(lon =
浏览 2提问于2022-02-08得票数 0
回答已采纳
我有多个多边形和一个光栅组成了我的地图。我想“放大”到这些多边形中的每一个，以便更好地可视化每个多边形的光栅信息。我想我最终可以使用facet_grid和ggdraw的一些组合来绘制最终的地图，但是现在我需要帮助创建地图的“缩放”部分。
示例：
#Raster layer
library(terra)
library(raster)
library(sf)
f <- system.file("ex/elev.tif", package="terra")
r <- rast(f)
#polygons
v <- vect(system.file
浏览 4提问于2021-08-20得票数 0
回答已采纳
我有两个多边形'.shp‘文件。我需要一个通过填写一个变量出现在地图上，另一个只出现在边界上，重叠第一个。
我已经用过“ggplot2”和“sf”了。
我绘制了一个地图('map1')，它是用'ggplot‘和'geom_sf’分层的多边形。
我使用“map1”中包含的变量( 'map1‘)作为“填充”。
现在，我需要添加(覆盖)另一层多边形的顶部('map2')。这必须是“透明填充”或“不填充”。只出现了边界的轮廓。
library(ggplot2); library(sf)
map1 <- st_read("m
浏览 3提问于2019-09-19得票数 4
回答已采纳
1回答
在R中将多边形转换为sf
、、、
使用这里的教程：https://www.r-spatial.org/r/2018/10/25/ggplot2-sf-2.html (在“状态(多边形数据)”一节中)，我尝试将多边形转换为sf对象。 library("maps")
states <- st_as_sf(map("state", plot = FALSE, fill = TRUE))
head(states) 我得到了错误消息：Error in as_mapper(.f, ...) : argument ".f" is missing, with no default 提前
浏览 10提问于2020-07-25得票数 2
回答已采纳
我有多边形，我想要计算它们之间重叠的百分比。其思想是，当一个多边形与另一个多边形相交时，可以根据一个多边形或另一个多边形的透视图(即最大值)计算百分比。因此，我想要制作一个脚本，在多边形之间生成百分比覆盖率--从一个多边形到另一个多边形的重叠百分比，并将所有的结果放在一个数据帧中。
下面是我目前的代码：
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回答
我设置了一个PostGIS数据库，我通过参数化SQL视图在GeoServer中添加了该数据库。我用传单显示这一层通过wms。
在我使用url "/geoserver/gwc/service/wms“而不是"/geoserver/wms”添加GeoWebCache之前，它工作得很好。我仍然可以看到我的多边形，当我在最小的缩放。但是当我放大时，我只看到一个红色多边形和一个绿色多边形的一半，如果我再次放大，我只看到红色多边形。您可以在下面的图像中看到这3种状态：
我想这是打瓷砖的一个问题:我得到了最小的瓷砖，还有一些红色多边形周围的瓷砖，以进一步放大，但出于某种原因，其他的
浏览 3提问于2021-01-26得票数 0
回答已采纳
1回答
我有5个多边形的形状文件，我有重力场的负载。我想看看这些多边形是否互相重叠。
我想识别所有重叠的多边形。我在精巧的手册中看到，is_valid工具可以检测到：
如果功能在1的意义上是“有效的”，则 object.is_valid返回True。
备注
有效性检验仅对多边形和MultiPolygons有意义。对于其他类型的几何图形总是返回True。
有效多边形不可能有任何重叠的外部或内部环。有效的MultiPolygon不能收集任何重叠的多边形。对无效功能的操作可能失败。
所以我在多边形上试了一下：
imprt geopandas as gpd
import shapely as sh
浏览 0提问于2020-09-24得票数 0
回答已采纳
我有3个多边形，一个被其他两个重叠。有没有办法“走到前面”？请展示一个编程示例或演示如何做到这一点。
浏览 2提问于2016-03-03得票数 0
回答已采纳
1回答
我正在尝试创造一种体验，我在地图上有几个详细的建筑物3D模型，通过矢量瓦片源突出邻近建筑物的建筑物轮廓。3D模型将是主要焦点，挤出的轮廓线将作为参考。我遇到的一个挑战是，我有一个全局建筑轮廓层，它有一个3D建筑的轮廓线，这并不完全匹配。此外，当挤出时，它最终会合并/重叠漂亮的3D模型。
我希望能够隐藏与3D模型重叠的各个轮廓线。我最初的想法是获取3D模型的边界框，然后使用新的within样式表达式，但看起来这只会过滤点和线，而不是多边形。建筑物轮廓多边形中没有我可以用来过滤的唯一信息。
我知道我可以监控地图移动，查询渲染的特征，并手动检测相交的多边形，但由于轮廓线上没有唯一的属性，我不能过滤或
浏览 0提问于2021-03-13得票数 0
3回答
我有一个多边形的列表，在这个列表中，一些多边形是重叠的，或者与其他多边形接触。
我的任务是合并所有相互重叠或接触的多边形。我有一个union方法可以做到这一点。
执行此操作的最有效方法是什么？我现在能想到的就是遍历多边形列表，检查合并列表，看看这个多边形是否已经属于合并列表中的一个多边形，如果是，合并它们，如果不是，将这个多边形添加到合并列表的末尾。
再次重复上述步骤几次，以确保所有多边形都已正确组合。
这种方法看起来很不雅观；有没有更好的方法呢？
浏览 3提问于2010-08-02得票数 4
1回答
我正在尝试使用Tilemap在2D场景中实现A*寻路算法。为了做到这一点，我需要计算网格中的哪些块是可以通过的。我有一个TilemapCollider2D连接到我的Tilemap上有墙的层，但我正在努力将对撞机的边界与每个节点的边界网格化(没有双关语)。
我的Tilemap中的每个瓦片都被细分为4个节点，用于A*算法，以便提供更精细的移动粒度。我需要能够获取每个探路器节点的边界，并确定它是否包含对撞器区域。我的第一个想法是使用TilemapCollider2D.OverlapPoint()检查每个节点的左上角、右上角、左下角和右下角的点，看看这些点是否与对撞机重叠，但后来我意识到这样做的问题是
浏览 4提问于2018-02-18得票数 0
3回答
目标是通过像多边形画笔一样在屏幕上拖动鼠标，创建画笔的和鼠标的路径，从而进行简单的矢量图像编辑。新的多边形将从任何先前存在的不同颜色的多边形中减去，并与任何相同颜色的现有多边形合并。
我们的计划是将每次鼠标移动作为从鼠标的先前位置到其当前位置的线段，计算该线段上的Minkowski和，然后使用更新现有多边形以包括该Minkowski和。
由于Weiler-Atherton似乎可能会导致UI延迟，如果每次鼠标移动都运行，我们计划延迟这一步，将其放入另一个线程，该线程可以花时间赶上最新的鼠标移动，或者保存所有Weiler-Atherton计算，直到绘图完成，然后在保存时作为批量操作执行。我们担心这
浏览 0提问于2017-09-04得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
