地理信息系统

我有一个小脚本，用于从GIS软件当前不支持的数据格式生成GeoJSON数据。 我想验证输出，但是当我使用ogr2​​ogr读取输出时，它没有为我提供有关哪些方面，行或字符包含无效数据的任何有用信息。 是否在某处有一个不错的GeoJSON验证器，或者使用OGR验证GeoJSON数据的好方法？ 编辑：我不是在寻找JSON验证程序或GeoJSON规范。我已经看到并使用了两者，但是OGR仍然返回错误。 谢谢！

16 geojson 

有没有一种方法可以计算到QGIS中到最近点的距离？在ArcGIS中，有一个称为Euclidean Distance的函数，该函数可以获取输入矢量或栅格，并生成距离栅格。 QGIS中有类似的功能吗？我在文档中找不到一个。我确实研究了RasterQGIS中的功能以及MMGIS中的Hub Distance计算（对于非光栅解决方案），但是我无法弄清楚如何使用这两种方法。如果向量方法是唯一的解决方案，那也可以。

16 qgis  raster  distance  point 

从几何字段获取SRID的功能是什么。 我已将一个shapefile插入定义了Geometry字段的SQL Server表中。我知道shapefile的投影为Swreff99_1800（EPSG：3011）。 现在我的问题是如何通过使用几何字段的某些OGC方法获得SRID：3011？

16 geometry  sql-server  srid 

我是Python的新手，我想知道是否存在一种快速方法来在ArcGIS 10中使用Python逐个像素获取栅格像元值以及坐标（每个像素中心的地图XY坐标）？ 为了进一步描述这一点，我需要获取第一个像素的地图X，地图Y和像元值，并将这三个值分配给三个变量，并对其余的其他像素重复此步骤（遍历整个栅格）。 我想我需要更多地描述我的问题。问题是，我需要获取第一个栅格的像素的XY位置，并获取与该XY位置相对应的其他几个栅格的像元值。此过程应遍历第一个栅格的每个像素，而无需创建任何中间点shapefile，因为这将真的非常耗时，因为我必须处理具有近80亿像素的栅格。另外，我需要在ArcGIS 10中使用Python进行此操作。 @JamesS：非常感谢您的建议。是的，这将适用于一个栅格，但是我还需要收集其他几个栅格的像元值。问题是，在获取第一栅格的第一个像素的X和Y坐标之后，我需要获取第二栅格的像元值，该值对应于第一栅格的X，Y位置，然后是第三栅格，依此类推。因此，我认为在遍历第一个栅格时，应该同时完成像素的X和Y位置以及与该位置相对应的其他栅格的像元值的获取，但我不确定。这可以通过将第一个栅格转换为点shapefile并在ArcGIS 10中执行将多值提取为点功能来完成，但是我 @hmfly：谢谢，是的，如果我可以获取数组的已知行和列值的坐标，则此方法（RastertoNumpyarray）将起作用。 @whuber：我不想执行任何计算，我要做的就是将XY坐标和单元格值写入文本文件，仅此而已

16 arcgis-10.0  python  raster 

关于所有数据中80％如何具有空间成分的报价在各地无处不在。GIS类简介，产品描述，讲座等。有人知道谁先说的吗？或者链接到原始文章会很好。

16 history 

如果以下问题有些愚蠢，我深表歉意，但是我对这整个GIS只是非常陌生。 我正在尝试使用python中的gdal将一些投影的geoTiff图像转换为WGS84。我发现了一篇帖子，概述了使用类似于以下内容的方法在投影的GeoTiff中转换点的过程： from osgeo import osr, gdal # get the existing coordinate system ds = gdal.Open('path/to/file') old_cs= osr.SpatialReference() old_cs.ImportFromWkt(ds.GetProjectionRef()) # create the new coordinate system wgs84_wkt = """ GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.01745329251994328, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4326"]]""" new_cs = osr.SpatialReference() new_cs .ImportFromWkt(wgs84_wkt) # create a transform object …

16 python  gdal  coordinate-system  wgs84 

新手可以在此处进行堆栈交换。 我目前正在使用此处找到的自动瓷砖切割器来制作为我的地图创建自定义地图类型所需的标准256x256瓷砖。使用地图拼贴切割器脚本感觉很笨拙，因此我很快用完了计算机上所有4GB的内存以打开photoshop并使用它。这样，在我的图像的顶部和底部留有空白，以使其正确放置非常麻烦。 有谁知道其他可用于制作Google地图图块的工具？

16 google-maps  tiles 

我正在将代码从sp包迁移到较新的sf包。我之前的代码中有一个多边形SpatialDataFrame（censimentoMap）和一个SpatialPointDataFrame（indirizzi.sp），并使用下面的指令获取了放置在其中的每个点的多边形单元格ID（“ Cell110”）： points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) 实际上，我创建了两个SF对象： shape_sf <- st_read(dsn = shape_dsn) shape_sf <- st_transform(x=shape_sf, crs=crs_string) 和 indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string) 而且我正在寻找上述说明的SF…可能是： ids<-sapply(st_intersects(x=indirizzi_sf,y=shshape_sfpeCrif), function(z) if (length(z)==0) NA_integer_ else z[1]) cell_ids <- shape_sf[ids,"Cell110"]

16 r  overlay  sp  sf 

是否可以将ArcGIS中的要素图层直接转换为GeoPackage？ 在ArcGIS Desktop 10.5的ArcToolbox中，在“转换工具/到地理包”下，向我显示的唯一选项是“将栅格添加到地理包”。没有“添加功能”之类的东西。 在大约3年前的这篇博客文章中，ESRI声称“在10.2.1或10.2.2的ArcGIS Desktop上，您可以创建一个空的GeoPackage并通过将要素数据复制到其中来填充GeoPackage”。但是，我看不到。在ArcCatalog中，没有创建GeoPackage文件的选项。

16 arcgis-desktop  convert  arccatalog  geopackage  arcgis-10.5 

我正在进行多边形（边界）清理，想知道是否可以将节点/顶点添加到一组选定的多边形甚至单个多边形中。 顺便说一句，QGIS拓扑编辑工具是令人难以置信的-非常简单但非常准确，并且从手工工作中省去了很多工作。 QGIS提供了删除选定节点的功能-惊人的是，在清理多余节点时，拓扑得以保留甚至得到改善。好东西！

16 qgis  editing  vertices 

我在postgis中有一个线串表和一个点​​表。 我知道最接近任何给定点的线。我需要知道的是那条线的哪一侧。我想我必须通过创建从给定点到该线的垂直线（该线上的最接近点）然后比较坐标来做到这一点，但是我不知道该怎么做，如果这是正确的方法，因为线改变了方向。 我做了一张照片来说明我的任务。 该线本身为黑色，其方向以绿色箭头显示。我需要在点表中添加一个“侧”列，以便红点的值应为“ right”，蓝点的值应为“ left”。 有人可以给出一个计算点的“边”值的SQL代码示例吗？

16 postgis  point  linestring 

如何在QGIS中获取栅格中像元值计数的摘要？我知道QGIS就是这样做的，以便在层的属性中提供摘要统计信息和直方图。我不需要摘要统计信息或直方图图像，而无需计算单个值。 似乎曾经尝试过为此添加特殊功能（功能＃4321，http://hub.qgis.org/issues/4321）。 我是否真的缺少明显的东西？

16 qgis  raster 

Python的geopy软件包具有两种距离测量技术：Great Circle和Vincenty的公式。 >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) - https://goo.gl/maps/TQwDd >>> p2 = (31.8300000,35.0708167) # (lat, lon) - https://goo.gl/maps/lHrrg >>> vincenty(p1, p2).meters 429.16765838976664 >>> great_circle(p3, p4).meters 428.4088367903001 有什么区别？首选哪种距离测量？

16 python  distance  great-circle  geopy  vincenty-formulae 

ArcGIS 10.2中提供了哪些方法来随机选择点集。例如，在所附的屏幕截图中，我有兴趣保留20％的选定点并删除其余的点。

16 arcgis-desktop  arcpy  editing  select  random 

我正在寻找有关为用Python编写的QGIS插件编写自动测试的建议。 过去，我使用PyUnit（unittest模块）编写了针对Python脚本的测试，但从未针对具有GUI的应用程序进行过测试。我已经找到了一个页面，描述了如何使用PyQt4.QTest对Qt小部件进行单元测试（http://www.voom.net/pyqt-qtest-example），但是我正在努力查看如何使用它带有设计为从QGIS内部运行的小部件。 PyQGIS文档中关于“测试”的部分明显不存在。 到目前为止，我有： 将实际的数据处理保存在隔离的模块或功能中，并为这些模块或功能编写单元测试； 使用QTest对UI进行基本测试； 当从QGIS内部使用插件时，请祈祷一切融合在一起。 有没有更好的办法？

16 pyqgis  qgis-plugins  unit-testing