Questions tagged «python»

我正在寻找针对空间统计问题的另一种更优雅的解决方案。原始数据由每个单独树的xy坐标组成（即，转换为点.shp文件）。尽管在此示例中未使用，但是每棵树还具有代表顶冠直径的相应多边形（即.shp）。左侧的两个图像显示了从单个树位置的点.shp文件得出的景观尺度内核密度估计（KDE），一个是1989年，另一个是2009年。右图显示了两个KDE之间的差异。仅显示平均值的+/- 2标准偏差的值。使用Arc的栅格计算器执行在右侧图像上生成栅格叠加所必需的简单计算（2009 KDE-1989 KDE）。 有没有更合适的方法以统计或图形方式分析树木密度或冠层面积随时间的变化？有了这些数据，您将如何评估地理空间环境中1989年和2009年树木数据之间的变化？鼓励使用ArcGIS，Python，R，Erdas和ENVI中的解决方案。

15 arcgis-desktop  python  r  remote-sensing  spatial-statistics 

我正在尝试使用Shapely（给定的多边形A，B和C）捕获下面指示的所有非重叠多边形。而且，我希望这样做而无需迭代，测试相交等。这个问题的公认答案表达了PostGIS方法，但似乎“工会”对不同的人而言意味着不同的事情。

15 python  polygon-creation  shapely 

有谁知道纯Python（没有C或Java依赖项）库来执行诸如缓冲区，联合，相交等几何运算？我正在寻找与Google AppEngine兼容的产品。

15 python  open-source-gis 

我在Python中使用GDAL导入shapefile。我想知道文件的字段名称，我当前的方式是： fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field) 但是通过这种方式，我得到了第一层的功能。这是否意味着不同的图层可能具有不同的功能？ 如果不是，是否可以立即获取字段名称，而不是深入了解？如果是，是否有更简单的方法来获取字段名称？

15 python  shapefile  gdal 

我已经开始学习如何在python中处理LAS数据，并希望了解其他人如何处理LAS文件。我想阅读要点（我使用的是numpy数组），并将1类和2类（未分类和分类）过滤到单独的数组中。我有以下代码，但似乎无法过滤点。 # Import modules from liblas import file import numpy as np if __name__=="__main__": '''Read LAS file and create an array to hold X, Y, Z values''' # Get file las_file = r"E:\Testing\ground_filtered.las" # Read file f = file.File(las_file, mode='r') # Get number of points from header num_points = int(f.__len__()) …

15 python  lidar  numpy  array  liblas 

尝试在Windows XP上为Python 2.7.3导入gdal： >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <mo _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swit_helper _mod = imp.load_module('_gdal', fp, pathname, description) ImportError: DLL …

15 python  gdal  windows-xp 

是否有可能在python中按功能分割一个shapefile？（最好是将临时生成的矢量对象临时保存到内存而不是磁盘的解决方案）。 原因：我想将gdal rasterizeLayer函数与shapefile的几个不同子集一起使用。该函数需要一个osgeo.ogr.Layer对象。 mkay，我尝试了一下，它可能如下工作。您可以按以下方式获取每个要素的gdal图层对象的几何形状。 # Load shape into gdal shapefile=str(vectorPath) layer_source = ogr.Open(shapefile) lyr = layer_source.GetLayer(0) for i in range(0,lyr.GetFeatureCount()): feat = lyr.GetFeature(i) ge = feat.geometry() 现在，我只需要知道如何基于此几何形状创建osgeo.ogr.layer对象。 为了澄清。我需要使用纯Ogr / gdal代码编写的函数！这似乎也引起了其他人的兴趣，我仍然希望没有任何辅助模块的解决方案（尽管来自此处的任何解决方案都将在免费的可用qgis插件中使用）。

15 python  gdal  rasterization  osgeo  split-by-attribute 

给定的数据点具有经度，纬度和该点的第​​三个属性值。如何基于属性值将点聚类为组（地理子区域）？我在Google上搜索后发现，这个问题似乎被称为“空间约束聚类”或“区域化”。但是，我对处理地理数据并不熟悉，也不清楚如何使用哪种算法，以及哪种python / R软件包适合此任务。 为了更直观地了解我想要的东西，假设我的数据散布图如下： 因此，每个点都是一个点，x是经度，y是纬度，并且颜色图显示值是大还是小。我想根据位置和值的相似性将这些点划分为子区域/组/集群。如下所示（这并不是我想要的，只是为了显示一个直观的想法。）： 那么我该如何实现呢？

15 python  r  clustering  grouping 

我一直在尝试创建一个功能，该功能基本上与QGIS的“溶解”功能相同。我认为这将超级容易，但显然并非如此。因此，从我收集到的信息来看，在这里使用fiona配合体形应该是最好的选择。我刚开始弄乱矢量文件，所以这个世界对我和python来说都是新的。 对于这些示例，我使用的是位于此处http://tinyurl.com/odfbanu的County shapefile， 所以这里是我收集的一些代码，但找不到使它们协同工作的方法 目前，我最好的方法是基于以下方法：https : //sgillies.net/2009/01/27/a-more-perfect-union-continued.html。它工作正常，并且我获得了52个状态的列表作为Shapely几何形状。如果有更直接的方法来完成此部分，请随时发表评论。 from osgeo import ogr from shapely.wkb import loads from numpy import asarray from shapely.ops import cascaded_union ds = ogr.Open('counties.shp') layer = ds.GetLayer(0) #create a list of unique states identifier to be able #to loop through them later STATEFP_list = [] for i in …

15 python  shapely  fiona  geopandas 

我有一个数据框列表，所有的结构都相同，并且都是多边形。 如何使用geopandas将它们全部添加在一起？

15 python  geopandas 

要使用rasterio在栅格中的某个点上获得单个像素值，请参见此处的示例：https : //github.com/mapbox/rasterio/pull/275 但是，在rasterio（而不是cli）中是否有直接API可用于在栅格中的单个点提取值？ -编辑 with rasterio.drivers(): # Read raster bands directly to Numpy arrays. # with rasterio.open('C:\\Users\\rit\\38ERP.tif') as src: x = (src.bounds.left + src.bounds.right) / 2.0 y = (src.bounds.bottom + src.bounds.top) / 2.0 vals = src.sample((x, y)) for val in vals: print list(val)

14 python  numpy  rasterio 

我在前端有传单的django应用程序，我需要计算geojson对象的bbox，以便可以将坐标传递给 map.fitBounds() 我尝试了一些库，还尝试了getBoundsLeaflet中某些对象的方法，例如FeatureGroup，但是它抱怨未定义边界。 谁能指出我的简单解决方案？

14 python  leaflet  geojson  django 

我一直在研究使用R来挖掘Twitter的数据，但是我还没有真正找到问题的答案或不错的教程。 我感兴趣的是在特定时间范围内从具有特定标签的Twitter推文中提取信息，并在QGIS或ArcMap中在地图上绘制这些推文的位置。 我知道推文可以绑定地理位置，但是首先如何提取这些信息？

14 python  point  twitter 

我一直在访问和重新访问PyQGIS Cookbook中有关几何处理的页面：http ://documentation.qgis.org/2.0/en/docs/pyqgis_developer_cookbook/geometry.html，但似乎无法弄清楚如何获取从Python控制台绘制的多边形。有人可以帮忙吗？

14 python  geometry  qgis-2 

我正在开始研究，希望将有一段漫长的遥感专家生涯。我目前正在与ArcGIS合作以开发某些应用程序，并正在为其他应用程序学习ENVI。我已经意识到，学习编程语言是当务之急，必须面对IDL和Python之间的选择。我很想听听社区推荐哪种编程语言来遥感地表过程。

14 python  remote-sensing  software-recommendations  idl