地理信息系统

我想使用“ QGIS栅格样式”对话框来显示修改后的栅格值。具体来说，我想查看使用栅格值平方根的渐变。 我知道可以使用栅格计算器来转换栅格值，但是-如果可能的话-我想跳过该步骤并在样式化过程中应用数学公式，并避免创建原始栅格的修改后的副本。 对于矢量文件，使用颜色/宽度/等数据定义的值，可以进行这种类型的操作。可以使用栅格进行此操作吗？

10 qgis  raster  style 

新手在这里苦苦于GIS。我正在尝试使用在其县网站（县）网站上找到的shapefile来为密尔瓦基市（Milwuakee）绘制地图。我在这里取得了一些成功。我的代码是给： from pyproj import Proj, transform # wisconsing EPSG:32054 # epsg:4326 is for the entire world, wgs 84...not obvious inProj = Proj(init='epsg:32054') outProj = Proj(init='epsg:4326') x1,y1 = 2560131.496875003, 406816.434375003 x2,y2 = transform(inProj,outProj,x1,y1) print(x2,y2) 与输出， -65.70220967836329 43.08590211722421 问题是这是错误的。密尔沃基的lon / lat为-87.863984和42.920816。 其次，我如何以编程方式对整个shapefile执行此操作。我想将其绘制到底图中。当我尝试遵循此线程时，出现错误代码： with fiona.open("ward2012/ward.shp") as shp: ori = Proj(init='epsg:32054' ), dest= Proj(init='EPSG:4326',preserve_units=True) …

10 python  shapefile  fiona  pyproj 

我正在尝试创建要打印的地图，并且多边形的大小差异很大。由于这将是参考地图，因此我需要标记每个多边形，但是某些名称将不适合该多边形。在这些情况下，我想改用功能ID，然后让属性表将ID链接到名称。 我想做的是具有某种表达式，可以用来选择在地图上使用的标签以及在属性表中包括的行。这样的isLabelled事情将能够考虑地图比例，标签长度，字体大小等。

10 qgis  labeling  print-composer  qgis-custom-function 

我想在模型中创建一个简单的if / else函数，我想知道是否可行。 我希望这件事做到这一点，如果选择了“细分层”的输入，则它不必执行“细分（均移）”，而将“细分层”用作下一个流程的输入。 “细分（平均转换）”的输入如下 尝试实现约瑟夫的答案后进行更新：我使用以下代码，并收到此错误 更新2，替换参数后，在第28行出现错误，这是“ else：”立场 ##Example=name ##BGREN=optional raster ##Segment_Layer=optional vector ##output=output vector if BGREN is not None: processing.runalg('otb:segmentationmeanshift', \ -BGREN\ -0\ -70\ -70\ -0.1\ -100\ -100\ -0\ -0\ -0\ -'True'\ -'False'\ -1\ -0.1\ -'Layer'\ -'DN'\ -1024\ -1\ -0\ -output) else: output=Segment_Layer 更新它的工作原理:)经过几次尝试，并在Joseph的帮助下，此代码现在像个魅力一样工作了，最后的技巧是用'None'替换两个可选值。 ##Example=name ##BGREN=optional raster ##Segment_Layer=optional vector …

10 qgis  python  pyqgis  qgis-modeler  if-else 

我在PostGIS上有一个多边形层，当我尝试使用QGIS添加特征时，在保存编辑时会出错。 这是错误通知“无结果缓冲区”。 我在这里发现了类似的情况：添加功能时出现PostGIS错误：没有结果缓冲区。但这并没有给我解决的答案。有任何想法吗 ？

10 qgis  postgis 

我正在尝试使用开源python工具从坐标列表中创建多边形shapefile。以下是我到目前为止从Python GDAL / OGR Cookbook和GIS SE答案中获得的信息。 还有一个类似的问题Python：如何从X，Y坐标列表创建多边形Shapefile？，尽管这个问题是指使用pyshp。但是，我对仅使用gdal / ogr Python工具创建多边形shapefile感兴趣。 import ogr def create_polygon(coords): ring = ogr.Geometry(ogr.wkbLinearRing) for coord in coords: ring.AddPoint(coord[0], coord[1]) # Create polygon poly = ogr.Geometry(ogr.wkbPolygon) poly.AddGeometry(ring) return poly.ExportToWkt() def write_shapefile(poly, out_shp): """ https://gis.stackexchange.com/a/52708/8104 """ # Now convert it to a shapefile with OGR driver = ogr.GetDriverByName('Esri …

10 python  gdal  polygon  ogr 

我正在使用Essen，并尝试遵循Anita Graser精彩的教程《照明轮廓》。我已经创建了azimuth函数并加载了该函数，但表达式选项卡中没有任何显示-甚至没有使用默认功能的Custom！表达式与关联文件一样存在azimuth.py于C :\Users\Max\.qgis2\python\expressions中.pyc。重新启动QGIS没有任何区别

10 qgis  field-calculator  qgis-custom-function 

我正在尝试创建一个显示与最近商店相对应的多边形的多边形图层。 对于一组30个商店位置，结果应为具有30个要素的面图层，每个商店一个。每个要素将代表一个区域，在该区域中对应的商店最近。例如，多边形12内的地址表示商店12最近。 我的商店位置图层和OSM街道图层具有“最大速度”列。生成的多边形层应基于OSM街道层，而不仅仅是线性区域。 这个想法是，给定一组固定的商店和道路，最近的商店应在地图上的任何给定点保持不变。因此，我希望最终的多边形层将无间隙。 这段 Maptitude软件的视频提供了我要完成的工作的完美示例（快进视频的1:55）。或者，请参见下面的图像作为示例： 对于每个医院，都有一个相应的区域，该区域内的任何东西都最靠近该区域内的医院。 如何使用QGIS复制此内容？

10 qgis  network-analyst  routing  network  proximity 

刚刚从OSGEO4W安装程序安装了QGIS 2.14.6 64Bit。处理扩展版本为2.12.99。从屏幕快照中可以看出，我无法像往常一样切换到“高级界面”以访问该Polygonize工具（到多边形的线仍然不起作用，cp。QGIS 2.12到多边形的线无法正常工作）： 在另一台计算机上并行运行的同一安装具有处理插件版本2.10.3，在这里我可以切换到高级界面并访问该Polygonize工具： 添加/更新： 完全令人困惑：在独立安装程序安装的2.14.5中，处理插件的版本也为2.12.99，但是在这里我可以访问该Polyonize工具： 到目前为止的研究工作：我意识到在情况2中，处理插件安装在username\.qgis2\python\plugins\核心插件的异常路径下。显然（或也许？）问题是在处理成为核心插件之前先安装了QGIS，因此此处的插件似乎覆盖了核心插件。这也将解释为什么更新2.下的插件不受更新的影响（由于处理是核心插件，因此在存储库中不可用）。解决方案只是卸载插件，重新启动QGIS和核心插件版本2.12.99。显着：Polygonize工具仍然可用（与情况1相反）。到目前为止，一切都很好。 总结一下： 在这三种情况下，处理插件的版本均为2.12.99： OSGeo4W安装程序，无Polygonize工具 OSGeo4W安装程序，提供Polygonize工具 独立安装程序，提供Polygonize工具 那么情况1可能是什么问题呢？我可以错过OSGeo4W安装程序中的某些内容吗？

10 qgis  qgis-processing  polygonize  user-interface 

有人可以帮助我创建QGIS处理脚本，以将序列添加到PostGIS中现有的唯一标识符列（类型：整数）吗？ 这将非常有帮助，例如，作为错误＃6798的解决方法。不幸的是，我没有任何Python经验。 CREATE SEQUENCE /*input_schema*/./*input_table*/_/*uic*/_seq OWNED BY /*input_schema*/./*input_table*/./*uic*/; SELECT SETVAL('/*input_schema*/./*input_table*/_/*uic*/_seq', (SELECT MAX(/*uic*/) FROM /*input_schema*/./*input_table*/)); ALTER TABLE /*input_schema*/./*input_table*/ ALTER COLUMN /*uic*/ SET DEFAULT nextval('/*input_schema*/./*input_table*/_/*uic*/_seq'::regclass);

10 postgis  pyqgis  qgis-processing  serial 

我想编写一些使用Qgis处理工具箱的独立脚本。 我已经读了一些线程（例如，在这里和这里），但是找不到有效的解决方案。 在Ubuntu Xenial 16.04 LTS上使用Qgis 2.16.1 我的脚本的import部分如下所示： # Python modules import sys import time import os # Qgis modules from qgis.core import * import qgis.utils from PyQt4.QtCore import QFileInfo, QSettings 有人知道我无法导入处理模块吗？ 通过简单的导入处理，我得到了： Original exception was: Traceback (most recent call last): File "/home/steph/Documents/Projets/20141227-CIM_Bishkek/Scripts/python/00-projets/20160811-AnalysesUAVs/20160811-UAVAnalyse.py", line 36, in <module> import processing File …

10 qgis  pyqgis  qgis-processing 

我正在尝试通过使用虚拟层来使QGIS中的多边形几何图形相交： SELECT sbqry.rowid AS gid, sbqry.geom FROM (SELECT ST_Intersection(land_parcels.geometry, flood_zone.geometry) AS geom FROM land_parcels, flood_zone WHERE ST_Intersects(land_parcels.geometry, flood_zone.geometry) AND NOT ST_Touches(land_parcels.geometry, flood_zone.geometry)) AS sbqry; 不幸的是，sbqry.rowid AS gid返回NULL而不是自动递增值。 有人知道如何创建名为“ gid”的唯一标识符列吗？据我所知，虚拟层是基于SQLite / Spatialite的。

10 qgis  spatialite  sqlite  unique-id  virtual-layer 

我正在使用ArcGIS Pro，无法找到撤消字段计算的方法。 使用ArcGIS 10.x for Desktop，当编辑会话处于活动状态时，字段计算可以撤消，但是Pro似乎处于永久性的编辑会话中，其中不包含撤消字段计算的功能。 是否可以在ArcGIS Pro中撤消字段计算？

10 field-calculator  arcgis-pro  undo-redo 

我需要根据要素类中的组自动增加字段。我在给定的多边形内有8个图，并且需要为每个多边形内的每个图集从1-8分配一个ID。多边形将具有其自己的唯一ID号，该ID号可用于对图进行分组。 我认为这是对它的修改： rec=0 def autoIncrement(): global rec pStart = 1 pInterval = 1 if (rec == 0): rec = pStart else: rec = rec + pInterval return rec

10 arcgis-desktop  field-calculator  python-parser 

我发现了一种有趣的比例符号映射和标记方式，我认为这是由Adobe Illustrator制作的，但我无法确定它是否可以在QGIS中完成。我最感兴趣的是文本（属性）在圆圈周围的显示方式。我真的不知道如何在QGIS中像这样制作这样的主题地图标签。 这是图片：

10 qgis  labeling  symbology  thematic-map