Questions tagged «geoprocessing»

我似乎对几乎所有的python地理处理都使用ESRI的Arcpy站点软件包。值得ESRI赞扬的是，这些工具令人难以置信，可帮助您完成很多工作。但是，我也想在ESRI Arcpy域之外创建地理处理脚本。例如，如果我想将栅格裁剪为多边形，我将从ESRI的以下脚本开始： # Import system modules import arcpy from arcpy import env from arcpy.sa import * # Set environment settings env.workspace = "C:/sapyexamples/data" # Set local variables inRaster = "elevation" inMaskData = "mask.shp" # Check out the ArcGIS Spatial Analyst extension license arcpy.CheckOutExtension("Spatial") # Execute ExtractByMask outExtractByMask = ExtractByMask(inRaster, inMaskData) …

69 python  geoprocessing 

在我看来，这似乎是一个简单的问题（可能是），但我似乎找不到能给我答案的例子。使用PostGIS，我只想选择不在多边形范围内的点。就我所知，这最终是ST_Intersects的逆函数。 示例：我有一个taxlot层和一个地址点层。我假设我应该使用ST_Intersects，但是如何告诉它进行反向选择？我以为可能在下面的代码前添加一个NOT语句，但这没有用。 CREATE table t_intersect AS SELECT hp.gid, hp.st_address, hp.city, hp.st_num, hp.the_geom FROM public.parcel as par, public.housepoints as hp WHERE ST_Intersects(hp.the_geom,par.the_geom);

41 postgis  geoprocessing  sql  intersection 

我有Python代码，旨在通过以下工作流程获取点shapefile： 合并点 积分点，使彼此之间1 m之内的任何点成为一个点 创建要素图层，在其中选择z <10的点 缓冲点 多边形到栅格的1m分辨率 重新分类，其中1-9 = 1；NoData = 0 每个shapefile大约有250,000至350,000个点，覆盖〜5x7 km。用作输入的点数据表示树的位置。每个点（即树）都有一个关联的“ z”值，该值代表冠半径，并在缓冲过程中使用。我的目的是在单独的过程中使用最终的二进制输出来生成描述遮篷的栅格。 我对四个shapefile进行了测试，生成了700MB的栅格，耗时35分钟（i5处理器和8GB RAM）。鉴于我将需要在3500个shapefile上运行此过程，因此，我希望获得有关简化该过程的任何建议（请参阅随附的代码）。一般来说，处理地理数据的最佳方法是什么？更具体地说，对代码或工作流程是否有任何有助于提高效率的调整？ 编辑： 地理处理任务的时间（占总数的百分比）： 合并= 7.6％ 积分= 7.1％ Lyr = 0的特征 缓冲= 8.8％ 多边形转栅格= 74.8％ 重新分类= 1.6％ # Import arcpy module import arcpy # Check out any necessary licenses arcpy.CheckOutExtension("spatial") # Script arguments temp4 …

34 arcpy  geoprocessing  big-data 

我想使用一个不是基于arcpy的Python脚本来执行诸如通过属性查询shapefile，通过选择创建新图层，计算多边形的面积并将多边形转换为点之类的事情。 有没有使用其他Python模块或库的代码示例？我可以使用arcpy轻松做到这一点，但我想探索其他选择。

33 python  geoprocessing 

假设我将一个Shapefile放在一起，所有要素的顶点都移动了恒定量。转移所有特征的最简单方法是什么（因此（x，y），其顶点位置）进行任意偏移？我有很多文件要应用此更正，因此首选Bash / OGR答案:) 最后，我最终使用了Spatialite，因为它具有不错的功能ShiftCoords。但是，该线程非常有用！谢谢大家！

33 python  geoprocessing  open-source-gis  fiona 

这篇文章是社区维基。编辑现有答案以改善此职位。它当前不接受新的答案。 我正在寻找现有的解决方案（软件，库，方法等），以实现在那里以及那里呈现的地理数据的一般化/简化。 您知道并使用一些工具来简化您的地理数据吗？ 我已经看过这些： maptimize以简化KML标记层的表示， ESRI的Modelbuilder 激光扫描的半径清晰度 mapshaper来计算几何上的在线过滤算法， 通过轴系统进行轴伸 在上概括和多重代表网站ICA佣金这个主题一个巨大的信息源。许多论文提出了与此主题相关的高级方法和用例。 的非常著名的道格拉斯＆普克滤波算法通常用于几何简化和存在于几乎所有的GIS软件。但是，其效果是使几何结构非常骨/角，而不真正简化几何形状。例如，它不适用于必须平滑且没有角度的代表道路或河流的线。

25 geoprocessing  simplify  generalization 

假设我的局域网上有50台计算机。每台计算机都有一个地理数据库，用于存放美国特定州的所有宗地多边形。 我想编写一个地理处理任务，查找所有价值超过x $ /英亩的地块，该地块在另一个价值小于z $ /英亩的地块的y英尺内。 我想制定和运行此查询，而又不知道或不关心数据分布在50台计算机上。请记住边界条件：我还希望查询返回一种情况，即一种状态下的昂贵包裹接近另一种状态下的廉价包裹。 是否有支持这种分布式地理处理的体系结构？ 可以抽象地描述该体系结构，也可以将其描述为Azure或Amazon Web Services的特定实现。或者，最好是作为典型的办公室，在该办公室中，计算机在夜间拥有大量ArcGIS桌面许可证而处于闲置状态。

24 geoprocessing  cloud-gis  amazon-web-services 

我一直在使用FUSION和命令行FUSION Lidar Toolkit（LTK）来处理LiDAR数据。广泛的Google搜索（“ Lidar Python”）产生了libLAS和pyLAS作为Python LiDAR库，但是，这些库似乎仅提供对LAS数据的读写访问。 除了来自点云的冠层曲面模型之外，我对创建强度和密度图像特别感兴趣。Python中是否有一套公认的工具可以完成FUSION LTK能够完成的同类任务？

22 python  geoprocessing  lidar  las  fusion 

GIS的核心概念是回答有关数据集的问题。从数据库的角度来看；具有空间扩展的SQL是询问此类问题的一种方式。可以通过基于机器可读文本的形式来表达问题吗？不同方法的好处是什么？

22 geoprocessing  analysis 

许多其他领域存在视觉/“空间”问题，可以使用GIS工具来解决。我知道我已经使用GIS工具来完成奇怪的非地理任务，并且想知道这会发生多少，以及是否有人在此方面有很好的例子。 我感兴趣的一种例子： 使用地理信息技术量化海洋叠层岩沉积颗粒内的内孔钻探过程的空间结构 [PDF] 总结方法： 用显微镜拍摄了Stromatolite（一种岩石）的岩屑图像。这些图像被导出为BMP，然后用Erdas Imagine进行监督分类。将分类导出到ArcView 3.x中，进行清理并将合适的分类转换为多边形。然后他们做了一些空间选择和一些密度映射。 论文中的数字： 如您所见，我正在寻找使用非地理数据但仍然具有空间元素的问题。GIS方面可以是任何东西：制图，可视化，分析等。 科学论文/研究特别感兴趣，但不是强制性的。

18 geoprocessing  cartography  analysis  visualisation  gis-principle 

我目前正在使用LANDSAT数据参加遥感课程。 我课程的一部分内容是基于辐射值对场景进行标准化。我认为由于这是一项繁琐的任务，因此有人已经想出了一种自动化的方法。 是否有ESRI或其他工具可以标准化LANDSAT场景？

17 geoprocessing  grass  remote-sensing 

关于如何处理GIS中的“时间”属性，是否有任何方法或想法？例如，有两种方法可以指定一个对象位于给定时间间隔内的某个位置，而又可以指定它位于另一个时间间隔内的另一个位置，比如说政治边界发生了变化。

16 geoprocessing  data  time  modelling  spatio-temporal-data 

这应该是另一个问题的显而易见的先驱（我没有问过）： 如何在PostGIS中创建蜘蛛图（集线）？ 如果我不知道A层中的一个点（商店）和B层中的一个点（客户）之间的关系，我通常要说“最近的商店为客户1提供服务”。虽然我意识到这个事实可能并非如此，但它可能是一个不错的替代品。 使用PostGIS，最有效的方法是将A层（商店）中最近的点的ID分配给B层（客户）中的每个点。我正在寻找的输出如下所示。 Customer | Store 1 | A 2 | A 3 | B 4 | C

15 postgis  geoprocessing  business  nearest-neighbor 

我们拥有esri个人地理数据库，需要将其导入PostGIS。 最终，我们希望允许人们将MDB上传到Web服务器上，并将其导入到PostGIS中。 在Windows上，使用FW工具，以下工作正常： ogr2ogr -f "PostgreSQL" PG:"host=10.10.10.10 user=someuser dbname=poi password=somepassword port=5432" S:\GISData\Test.mdb -a_srs EPSG:26986 但是，我想使它自动化（并以编程方式进行），因此[我认为]在Linux上要容易得多。

15 postgis  open-source-gis  geoprocessing  esri-geodatabase 

在FME中，是否可以将多个读取器连接到单个变压器？我正在使用FME 2015，并希望将一堆要素类剪辑为一个要素，但是似乎找不到一种简便的方法来实现。 编辑：试图选择多个功能类型的输入数据集，并将它们全部链接到一个单一的变压器，而无需点击功能类型，并链接到变压器上的输入箭头单独的输出箭头。 编辑2：这是一个屏幕截图，希望这会有所帮助：

14 geoprocessing  fme  clip