标签[gis] - 堆栈内存溢出

如何在 R 中使用 Leaflet 根据数字变量绘制颜色线 - How to draw a line with color based on a numeric variable using Leaflet in R

这是怎么回事！我有一个数据框，其中每一行都有 object 的即时速度和坐标（纬度和经度）信息。我正在尝试画一条连接所有这些经纬度的线，这条线的颜色就是速度。任何人都可以帮助我吗？这是一些代码： ...

从 shp 文件中获取瑞士的救济 - Get relief of Switzerland from a shp-file

我想要通过瑞士的 map 绘制浮雕的官方来源。我下载了 https://cms.geo.admin.ch/ogd/topography/DHM25_BM_SHP.zip来使用 dhm25_p.shp-File 现在使用代码我会得到错误 ...

从 tiff 文件中解除瑞士 - relief of Switzerland from a tiff-file

我想要通过瑞士的 map 绘制浮雕的官方来源。在 swisstopo 的其中一页上，我找到了以下文件： https://data.geo.admin.ch/ch.swisstopo.pixelkarte-farbe-pk1000.noscale/swiss-map-raster1000_2022_ ...

我可以在 R 工作区 (RData) 文件中“保存”一个 SpatRaster 吗？ - Can I 'save' a SpatRaster in a R workspace (RData) file?

我注意到，如果我使用 terra package 从存储在本地硬盘上的文件加载栅格，然后稍后加载工作空间，则 SpatRaster object 没有任何与之关联的数据。有没有办法在保存和加载工作区时保留与 SpatRaster object 相关的所有信息？这是一个示例代码来说明这个问题： ...

使用组合键合并 pandas 中的数据框 - Merge dataframes in pandas with a combination of keys

我有两个数据框，我需要根据一个键（一个“事件编号”）将它们组合在一起。然而，关键是重复的，因为它们将被摄取的数据库需要特定的坐标格式。如何根据键的组合加入必要的列？例如，这两个表如下所示：事故编号纬度/经度 GPS坐标 AB123 纬度 32.123 AB123 长 120.123 CD ...

Leaflet - 根据以米为单位的半径计算缩放级别 - Leaflet - Calculate zoom level from radius in meters

我正在尝试找到一种方法来计算给定半径（以米为单位）的缩放级别。例如，我想在point x,y内显示100km radius ，如何计算该点的适当缩放级别？ ...

如何在 r 中的光栅 map 之上正确组合多边形 map？ - How to properly combine a polygon map on top of a raster map in r?

我想在使用库 ggmap 获得的光栅 map 之上组合美国一些州的形状文件ggmap 。这是我的代码这就是我得到的：还收到以下警告，我不确定这是什么意思：如您所见，US state shapefile 和栅格 map 之间存在不匹配。shapefile map us.state.long ...

从 lapply 的结果中提取正确的列 R - Extract the Correct Column from the Result of lapply in R

我有一个测深栅格和一个边界网格化的 shapefile。我无法隔离 lapply 生成的列表的“值”列。目的是使用栅格平均值列作为 shapefile 中的新列最后一行应该分配 Mean_Raster_Values 列表 object 中的整个第三列，但 [[3]] 只提供第三行如何从 Me ...

如何在知道单个坐标、角度和距离的情况下创建正方形多边形？ - How to create a square polygon with knowing single coordinate, angle and distance?

如何使用 javascript 创建一个只知道单点坐标、每边距离 (3km) 和 90 度角的多边形（正方形）。如果我将拥有所有 4 个点的坐标，那么我可以绘制一个矩形。 . 任何建议都会有很大帮助。 ...

在地理数据框中分配新几何时出错 - error while assigning new geometry in geodataframe

我在尝试在地理数据框中分配新几何体时遇到了一些问题。我的原始地理数据框几何是 LINESTRING Z，我确实将其转换为 LINESTRING，但我无法评估新几何。我的代码：我得到：Value should be either a BaseGeometry or None, got LINES ...

将网格剪辑到 shapefile 但保持单元格完整 - Clip a grid to a shapefile but keep cells whole

我有两个 shapefile：一个六角网格，一个边界 Shapefile：我正在尝试将 Grid 剪辑到 Shapefile 但无论我如何尝试st_intersect() 、 st_overlap() 、 st_contains()等，结果始终相同：我想要的是剪辑网格并保持剩余的单元格完整，如下 ...

在 Openstreetmap 中计算每条道路的分数会产生意想不到的结果。我错过了什么？ - Calculating a score for each road in Openstreetmap produces unexpected result. What am I missing?

我有一个安装了 postgis 扩展的 Postgres 数据库，并填充了 open street map 数据。使用以下 SQL 语句：我计算了一个“绿色”分数。我的目标是为每个 osm_id 创建一个“绿色”分数。意思是; 有多少道路靠近水域、福雷斯特或类似的地方。例如，被公园包围的 ...

堆叠栅格并计算每个像素的最大值，然后在 R 中的 rest 个图层中保留该值 - Stack rasters and calculate each pixel's max value then keep that value for rest of the layers in R

我每天堆叠 181 个 ndvi 光栅。我想为每个像素找到峰值，然后将峰值分配给 rest。原始数据曲线：这就是我要的：此结果适用于单个像素。我已经在 excel 中完成了，但为了方便起见，我想在 R 中创建一个代码。现在每个像素在不同的日期达到峰值（因此堆栈中的不同层）。这是我试过 ...

优化代码以将大向量拆分为较小的文件？ - Optimise code to split large vector into smaller files?

在我下面的玩具/示例代码中，我制作了一个覆盖世界的网格，然后使用此网格将一个大型复杂的全球数据集拆分为每个网格单元格的一个文件。在我的实际工作中，这是一个瓶颈，需要很长时间。我会很感激一些关于优化它的想法和想法。我通过利用并行处理取得了一些成功，但我也认为它可以做得“更聪明”。 ...

我应该在哪里放置 layerList.multipleSelectionEnabled=false; ArcGIS Maps 中的属性 SDK for JavaScript 4.x？ - Where should I place the layerList.multipleSelectionEnabled=false; property in ArcGIS Maps SDK for JavaScript 4.x?

我正在尝试创建一个 web 应用程序，允许用户探索与气候安全社区相关的图层。当我让我的一些同事尝试使用 ArcGIS Online 的开箱即用模板时，他们很难理解他们需要关闭他们想要查看的图层上方的图层才能查看它。我想编写一个 map 代码，它只允许用户一次 select 一层，所以这不会发生 ...

Shapely Buffering，没有按预期工作 - Shapely Buffering, not working as expected

为什么缓冲我的一个几何图形时会出现意想不到的漏洞？通过删除一对坐标，它没有孔。当使用 Sedona 的ST_Buffer时，这种情况在更多情况下发生。 ...

R：将使用纬度/经度定义的 sf 多边形转换为 UTM - R: Convert sf polygon defined using lat/long to UTM

我有以下多边形，使用纬度/经度定义：当我使用 Tmap 在 map 上 plot 时，它出现在不列颠哥伦比亚省海岸外的太平洋中，而不是育空地区中部：我猜问题出在使用纬度/经度而不是 UTM，因为我使用 UTM 定义了其他多边形，这些多边形确实出现在它们（以及上面定义的多边形）应该出现的位置。 ...

在 Python 中使用循环打开和操作一个文件夹中的几个.tiff 文件？确保文件作为单独的变量导入 - Open and manipulate several .tiff files from one folder using a loop in Python? Ensuring the files are imported as separate variables

我有二十个.tif 文件，其中包含来自全局映射 model 的数据。我正在将文件导入 python，这样我就可以旋转它们，并使 .tif 文件以大西洋为中心而不是以太平洋为中心。我可以使用下面的代码旋转单个文件，但是，我想导入所有 20 个文件，旋转它们，然后使用循环将它们全部保存在新位置，而不是 ...

Openlayers GIS - 地理参考 DEM(GeoTIFF) 在 Map (OpenStreetMap) 中定位不正确 - Openlayers GIS - Georeferenced DEM(GeoTIFF) not positioned correctly in Map (OpenStreetMap)

使用 Openlayers，我在 OL Map 的 OSM（作为 WebGLTile 层）之上以GeoTIFF格式（使用 WebGLTileLayer）可视化 GIS 栅格（DEM、正射镶嵌）。我遇到的问题是，具有地理配准的 GIS 栅格未定位在 OSM 基础 map 上的正确地理 positio ...

通过 Apache Beam 读取 GPKG 格式文件 - Reading GPKG format file over Apache Beam

我需要通过 apache beam（数据流运行器）解析gpgk扩展文件并将其加载到 Bigquery 表。我可以看到 beam 具有名为Geobeam的功能，但我看不到有关加载gpgk文件的参考。 Q1：哪个Beam库可以帮我加载geopakage文件？ Q2：作为替代解决方案，我正在尝试将g ...