地理信息系统

我想在Raspberry Pi上运行GDAL QGIS，以充当负担得起的GIS处理引擎。该计划是购买大约50 Pis。 是否有人在Raspberry Pi上安装和运行GDAL或仅在RASPBIAN环境中安装GDAL方面有任何链接或信息？ 最好的是，Raspbian是我的首选，它是基于Debian的，但是ARCH，RASPBMC，PIDORA，RISC OS和OPENELEC是开箱即用的Linux风格。 希望使用Python。

11 gdal  linux  raspberry-pi 

我想按照本博文解释的方式设置PyCharm：为PyQGIS和Qt设置PyCharm。 我按照指示进行了所有操作，但是在添加项目解释器时，出现了以下错误消息： 有人对此有解决方案吗？

11 pyqgis  installation  pycharm 

这可能是一个幼稚的问题，但作为QGIS的新用户，我一直在努力。 我有一个非常大的shapefile（275,000个点，但是可以将其分成大约10个子区域，以便进行更快的处理）。 我想识别200米之内没有其他点的所有点，然后在文件的字段中用值“ unique”对每个点进行编码。 然后，对于属于本地集群的所有其他点，我想将其编码为“集群”。 实现这一目标后，我想为每个群集随机选择一个，以保留在数据集中，而丢弃其他群集。 目前，我无法完成第1步，因此欢迎您提供任何帮助。

11 qgis  clustering  spatial-query  map-point 

对于在ArcSDE中使用存储在ArcSDE地理数据库中的超大数据集的工作，哪一种方法能提供最佳性能？使用复杂的SQL查询标记特定功能的多个标签类，还是使用将相同的复杂SQL查询设置为每一层的定义查询的同一层的多个副本？

11 arcgis-desktop  labeling  performance  definition-query 

我正在尝试迭代模型构建器中的一系列要素类（多边形），以将一系列多边形转换为栅格，但是，输出文件名存在问题。 将“迭代要素类”工具插入模型构建器窗口，并将存储所有多边形的gdb链接为输入后，输出要素（绿色气泡）将自动成为我的第一个多边形的名称。结果，我尝试运行的模型似乎重命名为该多边形，而不是转到列表中的下一个多边形。创建了输出栅格，但是将覆盖该文件名，而不是生成一个具有与后续多边形匹配的新名称的新栅格。 我究竟做错了什么？

11 arcgis-desktop  modelbuilder 

这是我所拥有的： ArcGIS Desktop Advanced（10.2）和MS SQL Server 2008 R2 这是我的问题： 我需要能够连接到ms sql服务器并导入，导出，创建和编辑空间数据。我可以连接到数据库并使用“要素类到地理数据库”工具将特征导出到数据库，但是一旦将数据带回到ArcMap会话中就无法编辑数据。我尝试使用高级桌面许可证附带的“创建企业级地理数据库”工具，但它会索要我没有的授权文件（假设用于ArcGIS for Server）。我已经查看了ArcGIS for Server的价格，这是不可行的选择。 是否可以使用当前拥有的许可证来实现我所需的替代方案？ arcSDE在哪里适合？

11 arcgis-desktop  arcgis-server  arcgis-10.2  sql-server 

我正在为陡坡上的水流建模以进行侵蚀研究。您可以从屏幕快照中看到基于ArcGIS的流量分析的结果。为了生成这些流线，我生成了一个无凹陷的DEM，计算了流向，最后生成了流累积。用ArcScan对栅格流积积进行矢量化处理。 尽管此方法在显示水路方面做得很好，但它并不能提供任何可视的量级或累积量度。我如何更好地表示陡峭表面上的水流，以便表示大小和/或堆积？

11 arcgis-desktop  cartography  hydrology  hydrography 

我正在执行一项任务，要求我在海岸线上每1000公里获取一次采样点，并且遇到了南极洲的问题。据我所知，在函数中使用几何图形似乎是一个问题，而实际上我认为应该将地理信息用于此操作。 使用这个非常相似的问题中的函数，我可以产生如下所示的结果：。 如您所见，ST_AddMeasure()并且ST_LocateAlong()似乎没有球形地处理几何，这导致许多点位于南极上。甚至沿日期线（左侧）在剪辑上添加了一个点。根据这两个函数的文档，只能使用几何。 可以在此处找到用于生成多边形和点的代码，但这是用于生成点的SQL： CREATE TABLE atest AS WITH line AS (SELECT id, ST_ExteriorRing((ST_Dump(geom)).geom) AS geom FROM line_sample_test), linemeasure AS (SELECT ST_AddMeasure(line.geom, 0, (ST_Length(line.geom))::int) AS linem, generate_series(0, (ST_Length(line.geom))::int, 10) AS i FROM line), geometries AS ( SELECT i, ST_LocateAlong(linem, i) AS geom FROM linemeasure) SELECT * from geometries; 如何在此海岸线上每1000公里产生一个点？

11 postgis  point  spherical-geometry  antarctica 

我正在尝试使用Python从栅格中提取自然中断值。 预期的过程进行如下：将Arcpy栅格对象转换为NumPy Array（使用RasterToNumPyArray）以自然中断值（使用PySAL Natural Breaks函数）。 import arcpy, pysal from pysal.esda.mapclassify import Natural_Breaks as nb # code to create greenIndex arcpy Raster object here greenArray = arcpy.RasterToNumPyArray(greenIndex) breaks = nb(greenArray,k=2,initial=20) 该代码返回错误“ ValueError：矩阵必须为二维”。 据我所知，greenArray是一个二维数组。

11 arcgis-desktop  arcpy  numpy  pysal  valueerror 

我想使用具有多个不规则多边形的多边形要素类，并让Arc在每个多边形内部绘制平行线。理想情况下，对于Arc来说，找出多边形最长边的角度并在该边上绘制平行线会很好，但是为了简单起见，我想如果我可以为所有平行线输入一个角度，那会更容易。 因此，线的角度，线之间的宽度，最小/最大长度以及从多边形侧面开始的缓冲区宽度是我的基本标准。 如果有帮助，请附上图像。

11 arcgis-desktop  arcpy 

我有一个包含数以千万计的点的栅格，这些点都定义了X，Y和Z数据。我需要找出一种简单的方法来提取和创建新的栅格或矢量，该栅格或矢量的所有点之间的距离均为720 m（X，Y），彼此之间的高度差为120 m（Z）。 我对SQL和Python具有零知识。我一直在尝试在VBA上执行此操作，并提出了一些算法，但是处理时间不合理且不切实际。我确信必须有一种简单的GIS方法来完成此任务，但似乎找不到它。 我正在使用ArcMap。

11 arcgis-desktop  vba 

R具有图像分割的能力，尽管我遇到的所有示例都使用单个波段进行分割（示例）。我对将R的随机森林图像分类能力与面向对象的分割方法结合起来很感兴趣。 R对于适合基于自然资源的分析的多光谱图像分割具有什么功能？或如何链接单个波段分割的结果以进行进一步分析。

11 r  remote-sensing  digital-image-processing  image-segmentation  random-forest 

我使用GeoServer在此手册下创建ImageMosaic存储的位置。当我创建一个新的数据存储时，所有*.tif文件都被索引到数据库（带PostGIS的PostgreSQL）中。当我将新的geotiff放入我的数据文件夹，或在数据库中重新索引所有geotiff文件时，如何在数据库中创建新的索引条目？ 我的配置文件： datastore.properties： SPI=org.geotools.data.postgis.PostgisNGDataStoreFactory host=localhost port=5432 database=GeoServerRaster schema=public user=geoserver passwd=******* Loose\ bbox=true Estimated\ extends=false validate\ connections=true Connection\ timeout=10 preparedStatements=true indexer.properties： Caching=false TimeAttribute=ingestion ElevationAttribute=elevation Schema=*the_geom:Polygon,location:String,ingestion:java.util.Date,elevation:Integer PropertyCollectors=TimestampFileNameExtractorSPI[timeregex](ingestion),DoubleFileNameExtractorSPI[elevationregex](elevation) timeregex.properties： regex=[0-9]{8}T[0-9]{6}Z(\?!.\*[0-9]{8}T[0-9]{6}Z.\*) levationregex.properties： regex=(?<=_)(\\d{1,5})(?=_) 目录中的文件： tc_10000_20140807T033115Z.tif tc_10000_20140807T040046Z.tif tc_10000_20140807T043115Z.tif ...

11 postgis  geoserver  geotiff-tiff  image-mosaic 

我有一个带有管理边界及其几何的postgres数据库。这些边界中的每一个都有一个标识号。 我要实现的目标： 首先，我想选择所有以特定邮政编码开头的行。例如： SELECT * FROM "post" WHERE "post"."ident" LIKE '101%' 所有相邻的几何应合并，不合并的几何应返回。 我已经尝试过以下语句： 我从这个答案中得到这一点： 使用PostGIS将许多小多边形连接起来以形成更大的多边形？ with t as ( select (st_dump(geom)).geom from "post" where "post"."ident" LIKE '593%' ) select ST_UNION(geom) from t; 但是此语句仅返回那些相邻但不包括在内的多边形 with t as ( select (st_dump(geom)).geom from "post" where "post"."ident" LIKE '593%' ) select ((ST_Dump(ST_UNION(geom)))) from t; …

11 postgresql  polygon  geojson  merge 

使用Python将Point插入PostGIS数据库的正确方法是什么？

11 python  postgis  psycopg2