地理信息系统

为了在子文件夹上递归迭代并将所有shapefile合并为一个脚本，基本脚本是： #!/bin/bash consolidated_file="./consolidated.shp" for i in $(find . -name '*.shp'); do if [ ! -f "$consolidated_file" ]; then # first file - create the consolidated output file ogr2ogr -f "ESRI Shapefile" $consolidated_file $i else # update the output file with new file content ogr2ogr -f "ESRI Shapefile" -update -append $consolidated_file …

11 gdal  ogr2ogr  merge  bash  mergeshapes 

我有一堆线串的geopandas数据框，其中的一些数据与每个顶点/点相关： Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) 这些点将根据ID转换为线串，并按时间排序。 我想在满足某些条件的地方断开线串。此时，点之间的距离大于某个值。将来，数据字段的功能可能会有所价值。例如，当Speed超过5 kph时拆分线串。 当前的问题是，某些轨道是由具有重复ID的点形成的，因此线串在很长的距离上来回跳跃，我希望有一个阈值来打破这些线。 关于构建此方法的正确方法或可能有用的库/方法的任何想法？ 该数据帧有超过15万条磁道，每个磁道中有许多点，因此效率会很高。 这是轨道DF的示例： ID geometry 204235000 LINESTRING (37.62001 -28.99535, 37.62015 -28.9... 205400000 LINESTRING (3.807816666666666 -18.083181666666... 207138000 LINESTRING (22.73206 -34.97915833333333, 22.73... 209016000 LINESTRING (8.447673333333331 -23.522783333333... 这是DF点的示例。共有18列，包括日期时间，时间点（经度，纬度），速度，大小等： Index Heading Latitude Longitude ID 20 92.8 -35.946802 13.089695 210725000 21 …

11 python  shapely  geopandas 

我正在寻找使用GDAL（OSGeo4W）将16位无符号tiff图像缩放为8位图像。但是，我想限制像素值，以使它们保持在8位图像的范围内。我的意思是，我想确保16位图像中的视觉比例在8位图像中表示（0-255）。使用gdal_translate将16位转换为8位会剪切像素值，并且不会生成RGB图像。我不确定从gdal_translate中调用什么函数和设置，或者是否有更好的选择。 GDALINFO C:\>gdalinfo C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Driver: GTiff/GeoTIFF Files: C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Size is 17310, 11310 Coordinate System is `' Metadata: TIFFTAG_DATETIME=2015:05:03 12:27:13 TIFFTAG_IMAGEDESCRIPTION=UltraCam-Lvl03 -------------- CAM_ID: UC-SXp-1-50215465 [3] IMG_NO: 4761 RECORD_GUID: 8fc3d7a9-d2e5-40e5-babf-939f1f803dcf IMG_GUID: C850D0A4-39CA-4BC8-ACD7-81A363D78C6E FILE_GUID: 86C4F8E4-A2FD-43FB-B3B1-4C4ED8DF0E42 LICENSE_ID: 820342819 SOFTWARE: UltraCam Aerial Radiometry Core 12.6.1408.2501 PIXEL_SIZE_WIDTH: 6 [micron] PIXEL_SIZE_HEIGHT: 6 [micron] APERTURE: F_8 EXPOSURE_TIME: 0.002000000000 …

11 gdal  imagery  gdal-translate 

使用QGIS 2.16.2，我可以使用以下表达式从shapefile属性表中成功选择记录： left("start_time", 10) = '2015-08-01' 请注意，start_time是一个字符串字段。 但是，当我将相同的表达式复制并粘贴到shapefile的“图层属性”>“常规”>“特征过滤器”>“查询生成器”中时，出现以下错误： 我希望在QGIS中的任何地方都可以使用相同的表达式。我究竟做错了什么？

11 qgis  query  filter  expression 

是否有工具可以从较大的GTFS提要中提取较小的区域？据我所知， Google的Transitfeed库包含“合并”工具，但不包含“提取/剪切”工具。我发现的其他库只能读取和绘制数据，但不支持操纵。 背景信息：Transitfeed Feedvalidator显示了我的GTFS文件的很多错误。我想提取我感兴趣的区域以评估本地质量。 Twitter上的建议向我介绍了FME，它既具有GTFS读取器又具有写入器。我的简单工作流程包括一个GTFS阅读器（左下），Shapefile阅读器（左上角为感兴趣的区域），一个裁剪器和一个GTFS书写器作为输出（现在只写停靠点和形状文件，但最终我还是想要所有文件）： 到目前为止，问题在于输出形状内容（橙色点）不包含位于黄色区域内的所有输入点（绿色）： （数据来源：https : //transitfeeds.com/p/helsinki-regional-transport/735/latest/download） 使用“ Create Aggregates” =“ No”，剪切的形状看起来很完整 但是停站肯定不是

11 fme  general-transit-feed-specification 

是否可以同时拖动QGIS中某些平行线的多个顶点？ 我有一个平行线的线层。在某些区域中，我必须移动所有平行线的顶点，以使它们仍然彼此保持距离。 使用节点工具，我只能选择要拖动的一条线的顶点，而不能同时选择多条线的顶点： 在“多行图层”中，可以同时拖动所有行，但这不是一种选择，因为我必须保留单行的所有属性。 还有其他选项可以同时拖动所有行吗？ Edit1：刚刚找到了CadNodeTool（https://github.com/wonder-sk/CadNodeTool）。使用此插件，至少已经可以同时选择多条线的顶点了。也许可以扩展此插件以移动平行线的顶点。

11 qgis  pyqgis  editing  vertices 

我正在从Amazon S3下载Sentinel-2磁贴，但它们似乎没有100％对齐。右边的图像是添加到ArcGis Desktop的波段08 Sentinel-2（我所做的不过是将08.jp2添加到地图上），左边的是esri底图照片。Sentinel中的偏移量约为35-40米。我该如何纠正它们？ 这是下载的磁贴。等级1C，前哨2： 解： 从HDunn的答案，我做到了： 我使用道路矢量层添加了五个点对，并计算了x和y坐标的差。我在Shift工具中将平均差用作xvalue和y值。结果：

11 arcgis-desktop  remote-sensing  sentinel-2 

我有 NE lat lng和SW lat lng。我的目标很简单，我使用的是RStudio，我想从上述两个纬度和经度创建一个边界框，然后生成多边形，并想向正方形多边形添加随机点。我在CSV文件中包含NELat / Lng和SW-Lat / Lng，可以将其导入到data.frame中。 我只需要一点指导就可以开始使用，对于R语言来说我是非常陌生的人，它与我以前所使用的完全不同。 到目前为止，我有这个： coords = cbind(78.46801, 19.53407) coordsmax = cbind(78.83157, 19.74557 ) sp = SpatialPoints(coordsmax) sp2 = SpatialPoints(coords) r1 = rbind(coords, coordsmax[1, ]) # join P1 = Polygon(r1) Ps1 = Polygons(list(P1), ID = "a") plot(Ps1) 我只是在其他地方看到了一些示例，但无法绘制4个角的多边形。

11 r  rgdal 

如何使用细线显示水体来重现旧地形图的外观？我查看了填充选项，找到的最接近的是形状爆裂填充，但它使用了颜色渐变而不是线条。 此图显示了我要创建的效果： 图片来源Wikipedia Commons 我认为这种样式的名称称为“水衬线”，之所以使用该名称，是因为平版印刷版无法打印纯色区域。 有没有办法在QGIS中创建这种样式？

11 qgis  polygon  style 

我试图从ArcGIS Pro 1.2.0导出PNG图像文件，且背景设置为透明。在ArcMap 10.3中，这很简单，在导出过程中有一个用于将透明颜色设置为白色的按钮。 当我在ArcGIS Pro中“共享”图像时，看不到用于设置背景色或透明色的相同选项。此按钮是否已移至其他位置？还是不再有选择了？ 目的是找到与上面的屏幕截图类似的设置。对于这个项目，我有一个底图作为独立的PDF文件。我需要在该文件的顶部覆盖一个标注框-这要求标注框必须具有透明性，而不是白框（例如在裁剪时会发生什么）。

11 arcgis-pro  transparency  png 

我正在运行Ubuntu 14.4，并且已经安装了python虚拟环境。然后我pip安装了几个python库（gdal，lxml，psycopg2等） 同样在我的系统中，我有用于QGIS的模块。因此，当我输入python时，我可以将该模块导入为： import qgis.core 但是，当我尝试通过virtualenv执行相同操作时，出现错误：找不到模块。 创建virtualenv时，我使用了以下选项： virtualenv --no-site-packages env 这限制了我的环境不受系统环境的影响。问题是如何在我的virtualenv中访问QGIS模块？

11 qgis  python 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 3年前关闭。 我经常看地图，想知道何时创建。 我如何确定未日期的世界地图的年龄？ 如果我有世界政治边界的矢量文件，但是不确定日期，那么我可以使用哪些资源来确定日期？ 网络上是否有地图服务允许我输入日期并获取该日期的世界边界地图？ 假设存在，我将使用什么过程来确定匹配的接近程度？ 是否有特定于地图日期的基于规则的系统？

11 time 

使用以下代码，我可以将100％透明度设置为栅格值0： map=None for lyr in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr.name() == "some_layer": map= lyr tr=None tr=QgsRasterTransparency() tr.initializeTransparentPixelList(0) map.renderer().setRasterTransparency(tr) map.triggerRepaint() 如何定义两个值（0和8）为100％透明？ 当我用栅格值8重复最后一个块的命令时，第一个值（0）从图层透明度中删除。

11 qgis  python  pyqgis  transparency 

在最近的培训课程中，有人问我QGIS是否可以自动计算使用地图集生成器创建的地图册的下一页/上一页以及上一页/下一页。如果您知道网格的宽度和高度，那么我设法为常规网格制定了一个相当合理的标签表达式。 但是随后我们开始考虑一些现实的示例，在这些示例中我们不想绘制不包含我们感兴趣的地区的页面，例如我的家乡之一： 因此，今天下午我在一个python脚本上玩耍，计算出每个网格单元我感兴趣的4个邻居，并将这些值添加到我的网格中（这很大程度上基于Ujaval Gandhi的教程）： for f in feature_dict.values(): print 'Working on %s' % f[_NAME_FIELD] geom = f.geometry() # Find all features that intersect the bounding box of the current feature. # We use spatial index to find the features intersecting the bounding box # of the current feature. This will …

11 qgis  pyqgis  print-composer  atlas 

我有一个QGIS（2.14.3-Essen）项目，可以在办公室和家中使用。我保留该项目的所有内容，只有两层。 剩下的两层非常大-DEM和山体阴影总计超过20 GB-它们分别存储在工作场所和家庭中的不同路径中。 我想忽略Handle Bad Layers对话，并保留这些图层的折断路径。我常常不在乎他们是坏人。我想对项目的其他部分进行更改。 有没有办法忽略或跳过“处理不良图层”对话框，以便在不更新图层的情况下不会丢失对图层的引用？

11 qgis  layers