地理信息系统

我在QGIS中有一个投影线层（Shapefile）。我想为每个线段创建一个没有线端盖的多边形。生成的多边形不应分解为单个多边形。 我尝试了fTools缓冲工具。这样可以创建一个缓冲区而不会溶解，但无法定义线路端盖。 我也尝试过GRASS的v.buffer。这可以创建具有定义的行尾的缓冲区，但是我不能禁用此溶解。 是否有可能同时存在：未溶解且管线末端没有瓶盖？结果应该是每个线段都带有多边形（缓冲区）的Shapefile。

10 qgis  buffer  dissolve 

我刚刚发现了惊人的汉森（Hansen）关于森林丧失和再生长的分类，该分类可在http://www.globalforestwatch.org/上找到，该分类发表于2013年《科学》杂志上，网址为： Hansen，MC，Potapov，P.V，Moore，R 。，Hancher，M.，Turubanova，SA和Tyukavina，A.（2013年）。21世纪森林覆盖率变化的高分辨率全球地图。Science，342（6160）（11月15日），850–854。doi：DOI：10.1126 / science.1244693。 但是，我在本文/网站中找不到确切的方法来再现这种地图，因此汉森使用了哪种分类？ 我唯一能找到的是，有监督的学习算法用于识别树的覆盖率，但这是一个广义的术语。 如果可能的话，我想使用相同的方法（但要在90年内应用），所以在我选择的地区进行汉森分类之前。

10 raster  remote-sensing  classification  landsat  forest-ecology 

我相信我应该能够从指定点开始产生步行和骑自行车等时线，所以刚刚安装了QGIS。在使用ArcMap之前，我已经很容易做到这一点，但是，我确实在QGIS方面苦苦挣扎。 从我已阅读的内容中，我应该使用PGRouting和OSM2PGrouting映射（已安装）。我以前从未真正使用过SQL，并寻求做这些等时线的虚拟指南。我已经阅读了各种方法，其中最详细地介绍了我所苦苦挣扎的SQL。

10 qgis  pgrouting  osm2pgrouting 

如何从GeoTiff文件中获取投影坐标以及这些坐标处的实际像素值，然后将其保存到numpy数组中？我有arsenci020l.tif文件，其坐标以米为单位。以下是我在上面运行的gdalinfo的摘要输出。 ~$ gdalinfo arsenci020l.tif Driver: GTiff/GeoTIFF Files: arsenci020l.tif arsenci020l.tfw Size is 10366, 7273 Coordinate System is: PROJCS["Lambert Azimuthal Equal Area projection with arbitrary plane grid; projection center 100.0 degrees W, 45.0 degrees N", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Lambert_Azimuthal_Equal_Area"], PARAMETER["latitude_of_center",45], PARAMETER["longitude_of_center",-100], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] …

10 python  gdal  geotiff-tiff  numpy 

我一直在阅读ArcGIS 10.3更新的“ ArcGIS 新增功能 ”，并且遇到了Python和ArcPy部分。它指出： ArcGIS 10.3已升级为包含Python 2.7.8。 现在这很好，任何2.x版本的python我的问题都为零。我担心的是Esri最终会将ArcGIS推向Python 3x。我希望它们不会在不久的将来出现，但同时需要在某个时候实现Python 3x。 我在python中为各种GIS和非GIS应用程序进行了大量编码。我始终在ArcGIS版本中使用最新版本。如果他们升级到3倍，我将花一个月的时间为我的应用程序进行2到3倍的升级。 我只是想知道是否有人对此主题有任何内幕信息。 Esri在不久的将来是否有任何计划迁移到Python 3，还是会继续使用Python 2x？

10 arcgis-desktop  arcgis-pro  python-3  arcgis-10.3 

我的工作任务之一是将包裹分成几组。这些小组将由代理商用来与业主交谈。目的是通过将彼此靠近的宗地分组在一起，以及将宗地划分为相等的数量，以便使工作平均分配，来简化代理的工作。特工的数量可以从一对波动到10+。 目前，我会手动执行此任务，但如果可能的话，希望自动执行该过程。我探索了各种ArcGIS工具，但似乎没有一个适合我的需求。我尝试了一个使用near_analysis和选择多边形的脚本（在python中），但它是随机的，要永久完成一个半正确的结果，比起从头开始手动完成所有工作要花费更长的时间进行修复。 有没有可靠的方法可以自动执行此任务？ 结果示例（希望没有出现黄色的除法）：

10 arcpy  arcgis-10.2 

我想使用OpenStreetMap数据以Mapbox矢量图块格式（MVT）格式渲染自己的矢量图块（最终与Mapbox的热Mapbox GL JS客户端渲染器一起使用）。 TileStache + Mapnik + mapnik-vector-tile + PostGIS是这样做的方式吗？如何配置Mapnik渲染MVT切片？

10 mapnik  mapbox  vector-tiles 

我创建了一个Python工具箱工具来对字段重新排序，并使用重新排序的字段创建新的要素类。该工具运行良好，并且我可以使用值表让用户按他们选择的顺序排列字段，也可以为每个字段填写排名值。但是，此工具的烦人之处在于，在重新排序之前，所有字段都必须一次添加到值表中。 我试图将其设置为默认情况下将所有字段引入值表，并且可以在重新排序之前删除所有不需要的字段。以前有人做过这样的事吗？我试图在UpdateParameters方法中实现这一点。这是我正在尝试的代码： import arcpy import os class Toolbox(object): def __init__(self): """Define the toolbox (the name of the toolbox is the name of the .pyt file).""" self.label = "Reorder Fields" self.alias = "Reorder Fields" # List of tool classes associated with this toolbox self.tools = [ReorderFields] class ReorderFields(object): def __init__(self): """Define …

10 arcpy  arcgis-10.2  python-toolbox  tool-validation  value-table 

是否可以直接通过管道传输数据到ogr2ogr？我正在尝试对curlGeoJSON文件进行处理，并将其通过管道转换ogr2ogr为另一种格式。 在不使用管道的情况下，我能够通过首先编写文件然后使用以下命令对其进行转换来实现此目的，但是我很好奇是否可以跳过该步骤。 curl "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson" -o "countries.geojson" && ogr2ogr -f "KML" countries.kml countries.geojson 我也尝试用curl输出创建变量而不是管道输出，但ogr2ogr无法读取它： DATA=$(curl -s "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson") && ogr2ogr -f "KML" countries.kml $DATA

10 gdal  ogr2ogr 

即使对于最简单的查询，我也无法在PostgreSQL 9.3.5上运行PostGIS 2.1来使用空间索引。在整个数据集为800万点（人口数电网从这里）。该表创建为 CREATE TABLE points ( population DOUBLE PRECISION NOT NULL, location GEOGRAPHY(4326, POINT) NOT NULL ) CREATE INDEX points_gix ON points USING GIST(location); 查询很简单 SELECT SUM(population) FROM points WHERE ST_Distance( location, ST_GeographyFromText('SRID=4326; POINT(0 0)') ) < 1000 PostgreSQL总是使用Seq扫描，我尝试了10000点的子集-仍然是Seq扫描。有任何想法吗？

10 postgis  postgresql  spatial-index  st-distance  st-dwithin 

我正在使用QGIS 1.8和2.6，无法使用任何“横截面轮廓”工具来提供椭圆（“适当”）距离。所有距离均使用太长/远的地理坐标。 例如，从卑诗省温哥华到基洛纳的直线距离约为270公里。我的横断面都说是420公里。 QGIS距离测量工具允许您使用椭球体（我将其设置为WGS 84）并给我271公里。 谁能阐明如何使用任何轮廓工具来给我“适当的”距离？

10 qgis  qgis-plugins  distance  ellipsoid  cross-section 

我试图根据给定坐标（起点和终点）绘制一条线。 Googled，发现了很少的示例，但似乎都不起作用，可能是因为它们适用于OL2，所以这是我的最后选择。 坐标位于标记数组中 <!doctype html> <html lang="en"> <head> <link rel="stylesheet" href="css/ol.css" type="text/css"> <style> .map { height: 100%; width: 100%; } </style> <script src="build/ol.js" type="text/javascript"></script> </head> <body> <div id="map" class="map"></div> <script type="text/javascript"> // inicijalizacija mape var map = new ol.Map({ target: 'map', layers: [ new ol.layer.Tile({ source: new ol.source.MapQuest({layer: 'osm'}) // …

10 openlayers 

阅读完在QGIS中启用多线程处理后？，我想知道QGIS 2.6在运行各种过程时是否将包括此功能。我在网上检查了一下，发现QGIS 2.4引入了多线程渲染（我认为它在2.6中将再次可用）。现在，可以非常快速地重绘地图，尤其是在处理大型数据集时。 我读到将多线程包含到核心功能中是相当困难的，建议用户操纵Python代码以实现此功能，例如在PyQGIS中的并行GIS操作中？。 我还检查了功能请求，但Tim Sutton已将其关闭了9个月，最新评论为： “我要结束这一点-Martin Dobias在分支机构中有一个实现，它将在QGIS 2.2之后合并。” QGIS 2.6中的某些功能是否会利用多线程（或者将再次专注于呈现），如果没有，在QGIS 3.0中会有所改变吗？

10 qgis  qgis-processing  multithreading  qgis-3 

我阅读了有关30米（SRTM 30米）的数字高程模型的新版本的信息，但是我无法下载。 一些来源： http://glcfapp.glcf.umd.edu:8080/esdi/index.jsp http://asterweb.jpl.nasa.gov/data.asp http://www2.jpl.nasa.gov/srtm/

10 dem  remote-sensing  srtm 

我正在尝试使用2013年以来的Landsat 8卫星图像对蒙古北部的某些地区进行分类。由于记录是在冬季进行的，因此在采集时太阳非常低。因此，山上有很长很深的阴影。 我可以使用DEM识别此阴影，如以下问题中所述：如何从Landsat图像中删除诸如阴影之类的地形效果 如何对修剪的阴影区域执行监督分类？有可能加强这一领域吗？我尝试了一些频段比率，但我不知道哪个最适合我的任务。 在此图像中，您可以看到暗影中有一些植被区域，但是无法对它们进行分类。

10 dem  remote-sensing  classification  landsat-8  shadow