地理信息系统

如何找到保证在PostGIS中给定多边形内的点？ 我知道这个ST_Centroid功能。但是，质心并不总是在多边形内，请参见下文： 此外，我想避免使用一个在多边形边界上的点，而是想要一个在边界内的点（而不是在甜甜圈形多边形的孔内）。

22 postgis  polygon  point  geometry  point-in-polygon 

我已经使用了R与QGIS 的集成，并发现它是非常实用的工具。看过R工具箱的9 扩展和10的扩展后，它们似乎都只提供点聚类来演示概念验证。 R中可用的其他分析方法是否可以通过python脚本在ArcGIS Desktop中实现（例如，实现spatstat包的特定部分）？

22 arcgis-desktop  r  spatstat 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 4个月前关闭。 我正在创建一个目录，该目录使用位置对所有列表进行分类。我想提供用户可以选择的位置列表，而不是依靠用户输入。我已经尝试过GeoBytes提供的GeoWorldMap，但存在某些问题。例如，它列出了275个国家。除此之外，他们的数据库会定期更新，这意味着，如果我维护本地数据库，则需要不时更新它，以确保它是最新的。 我想知道开发人员为此使用了什么。他们是否使用数据库而不是数据库来获取位置列表？ 恐怕手动用户输入不是一个开放选项，因为按位置进行分类将是必不可少的，因此我不能信任或依赖用户手动输入的内容。 我见过很少有网站这样做，所以我确信那里有一个实用的解决方案。

22 data  geolocation 

我是Web映射的新手，因此这个问题可能会被误导。 我希望有一个使用OpenLayers（甚至更简单的Google Maps）的非常简单的Web界面，它可以显示PostGIS中当前的许多点和多边形表。 我希望这尽可能简单（且维护成本低）： 我不想创建图块（除非有充分的理由），也不希望拥有庞大，复杂的中间件层。 性能不需要是惊人的，因为这是供内部使用的，因此不需要扩展到几个并发用户。 矢量层的样式至关重要。 开源工具是首选，因为我对此几乎没有预算。 推荐的方法是什么？

22 postgis  openlayers  web-mapping  software-recommendations 

建议的简化几何形状的方法是什么？记住投影，例如简化状态几何。 我听说过要转换为“等距”投影，该投影可以简化而不会失真，然后再转换回您选择的投影。

22 coordinate-system  simplify  generalization 

我正在尝试使用QGIS Atlas工具来构建学生的“散点图”，即。显示215所学校的学校和该学校的学生。 学校层是覆盖层，并受SCHNUM（学校编号）值的驱动。这工作正常，即正确驱动页面。 但是，如何仅显示来自其他层SCHNUM的要素（在这种情况下为学生点层），这些要素与coverage层中当前要素的对应值匹配？ 我还想将覆盖层（学校）切换为“隐藏”，并使用相同的技术仅显示来自重复学校层的当前学校。 这里有一个演示文稿，暗示了使用的图层渲染规则设置$atlasfeatureid，但是我不知道它是如何工作的。 我还找到了这个错误跟踪器页面，该页面显示了此工作方式，但是同样，我也无法将其用于我的数据：

22 qgis  atlas 

我在QGIS中创建了2个多边形。在R中使用它们，多边形将自动变为SpatialPolygonsDataFrame（SPDF）。我想将它们合并为单个SPDF（这在ArcGis中使用Tool Merge极为简单）。我确信应该有一种简单的方法来完成R中的操作，但是我找不到方法。该合并功能似乎只合并data.frames，聚合函数溶解多个多边形为一个SHP，gIntersect（由打字加入功能）返回逻辑值，而不是在所有的SPDF。 数据可在此处获取：http : //ulozto.cz/xpoo5jfL/ab-zip library(sp) library(raster) library(rgeos) library(spatstat) library(rgdal) library(maptools) setwd("C:/...") a<-readOGR(dsn=getwd(), layer="pol.a") b<- readOGR(dsn=getwd(), layer="pol.b") ab<-merge(a, b) # what tool if not "merge" to use??

22 r  spatial-analyst  merge  maptools 

我花了整整两天的时间来了解中国臭名昭著的GPS转换的本质。筛选出噪音后，有两个重复出现的症状： GPS坐标（WGS-84）在中国地图上的绘制效果不佳，偏移范围为50至1000米。 中国的Google和Apple Maps的街道地图和卫星图像之间存在严重的偏差（现场例子） 各种 来源提出矛盾的主张： 出于安全考虑，冷战时期的法规要求将随机偏移量引入GPS芯片 保护主义措施旨在青睐中国制造的GPS设备 中国只是使用了不同的投影系统，因此偏移量不是随机的 中国制造的GPS芯片会抵消返回的坐标，以匹配必须遵循GCJ-02基准的批准的中国地图。 实际的故事是什么？这看起来像是应该在Wikipedia页面上普及的问题，但是我还没有找到任何权威的方法。 建议的两个可能重复的问题并不能完全回答原因。而是描述症状。

22 gps  wgs84  china 

我已将相应的“本地企业合作伙伴”表数据连接到包含UK Local Authority的多边形shapefile中。 如屏幕截图所示，如何溶解或合并LEP_NAME的“本地权限”多边形？ 理想情况下，我想保留带有合并后的多边形的“本地权限”属性。我的脚本知识为零，因此我正在寻找基于工具的解决方案。任何帮助或指针，不胜感激。 附加说明：尝试使用溶解工具时，LEP_NAME字段不可用，如上所示。LEP属性数据是在Excel中创建的，另存为CSV格式，并已加入到Local Authority多边形shapefile中。如其他屏幕快照所示，有重复的IDENTIFI0字段-这可能是问题的根源吗？

22 qgis  polygon  fields-attributes  merge  dissolve 

是否有可以在GIS环境中使用的虚构世界（例如，中土世界，Discworld等）？ 例如，MBTiles存档和shapefile或可与OpenLayers / Polymaps / Leaflet一起使用的任何其他格式？

22 data 

在QGIS中，我有两行。端点彼此靠近，但不完全对齐。我想将两条线合并为一条，将两个端点合并为中点。 我最接近的功能是“合并所选功能”。这给了我一个功能，但是它由两条不相交的线组成-端点之间没有线。 我正在寻找行为与Illustrator的“连接路径”功能类似的对象（“对象”>“路径”>“连接”）。 QGIS会这样做吗？

22 qgis  polyline-creation  editing 

我正在寻找一个函数来分解表中多边形要素之间的共享边界。ST_UNION（）几乎可以满足我的要求，但是它会根据图层中的所有多边形创建一个多面体，无论它们是否共享公共边界。我宁愿只溶解彼此接触的多边形之间的边界。我认为，应该有某种使用ST_TOUCHES（）的方法，但是对于溶解函数的需求似乎太普遍了，如果没有内置函数可以实现这一点，我会感到惊讶。 用例如下所示：我下载了一个欧洲大国的Corine Landcover数据，我想消除不同森林类型之间的界限（一张桌子中约有75,000个多边形）。我尝试了ST_UNION，但由于出现“内存不足”错误而失败了（不过30,000个多边形确实起作用了）： create table corine00 as select st_union(the_geom) as the_geom, sum(area_ha) as area_ha, substr(code_00,1,2) as code_00 from clc00_c31_forests group by substr(code_00,1,2) 注意：所有森林代码均以“ 31”开头，而我使用的是PostGIS 1.4，GEOS版本：3.2.0-CAPI-1.6.0

22 postgis  dissolve 

我正在创建虚拟地图，因此我需要创建许多点，线和确定的多边形。稍后，我将数据导出为geojson。但是在此之前，我必须离开并给每个元素一个唯一的ID。 我不需要特殊的排序，例如最大的多边形得到的ID最小。我只是需要所有多边形 与一个ID结尾，没有 做它手动像我现在要做的。 如果有人知道该怎么做会很棒。

22 qgis  polygon  unique-id 

我是GeoJSON的新手。如图所示，我有一个GeoJSON功能集，并希望将其存储在postgres表（测试表）中。我的postgres表有一个序列号和几何列。 { "type": "FeatureCollection", "features": [ { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ 2565453.1826721914, -3835048.659760314 ] } }, { "type": "Feature", "geometry": { "type": "LineString", "coordinates": [ [ 2727584.7219710173, -3713449.1942418693 ], [ 2732476.691781269, -3992291.473426192 ] ] } }, { "type": "Feature", "geometry": { "type": "Polygon", "coordinates": [ [ …

21 postgis  geojson  geometry 

我忘记了特定GIS原理的名称，需要复习... 场景：假设有一个代表某些管理区域的多边形。为了简单起见，假设它是一个正方形。我们还假设该广场内某处有9栋房屋，但我们不知道它们在哪里。如果将这个正方形分成9个相等的部分（例如井字游戏板），我想计算每个正方形内的房屋数量。 现在，很明显，分成9个小正方形的9栋房子就是每平方1个房子。但是，我知道所有九座房屋都可以在任何地方。它们全部九个都可以在右上角这样： 或者，最下面一行可能每个都有三个，而最上面的六个可能都不像这样： 如果我这样假设每平方一间房子， 这个错误的假设被称为什么？这是可修改的地域单位问题（MAUP）吗？

21 gis-principle