地理信息系统

制图师，地理学家和GIS专业人员的问答

1
在R中使用writeOGR（）将多层写入GeoPackage？
我试图在R中的同一GeoPackage中写入多个图层，但出现错误Creation of output file failed。我试图搜索使用RGDAL读写.gpkg文件的文档，例如，弄清楚是否writeOGR()实际上支持多层，但收效甚微。如果可能的话，这怎么可能呢？最小的工作示例： library(sp) library(maptools) library(rgdal) data(wrld_simpl) norway <- wrld_simpl[wrld_simpl$NAME == "Norway", ] sweden <- wrld_simpl[wrld_simpl$NAME == "Sweden", ] file <- tempfile("scandinavia", fileext = c(".gpkg")) writeOGR(norway, dsn = file, layer = "norway", driver = "GPKG") writeOGR(sweden, dsn = file, layer = "sweden", driver = "GPKG") ogrListLayers(file) 显然，有一个可以完成技巧的ogr2ogr shell命令（hat …
11 gdal  r  rgdal  geopackage  sf 
2
QGIS打印作曲家pdf导出中的可单击HTML链接？
根据此处的 qgis帮助文件： '标签可以解释为html代码：选中“渲染为HTML”复选框。现在，您可以插入网址，链接到网页的可点击图片或更复杂的内容。” 但是，我的标签内容如下： <a href="http://google.com">google</a> 产生一个以链接形式显示为“ google”的pdf文件，该pdf文件为蓝色文本（如链接），但不可点击。 输出中是否可以包含可点击的链接？如果是这样，我在做什么错？
4
尝试使用python从字段中提取唯一值列表
我在FGDB的许多表中都有许多列，我需要在其中提取每列的唯一值。 例如：值可能是[1,2,2,2,3,4]，而我正在尝试返回[1,2,3,4] 我可以在ARCGIS中通过其他多种方式来完成这项工作，但是我正在努力扩展自己。 我在网络上找到了一条我认为可以完成工作的python，但是我正在努力使其运行（我不断收到无效的语法错误，因为我一直在第3行中收到语法错误），这无疑是一个非常简单的用户错误。 下面的代码段 import arcpy def unique_values(r'N:\GISProjects\Landuse\Plant_Biosecurity_Project\ArcGIS_Online.gdb\Holdings_Property_Merge' , 'LU_ALUMMaj'): with arcpy.da.SearchCursor(table, [field]) as cursor: return sorted({row[0] for row in cursor}) 这是我从崇高的文字中得到的错误消息： File "C:\Users\hawkinle\Desktop\STDTAS\Unique_Data.py", line 3 def unique_values(r'N:\GISProjects\Steve_Eastwood_Landuse\Plant_Biosecurity_Project\ArcGIS_Online.gdb\Holdings_Property_Merge' , 'LU_ALUMMaj'): ^ SyntaxError: invalid syntax [Finished in 0.1s with exit code 1] 原始问题的更新 我现在使用下面提供的答案更新了我的代码，但是正在收到一个辅助错误。 新代码段： import arcpy def unique_values(table , …
2
在PyQGIS中检索可用的PostGIS连接
我可以在PyQGIS中检索到PostGIS数据库的可用连接吗？我想提供一个可用的数据库连接列表，然后在我的插件的用户界面中提供一个表列表。 我检查了食谱，但找不到进一步的方法。
11 qgis  postgis  pyqgis 
1
在R传单中包含风数据的交互式动态地图？
我有一些这样的数据： longitude latitude speed direction 1 6.10722222 46.23639 4 360 2 6.95416667 43.54694 4 360 3 7.21472222 43.66556 13 330 4 4.01666667 48.32167 7 290 5 2.30833333 43.21611 14 290 6 2.48305556 44.40806 13 320 7 5.21500000 43.43694 19 330 8 4.92361111 43.52278 32 320 9 5.10805556 43.60306 26 330 …
3
使用ST_Difference删除重叠的功能？
我正在尝试使用PostGis 2.1（和Postgres SQL 9.3）使用ST_Difference创建一组多边形（processing.trimmedparcelsnew），这些多边形不包含另一组多边形（test.single_geometry_1）所覆盖的任何区域。这是我的查询： CREATE TABLE processing.trimmedparcelsnew AS SELECT orig.id, ST_Difference(orig.geom, cont.geom) AS difference FROM test.single_geometry_1 cont, test.multi_geometry_1 orig; 但是生成的多边形尚未修剪，相反，它们似乎在与另一层相交的地方被拆分了。我试过只运行选择而不将结果放在表中以及其他我能想到的东西，但是我似乎无法使该功能正常工作。 我已附上结果图 评论后，我尝试添加WHERE子句。我希望没有交集的宗地，而其他宗地的相交区域被删除（图层test.single_geometry表示要从宗地中删除的污染）。我尝试了一个相交，但是我当然想要非相交，所以我现在尝试分离。我也尝试过将折纸添加到表中，但是ST_Difference的文档（http://postgis.net/docs/ST_Difference.html）确实说它返回了我需要的精确几何图形（表示几何图形A的那部分，不会与几何B）相交，所以我对为什么要在表格中使用原始多边形感到困惑。无论如何，这是我修改后的代码： CREATE TABLE processing.trimmedparcelsnew AS SELECT orig.id, ST_Difference(orig.geom, cont.geom) AS difference, orig.geom AS geom FROM test.single_geometry_1 cont, test.multi_geometry_1 orig WHERE ST_Disjoint(orig.geom, cont.geom); 根据dbaston的回答，我现在尝试： CREATE TABLE processing.parcels_trimmed AS SELECT id, …
3
如何在QGIS中查看GTFS数据？
我想使用QGIS创建我的地图。关于公交路线的gtf是一种数据馈送。我可以在Google Earth上查看它（尽管比较笨拙），但是当我尝试将.kml导入QGIS时，它刚用完RAM（持续30分钟），并使我的地图崩溃了。如果我在18xx图层中导入10，则该路径将重复10次（“ Google Earth”可能会显示所有路径，尽管比较费时）。 我的程序 我正在尝试一种新方法，所以... 我所做的是在TfNSW网站（Google tfnsw数据和登录名）上下载.txt，然后使用bdferris称为Transitfeed的.py脚本将.txt转换为.kml。最终将.kml导入为矢量层，但是我之前曾在Google Earth中查看过，所以我可以。 我打算将shape.txt转换为路线形状。 问题 我在考虑各种可能性和解决方案，例如： .kml错误 gtfs中的.txt错误 从.txt到.kml的错误转换 我在想如果前两种可能性都有效，那么（我的.kml）解决方案合适的验证方法是什么？否则，我想要将gtfs .txt文件查看到QGIS中的低成本解决方案。
5
免费雷达（SAR）卫星数据[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 3年前关闭。 我正在寻找雷达（SAR）卫星图像的免费（非商业用途）来源。有人知道来源吗？
3
PostGIS中的几何SRID问题
我在PostGIS中选择了我的桌子（实地） id_borne serial NOT NULL, num_borne character varying, shape_borne geometry 我成功使用更改了表的SRID select UpdateGeometrySRID('public', 'borne', 'shape_borne', 26191) ; 旧的SRID为0。 问题是，当我尝试向表中添加新元素时，borne出现此消息错误： geometry srid (0) does not match column srid (26191) 那是我的问题！SRID已更新，但是当我尝试添加新点时，出现此消息。
3
是否可以在QGIS中旋转画布
画布是否可以旋转，使北方指向投影定义以外的其他方向？这个问题针对作曲家，而不是画布。错误跟踪数据库中有一个相关的5年旧条目，由于不必要而将其关闭。这是当前状态，还是有办法在较新版本中旋转画布？ 创建新的投影确实不是答案
11 qgis  qgis-2 
1
如何在QGIS 2.2.0中的COMPOSER中将North Arrow和Scale放置在地图上
我是QGIS的新手，并且在此处询问和回答的有关North Arrow的问题似乎并不能帮助我解决以下问题。 在QGIS 2.2.0 Valmiera的主映射窗格中，我可以添加比例尺和向北箭头，但是一旦我进入Composer来打印地图，这两件事就会从我在Composer中打开的地图中消失。 Composer可以让我添加一个比例尺，但是它显示了我将在其上打印的纸张上的实际比例，而不是所显示地图的地球表面的实际比例。我什至尝试刷新地图，但是没有用。 底线是：如何在Composer中获得缩放比例的东西（对不起！）（在缩放时放大和缩小值），以及如何在其中获得北向箭头。
7
如何对Google SpreadSheet进行地理编码？
我在Google文档中有一个电子表格，该电子表格大约有20-30行，包括以下6列： ID, Name, Address, Cell_No, Latitude, Longitude 最后两列为空。 对于非编程人员的新手用户来说，对地址进行地理编码并使用地理编码的结果填充纬度和经度字段的简便方法是什么？ 我只需要填充纬度和经度字段。我不需要在地图上看到任何数据。
3
如何从QGIS处理加载内存输出？
许多处理算法可以选择将输出保存为内存。如果从工具箱运行，则可以正常工作，因为您可以选中“运行算法后打开输出文件”框。 查看QGIS的源代码，似乎在函数中定义了用于加载内存层的代码Postprocessing.handleAlgorithmResults。该函数使用对alg.outputs列表的访问权，并执行以下操作： for out in alg.outputs: progress.setPercentage(100 * i / float(len(alg.outputs))) if out.hidden or not out.open: continue if isinstance(out, (OutputRaster, OutputVector, OutputTable)): try: if out.value.startswith('memory:'): layer = out.memoryLayer # NOTE!! QgsMapLayerRegistry.instance().addMapLayers([layer]) # NOTE!! else: # ... 当您从控制台运行处理算法时，是否可以在不访问该对象的情况下加载图层？我可以跑 processing.runalg("qgis:intersection", layer1, layer2, "memory:") 甚至 processing.runalg("qgis:intersection", layer1, layer2, "memory:myLayerName") 但是，我找不到一种方法来获取结果输出。
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.