地理信息系统

有没有办法在QGIS中获得矢量层的范围（bbox）？ 我看到可以更新范围，但是我正在寻找范围的实际坐标。

11 qgis  layers  vector  extents 

我正在使用shapely非常简单地计算折线的长度： from shapely.geometry import LineString ... xy_list = [map(float,e) for e in xy_intm] line = LineString(xy_list) s = '%s,%s,%s' % (fr,to,line.length) 我的坐标在WGS84中。我似乎找不到有关shapely的length属性的任何信息。长度属性的单位是什么？有没有简单的方法可以转换成公里或米？

11 shapely  length  units 

我正在开发将GeoJSON渲染到OpenLayers地图中的应用程序。 这个想法是关于不从外部服务（OSM或Bing）下载PNG图像，而是在GeoJSON文件中具有地理信息并实时呈现它们。根据我的研究，这种技术似乎被称为'vector tile'。我看过一些家伙的网站，例如 http://openstreetmap.us/~migurski/vector-datasource/ 和 http://mike.teczno.com/notes/postgreslessness-mapnik-vectiles.html（相同家伙）。 他的工作似乎是我想要实现的目标，但我并没有真正了解它的工作方式。 据我了解，由于我使用的是OpenLayers，因此客户端会使用GeoJSON格式检索地图数据，然后使用OpenLayers的GeoJSON解析器将其转换为矢量。但这只是我不了解的这项技术中使用的术语。 在他的博客中，他描述了将有一个服务器提供称为'MapNik Vector Tile'格式的数据。矢量图块和GeoJSON有什么区别？它仅表示与边界框对应的平铺的GeoJSON吗？还是新的图像类型？ 如果看一下他的示例，您可以观察到道路上的标签与道路的笔触对齐。我尝试通过在OpenLayers中更改Renderer代码在OpenLayers上执行相同操作，但我只能旋转标签，而不能使它们沿道路弯曲（例如，当道路有拐角时，地图以“ L”形显示标签。* * + 渲染效果非常好！ 看起来Mike（博客的作者）除了使用OpenLayers以外，还使用了其他地图，但我不知道他是如何做到的。有很多问题要问，但是我想首先理解以上内容。

11 openlayers-2  openstreetmap  vector  vector-tiles 

我有多个多边形组成，有几十万的地图点非常小多边形上内它。 我需要较大的多边形，但不需要点较小的多边形。如何在保留较大多边形的同时选择并删除所有点的小多边形？（我显然不想手动执行此操作-有太多的要点！） 我对QGIS完全陌生，并且在Google上搜索过很多，没有运气。可能是因为我不确定Google搜索中使用的正确术语是什么！ 我正在使用的文件是.dbf，.prj，.qpj，.shp和.shx。我只是将.shp文件拖到QGIS上，地图显示为一层。基本上，英国的地图分为多个区域（我想要），并以点数覆盖（不需要）。 这就是近距离的样子。我进行了一次观察，看看他们是否可以进入周围的环境，但是没有运气：

11 qgis  polygon  point 

我有几百个形状（polygons和multipolygons），每个形状都由数以万计的点组成，这些点正试图进入SQL 2008。 不幸的是，我尝试导入的形状是“右旋”的（每个形状的周长围绕其包含的点顺时针绘制）。SQL Server至少对于geography类型采用“左手”形状（围绕内部逆时针旋转）。这意味着SQL假定我要选择除形状之外的整个地球。有人将其描述为“由内而外”的形状。 从MSDN来看，令人沮丧的是没有说要使用哪种环定向： 如果使用geography数据类型存储空间实例，则必须指定环的方向并准确描述实例的位置。 如果您在SQL 2008中使用错误的环向，它将崩溃并显示以下错误（强调我的意思）： 用户定义的例程或聚合“地理”的执行过程中发生了.NET Framework错误：Microsoft.SqlServer.Types.GLArgumentException：24205：指定的输入不代表有效的地理实例，因为它超出了一个半球。每个地理实例必须适合一个半球。导致此错误的常见原因是多边形的环方向错误。 导入形状geometry而不是导入形状geography效果很好，但我想使用geography。 在SQL 2012中，解决此问题显得微不足道，但我与2008紧密相关。 我应该如何转换形状？

11 sql  sql-server-spatial  geography-data-type 

我正在尝试使用gdalwarp将对地静止映射为Lambert保形。我的输入数据在netcdf中，并且在地理坐标（度）中，我想将重新映射的数据输出到netcdf。我为输入的netcdf数据创建了一个相应的vrt文件。Gdalwarp将输出netcdf文件，但输出数据全为零，并且我收到以下错误： Creating output file that is 5120P x 5120L. Processing input file netcdf.vrt. ERROR 1: Too many points (441 out of 441) failed to transform, unable to compute output bounds. Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping. 0...10...20...30...40...50...60...70...80...90...100 - done. 我尝试了以下命令： /usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

在以GeoJSON规范说“的默认CRS是地理坐标参考系，用WGS84数据，并用经度和十进制度的纬度单位”。似乎没有指定点的海拔/高度/ z单位。 GeoJSON是否有默认的海拔高度单位？如果没有，是否存在通用约定或指定单位的方法？

11 geojson 

长期以来，我一直在使用基于Windows的软件来可视化和分析用于森林应用的LiDAR数据集。最近，我开始将所有工作移至Ubuntu平台，但我仍在寻找可与LiDAR数据集配合使用的优秀软件。 与Windows有关，我体验过的最好的免费软件是： 仅用于可视化： a）Quick Terrain Reader能够打开预建的数字高程模型（DEM）和点云，并允许用户以快速直观的方式自由地在地形中移动。 b）FugroViewer是一款功能强大，易于使用的免费软件，旨在帮助用户充分利用其地理空间数据。 c）PointVue LE是一个免费的3-D LIDAR可视化工具，可用于以ASPRS LAS版本1.1格式可视化LIDAR数据。 用于分析和可视化： a）FUSION / LDV提供对LiDAR和terrain数据集的快速，高效和灵活访问。 自从我迁移到Ubuntu以来，很难替代我以前使用的软件。到目前为止，我发现了一些替代方法，但不如Windows更好： 为了可视化： a）LiDAR可视化被实现为核心外的多分辨率点云渲染器。渲染器能够使用固定大小的内存缓存以交互式帧速率可视化我们目前拥有的最大LiDAR扫描，其中包含多达113亿（11.3 * 109）个采样点。 b）LAG是用于可视化，检查和分类LiDAR点云的工具。目前，它支持LAS和ASCII文件格式。 不幸的是，我找不到像FUSION这样的要处理的东西。 是否有人体验过与LiDAR文件分析相关的Linux优秀免费软件？ 不要只列出其他软件，而要解释为什么使用并推荐它！

11 geoprocessing  lidar  visualisation  linux  fusion 

我主要使用QGIS并有一个大型项目要完成。完成后，我需要将数据提交给ArcGIS用户。数据将主要由点和多边形形状文件组成。 如果我知道要把它交给使用ArcGIS Desktop的人，如何从头开始设置QGIS项目？

11 qgis  arcgis-desktop 

我一直在尝试将我的传单地图的高度更改为引导程序内的某个百分比，但是每次我做时地图都不会绘制。因此，我总是必须恢复为px值。我很确定这是我是CSS新手，因此我不希望使用的简单设置。这是我的CSS。 <style type="text/css"> body { padding-top: 60px; padding-bottom: 40px; } #map { height: 75%; } </style>

11 web-mapping  leaflet  bootstrap-framework 

使用ArcGIS 10.1 for Desktop，我试图在两个大型数据集上运行“身份”功能，但是每次都会失败。我认为这不是缺乏硬件，因为我可以在同一数据集上运行其他进程。 我已经更新了软件，并添加了64位地理处理加载项，检查并修复了几何，并完成了KB29559中的修复（将ARCTMPDIR env变量设置为可写的本地文件夹）。 Executing: Identity DudleyStamp OA2011 "D:\Case Studies\National2\National.gdb\DudleyStampOA11" ALL "1 Meters" NO_RELATIONSHIPS Start Time: Wed May 08 21:11:43 2013 Reading Features... Processing Tiles... Assembling Tile Features... ERROR 999999: Error executing function. The table was not found. The table was not found. [DudleyStampOA11] The table was not …

11 arcgis-desktop  arcgis-10.1  identity  error-999999 

我创建了一个Python加载项按钮，通过将一个要素类属性复制到另一个要素类来帮助加快我的同事的工作流程。它使用arcpy.UpdateCursor函数更新目标要素类中的一行。现在，无论编辑模式如何，都可以运行此按钮脚本。显然，当脚本在编辑会话中运行时，用户可以选择停止编辑而不保存更改，但是当脚本在编辑会话之外运行时，情况并非如此。 如果ArcMap当前不在编辑会话中，如何向脚本添加检查以停止脚本运行？ 这关系到ArcMap 10和10.1 我还想与其他ArcMap用户进行核对，以验证在没有进行编辑会话的情况下，通常不允许对表进行更新。 那么该脚本如何在编辑会话之外运行？ 该脚本还提出了另一个问题，即当我从列表更新第二要素类表时，ArcMap执行的看似偶然的选择顺序正好对我有用，但这又是一天了。 这是现在可以正常运行的脚本（没有任何10.1编辑器实现）： 如何添加检查以确保用户处于编辑会话中？ def onClick(self): #Reference mxd mxd = arcpy.mapping.MapDocument("CURRENT") #Reference the main Data frame mm = arcpy.mapping.ListDataFrames(mxd, "MainMap")[0] #Reference the Water System Valve feature class waterValves = arcpy.mapping.ListLayers(mxd, "Water System Valve", mm)[0] #Reference the fire hydrant feature class fireHydrants = arcpy.mapping.ListLayers(mxd, "Water Hydrant", …

11 arcpy  arcmap  editing 

我想知道是否有人可以澄清使用ArcGIS IDW空间插值时是否应将点数据投影到等距投影？ 我正在研究来自北美西部的数据集，其跨度约为30度。数据当前为纬度/经度（NAD83）。 在插值过程中计算到样本点的距离时，ArcGIS是“即时投影”还是以某种方式调整了纬度？或者我应该在保留距离的投影中提供所有内容吗？

11 arcgis-desktop  coordinate-system  interpolation 

有没有办法从python控制台调用QQGIS中的clip函数？可在矢量菜单的地理处理工具下找到它。

11 pyqgis  clip 

当我在另一个Postgres图层中输入新元素（功能）时，可以通过两种方式进行操作： 绘制新元素（带有“添加功能”），我很少这样做或 复制（或剪切）另一个Postgre层（Source层）中的某些元素，并将其粘贴到我经常执行的目标层中 在第一个示例中，保存编辑工作正常，因为此层从postgre数据库序列获取gid * nextval（'layer_name_gid_seq':: regclass）* 在第二个示例中，在保存编辑时出现错误，因为在将元素从源层复制到目标层时，qgis 从源层复制 了元素的gid。尝试保存编辑内容时，将返回以下错误： 无法提交对“ Cjevovodi”层的更改。 错误：错误：未添加1个功能。 提供程序错误： 添加功能时出现PostGIS错误：错误：重复的键值违反了唯一约束“ cjevovodi_okill_pkey” 详细信息：键（gid）=（5）已存在。 我试图在字段gid中复制* nextval（'layer_name_gid_seq':: regclass）*，但是由于字段被定义为数字，因此无法将该序列粘贴到字段gid中。 有谁知道从源层复制元素的简单方法（使用现有的gid）分配新的gid？ 谢谢！

11 qgis  postgis