地理信息系统

我有大量的多边形代表美国各县的边界，我需要将它们合并成尽可能少的多边形。 是否可以使用非商业软件或价格合理的软件来做到这一点？ PostGIS不可用，Java也不可用，但是C可用。

16 ogr2ogr  aggregation 

我想使用一系列单边界多边形在土壤多边形的R中创建一个真正的GIS剪辑，但是我找不到能够正确执行此操作的R函数。它应该像clipESRI的ArcMap中的功能一样工作。我已经尝试过该over方法在sp包中，但它似乎不适用于多边形而不是多边形。 一种建议是使用以下代码将gIntersectionin rgeos包用作剪辑： #------------------------------------ library(rgeos) library(maptools) #Read layers as SpatialPolygonsDataFrame (both the same Albers projection) Soils_poly = readShapePoly("Soils_polygons") #Note - Has 400 polygons clipper_poly = readShapePoly("clipper_polygon") #Note - Has 1 polygon #Try gintersection as clip Clipped_polys = gIntersection(Clipper_Tile_poly, Soils_poly) #----------------------------------- 这需要5分钟才能运行（速度太慢），并且会出现以下错误： RGEOSBinTopoFunc中的错误（spgeom1，spgeom2，byid，id，drop_not_poly，“ rgeos_intersection”）：TopologyException：在-721459.77681285271 2009506.5980877089中找不到传出的dirEdge 我还尝试了以下代码来检查是否重叠： gIntersects(Clipper_Tile_poly, Soils_poly) 结果为TRUE。clipESRI ArcMap中的函数对此数据可以正常工作。 有人知道R函数可以使用空间多边形对空间多边形进行适当的剪辑吗？

16 r  polygon  layers  vector  clip 

我对坐标系统不太了解... 在我的办公室里，我们用来处理来自考古现场的空间数据。每个站点都有自己的xyz坐标系（GCS）。三个简单的正交笛卡尔轴。在过去的几年中，我们一直在通过GIS软件（ArcGIS）来管理此空间数据，而无需使用特定的坐标系（只需将其保留为“未定义”即可） 我想知道是否存在设计用于使用简单笛卡尔正交轴处理此类数据集的任何GCS，而没有典型GCS的网格失真。另外，我想知道这个系统是否适合在在线地图应用程序中使用。 顺便说一下，我们管理2D（ArcMap）和3D（ArcScene）环境，并使用“ mm”作为长度基础单位。 如果不存在这样的东西，也许有人知道如何创建它。

16 arcgis-desktop  coordinate-system  epsg  srid 

我想知道是什么触发了pyQGIS中文件锁定的释放？ 我试图通过调用删除一些数据源（临时使用）QgsVectorFileWriter.deleteShapeFile，但是必须退出QGIS才能这样做。我已经将源代码加载到QgsVectorLayer对象中。在删除源之前，是否必须对所有这些对象及其引用进行垃圾回收？有没有办法强制执行此操作？ 我设法创建了一个失败的最小代码示例。在运行之前，请确保temp目录为空。 from qgis.core import * import processing, os, gc project_temp_dir = "C:/Path/To/My/Dir/" layer1_path = project_temp_dir + "layer1.shp" layer2_path = project_temp_dir + "layer2.shp" input_layer = QgsMapLayerRegistry.instance().mapLayersByName('in_layer')[0] if not input_layer.isValid(): raise Exception("Failed to grab input layer") # Create layer 1 err = QgsVectorFileWriter.writeAsVectorFormat(input_layer, layer1_path, "utf-8", input_layer.crs()) if err != QgsVectorFileWriter.NoError: …

16 pyqgis  lock 

早期的问答（Esri的新“专业用ArcGIS”和“ ArcGIS Pro”是什么？）使人确信ArcPy代码将在ArcGIS Pro中具有重用价值： ArcGIS Professional可扩展性。这是我们目前对ArcGIS Pro可扩展性的想法。您可以在ArcGIS中编写和运行Python脚本，以调用地理处理工具，并使用ArcGIS Python API ArcPy中提供的详尽的脚本功能套件来自动执行GIS任务。 但是，它没有提及ArcGIS Pro是否将使用（等同于）现有ModelBuilder或对其进行重大改进以解决其他人和我以前提出的增强要求的GUI，来使模型构建能够使工作流程自动化。 有人听说过ArcGIS Pro中的ModelBuilder（或至少是模型构建）的道路吗？

16 arcgis-desktop  modelbuilder  arcgis-pro  arcgis-10.3 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 9个月前关闭。 除ArcPy之外，是否有任何可以与shapefile进行地理处理的python库（例如缓冲区/相交）？

16 python  buffer  software-recommendations 

我正在寻找生成原产地（OD）成本矩阵所需的步骤。 我正在寻找的是道路距离矩阵。距离矩阵工具仅产生欧氏距离。QGIS Roadgraph插件产生从一个起点到N个目的地的最短道路距离。我需要的是一个文件，其中包含从所有起点到N个最近终点的距离。我有很多渊源，因此无法使用Roadgraph重复分析。再次感谢。 我正在寻找从头开始的逐步说明（即，使用我的3个shapefile：起点，终点，道路网络）。

16 qgis  routing  network  cost-path 

我在为传单创建自定义控件方面获得了一些帮助，但是在升级到最新版本的Leaflet所使用的CSS样式方面仍然遇到困难。基本上，我所需要做的只是在白框中显示一个图标，并在其周围添加阴影。 到目前为止，这是我一直在努力的工作，请随意看一下： http://codepen.io/DrYSG/pen/zJsiG

16 leaflet 

如何使用python在Qgis中按特征（类似于arcobjects中的Iqueryfilter）的属性获取特征？除了获取所有功能并手动对其进行过滤之外，是否可以使用where子句将其过滤掉？ 示例：我有一个名为“县”的字段。它具有五万多个特征。即，由于耗时，无法获取所有特征并对其进行过滤。所以我可以在arcobjects中使用iqueryfilter.whereclause ='Counties = Norwich'来查询它。我在PyQgis中需要类似的东西。

16 qgis  pyqgis 

我对学习如何利用NumPy数组优化地理处理感兴趣。我的大部分工作都涉及“大数据”，其中地理处理通常需要几天的时间才能完成某些任务。不用说，我对优化这些例程非常感兴趣。ArcGIS 10.1具有许多可通过arcpy访问的NumPy函数，包括： NumPyArrayToFeatureClass（arcpy.da） RasterToNumPyArray（arcpy） TableToNumPyArray（arcpy.da） 出于示例目的，假设我想利用NumPy数组优化以下处理密集型工作流程： 这里的总体思想是，有大量的基于矢量的点在矢量和基于栅格的操作中移动，从而生成二进制整数栅格数据集。 如何合并NumPy数组以优化此类工作流程？

16 arcgis-desktop  python  arcpy  numpy  big-data 

我有一个简单的问题：我想计算一组多边形中的点数。 我已经有一个SQL，但它只返回实际包含点的多边形的gid。 我的表格：包含19.000行的多边形层和包含450行的点层。 下面的SQL select grid.gid, count(*) AS totale FROM grid, kioskdhd3 WHERE st_contains(grid.geom,kioskdhd3.geom) GROUP BY grid.gid; 仅返回一些实际包含点的320个多边形。但是我想返回所有多边形，甚至以为点数为0。 当然，这与我的WHERE条款有关。我必须在哪里放置我的st_contains？ 谢谢卡斯滕

16 postgis  point-in-polygon 

在Leaflet或OpenLayers中是否可以将两个地图与一个在两个之间滑动的条形图进行比较？ 示例包括： http://www.esri.com/services/disaster-response/hurricanes/hurricane-sandy-the-aftermap http://www.arcgis.com/apps/StorytellingSwipe/index.html?appid=c62112f6ba834227978d8f1cf2b6f5af http://www.abc.net.au/news/specials/qld-floods/（不是地图） 我知道这是可能的，但我不知道该怎么做。任何人有任何想法吗？例子？

16 openlayers-2  web-mapping  leaflet  swipe 

我有一个ArcGIS 10.1 Server实例，该实例在Internet上公开了受保护的地图服务。我需要编写一个客户端应用程序（我目前正在使用ArcGIS Javascript API的3.3版本来构建）的代码，使用户可以查看那些受保护的Web服务： 我认为这个在线ESRI示例是一个好的开始。 我的意愿是不提示用户进行EACH mapservice身份验证，因为我已经知道所有mapservices都属于她/他，因此她/他可以使用相同的用户名和密码进行访问。在我的想法中，凭据提示仅应显示一次，因此JS代码应通过某种RESTful登录调用将凭据提供给每个mapservice。在我看来，ArcGIS Server REST API似乎没有提供这样的调用...也许我错了。 因此，这种使用ArcGIS Server登录到安全地图服务的“ RESTful”方式是否可行（从而可以通过编程方式访问安全服务）？如果是这样，你们中的任何人都可以提供示例或指向Web资源的链接对此进行解释吗？

16 arcgis-10.1  arcgis-rest-api  security  authentication  arcgis-javascript-api 

有什么网站可以让我从Google Road网络获取shapefile吗？据我所知，它不是由Google分发的。因此，任何人都可以建议我如何创建此路网以用于自定义用途？

16 shapefile  google-maps  download 

有人可以提供什么是“非节点相交”的精确定义吗？有各种各样的GIS平台和工具包对此文本有错误，并且针对此的Google搜索结果往往属于“降低精度，所有问题都消失了”的类别，而没有任何真正的错误原因的解释。给定两个有效的几何形状，为什么会发生此错误？ 参考文档将是最有帮助的。

16 postgis  nettopologysuite