geospatial相关内容

在R中通过世界贴图裁剪空间多边形

这是我第一次在R中进行任何类型的空间数据可视化，我被某个特定的问题困住了。我想根据世界地图剪裁一个空间多边形(由一系列经度/经度坐标指定)，以便删除与地图多边形重叠的任何多边形部分。以下面代码中的内容为例，我希望剪裁矩形空间多边形，以便只保留该多边形的海洋部分。 我已经找到了如何保持两个空间多边形之间的交集的示例，但我想要做相反的事情。也许有一种方法可以定义交点，然后从我想要裁剪的多边形中减 ..
发布时间：2022-07-21 09:30:40 其他开发

一组坐标中的边界点

我有一组很长的点，我想从这些点中提取形成边界的点，我使用了凸包，但对于我的包来说，这是不够的，因为Transporthull只返回形成多边形的最远的点，在所有点都适合的地方，我需要形成透射器的所有点，类似于我附加的图像。我能做什么呢？是否有某种包可以随时使用，而不是实现任何空间算法？ 谢谢 推荐答案 您必须使用用于凸面的包。下面是一个例子： import alphashape i ..
发布时间：2022-07-14 16:49:42 Python

从R中的netcdf创建栅格最准确的方法是什么？

我处理netCDF数据已多年。NetCDF用于空气污染物数据，纬度和经度作为单独的变量提供，而不是原始网格的一部分。 链接到日期：Sample Netcdf 这些netCDF文件提供2级二氧化氮数据，可从NASA EarthData门户网站下载。卫星是Sentinel-5P，仪器是Tropomi。 因此，在处理这些数据时，您必须为NO2、纬度和经度创建变量。我正在尝试创建栅格图层，然 ..
发布时间：2022-05-13 19:59:28 其他开发

对于给定点，查找具有负值R的最近格网像元

我在R中有一个栅格类型对象，它是一个经度/经度栅格，每个单元格都有一个值，即该单元格的深度或高程(从玩具示例中的marmapgetNOAA.bathy()，bathymap下载)。 然后我就有了一个具有最晚、最长坐标的点列表(points)。对于每个点，我希望找到距离值为负(深度)的点最近的栅格单元。如果有多个单元格，那么我想随机选取其中一个。一旦我有了/一个最近的负网格单元格，我想要查找并保存 ..
发布时间：2022-05-13 19:55:34 其他开发

如何在R中为.shp多边形文件设置CRS

我对在R中使用地理空间数据非常陌生。我在我的研究区域中有一个由几个象限组成的Shapefile，我根据“QUIDID”对其进行了分类。我可以很好地使用gglot和sf包来显示我的数据，但是，我的数据似乎不能正确显示。在形状文件中检查任何空间数据信息时，似乎没有： > proj4string(quads) Error in (function (classes, fdef, mtable) ..
发布时间：2022-04-09 20:52:31 其他开发

R RASTER：：EXTRACT无法创建数据。帧

我希望使用RASTER：：EXTRACT和SF对象从栅格中提取值。然而，我遇到了一些困难。我希望创建一个数据框，其中包含单元格号(以派生XY值)、提取的单元格值和ID字段(如果我从多行提取)。我通过以下几点做到了这一点： dtm_values ..
发布时间：2022-04-09 20:49:34 其他开发

SF&；dplyr：按相同坐标分组不起作用

我有一个从S2卫星栅格(10x10米)提取的数据集，其中有12个值(ras.df.ll)，但6个在一个块(T21JYG)中，第二个在另一个(T21JYG)中。 我想计算瓷砖之间相同(x，y坐标)的平均值，但没有成功。我找不到任何方法来识别第一个块中的第一行与第二个块中的第一行的坐标相同，只是我的数据集的末尾。在我的示例中： library(sf) library(sfheaders) lib ..
发布时间：2022-04-09 20:48:25 其他开发

两个地理序列的指数是不同的--理解指数

我正在使用GeoPandas，我有两个GeoDataframes和相同的CRS。其中一个包含带有多边形geometry的geometry列，另一个包含带有点geometry的列。我想检查哪些点在多边形内部。 我天真地尝试 shape.contains(points) 这给了我 > The indices of the two GeoSeries are different ..
发布时间：2022-04-09 20:42:32 Python

Bing地图地理编码器的结果不一致

我编写了一个应用程序，用于解析CSV文件、提取地址并使用Bing Maps rest API对其进行地理编码。我遇到的问题是，如果我多次运行同一文件，我会得到不同的结果。在第一次尝试时，第2行和第6行可能不返回匹配项。在下面的试验中，第2行和第6行将返回匹配项，而第4行可能不会。我使用Fiddler检查了HTTP流量。我的请求格式正确，每个响应都返回HTTP状态代码200和有效的JSON，只有一些 ..
发布时间：2022-04-09 20:40:03 其他开发

绘制R区粮农组织捕鱼区地图

我想根据一个数据集(在我的例子中，是鲨鱼物种的长度数据)，用R标出粮农组织捕鱼区的颜色。 我更喜欢在ggplot中做一个全息地图，但其他类型的地图也是可以的。最坏的情况是粮农组织区域的基本地图，我可以在其中添加泡沫。即使只是一张现有的粮农组织地区基本地图也是很棒的。欢迎提出任何建议！ 推荐答案 我转到this page并单击以查找this link以检索GeoJSON文件： ..
发布时间：2022-04-09 20:38:38 其他开发

用geom_sf绘制基于点的核图

有人知道如何使用ggplot2和sf绘制基于事件发生的内核密度图吗？ 例如，考虑来自sp包的meuse数据集(假设每个点都是一个事件发生) library(sf) library(ggplot2) # data data(meuse, package = "sp") # as_sf meuse_sf ..
发布时间：2022-04-09 20:36:45 其他开发

处理空间数据：如何在不替换的情况下找到最近的点邻居？

我目前正在处理一些森林调查数据。 这些数据是在样地上收集的，样地的位置可以作为点数据(空间数据)。 我有两个数据集： 数据集dat.1，物种A的n个样地 数据集data2，物种B的k个样地 使用n<；k 我想要做的是将数据1的每个点与一个数据点2进行匹配。结果应该是n对点。因此，应从dat2中选择k个曲线图中的n个。 匹配条件为： 一对点之间的空间距离应尽 ..
发布时间：2022-04-09 20:35:28 其他开发

Db地理多边形获取点

我在SQL Server 2012数据库上以Sys.Geography类型持久保存了一个Polygon。如何获取该多边形的所有点？ 我正在考虑使用AsText()方法并分析字符串，但或许还有更好的选择？ 推荐答案 找到了方法，下面是一个扩展方法： public static IEnumerable GetPointsFromPo ..
发布时间：2022-04-09 20:33:31 其他开发

缓冲区重叠时计数点数

我在下面包含了我的所有代码和指向示例数据的链接。 简要说明：我的缓冲区重叠；我想计算离学校一定距离内的商店数量。 我特别想知道离学校1000米以内有多少家商店，离学校2000米以内有多少家商店，因为我想比较一下差异。当然，其中一些学校的缓冲是重叠的。因此，一家商店可能距离A校1500米，但距离B校只有750米。因此，它被算作距离一所学校1000米以内，应该只被计算在B校的1000米内， ..
发布时间：2022-04-09 20:32:02 其他开发

为什么使用st_intersects而不是st_intersects？

st_intersection与st_intersects相比非常慢。那么，为什么不使用后者而不是前者呢？下面是一个带有小玩具数据集的示例，但是执行时间上的差异是巨大的，因为我的实际集只有62,020个点，与实际的地理区域边界相交。我有24 GB的RAM，st_intersects代码需要几秒钟，而st_intersection代码需要15分钟以上(可能更多，我还没有耐心等待...)。st_int ..
发布时间：2022-04-09 20:28:44 其他开发

R|SF：我有点，每个点周围有2个缓冲区。如果较大的缓冲区重叠(但较小的缓冲区不重叠)，如何将这些点合并为单个单元？

我目前在地图上有点(学校)。每个人在点(学校)周围都有两个缓冲区。一个是450米，一个是250米。如果点重叠，我希望将它们视为单个单元(因为否则事情会变得复杂)，但我希望它们保持其覆盖的几何图形/面积。 因此，在这里给出的示例地图上，我希望将排名前三的学校/点合并为一个单位。&我希望它们保留其覆盖的面积，但仅将R作为一个单位计算。如果我使用"；st_Union"；函数，则必须 ..
发布时间：2022-03-23 17:14:02 其他开发