Questions tagged «python»

在ArcGIS 10中，Python已成为ESRI的事实上的脚本语言。虽然领先的开源软件（例如QGIS，GeoServer，PostGIS）都支持Python。因此，对于GIS开发人员和用户来说，了解/学习Python变得很重要。 谁能列出GIS中至关重要的Python这样的工具/模块/插件吗？

31 python  software-recommendations 

Kyngchaos安装程序未找到python3。我使用HomeBrew安装了python3，但位于/usr/local/Cellar/python/3.6.4_4/Frameworks/Python.framework/Versions/3.6/bin/python3.6 如何告诉QGIS安装程序识别已安装python3？ 我以前曾使用HomeBrew安装GDAL，因此尝试重新安装（或在可能结束的地方安装其他版本），认为这可能有助于解决这种情况，但QGIS3安装程序仍然找不到python 3。 安装错误是在安装程序打开后弹出的错误消息，提示“ python3 QGIS需要Python 3.6”。 对于非Mac用户，大多数人似乎都依赖Kyngchaos提供QGIS安装程序，而他刚刚发布了v3安装程序。 更新：在下面回答了我自己的问题。除非您有足够的知识将QGIS安装程序/ QGIS指向正确的版本（我当然不是），否则不要使用HomeBrew，而要使用自述文件中链接到的python3安装程序版本。

30 qgis  python  installation  mac  qgis-3 

有没有办法使用gdal的Python绑定从栅格文件中获取角坐标（以纬度/经度为单位）？ 网上的一些搜索使我确信没有这种搜索，因此我通过解析gdalinfo输出进行了开发，它有些基本，但是我认为这可能为那些不太熟悉python的人节省一些时间。仅当gdalinfo包含地理坐标和角坐标时，它也才有效，我认为并非总是如此。 这是我的解决方法，请问有人有更好的解决方案吗？ 在适当的栅格上的gdalinfo会在输出中途产生如下所示的结果： Corner Coordinates: Upper Left ( -18449.521, -256913.934) (137d 7'21.93"E, 4d20'3.46"S) Lower Left ( -18449.521, -345509.683) (137d 7'19.32"E, 5d49'44.25"S) Upper Right ( 18407.241, -256913.934) (137d44'46.82"E, 4d20'3.46"S) Lower Right ( 18407.241, -345509.683) (137d44'49.42"E, 5d49'44.25"S) Center ( -21.140, -301211.809) (137d26'4.37"E, 5d 4'53.85"S) 此代码将在gdalinfo看起来像这样的文件上工作。我相信有时坐标将以度和小数表示，而不是度，分和秒。为这种情况调整代码应该是微不足道的。 import numpy as np import subprocess …

30 python  gdal 

使用envi图像处理和分类工具，您可以从具有某些光谱值的图像中获取屋顶，然后可以将其转换为应用程序的一些矢量数据。 在Intel开发的具有OpenCV的 python中（具有许多对象检测算法），您可以从图像中检测面部。 OpenCV示例： 我的问题是，我们可以使用python开源工具从协调或不协调的卫星图像中检测屋顶或任何东西吗？ 卫星图像示例：

30 python  remote-sensing  satellite  envi 

我想发展我的GIS编程技能，从哪里开始呢？人们说学习C ++或Python，但是我在哪里可以在地理环境中学习？教程或其他任何东西都将非常有用，有关使用的语言/程序的任何信息也将非常有帮助。

30 python  c++  references 

我是GIS的新手。 我有一些代码将火星的红外图像转换成热惯性图，然后将其存储为2D numpy数组。我一直将这些地图另存为hdf5文件，但我真的很想将它们另存为栅格图像，以便可以在QGIS中对其进行处理。我经历了多次搜索，找到了如何执行此操作的方法，但是没有运气。我尝试按照http://www.gis.usu.edu/~chrisg/python/上的教程中的说明进行操作，但是当我将其示例代码生成的文件导入QGIS时，它们以纯灰色框打开。我觉得如果有人可以将最简单的过程建议为我想做的事的简化示例，那么我也许可以取得一些进展。我有QGIS和GDAL，很高兴安装任何人都可以推荐的其他框架。我使用Mac OS 10.7。 因此，例如，如果我有一个numpy的热惯性数组，看起来像： TI = ( (0.1, 0.2, 0.3, 0.4), (0.2, 0.3, 0.4, 0.5), (0.3, 0.4, 0.5, 0.6), (0.4, 0.5, 0.6, 0.7) ) 对于每个像素，我都有纬度和经度： lat = ( (10.0, 10.0, 10.0, 10.0), ( 9.5, 9.5, 9.5, 9.5), ( 9.0, 9.0, 9.0, 9.0), ( 8.5, 8.5, 8.5, 8.5) ) …

30 python  raster  gdal  numpy 

我具有多边形特征，并希望能够在其中生成点。我需要一项分类任务。 生成随机点直到多边形中的某个点不起作用，因为这确实花费了很多时间。

30 python  qgis  polygon  ogr 

我正在考虑通过ESRI 进行python地理处理课程来提高自己的编程技能。在我付高额的注册费之前，我想请您对24小时讲师指导的python编程的价值提出意见。我主要对学习以下内容感兴趣： 对于那些参加过此类课程的人，您是否拥有适用的技能或更多的理论知识？ 对于最近进入就业市场的那些人，您是否发现雇主在履历表上开这样的一类的工作看起来不错？ 对于雇主来说，像这样的专业课程真的很重要吗，还是大多数GIS分析人员/专家都自己学习python编程？

29 python  education 

我正在寻找有关如何提高我的python代码效率的一些建议。通常，效率对我来说并不重要，但是我现在正在处理一个美国地点超过150万点的文本文件。使用给定的设置，在一点上运行操作大约需要5秒钟；我需要把这个数字降下来。 我正在使用三个不同的python GIS软件包对这些点进行一些不同的操作，并输出一个新的带分隔符的文本文件。 我使用OGR读取县边界shapefile并访问边界几何。 匀称检查某个点是否在这些县中的任何一个县内。 如果在1之内，则使用Python Shapefile库从边界.dbf中提取属性信息。 然后，我将两个来源的一些信息写入文本文件。 我怀疑效率低下在于存在2-3层的循环...不太清楚该怎么做。我特别希望与有经验的人一起使用这3个软件包中的任何一个，因为这是我第一次使用它们。 import os, csv from shapely.geometry import Point from shapely.geometry import Polygon from shapely.wkb import loads from osgeo import ogr import shapefile pointFile = "C:\\NSF_Stuff\\NLTK_Scripts\\Gazetteer_New\\NationalFile_20110404.txt" shapeFolder = "C:\NSF_Stuff\NLTK_Scripts\Gazetteer_New" #historicBounds = "C:\\NSF_Stuff\\NLTK_Scripts\\Gazetteer_New\\US_Counties_1860s_NAD" historicBounds = "US_Counties_1860s_NAD" writeFile = "C:\\NSF_Stuff\\NLTK_Scripts\\Gazetteer_New\\NewNational_Gazet.txt" #opens the point file, reads …

29 python  shapefile  ogr  shapely  code-review 

我有一个DEM栅格，像素值在3000到-0.0003之间。我需要将所有值为0 且更小（<= 0）的像素设置为“ nodata”（用于以后的栅格计算并减小文件大小）。 对我来说，这似乎是栅格计算器的一项任务（像素<= 0的设置值到“ nodata”）或重新分类（像素<= 0的设置值到“ nodata”并保留所有其他值），但是我没有知识。

28 qgis  python  grass  saga 

左侧的像素表示树的位置及其关联的树冠半径（即像素值范围为2-5）。我想通过它们的冠状半径值来缓冲这些栅格像素。右边的图像是我希望仅使用栅格处理方法来完成的图像。 我最初会考虑在ArcGIS中使用圆形焦点总和，尽管邻域设置是一个固定值，但不会考虑可变大小的冠状半径。 用像素值“缓冲”像素的好方法是什么？

28 python  arcgis-desktop  raster  arcpy  buffer 

我试图在新计算机上运行Shapely，并尝试从源代码安装它，使用其网页上的安装程序并使用pip install shapely。我尝试从库中进行的每个引用都会from shapely import *引发错误WindowsError: [Error 126] The specified module could not be found。在尝试利用Shapely并弄清楚如何使其正常运行时，是否有人有类似的经验？ 引发的Python错误： WindowsError: [Error 126] The specified module could not be found 环境： Windows 7 64位 Python 2.7

26 python  python-2.7  shapely 

有人可以演示一种简单的方法将几何数据结构从形状写入shapefile吗？我对带有孔和线串的多边形特别感兴趣。远离arcpy也将是有益的（因此osgeo，pyshp等都会更好）。

26 python  shapefile  shapely 

我一直在尝试检查我在DEM栅格上的过滤器以进行模式识别，并且它总是导致缺少最后一行和最后一列（如..20）。我已经尝试了PIL库，图像加载。然后用numpy。输出是相同的。 我以为我的循环出了点问题，当检查数组中的值时（只是用ArcCatalog中的Identification选择像素），我意识到像素值未加载到数组中。 因此，只需打开，放入数组并保存数组中的图像即可： a=numpy.array(Image.open(inraster)) #raster is .tif Float32, size 561x253 newIm=Image.new(Im.mode, Im.size) Image.fromarray(a).save(outraster) 导致删除最后的行和列。抱歉，无法发布图片 任何人都可以帮助理解为什么？并建议一些解决方案？ 编辑： 因此，在伙计们的帮助下，我成功地将小型栅格加载到numpy数组中，但是当图像较大时，我开始出现错误。我想这是关于numpy数组的限制，因此数组会像这样自动重塑或变形……所以ex： Traceback (most recent call last): File "<pyshell#36>", line 1, in <module> ima=numpy.array(inDs.GetRasterBand(1).ReadAsArray()) File "C:\Python25\lib\site-packages\osgeo\gdal.py", line 835, in ReadAsArray buf_xsize, buf_ysize, buf_obj ) File "C:\Python25\lib\site-packages\osgeo\gdal_array.py", line 140, in BandReadAsArray ar = numpy.reshape(ar, [buf_ysize,buf_xsize]) File …

26 python  raster  gdal  numpy  floating-point 

我已经使用Osgeo安装程序安装了GDAL。如何以编程方式裁剪栅格图层和矢量图层？有可用的GDAL API可以帮助我吗？我正在使用Python。

26 python  gdal