地理信息系统

5

如果我在地理坐标（WGS84）中具有多边形，那么如何考虑到地球的曲率来测量每个在地球表面上的总面积？

12 coordinate-system area

6

我在运行ArcGIS Desktop 10的虚拟PC上配置了RAM磁盘，并将HKCU \ Software \ ESRI \ Output \ TempPath设置为RAM磁盘上的目录。很难感觉到这种配置有任何改进，我很好奇，是否有其他人在非虚拟计算机上完成过此测试以及结果如何。轶事的答案很好。另外，将缓存路径放在ram磁盘上可能不是一个好主意，但是有人这样做吗？（HKCU \软件\ ESRI \ Desktop10.0 \ Settings \ Cache路径）

12 arcgis-desktop arcgis-10.0 performance

4

获取QGIS 2.x python代码以在QGIS 3.x中工作？对象没有属性“ legendInterface”的问题

下面的代码在QGIS 2.x中运行良好，但是在QGIS 3.x中不起作用。 myDir = 'd:/work/output_folder/' layers = iface.legendInterface().layers() pipe = QgsRasterPipe() for layer in layers: extent = layer.extent() width, height = layer.width(), layer.height() renderer = layer.renderer() provider=layer.dataProvider() crs = layer.crs().toWkt() pipe.set(provider.clone()) pipe.set(renderer.clone()) opts = ["COMPRESS=LZW"] file_writer = QgsRasterFileWriter(myDir + layer.name() + ".tif") file_writer.setCreateOptions(opts) file_writer.writeRaster(pipe, width, height, extent, layer.crs()) 这是我运行代码时遇到的错误： …

12 pyqgis qgis-3 qgis-python-console

2

如何从R包SF组合SFC对象

使用R包sf，如何组合sfc对象？例如，给定以下代码，将如何创建包含来自和的几何的单个sfc对象？（应为2。）sfc12sfc1sfc2length(sfc12) library(sf) pt1 = st_point(c(0,1)) pt2 = st_point(c(1,1)) sfc1 = st_sfc(pt1) # An sfc object sfc2 = st_sfc(pt2) # Another sfc object # sfc12 = ? 一些无效的方法： sf_sfc(sfc1, sfc2) # Error in vapply(lst, class, rep("", 3)) : values must be length 3, # but FUN(X[[1]]) result is length 2 sfc1 …

12 r sf

3

根据QGIS中的线数据制作热图？

我想基于属性对线数据执行线密度分析。热图选项（内核密度分析）似乎仅适用于点数据。我试图将线数据转换为点，并创建一个热图，但是权重的属性在转换时丢失了。这是使用ArcGIS执行的线密度分析。我一直想对QGIS做同样的事情。是否可以在QGIS中执行类似的线密度分析？

12 qgis heat-map density

4

定性调色板-复印安全（可用于灰度）

我正在寻找可以在黑白打印机中打印的5类定性调色板。下图显示了彩色地图（上图：定性图，下图：发散图），以及每张图的灰度外观（右侧）。这两个调色板（Set3和Spectral）都属于ColorBrewer。 [UL]定性 Set3 [UR] Set3 灰度定性 [LL]发散 Spectral [LR] Spectral灰度差异我尤其会在右上方的图像中辨别组，而底部（光谱）则更好，即使不是完美的（我不能说是蓝色还是红色）。根据ColorBrewer所说，“光谱”配色方案（我在这里使用的方案）是唯一photocopy safe适用于5类调色板的。减少类（减少到3或4）可能是一种选择，但恐怕这不是此地图的预期目的。这可能与GIS软件无关，但是我基本上使用QGIS，它提供了ColorBrewer，cpt-city和其他调色板和色带。 @obrl_soil建议这是Viridis颜色集 @underdark建议这是两种类型的4 colors + white方法。左：橙红色+白色（顺序），右：紫橙色+白色（发散）

12 qgis symbology

1

更改默认的ArcGIS Home文件夹？

ArcGIS Desktop将主文件夹的默认值设置为类似C:\Users\<username>\Documents\ArcGIS- 根据Windows的版本，在My Documents或Documents文件夹中创建一个文件夹。我个人试图将GIS工作排除在Documents文件夹之外，通常将所有文件存储在名为的文件夹中D:\GIS。主要原因是因为我的工作是网络存储的Documents文件夹，其配额为5mb，对于存储任何类型的GIS数据而言，这太小了。Esri加载项默认也安装到该文件夹中，这使情况变得更糟。我一直在尝试寻找一种方法，D:\GIS以在不先保存MXD的情况下将ArcGIS默认设置为我的文件夹作为家庭位置。我知道我可以将MXD保存到该文件夹中，然后它将成为该MXD的主文件夹。我希望能够打开ArcMap并从该Home文件夹访问我的文件/文件夹，D:\GIS而不必去查找它，甚至不必将其作为文件夹连接。默认情况下，ArcMap中使用Home文件夹保存结果，存储新数据集以及访问基于文件的信息。该位置最初是用于向ArcMap添加新图层以及存储在编辑和地理处理期间创建的新数据集的位置。非空间数据也是如此。当您将页面布局导出为PDF时，默认位置将是地图的原始位置（地图文档的工作区文件夹）。使用主文件夹我发现了几篇文章（主要是在Esri的Geonet上）询问如何执行此操作，但是除了我上面提到的两个选项（实际上不是默认选项）之外，这些文章都没有解决方案。如何强制ArcGIS将自己的文件夹位置用作默认主文件夹，而不必先保存MXD？

12 arcgis-desktop default

3

在Python中高效获取多个多边形的交集

我想得到多个多边形的交集。使用Python的shapely程序包，我可以使用intersection函数找到两个多边形的交点。是否有类似的有效函数来获取多个多边形的交集？这是一个代码片段，以了解我的意思： from shapely.geometry import Point coord1 = ( 0,0 ) point1 = Point(coord1) circle1 = point1.buffer(1) coord2 = ( 1,1 ) point2 = Point(coord2) circle2 = point2.buffer(1) coord3 = ( 1,0 ) point3 = Point(coord3) circle3 = point3.buffer(1) 可以找到两个圆的交点circle1.intersection(circle2)。我可以找到所有三个圆的交点circle1.intersection(circle2).intersection(circle3)。但是，这种方法不适用于大量的多边形，因为它需要越来越多的代码。我想要一个接受任意数量的多边形并返回其交点的函数。

12 python polygon intersection polygon-creation shapely

3

使用QGIS将点触线属性返回到线层？

我正在使用QGIS 2.14.4-Essen。我有两层： points.shp包含带有YEAR列的点 lines.shp包含精确连接点与point.shp的线我想将point.shp的YEAR属性返回到lines.shp。每条线都是一个单独的线段，两端各有一个点（请参见下图）。我想得到的是从第一点开始的YEAR和从第二点开始的YEAR，回到每条线属性。例如：第1行触及YEAR = 2010的第一点和YEAR = 2011的第二点。我想将“ 2010-2011”之类的内容恢复到第1行属性中。结果应如下所示： id | points ----+----------------- 1 | 2010-2011 2 | 2011-2012 3 | 2012-2016 有没有一种使用QGIS来实现这一目标的方法？

12 qgis qgis-2 spatial-join

2

显示QGIS中线要素的定向方式？

我正在街道shp文件的QGIS上工作，该文件具有Left From，Left To，Right From，Right To字段字段，正在构建地址解析器。我需要知道这些特征的定位方式，我认为这与知道使这些线最先出现的两个（或多个）点中的哪一个相同。有没有一种方法可以自动旋转标记线样式上的三角形标记，以显示线的方向？

12 qgis style line markers

4

在地图上显示天线符号：点符号或要素（多边形）

我想在地图上显示一个蜂窝网络。输入数据是一个.csv文件，其中每个字符串都是一个蜂窝扇区。这些属性是：扇区ID，其坐标，其方位角和天线波束宽度的角度。天线波束宽度值在30到360度之间。360天线波束宽度意味着它必须在地图上显示为圆形。具有其他波束宽度的天线必须显示为具有适当孔径角的扇区。是否可以仅使用符号显示天线？我知道如何创建自己的SVG符号，并希望找到如何根据方位旋转它。但是，是否有任何方法可以根据其30至360度的属性值来应用可变的天线波束宽度？我认为符号是绘制天线的最佳方法，因为如果在QGIS中可以根据视图比例尺在地图上进行动态可视化，则这些符号是可行的。当然，可以通过将适当的多边形绘制为图层要素来解决该任务，但这将是一种解决方法。

12 qgis symbology style visualisation

5

如何自动填充多边形中的孔？

QGIS中是否可以通过设置参数（例如最大面积，边界线之间的距离）来自动填充多边形中的孔？我想删除/填充非常“薄” /小的孔。

12 qgis polygon

4

我在postgis安装中找不到shp2pgsql

我已经按照本指南安装了postgis 。一切都OK。 SELECT postgis_full_version(); "POSTGIS="2.1.7 r13414" GEOS="3.4.2-CAPI-1.8.2 r3921" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.10.1, released 2013/08/26" LIBXML="2.9.1" LIBJSON="UNKNOWN" TOPOLOGY RASTER" 但是我没有看到shp2pgsql，/usr/lib/postgresql/9.4/bin如您所见： foo@bar:/usr/lib/postgresql/9.4/bin$ ls clusterdb initdb pg_dump pg_standby psql createdb oid2name pg_dumpall pg_test_fsync reindexdb createlang pg_archivecleanup pg_isready pg_test_timing vacuumdb createuser pg_basebackup pg_receivexlog pg_upgrade vacuumlo dropdb pgbench pg_recvlogical pg_xlogdump droplang pg_controldata …

12 postgis postgresql ubuntu shp2pgsql

2

将向量图块源添加到QGIS

是否可以将矢量图块源（使用MapBox的协议缓冲区格式或GeoJSON或其他方式）添加到QGIS？如果是这样，怎么办？（我的问题有点类似于可以将新的Mapbox矢量图块用作QGIS底图吗？但不限于MapBox源，而不仅仅是底图或OpenLayers。）

12 qgis mapbox vector-tiles

2

使用键盘停止具有无限循环的PyQGIS脚本？

因此，我编写了一个脚本，该脚本在一定时间间隔（使用time.sleep（5）每5秒一次）中使用“ while True：”一次又一次地执行我想要的操作。到目前为止，还不错，但是当我想要停止它时，我只是做不到。我尝试了Control + C，Control + Break，Escape，但实际上它忽略了我的键盘。停止它的唯一方法是关闭QGIS。有任何想法吗？此外，当脚本命中time.sleep（5）时，QGIS会滞后5秒钟，并且无法进行平移，但是我认为这是正常的。这是我的脚本： from PyQt4.QtGui import * from PyQt4.QtCore import * from qgis.core import * from qgis.utils import iface import time while True: def change_color(): active_layer = iface.activeLayer() pipeline=[] txt=open('C:/users/stelios/desktop/project/Sensor.txt','r') for line in txt.readlines(): pipeline.append(line.split()) print pipeline pipeline2=[] for label,color in pipeline: if "0" …

12 pyqgis loop