Spatial

茫茫宇宙，寻觅星辰

Blender加载天地图及建筑白膜

思路 Blender中的插件BlenderGIS可以加载Google、ESRI、BING等数据，通过修改代码，增加对天地图的支持，然后通过获取OSM建筑数据和OpenTopography SRTM 30m的数据获取建筑高程，最后生成白膜数据。 Blender 下载建议下载版本为 2.8.3。可以直接从官网下载，点击历史版本，可以自由选择需要的版本。也可以通过分享的百度网盘下载。 Ble...

Posted by WT on July 25, 2023

建筑白膜数据赋高程值

思路建筑白膜数据为多边形面数据，赋高程值有两种方法：方法一：将面转点，根据点提取DEM数据中的高程数据，将提取的值赋给对应的多边形方法二：计算每个多面性区域的DEM数据的均值，将该值赋给对应的多边形。方法一 1、通过 Add Field增加建筑白膜数据的高程字段Z。 2、通过Feature To Points 工具将建筑白膜数据转换为点，选中 Inside选项，保证点在...

Posted by WT on May 9, 2023

MacOS下使用Anaconda和Pycharm配置R环境

安装前提 1、MacOS 10.15+(Intel) 2、Anaconda 3、PyCharm 2022.3 Anaconda 安装R环境 Pycharm下的R插件 R Language for IntelliJ 只支持R版本大于3.4的，Anaconda默认安装3.3版本的R，所以使用conda强制安装高版本的R。通过以下命令可以在Anaconda中创建r4-base的虚拟环境，并安...

Posted by WT on May 4, 2023

水动力模型堤防局部更新方法

二维模型分块建模 IFMS二维模型建模过程中，可能会遇到同一个模型被不同的人员进行修改，最后合并时无法准确的合并的问题。情景1：在公共模型的基础上，人员A仅更新了区域1的堤防数据（边元），人员B仅更新了区域2的网格高程数据（单元）。此时，以人员B的模型为基础，将A修改的边元数据更新至B的模型。 1、边元更新数据导出边元默认的高程值是-999,更新后的值>0，通过右键边元图层，选择...

Posted by WT on April 21, 2023

Mac结束进程的方法

Linux / MacOS 结束运行中和挂起状态下的进程 jobs # 列出所有作业 jobs -r # 列出运行中的作业 jobs -s # 列出挂起的作业执行以上任意一条命令，再执行下面的命令 kill %1 其中 %1 是从 jobs 命令中读取的作业数。还可以使用 jobs -l 列出 pid，并像往常一样使用 kill 命令。 DEMO 1、通过命令行 ./pg_ti...

Posted by WT on August 3, 2022

Mac内置Web服务器设置及Safari浏览器跨域

Mac内置Web服务器设置 Mac OS 系统自带了Apache服务器，启动和关闭的方法如下： $ sudo apachectl start // 启动。此默认情况下对应电脑路径为“LibraryWebServer/Document/” $ sudo apachectl stop // 关闭。 $ sudo apachectl restart // 重新启...

Posted by WT on June 23, 2022

ArcGIS分区域统计栅格数据

上一篇文章提到了Geohash网格生成，本文讲解如何基于ArcGIS利用Geohash 5级网格统计京津冀地区PM2.5的浓度。数据 1、京津冀地区PM2.5浓度，分辨率为0.01度（约1km） 2、京津冀Geohash 5 网格步骤 1 打开ArcMap的工具箱，找到 Spatial Analyst Tools–>Zonal–> Zonal Statistics...

Posted by WT on June 16, 2022

QGIS和Python生成中国区域Geohash网格

步骤 1 获取中国地图，利用QGIS绘制包含中国全境的简单多边形。 2 安装Python库 3 按需求生成Geohash网格 from polygeohasher import polygeohasher import geopandas as gpd import matplotlib.pyplot as plt %matplotlib inline gdf=gpd.read_fi...

Posted by WT on March 30, 2022

Cesium加载TMS切片

Cesium 可以加载来自ArcGIS、BingMap、GoogleMap、OSM等各种数据源，本文主要介绍加载TMS切片的功能。 TMS的切片有不同的变种。我们使用QGIS自带的GDAL2TILES生成TMS切片。步骤 1 生成切片打开QGIS的工具箱，目录GDAL–Raster miscellaneous–gdal2tiles 打开生成切片的工具后，加载要切图的TIF文件，T...

Posted by WT on December 31, 2021

DEM数据合并

步骤： 1 DEM数据下载在地理空间数据云可以免费下载30m的DEM数据，一共有三个版本，注意区别，建议使用GDEMV3 30M 分辨率（2019年发布的修正数据）。 2 ArcGIS 镶嵌合并DEM 点击toolbox-数据管理工具-栅格-栅格数据集-镶嵌至新格栅注意事项：输出的文件要有后缀名，否则其他软件无法识别，建议加.tif；像素类型要选择合适的范围，如果选8位，最大值为2...

Posted by WT on December 9, 2021