Spatial

茫茫宇宙，寻觅星辰

GIS常见操作

ArcGIS 多层标注使用下一行的代码可以实现图片中的效果 “"+ [水深] +"” + vbnewline +”"+ [水位] +"”+ vbnewline+ [单元编号] UE 安装 Python 环境：3.9.x 使用这种方法安装，不能直接使用Anaconda中的opencv库，否则会显示dll加载过程中有问题 pip install opencv-python

Posted by WT on March 7, 2024

长水面数据局部更新思路

问题不同时期的河流水面宽度不同，不同来源的影像纠正后位置偏离也不同，长水面矢量数据可能不会和底图完全匹配，但如果差别特别大，则需要更新。思路蒋总写代码cellclip.py前已经将数据区域划分为0.1°*0.1°的小方格，然后根据小方格内的三角形切片生成一个单独的文件，如果哪个网格内的数据出现问题，则重新勾画该网格及接边格区域内的一小块河道（保证新旧河道平滑过渡）即可。然后将这...

Posted by WT on January 25, 2024

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