Spatial

茫茫宇宙，寻觅星辰

GIS 四色图与五色图制图

乡镇绘图四川省、市、县绘图时，由于数量少，GIS自带的色带渲染基本能解决这个问题。但绘制乡镇时，全省乡镇数量较多，想通过四色图来绘制乡镇的颜色。 QGIS 四色图 Topological coloring 将数据转换为投影坐标系，方便设置不同颜色的距离。然后打开工具箱，找到对应的工具，输入数据，最小值选4，其他可以选择默认参数，也可以自定义。注意：最小值是4，如果一个多边形周边的邻...

Posted by WT on January 18, 2024

使用R绘制Chord弦图

需求想绘制一个不同土地覆盖类型直接相互转换的图，找到了桑吉图和弦图，最后确定用弦图。搜了很久，发现Python没有R绘制的好。步骤按照参考文献1 写了代码，最后一步出图时显示不出来.于是找到了参考文献2，二者结合，把上图绘制出来了。后续更详细和复杂的绘制可以参考参考文献中的其他文献。 > install.packages('sankeywheel') > ins...

Posted by WT on January 11, 2024

倾斜摄影DSM和不同来源的DEM差值比较

问题倾斜摄影获得的三维模型和水动力模型所使用的DEM来源不同，导致UE在演示水流演进的时候，部分区域水流呈现倾斜状。数据处理步骤 1、通过DasViewer打开倾斜摄影模型，导出DOM时勾选DSM，同时获得DOM和DSM 2、根据DOM的范围获得倾斜模型的矢量范围，用此范围去裁剪水动力模型使用的DEM，然后将裁剪后的DEM重采样到倾斜摄影模型输出的DSM的分辨率，然后使...

Posted by WT on January 4, 2024

Python 中的Rasterio与Postgresql冲突

问题 Anaconda中的Rasterio使用过程中无法获取遥感影像的crs。具体错误如下： CRSError: The EPSG code is unknown. PROJ: proj_create_from_database: SQLite error on SELECT name, coordinate_system_auth_name, coordinate_system_cod...

Posted by WT on November 10, 2023

Open Data Cube 说明

本文记录 2023年Open Data Cube 安装过程及说明 datacube 1.8.17.dev2+g0cfe300d.d20231027 1 Install the Open Data Cube Software ODC 核心安装 ODC8.16 以后的版本要求python 3.10 建议直接安装Ubuntu 22，它自带的版本是Pyhon 3.10.12 按照Ubunt...

Posted by WT on October 23, 2023

统计下载的卫星影像的时相和景数

从ESA的网站下载了Sentinel 2的数据，在撰写报告时会有统计不同时相的卫星影像景数的问题，Python结合Excel可以解决这个问题步骤： 1 利用Python批量获取文件名，保存至csv文件 import os import csv filelist=[] filelist1=[] def file_name(file_dir): for root, dirs, ...

Posted by WT on September 13, 2023

利用CNBH-10m数据赋值建筑高度

思路复旦大学生命科学学院GC3S团队生产发布了中国首套首套10米分辨率的建筑高度栅格数据（CNBH-10m），前一篇文章介绍了利用QGIS、ArcGIS与ENVI提取建筑数据，也介绍了建筑白膜数据赋高程值,二者结合，就可以实现建筑白膜赋高度数据。问题与DEM不同，建筑高度数据只有在建筑存在的地方有值，其他位置是NoData，如下图所示：绿色的框中的建筑多边形的中心点处建筑高度数据为空...

Posted by WT on August 28, 2023

利用QGIS、ArcGIS与ENVI提取建筑数据

思路利用QGIS及其插件下载天地图带建筑轮廓的电子地图，使用ENVI识别这些数据并导出为SHP格式，使用ArcGIS将锯齿状的数据修正为正常的建筑物边界数据。申请天地图Key 登录天地图的网站,点击登录按钮登录天地图账号，如果没有账号可以点击注册，注册一个账号，登录账号,点击账号–>开发资源–>控制台—创建新应用，弹出对话框，填写便于识别的应用名称，从下拉框中选择行业类别及...

Posted by WT on August 8, 2023

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