Questions tagged «de-9im»

我正在进行一个环境流行病学项目，该项目中我有接触点（约2000头工业猪操作-IHO）。这些IHO喷洒在附近的田地上，但粪便中的水滴和气味可能传播数英里。因此，这些点曝光获得了3mi的缓冲区，我想知道每个NC人口普查区块的IHO曝光数量（各种-粪便量，猪的数量之和；最简单的，只是重叠曝光缓冲区的数量） （〜200,000）。排除人口普查区（蓝色）是（1）在人口最多的前5个城市中的任何事物，以及（2）不与县内有IHO接壤的县（注：这是使用gRelate函数和DE-9IM代码完成的-非常光滑！）。见下图可见 最后一步是将缓冲的曝光量表示汇总到每个人口普查区块。这就是我感到难过的地方。 到目前为止，我在sp程序包中使用％over％函数已经过得很愉快，但是从过分的插图中了解到poly-poly和poly-line over是在rgeos中实现的。该小插图仅涵盖线型多边形和自引用多边形，而不涉及聚合，因此，对于具有函数聚合（例如求和或均值）的多边形的选择，我有些困惑。 对于测试用例，请考虑以下使用世界国家边界文件的冗长代码段。这应该能够被复制出来并按原样运行，因为我为这些点使用了随机种子，并且由于我正在代码中下载和解压缩世界文件。 首先，我们创建100个点，然后将over函数与fn参数一起使用以在数据框中添加元素。这里有很多要点，但请看一下澳大利亚：3分，标号3。到目前为止，一切都很好。 现在，我们变换几何形状，以便我们可以创建缓冲区，向后变换并映射这些缓冲区。（包括在以前的地图中，因为我仅限于两个链接。）我们想知道每个国家/地区重叠了多少缓冲区-以澳大利亚为例，这是4。虽然要通过over函数来实现。在最后的代码行中看到我的一团糟。 编辑：请注意，关于r-sis-geo的评论者提到了聚合函数-也引用了堆栈交换问题63577-因此，解决/流程可能通过该函数进行，但是我不明白为什么我需要去当结束似乎具有其他空间对象的功能时聚合。 require(maptools) require(sp) require(rgdal) require(rgeos) download.file("http://thematicmapping.org/downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip", destfile="world.zip") unzip("world.zip") world.map = readOGR(dsn=".", "TM_WORLD_BORDERS_SIMPL-0.3", stringsAsFactors = F) orig.world.map = world.map #hold the object, since I'm going to mess with it. #Let's create 500 random lat/long points with a single value in the data frame: …

12 r  point-in-polygon  sp  de-9im 

考虑以下简单情况，其中三个边在一个节点处连接： 我想对A和B之间的关系进行简洁明了的描述，以使其与A和C之间的关系区别开。到B，但是A不相邻？到C。” 但这并不是真正的邻接。 用另一种方式说：假设您站在节点上并且正对着A。您开始按顺时针方向旋转。下一个优势是B，而不是C。 有没有一种方法可以比我上面写的更简洁，正式或正确的方式描述A和B之间的这种关系？ 它必须是有方向的（这种关系从A沿顺时针方向存在，而另一种沿逆时针方向存在）。并且必须扩展到节点上连接三个以上边缘的情况。也许与路由有关？（我在道路网络的背景下正在考虑这一点。） 我已经尝试过但还没有走很远的两种方法： 类似于9IM的拓扑参考：我研究过DE-9IM，即使我不是数学家，我仍然可以从图表和术语中看出它并不涵盖这种类型的关系。我在ESRI帮助或Oracle帮助的拓扑描述中都找不到它。（也许那里有东西，但是我还没找到呢！） 面孔：我一直在考虑这样一个事实，即A的“北”侧的面孔也可能受B限制，但不受C限制。但是，正如您在此图中看到的那样，这并不总是正确的。想象一下，我的图是从道路网络中摘录的，其中A和C是主干道路，B是一条短死角道路。 我怀疑我要说的话可能没有一个术语。至少我希望能够以一种比上述简单的方式描述这种关系。这是一个与平台无关的问题。现在，我只是在寻找正确的词。稍后，我将尝试在shapefile中以python（pyqgis或arcpy）实现该概念，因此考虑到该终点的任何答案都将特别有趣，但并非必需。

11 vector  topology  terminology  de-9im 

我尝试了解DE-9IM交集矩阵的空间谓词。我不明白谓语之间的区别封面（各自CoveredBy）和包含（相应地，在）。 我试着用画几何JTS拓扑套件，他们总是Covers = True 和 Contains = True（或两者都是假），但我不管理，使他们不同，即看到的区别Covers和Contains。是否有人举例说明了差异（最好是WKT或图片）？

9 intersection  jts-topology-suite  de-9im