地理信息系统

是否可以env.workspace = "in_memory"使用ArcGIS Pro和ArcGIS 10.2.2（或10.3）体系结构在ArcPy中进行设置？ 我想做的是将快照倾点操作的输出写入内存而不是磁盘。 我意识到我可以写入磁盘然后将其带入内存，但这无济于事。我正在尝试优化一系列任务，以最大程度地减少将仅作为下一个过程的输入所需的过程输出写入物理介质。

11 arcgis-desktop  arcpy  arcgis-pro  in-memory 

我正在尝试学习使用Python GDAL绑定和numpy进行遥感图像处理的绳索。作为首次尝试，我正在读取Landsat8 geotiff文件，进行简单的操作并将结果写入新文件。下面的代码看起来工作正常，除了原始栅格转储在输出文件中，而不是转储在操纵的栅格中。 欢迎提出任何意见或建议，但要特别注意为何未在结果中显示操纵的栅格。 import os import gdal gdal.AllRegister() file = "c:\~\LC81980242015071LGN00.tiff" (fileRoot, fileExt) = os.path.splitext(file) outFileName = fileRoot + "_mod" + fileExt ds = gdal.Open(file) band = ds.GetRasterBand(1) arr = band.ReadAsArray() [cols, rows] = arr.shape arr_min = arr.Min() arr_max = arr.Max() arr_mean = int(arr.mean()) arr_out = numpy.where((arr < arr_mean), …

11 python  gdal  geotiff-tiff  landsat  numpy 

SpatiaLite中的空间层（表）包括各种支持对象和元数据，包括（至少）geometry_columns表中的触发器，索引和条目。我正在寻找（最好是）一口气处理所有必要更改的GUI，或者寻找所有必需更改的（备用）文档，以免破坏空间层。 表有： 输入geometry_columns.f_table_name。 有五个名为的触发器[prefix]_[table_name]_geometry，其中前缀表示这是事务回滚还是索引更新。 命名为三个空间索引 idx_[table_name]_geometry[_suffix] 我已经在两个应用程序QGIS DB Manager和SpatiaLite-GUI中进行了尝试。 QGIS DB Manager出现具有以下效果： geometry_columns用新表名正确更新 不重命名触发器。触发器定义已部分修改，以BEFORE [INSERT|UPDATE|DELETE]引用新的表名，但条件仍在中搜索旧名称geometry_columns。 不重命名空间索引。不确定这是否重要，因为触发器仍引用旧的索引名称。 当您在SpatiaLite-GUI中选择“维护”→“重命名”表时，您仅获得一条SQL ALTER TABLE语句的存根。这是简单易懂的SQL，甚至比QGIS DB Manager还少。如果您填写新的表名，该表将被重命名。其他效果： 表格未在中重命名geometry_columns.f_table_name，这意味着许多GIS不会将表格视为空间层。 不重命名触发器。触发器定义已部分修改，以BEFORE [INSERT|UPDATE|DELETE]引用新的表名，但条件仍在中搜索旧名称geometry_columns。 不重命名空间索引。不确定这是否重要，因为触发器仍引用旧的索引名称。 请注意，Spatialite-GUI确实提供了恢复几何列的选项，该列将向其添加正确的条目geometry_columns（但要求您指定SRID，几何类型和尺寸），构建或重建空间索引以及恢复触发器，但没有这些功能中的一部分会删除旧的表行，触发器或索引，从而导致其中包含很多（可能无害但令人讨厌）schmutz的数据库。

11 qgis  spatialite  spatialite-gui 

我使用“ R”脚本下载MODIS图像。但是，它不提供MOD16产品。我在Earthexplorer网站上都找不到MOD16产品。虽然可以从http://www.ntsg.umt.edu/project/mod16#data-product获得所需的产品，但是我找不到要下载的任何脚本。 有可用的脚本吗？或我需要单击每个图像进行下载。 例如，如果我想从一年中的每个月下载“ h02v08 ”图像，该怎么做？

11 qgis  r  remote-sensing  modis 

我知道，我的问题类似于该网站上的一些旧问题。 我有很多CSV文件（地理坐标）要导入到qgis（然后进行转换），通常的方法不是最好的方法（太长）。 我有将近500个CSV文件（wgs84坐标），这就是我想要做的： 一次将所有CSV文件导入QGIS 投射他们 再次将它们导出到CSV文件中，但坐标不同（转换为UTM33N） 我试图了解如何使用python控制台，但我没有继续:( 谁能向我解释如何逐步实现它？

11 qgis  python  pyqgis  export  csv 

我有一系列代表道路标志的点要素，每个点都有与之相关的每个标志的一张照片。 我想创建一个地图簿，将标志的位置显示为点，并在该点旁边显示照片。地图册中的每个页面上都会有多个标志位置和照片。 有没有简单的方法可以自动将数据导入到点功能旁边的数据视图中的MXD中？ 因此大约有1000个自动化需求。还是以某种方式用照片“标记”功能？

11 arcgis-desktop  data-driven-pages 

在ArcMap中，很容易将逗号作为测量的网格上的数千个分隔符包括在内，如下面的屏幕快照所示，但是我一直无法找到一种方法来将使用的分隔符从逗号改为空格。 似乎没有执行此操作的高级ArcMap设置，我也一直在寻找但未找到用于此目的的ArcGIS Idea。 有没有办法让我忽略呢？

11 arcgis-desktop  labeling  grids-graticules 

我需要计算多边形（一些真实的地理对象，如湖泊，城市，国家等）的面积和相交面积。多边形位于俄罗斯的加利福尼亚州，新西兰，瑞典的Anadyr 所有多边形都在WGS84中。 我使用GeoTool Java API所做的工作： 使用EPSG：3488，EPSG：NAD83（NSRS2007）/ California Albers和计算的面积和重叠面积投影所有多边形。 使用World_Mollweide和World_Eckert_IV进行了相同的操作 为来自加利福尼亚，新西兰等地的多边形选择了“ 局部特定的投影 ” 我认为＃3是最准确的结果，因为我选择了覆盖多边形区域的投影 结果： 与＃3相比，＃2表现最差 '＃1和＃3区域和交叉区域的差异小于0.1％ 为什么？对于来自瑞典的多边形，我选择了绝对错误的投影EPSG：3488（加利福尼亚），得到的面积和交集面积几乎相同？ UPD：看来我没有正确解释我的困惑。这是带有说明的示例输出 #area_from_new_zealand_1 EPSG_27200 area[11733479] CRS[World_Mollweide] area[11736023] diff[2544] [0.0%] EPSG_27200 area[11733479] CRS[World_Eckert_IV] area[11736033] diff[2554] [0.0%] EPSG_27200 area[11733479] CRS[EPSG:NAD83(NSRS2007) / California Albers] area[11736034] diff[2555] [0.0%] #area_from_new_zealand_2 EPSG_27200 area[2952725] CRS[World_Mollweide] area[2953281] diff[556] [0.0%] EPSG_27200 area[2952725] CRS[World_Eckert_IV] …

11 coordinate-system  geotools 

想象一下，从法律的角度来看，您有一个要匿名的GPS轨迹。你会怎么做？捕捉到最近的x距离并剥离时间是否足够？是否有国际认可的标准？有没有人已经写过算法来做到这一点？如果没有，我计划在不断发展的stplanr程序包中创建一个函数来执行此操作。 可重现的示例（使用来自@geospacedman的出色旋转功能）来自我自己的' Identifiable '数据： library(rgdal) library(tmap) downloader::download("https://www.openstreetmap.org/trace/1619756/data", "test.gpx") r <-readOGR(dsn = "test.gpx", layer = "tracks") r <- spTransform(r, CRS("+init=epsg:27700")) rproj <- rotateProj(rs, 90) # rotate projection for plotting r <- spTransform(r, rproj) rs <- rgeos::gSimplify(r, 1000) # snap to nearest km qtm(r) + qtm(rs, line.col = "red") + tm_layout(draw.frame = …

11 r  gps  gpx  simplify 

我们同时在一个较大的层上工作（10-20人）。我们中的一些人使用QGIS 2.8.1和其他2.6。使用2.8的用户会发生很多崩溃（QGIS）。2.6用户没有崩溃，但是添加功能真是令人难以置信的缓慢。Postgres 9.4服务器在带有Win7的虚拟机（Hyper V）上运行。也许是网络而不是QGIS？也许有人有一个主意。会很好 ：）

11 qgis  postgis  postgresql  layers 

已关闭。这个问题是基于观点的。它当前不接受答案。 想改善这个问题吗？更新问题，以便通过编辑此帖子以事实和引用的形式回答。 4年前关闭。 MapBox 最近宣布，他们将从CartoCSS转向完全可视化的界面。显然，基于UI的新样式表将基于JSON，因此尽管新界面是点击系统，但从事此项目的开发人员之一仍然认为，对于编码人员而言，最终它将比以前的系统更好： 新的编辑器https://t.co/i7uy2V8YGL对于编码人员来说甚至更好：GL样式可以在编辑器中修改为JSON，并在JS中动态修改—汤姆·麦克赖特（@tmcw），2015年7月29日 但是，值得注意的是（并且希望我拥有此权利）这些JSON样式表由WebGL呈现，而CartoCSS用于创建由Mapnik呈现的XML。因此，尽管其中一种可以替代另一种，但至少在某种程度上可以替代，但它们的功能却不尽相同。 因此，这使我想知道CartoCSS会留在哪里？ 目前，许多主要的地图图块样式都使用该语言，也许最值得注意的是OpenStreetMap的默认样式。CartoDB也使用它，正如他们在下面的推文中指出的那样，他们的开发人员之一计划暂时坚持使用它，但这似乎也暗示他们将在某个时候放弃它。 @Anonymaps @tmcw但是，我们坚持使用它的时间更长了！— Francisco Dans（@fdans）2015年7月29日 在最近的一个方面，有一个名为kosmtik的项目，我通过此线程找到了。该工具读取与TileMill相同的.mss（cartocss）和.mml（json）文件，并使用Mapnik进行渲染，因此它确实为TileMill提供了替代方法，后者不再得到积极开发。 这个问题的最终根源是，我在一两个月前在TileMill进行了一个相当复杂的项目，也许要完成它需要25％到50％。但是随着这一公告的发布，我开始怀疑我是否正在将自己推向一个无法维护/发展的媒介。 所以我的问题是，还有什么其他证据可以证明GIS社区对CartoCSS的使用水平（或缺乏）？

11 carto  mapbox  cartocss 

关闭。这个问题是题外话。它当前不接受答案。 想改善这个问题吗？ 更新问题，使其成为地理信息系统堆栈交换的主题。 2年前关闭。 墨卡托投影（1536）必须非常适合航海。它可以让人们找到方位，并保持这一恒定不变，就应该达到目标（如果不在最短的直线上）。但是人们在此之前如何导航？例如，麦哲伦环游地球。1520。

11 navigation  history  mercator 

我正在寻找一个简单的point shapefile，它指示世界上的城市在哪里。我似乎找不到一个可以免费满足我要求的产品。香港专业教育学院得到了这个：http：//techcenter.jefferson.kctcs.edu/data/，但我也想要较小的城市... 有什么想法我可以得到吗？

11 shapefile  point  download  data 

shapefile是否支持混合几何类型（即点和面都在同一文件中）？

11 shapefile  geotools 

我对Postgres和PostGIS很陌生。 我正在尝试对一些GPS点进行简单的地图匹配（如果需要，也可以进行路径推断）。我正在使用postgres（postgis）。 我的数据库中有一个表，其中包含我的GPS点： gpspoints（纬度，经度，时间，速度，方位角，几何形状，...） 我也将道路网络的shapefile导入到我的数据库中： 路线（gid，idrte，version，nomrte，norte，clsrte，geom等）-我没有链接的方位角，而我不知道如何使用postgis函数ST_azimuth计算它。我所拥有的只是道路网的shapefile，其中包含上面列举的列。 仅当链接的方向与GPS点的方位角（+或-15度）一致并检索投影位置时，我才想将每个GPS点关联到最近的链接（在该点周围20米的缓冲区中）。否则，我希望它在20米缓冲区内搜索具有可接受方位角的NEXT最近的链接！（就像图片一样！） 我希望将投影GPS点的新坐标以“ projectedLat”和“ projectedLong”的形式添加到表“ gpspoints”中。 （在下图中，实际点是使用方向显示的，而投影点没有任何方向标记）

11 postgis  gps  map-matching