优步h3能以六角形网格覆盖全球吗？
优步发布了h3，这是一个高效处理地理空间文件中大数据的框架。使用h3，我尝试获得如图所示的六边形网格位置。 (https://eng.uber.com/h3/) 我从以下代码中获得了六边形网格的位置。然后我把它画在一张二维地图上，看看它是否覆盖了整个地球。但是，我在边界(-90°、90°、-180°、180°)上得不到有效的六边形。而且它似乎并没有覆盖整个全球。 (hexagonal gri
..
gis相关内容
如何让一个cran包在任何操作系统下只下载一次数据？
CRAN policy limits R package size to 5 Mb，这对于地图等图形应用程序来说很少。处理包裹大小限制的方法有多种，所有这些方法都有其缺点。下面列出了备选方案。 我的问题是：如何使R包只下载一次数据文件(即在重启后将它们保存到R可以找到它们的位置)？该解决方案应适用于所有常见的CRAN平台。 我一直在开发一个mapping package for R，它应该
..
GeoDjango：查找半径内的对象
我当前正在尝试获取半径内包含的点的列表，但无法使其工作。到目前为止，我的视图代码如下： from django.contrib.gis.geos import Point from django.contrib.gis.measure import Distance class AreaInfoViewSet(viewsets.ViewSet): queryset = models
..
GeoPandas中的过度功能不起作用
我只想使用geopandas获得两个多边形区域的并集和交集。我定义： import geopandas as gpd from shapely.geometry import Polygon polys1 = gpd.GeoSeries([Polygon([(0,0), (2,0), (2,2), (0,2)]), Poly
..
在Python中根据地理数据框内的点创建地图边界
我有一个名为map的地理数据框，其中包含一个点列表和一个Closest_TrainStation_name列，其中包含距离该点最近的火车站的名称。 几何 Cost_Trainstation_Name 1 point(1，1) Station_1 2 点(10，10) Station_2 ... ... ... 是否可以创建包含每个组的边界多边形，如下图所示？其中每个面具有距离原始
..
在R中通过世界贴图裁剪空间多边形
这是我第一次在R中进行任何类型的空间数据可视化，我被某个特定的问题困住了。我想根据世界地图剪裁一个空间多边形(由一系列经度/经度坐标指定)，以便删除与地图多边形重叠的任何多边形部分。以下面代码中的内容为例，我希望剪裁矩形空间多边形，以便只保留该多边形的海洋部分。 我已经找到了如何保持两个空间多边形之间的交集的示例，但我想要做相反的事情。也许有一种方法可以定义交点，然后从我想要裁剪的多边形中减
..
使用R处理多个栅格文件时保留原始文件名
我有以下问题：我需要使用R包landscapemetrics中的同一函数处理多个栅格文件。基本上，我的栅格文件是国家地图的一部分，形状和大小都相同(即象限)。我想出了一个文件的代码，但我必须对600多个栅格做同样的事情。因此，手动操作是非常不合理的。我的代码中的步骤如下： # 1. I load "raster" and "landscapemetrics" packages: libra
..
使用OpenLayers 5在地图上添加旋转的卫星图像
我正在尝试使用OpenLayers 5在地图上添加卫星图像。 问题是我无法做到这一点，因为我刚刚找到了一个选项，可以在地图上添加经过图像范围(xmin、ymin、xmax、ymax)而不是边界框的图像。图像应适合边界框内。因此，图像被扭曲。 图像在JPG文件中(属性Feature.Properties.ICON)。示例：http://exampleserver.com/220/063/
..
LEAFLE+EPSG：31466 WMS层
我要使用epsg：31466投影创建传单地图。我已经包含了proj4库和proj4leaflet插件。我计划使用fp_plan层包括thisWMS服务。 我已尝试使用以下方法(fiddle)设置我的地图。 // Define CRS (EPSG:31466) var rs31466 = new L.Proj.CRS( 'EPSG:31466', '+proj=tmerc +la
..
为什么也具有LEFT JOIN的for循环中的Tmap函数不起作用
我有一个Shapefile“GRID”，基本上是美国和加拿大的网格地图。每个网格都有自己的特定ID(GRID_CODE)。我还有一个CSV文件，其中包含与各个网格相关的鸟类物种的相对丰富值。由于每个物种都有其唯一的分布，因此物种表仅包括所有可用网格代码的子集。 我正在尝试制作一张专题地图，显示每个物种的相对丰富度。为此，我首先使用leftJoin函数将CSV文件连接到网格Shapefile的
..
在RStudio中使用do.call合并多波段卫星影像栅格块时出错
我正在尝试将一些多波段卫星图像合并到R中(不幸的是，由于版权问题，我不能分享这些图像)。它们具有相同的大小和分辨率，但它们之间略有重叠。当我单独导入每个图像并合并它们时，没有问题： library(raster) raster1
..
使用PostGIS将点转换为面
我想使用PostGIS创建一个多边形表。表‘point’中的每一行都有三个点ID. 表‘point_location’具有点的位置信息。我用谷歌搜索了这个问题，但没有找到答案。以下代码有什么问题？ SELECT ST_GeomFromText('POLYGON((' || b.x || ' ' || b.y || ',' || c.x || ' ' || c.y || ',' || d.x
..
增加边框大小，但保持比率(Lon/Lat)
我有一个由下列值定义的Bbox： xmin: 11.555333537980914 ymin: 47.76067947037518 xmax: 11.995692579075694 ymax: 48.281587762758136 我希望增加此Bbox的大小，但保持该比例。 我尝试的一种方法是计算Bbox的中点，并计算半径值增加50%的新Bbox。 问题是：比率丢失了。 我如
..
R|SF：我有点，每个点周围有2个缓冲区。如果较大的缓冲区重叠(但较小的缓冲区不重叠)，如何将这些点合并为单个单元？
我目前在地图上有点(学校)。每个人在点(学校)周围都有两个缓冲区。一个是450米，一个是250米。如果点重叠，我希望将它们视为单个单元(因为否则事情会变得复杂)，但我希望它们保持其覆盖的几何图形/面积。 因此，在这里给出的示例地图上，我希望将排名前三的学校/点合并为一个单位。&我希望它们保留其覆盖的面积，但仅将R作为一个单位计算。如果我使用"；st_Union"；函数，则必须
..
如何将GeoDataFrame导入MySQL？
我有一个具有以下结构的GeoDataFrame： RangeIndex: 28 entries, 0 to 27 Data columns (total 3 columns): # Column Non-Null Count Dtype --- ------ -
..
Networkx读取军械调查-ITN综合传输网络？/读取GML文件
我遇到一些问题。 我想将军械勘测-ITN综合传输网导入networkx进行一些网络分析。 由于ITN是GML 2.1.2格式，我尝试使用 nx.read_gml 但是它警告我 ParseException: Expected "graph" (at char 0), (line:1, col:1) 然后，我花了一些时间查看我的文件的内容，发现它与其他人的示例中使用的
..
为什么terra：：cellSize()和raster：：area()会生成不同的栅格单元面积估计值？
我刚刚注意到terra::cellSize()生成的单元格面积估计值与raster::area()生成的单元格面积估计值不匹配。 首先，为什么这两种方法不能提供相同的答案？第二，哪个估计最准确？请参见下面的示例。 library(raster) #> Loading required package: sp library(terra) #> terra version 1.3.4 #
..
如何将世界划分为单元格(网格)
如何将世界划分为大小几乎相等单元格，以便以后的每个单元格都可以映射到不同的单元格？ 我非常确定我见过一个库这样做，将单元格标记为S1、S2等。 假设我们有62.356279，-99.422395，如何将其映射到名为&FR，23&q；的2 km*2 km单元格？ 谢谢！ 推荐答案 PostGIS3.1+ PostGIS 3.1引入了非常易于使用的网格生成器，即ST_S
..
从巨大的栅格文件中删除NA值
我当前正在尝试从一个巨大的光栅文件中删除NA值(1.9*10^7个观测值)。在这些栅格中，99.9%是NA值。我的目标是删除NA并创建一个包含所有非NA值的.csv文件。 我的尝试如下： # Load packages packs = c('raster', 'rgdal') sapply(packs, FUN = 'require', character.only = TRUE) xy
..
解决IF(obs<；=ei)2*PV ELSE 2*(1-PV)中的错误：APE包Moran&39；的I函数在R中需要TRUE/FALSE的位置缺少值。
开发人员！ 我遇到错误消息 IF(OBS<；=EI)2*PV ELSE 2*(1-PV)中出错：缺少值，其中 需要TRUE/FALSE 阻止我从APE包中获取Moran的I函数的值。以下是我所做的： library(ape) nrstp
..