天地图下载地表覆盖数据
选中所在区域网格,弹窗点击下载数据
zip包解压如下图:
行政区域数据下载
天地图数据中心提供了行政区划下载服务(阿里提供的需要坐标转换),大赞👍👍👍
打开QGIS加载矢量数据和栅格数据![在这里插入图片描述](https://img-blog.csdnimg.cn/direct/e1c75b8369e14ea3aa74d550430d3145.png)
栅格数据裁剪
工具箱 > 栅格 > 提取 > 按掩膜图层裁剪栅格
设置目标坐标系为EPSG:3857(其他投影坐标系也可,主要方便后续过滤面积小的区域)
裁剪完成
查看波段值
点击识别要素获取不同颜色对应的波段值如下图:
或查看下载数据压缩包中的xls文件如下图:
波段值映射表
代码 | 中文 | 波段值 |
---|---|---|
No Value | 无值 | 0 |
Cultivated Land | 耕地 | 10 |
Forest | 林地 | 20 |
Grass Land | 草地 | 30 |
Shrubland | 灌木地 | 40 |
Wetland | 湿地 | 50 |
Water Body | 水体 | 60 |
Tundra | 苔原 | 70 |
Artificial Surfaces | 人造地表 | 80 |
Bareland | 裸地 | 90 |
Permanent Snow and Ice | 冰川和永久积雪 | 100 |
Sea | 海水 | 255 |
栅格矢量化
工具箱 > 栅格 > 转换> 栅格矢量化(栅格转矢量)
结果:
矢量化图层右键 > 属性 > 符号化 > 分类 > 值选择 > DN > 点击分类 确认完成
随机颜色渲染的矢量数据
提取波段值为20的面(森林矢量面)
工具箱 > 矢量选择集 > 按属性提取
选择集字段选择 “DN” > 操作符选 “=” > 值选择 “20”
提取结果
截止目前数据提取部分已经完成。
当前提取的数据为真实的数据,但是过于零碎,根据业务需求考虑是否继续优化。
以后步骤是删除掉部分面积较小的面,完成大面中空并补洞操作:
计算面的面积
打开属性表 > 字段计算器 > area(@geomtry) > 点确认
删除面积较小的部分(当前设置小于100000)
点击编辑删除选中部分即可。
处理完成叠加影像图看效果
放大效果
可以看出面中间有部分中空,边缘锯齿明显;
矢量平滑(去锯齿)
工具箱搜平滑工具,根据实际需求设置参数,结果如图;![在这里插入图片描述](https://img-blog.csdnimg.cn/direct/6f9e79e8e8f84bd6a876abf635f42946.png)
去中空
修正几何图形
工具箱搜索修正几何图形
行政区面提取范围面
矢量 > 研究工具 > 提取图层范围
相交取反
矢量 > 地理处理工具 > 交集取反
输入图层:范围面
叠加图层:修正几何图层
转换单部件
矢量 > 几何图层工具 > 多部件转单部件
删除单部件中的多余的大面。