Questions tagged «python»

我希望能够编写一些未通过arcgisscripting或ArcPy 公开的内容。 如何从Python 访问ArcObjects？

132 python  arcobjects  comtypes 

谁能解释如何在Windows上使用Python安装GDAL / OGR？ 我使用Windows Vista，并且尝试按照网站上的信息进行操作，但似乎没有为我绑定exe文件。 有人可以描述该过程，包括指向我需要的文件/文件夹的链接吗？ 我现在尝试使用minGW运行gdal设置，但这也失败了：

92 python  gdal  installation  windows 

我似乎对几乎所有的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 

我意识到这个网站存在之前，我曾在Stack Overflow上问过这个问题。由于该站点更合适，因此我也会在这里提出要求。 我正在尝试从pip install gdal虚拟环境（Ubuntu）中的pip安装gdal 。失败，因为找不到cpl_port.h extensions/gdal_wrap.cpp:2853:22: fatal error: cpl_port.h: No such file or directory compilation terminated 但是GDAL已正确安装，并且头文件位于/usr/include/gdal/cpl_port.h。是否需要为GDAL设置一些环境变量，以便pip查找头文件？

68 python  gdal  ubuntu  installation 

这是等高线图，所有级别的多边形都可用。 让我们问问如何平滑多边形，使所有顶点保留在其精确位置？ 实际上，轮廓是在网格数据的顶部制作的，您可能会建议然后对网格数据进行平滑处理，因此最终的轮廓将更加平滑。请注意，这并不是我所希望的，因为诸如高斯滤波器之类的平滑功能将删除少量数据，并将更改第三个变量的范围，例如在我的应用程序中不允许的高度。 实际上，我正在寻找一段代码（最好在Python中），该代码可以使2D多边形（任何类型：凸形，凹形，自相交等）平滑化，并且相当容易（忘记代码页）并且准确。 仅供参考，ArcGIS中有一个功能可以完美地做到这一点，但是对于这个问题，我不是选择使用第三方商业应用程序。 1） Scipy.interpolate： 如您所见，样条曲线（红色）不令人满意！ 2） 这是使用此处给出的代码得到的结果。运行不正常！ 3） 对我来说，最好的解决方案应该是如下图所示，其中仅通过更改一个值即可逐渐平滑正方形。我希望有类似的概念可以平滑任何形式的多边形。 满足样条通过点的条件： 4） 这是我在数据上逐行在Python中实现“ whuber的想法” 的方法。由于结果不好，可能存在一些错误。 K = 2是一场灾难，因此对于k> = 4。 5） 我在有问题的位置上删除了一个点，现在所得到的样条曲线与凸纹样条相同。但这仍然是一个问题，为什么该方法不适用于所有情况？ 6） 可以通过以下方式很好地平滑whuber的数据（由矢量图形软件绘制），其中平滑添加了一个额外的点（与更新比较） 4）： 7） 有关某些标志性形状，请参见Python版本的Whuber代码的结果： 请注意，该方法似乎不适用于折线。对于角折线（轮廓），绿色是我想要的，但变成红色。这需要解决，因为等高线图总是多段线，尽管闭合的多段线可以像我的示例中那样视为多边形。还不是更新4中出现的问题尚未解决。 8）[我的最后一个] 这是最终的解决方案（不完美！）： 请记住，您将不得不对恒星所指向的区域执行某些操作。我的代码中可能存在错误，或者建议的方法需要进一步开发以考虑所有情况并提供所需的输出。

52 python  open-source-gis  contour  smoothing  spline 

我认为这是一个基本问题，但我似乎找不到或认识到解决方案。 该网站返回 Point: X: -11705274.6374 Y: 4826473.6922 以第一个键值000090为例进行搜索时。我想这是一个空间参考，我有点明白了。 我正在寻找有关如何使用Python将其转换为纬度和经度的说明或示例。

50 python  coordinate-system  convert 

您认为学习GIS目标的最佳书籍/网站是什么？ “最好”是指： 不太长（书） 易于理解（书籍/网站） 良好的实际例子（书/网站） 对于特定于学习如何使用Python的ArcPy网站包来自定义ArcGIS for Desktop的答案，有一个单独的问答：学习ArcPy有哪些资源？

49 python  references 

我正在研究整个景观中大量野生传粉媒介的计算模型。该模型本身已经完成，现在我正在努力进行后处理步骤。 我的GDAL授粉媒介供应栅格看起来像这样（较浅的颜色表示授粉媒介对像素的访问量较高）： 我有一个OGR形状文件，该文件的点代表景观上的样本位置： 我正在尝试对这些点下的像素进行一些分析，但是要这样做，我需要能够提取一个点下的像素的值。 是否可以通过Python仅使用OGR和GDAL提取点下的像素值？我宁愿避免通过读取整个栅格到内存中ReadAsArray()，因为我的输出栅格非常大（太大而无法容纳到内存中）。 我注意到这篇文章，这很相似，但是需要命令行调用。

45 python  gdal  ogr  pixel 

我想知道是否可以在没有获得ArcMap许可的情况下使用Python查看shapefile的内容。这种情况是，您不仅可以通过ESRI软件，还可以通过许多不同的应用程序创建shapefile。我想创建一个Python脚本，该脚本检查空间参考，要素类型，属性名称和定义以及shapefile中字段的内容，并将它们与一组可接受的值进行比较。即使组织没有任何ESRI许可证，我也希望此脚本能够运行。为此，您是否必须使用ArcPy或不使用ArcPy就能挖入shapefile？

40 python  shapefile 

是否可以使用某些gdal API gdal_translate从Python代码调用？我的意思不是简单地从文件系统执行gdal_translate.exe，而是以某种方式在代码中调用它，所以我不需要知道gdal_translate可执行文件所在的确切目录？

40 python  gdal 

我希望在不使用ArcGIS / ArcPy的情况下使用Python可视化地理数据并制作地图。 在互联网上，我发现了如何使用Python制作专题图： 这是我尝试过的一些代码： import shpUtils import matplotlib.pyplot as plt shpRecords = shpUtils.loadShapefile('C:\\Users\\shapefile.shp') for i in range(0,len(shpRecords)): x = [] y = [] for j in range(0,len(shpRecords[i]['shp_data']['parts'][0]['points'])): tempx = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['x']) tempy = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['y']) x.append(tempx) y.append(tempy) plt.fill(x,y) plt.axis('equal') plt.title("Testing") plt.show() 但是，当我运行它时，它会给我随机的颜色。 如果我想可视化shapefile的某一列，如何使用类似的代码实现呢？ 在上面提供的链接中，这很不清楚，他只讨论颜色的使用... 我是否可能需要其他模块来完成此任务，例如笛卡尔和PySAL？

37 python  shapefile  visualisation  matplotlib 

我已经在#qgis和#postgis之间的stackoverflow和irc上多次问过这个问题，我也尝试在postgis中对其进行编码或实现，但没有真正的答案。 我想使用编程（最好是python），从点层绘制一条线，使其在直线或多边形层的最近线上投影。 截至目前，我的大部分数据都采用ESRI的格式和邮政格式。但是，我宁愿远离postgis解决方案，因为我主要是shp + qgis用户。 一个理想的解决方案是用python或类似的库实现GDAL / OGR 从哪里开始使用GDAL / OGR库？有可能给出解决方案吗？ 我可以使用NetworkX进行最近邻居分析吗？ 这实际上可行吗？ 如果更简单，则这些点可以连接到线段终点，而不是投影点

37 python  postgis  vector  proximity 

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

33 python  geoprocessing 

我不清楚如何最好地使用in_memoryArcGIS中的工作区。我找到的关于该主题的最佳资源来自ESRI的“ 使用内存中的工作空间”。不幸的是，帮助部分的信息不是很多。另外，此网站上有一篇帖子，简要介绍了此处和此处的主题。 我的具体问题： in_memory在基于ArcGIS / arcpy的脚本中使用工作区的正确用法和语法是什么？ 是in_memory工作区中的相同，例如，创建使用层arcpy.MakeFeatureLayer_management()？ 是否有任何标准，例如in_memory在脚本末尾删除工作区？

33 arcgis-desktop  python  arcpy  in-memory 

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

33 python  geoprocessing  open-source-gis  fiona