地理信息系统

制图师,地理学家和GIS专业人员的问答

5
根据点将线划分为不重叠的子集
给定一张具有线几何形状的表,并在单独的表中捕捉到该线的一个或多个点,我想在每条线与一个点相交的位置处使用一个或多个相交点来分割每条线。 例如,一条直线L沿直线几何顺序依次具有三个相交点A,B和C。我想将L作为四个不同的几何返回:从L的起点到A,沿着L的从A到B,沿着L的从B到C,以及从C到L的终点。 过去,我经常使用shape来完成此任务,这是一个线性引用问题(http://sgillies.net/blog/1040/shapely-recipes/)。但是,在这种情况下这是不切实际的,因为它有数百万条线和点。相反,我正在寻找使用PostgreSQL / PostGIS的解决方案。 注意,点被限制在一条线上。此外,一个点可以有效地位于一条线的起点或终点,在这种情况下,无需拆分该线(除非存在其他与同一条线的起点或终点不一致的点)。子线需要保留其方向和属性,但是点要素的属性无关紧要。

2
在10.2上运行ArcGIS 10.1应用程序
我支持已在装有ArcGIS 10.1的PC上编译的.NET应用程序。该应用程序包括一个用户界面和一个ArcGIS桌面扩展,以及一个Windows Installer调用的安装程序程序集,用于在安装过程中注册该桌面扩展。 我的问题是该应用程序无法在装有ArcGIS 10.2的PC上运行,因为未找到10.1 ESRI程序集,因此它崩溃了。从v9.3 ESRI开始,似乎不再提供发布者策略文件,从而在运行时将程序重定向到最新版本。 通过将运行时程序集绑定添加到应用程序的.exe.config中,我设法重定向了程序以加载10.2程序集。但是,这不会重定向桌面扩展,也不会解决安装程序在安装应用程序时寻找10.1库的问题。 ESRI的解决方案似乎是每次发布新版本的ArcGIS时都重新编译应用程序。但是,这意味着我将需要拥有最终用户使用的每个ArcGIS版本的副本。每次需要更改应用程序代码以重新编译时,我都必须卸载并重新安装不同的版本。 有没有人有什么建议?是否可以在运行时动态加载库以避免寻找特定版本?还有其他解决方案吗?

5
栅格差异:如何检查图像是否具有相同的值?
是否可以检查给定的2个栅格图层是否具有相同的内容? 我们的公司共享存储量有一个问题:它太大了,需要三天以上的时间才能进行完整备份。初步调查显示,占用空间最大的罪魁祸首之一是应使用CCITT压缩将开/关栅格真正存储为1位图层。 该示例图像当前为2位(因此有3个可能的值),并保存为文件系统中11 MB的LZW压缩tiff。转换为1bit(所以有2个可能的值)并应用CCITT Group 4压缩后,我们将其压缩到1.3 MB,几乎节省了一个完整的数量级。 (这实际上是一个表现良好的公民,还有其他人存储为32位浮点数!) 这真是个好消息!但是,也有将近7,000张图像可以应用。编写脚本来压缩它们将很简单: for old_img in [list of images]: convert_to_1bit_and_compress(old_img) remove(old_img) replace_with_new(old_img, new_img) ...但是它缺少一项至关重要的测试:新压缩的版本内容相同吗? if raster_diff(old_img, new_img) == "Identical": remove(old_img) rename(new_img, old_img) 有没有一种工具或方法可以自动(取消)证明Image-A的内容与Image-B的内容具有相同的价值? 我可以使用ArcGIS 10.2和QGIS,但是除了可以避免手动检查所有这些图像以确保覆盖之前的正确性之外,其他所有功能都可以使用。错误地转换和覆盖确实具有超过开/关值的图像将是可怕的。其中大多数会花费数千美元来收集和生成。 更新:最大的违规者是32位浮点数,每边最大可达100,000px,因此约30GB未压缩。

1
为什么多个栅格合并的结果这么大?[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗?更新问题,使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 我尝试像这样合并14个geotiff: 每个geotiff约为50Mb。我在输出端需要一个Geotiff 我的工作流程: gdalbuildvrt -input_file_list list.txt test.vrt (我的列表中包含tif的名称) 然后 : gdal_translate -of Gtiff test.vrt test.tif Input file size is 79841, 59955 它可以工作,但是结果是13,3 Gb的geotiff!对于14个文件(每个50 Mb),我尝试了700 Mb(而不是13 Gb)的geotiff。 我知道gdal默认不会压缩,所以我尝试了以下命令: gdal_translate -of Gtiff -co COMPRESS=JPEG test.vrt test_compressed.tif 但是文件的“合并”对于JPEG压缩来说太大了: Input file size is 79841, 59955 0ERROR 1: JPEGPreEncode:Strip/tile too large for …


1
限制光栅包在R中使用的内存?
将raster软件包使用的内存限制R为5GB RAM 的最简单方法是什么? 理想情况下,我希望将超出此大小的所有栅格写入硬盘驱动器并从临时文件进行分析。 我正在使用Windows 7(64位)和R版本3.0.3。 这是我在具有16GB RAM的计算机上使用的一些代码,但是在我处理大型栅格时,它一直保持最大值: require(raster) memory.limit(size = 5000) rasterOptions(format="GTiff",overwrite=TRUE,datatype="INT1S", tmpdir="C:/Research/BIN",tmptime=1.1,progress="text",chunksize=1000, maxmemory=1000) rasterTmpFile("delete_me_") r <- raster("myraster.tif") r[r==0] <- NA 正是在脚本的此阶段,RAM使用量增加到全部16GB,然后停止工作并取消了该功能。 我需要在这里更改什么?
10 raster  r  memory 

2
USGS DEM中水平和垂直条带的来源是什么?
当处理通过National Map Viewer从National Elevation Dataset下载的30m和10m DEM数据时,我们注意到水平和垂直条纹不仅出现在生成的结果中,而且还只是原始DEM的分析性山影。有人知道来源吗?如果不是来源,也许如何去除这些工件?当使用DEM计算地形指数时,这些伪影变得非常明显。即使在发生凹陷填充后,这些伪影仍然保留。 下图显示了宾夕法尼亚州和科罗拉多州分水岭的30m和10m数据的条带化,以及完成的地形指数计算,显示了纽约州锡拉丘兹的分水岭的文物。 科罗拉多州-HUC8-10190004-10m 科罗拉多州-HUC8-10190004-30m 宾夕法尼亚州-HUC8-02040103-10m 宾夕法尼亚州-HUC8-02040103-30m 纽约州锡拉库扎的奥农达加溪流域的TI计算完成

2
R中的镶嵌栅格?
我正在尝试将多个栅格镶嵌到R中的单个大栅格中。使用在/programming/15287807/how-can-i-create-raster-mosaic-using-list-of-rasters中发布的脚本 但是,我收到了警告消息和错误消息。 rasters1 <- list.files("F:\\MOD15A2_LAI_1km\\MOD15A2_LAI_2009", pattern = "mod15a2.a2009001.*.005.*.img$", full.names = TRUE, recursive = TRUE) mos1 <-mosaic(rasters1, fun=mean) 它报告错误如下 Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ‘extent’ for signature ‘"character"’ 然后我尝试了另一个版本。 rasters1.mosaicargs <- rasters1 rasters1.mosaicargs$fun <- mean 但是这里有一些警告信息如下 Warning message: In rasters1.mosaicargs$fun <- …
10 raster  r  mosaic 

1
如何在绘制多边形时撤消误点击?
我试图在QGIS 2.2中绘制一些多边形,但我却迷失了,因为每次单击不正确时,我都无法找到一种方法来撤消该点并重新开始绘制。 我敢肯定有办法。有人知道吗?
10 qgis  editing 

1
是否有任何类似Google Maps的可缩放世界地图,但具有相等面积的投影?
每当我在全球或极地查看Google Maps时,投影都会很麻烦。是否存在具有可缩放世界地图的Web服务,该服务使用等面积投影-最好是伪矩形(例如Mollweide投影)?甚至《国家地理》在线地图也使用墨卡托。 另一种选择是Google Earth,它并没有特别扭曲地球,但是也没有给出全局视图。 另请参阅:为什么Web Mercator(辅助领域)成为Web地图标准? 编辑2018-08-24: Google Maps不再使用Web Mercator,但已移至看起来是正交投影的位置。

3
元数据和GeoJSON?
我的代理机构目前使用FGDC XML元数据维护和分发shapefile。 现在,我们已经将许多形状转换为GeoJSON,但担心如果我们的开发人员忘记了从何处获取数据,那么当节点应用程序中断时,他们将不知道应该归咎于谁。 是否有GeoJSON元数据规范?

4
从未分类的LAS文件确定裸露DEM?
我有LAS格式的数据,其中的RGB值是使用无人机从航空摄影测量中创建的。我正在尝试找到一种从点云中提取地球DEM的解决方案。 我已经尝试过SAGA,Fusion和MCC-LIDAR,但似乎他们需要对LAS文件进行分类(自然而然就没有)。有人可以向我指出正确的方向,并简要说明该过程吗? 通常,我一次需要处理大约100个铣削点(如果需要,可以平铺它们)。

4
开放数据门户解决方案的比较
这个问题的答案是社区的努力。编辑现有答案以改善此职位。它目前不接受新的答案或互动。 对于许多组织,例如地方/县/州政府,提供对GIS数据的公共访问是一项挑战。 仅知道什么数据可用,什么数据,谁拥有和维护以及从何处获取数据可能是最紧迫的挑战。因此,GIS数据清单是必需的。当然,创建一件事是一回事,让人们对其进行维护(并对数据负责)是另一回事。因此,从内部利益相关者那里获得支持也很重要。 另一个重要的考虑因素是解决方案对元数据的支持。如果无法描述数据,则几乎没有用。对于某些解决方案,可能会以有限数量的专有字段的形式对元数据提供有限的支持。完整和原始格式的元数据支持似乎很少。 通用GIS数据(例如地址点或地块数据)的大小,记录/属性的数量,字段长度或几何复杂性可能受到限制,这可能会使某些解决方案无法使用。通常,某些类型数据的绝对大小(例如LiDAR,航空影像)使得很难或不可能以原始格式在互联网上托管它们,而必须通过物理介质(例如硬盘)进行传递。充其量,这些大文件可能托管在FTP服务器或Dropbox / Box之类的东西上,但是这样做可能会非常昂贵。 另外,在许多实现中,数据的空间性质可能不被认为是重要的,或者根本没有被考虑,这使得在空间上寻找,查看和使用空间数据变得困难或不可能。诸如以任何有用的方式符号化数据之类的事情可能都不存在。 最后,保持数据最新至关重要。如果没有自动的方法来更新数据,那么对于不重要的数据量,整个事情就变得不可行了。 因此,主要的问题是: 有哪些可用的开放数据门户解决方案,它们如何进行比较? 附带问题: 在决定开放数据门户解决方案时,最重要的考虑因素是什么? 是否存在使用多个门户网站有益(且具有成本效益)的方案? 是否可以将空间和非空间数据成功集成到同一门户中? 是否有开放式数据门户网站实现的示例,您会认为它们是成功的或具有示范性? 网上列出的美国政府免费GIS数据清单以及Wikipedia有关地理门户的文章中列出了一些实施方案。 附录: 我们目前正在评估的门户解决方案列表: ArcGIS Online ArcGIS开放数据 Portal for ArcGIS Server(本地) 苏格拉底和蒙达拉 Google Maps Engine和图库 我们尚未考虑的其他门户网站: 数据中心 卡恩 DKAN 开放式门户 门户网站不是真正的,但是可以很好地用于提供平面文件: 投寄箱 框 FTP(本地) 使用托管的存储/服务器或本地滚动自己的 在我们的案例中,由于我们拥有ESRI ELA并具有基于ArcGIS的空间数据基础结构,因此很可能会选择至少一个ArcGIS解决方案。 就个人而言,我宁愿看到我们使用更多的开源软件,但似乎不太可能将我们的产品整合在一起,以与ESRI的产品竞争,至少在没有改变人们使用我们组织的GIS方式的情况下。 但是,我希望这个问题对任何希望改善对公共数据的访问的人来说都是有用的资源。

2
比例和Z因子对QGIS中的山坡阴影分析没有影响
我使用SRTM DEM在QGIS 2.2.0中创建了山体阴影。但是,由此产生的山影似乎在垂直方向上被夸大了。同样,更改比例和z因子的设置对所得的山影没有任何影响。我已经尝试了QGIS,Analysis-> DEM(地形模型)和Terrain analysis-> Hillshade中的两个选项。两种方法均得出相同的结果。 这是孤立的事件,错误还是我缺少某些东西?
10 qgis  scale  hillshade 

1
通过shapefile中的前两位数字自动汇总德国邮政编码
我有一个shapefile,其中包含所有以五位格式(例如80637)的德语邮政编码。对于我打算做的事情,这是很多细节的方法。我想按前两位数字(例如80)来汇总所有邮政编码区域。 我的方法是对QGIS属性表中的邮政编码进行数字排序,然后手动选择它们,然后将它们加入。这有点麻烦,我不得不在这段时间里做些麻烦。 有没有办法以编程方式执行此操作?我按照此说明尝试自动执行此过程,但是导入六分仪失败。

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.