案例!从天地图中提取全市的建筑物矢量轮廓-以苏州市为例

简介: 案例!从天地图中提取全市的建筑物矢量轮廓-以苏州市为例

制作方法

步骤一:下载

先从 QGIS 中下载了苏州市部分主城区的天地图图像,参考系为3857,空间分辨率为0.5米,共1.6G。

 

左:下载的天地图矢量图像;右边:放大后

 

(虽然叫天地图矢量图像,但本质是服务器端渲染使用矢量数据制作的切片图片,在客服端显示而已,并不是矢量。)

 

步骤二:灰度图化

将下载的图像添加到任意 GIS 软件中,比如 ArcGIS Pro 或 ArcMap,我们需要将其灰度图化,我也喜欢叫二值化。

在 GIS 软件的符号系统中选择拉伸,拉伸色带使用黑白色带,然后导出,导出设置选择渲染器导出,他就会根据当前栅格的渲染结果进行导出。

 

 

导出结果如下

 

密密麻麻的白点全是建筑物

 

步骤三:分类并导出

然后根据灰度图设置合适的阈值将建筑物和非建筑物分开,然后导出栅格(勾选渲染器)。

 

ArcGIS Pro 符号系统

 

再次提醒,这不是一篇详细的教程文章,因为之前已经写过了,所以步骤不会非常详细。如果卡在某些地方,可以翻阅一开头提到的教程,非常详细。

 

步骤四:栅格转面

将分类导出的栅格数据转为面矢量数据。

需要注意的是,根据工具中是否勾选简化面,结果会有差别。

 

 

下面左边没有勾选简化面,右边勾选了简化面,有较为明显的差别。

另外最后导出的矢量效果和下载栅格数据的空间分辨率也一定的关系,注意事项中会详细说明。

 

 

 

注意事项

 

栅格空间分辨率

一开始下载的栅格数据的空间分辨率对最后的结果有非常大的影响,下面左边是下载的1米空间分辨率的天地图矢量瓦片,右边是0.5米,仅从图片上就能明显看到清晰度的差异。

 

 

在处理后,得到的最终结果也有非常大的差别,左边是从1米空间分辨率的天地图中提取的,右边则是0.5米。可以明显看到左边的锯齿更多,不光滑。

 

 

所以提高下载栅格的空间分辨率能有效把关最后的提取轮廓矢量质量。

当然一味的高空间分辨率也是不可取的,因为原始栅格影像的切片级数已经到头了,再高也不会有清晰度的提升,只会增加数据存储大小。

 

简化面、平滑面

在最后从栅格转面的时候,勾选简化面能有效平滑矢量轮廓的锯齿,如果还是没有达到期望的话,可以尝试以下工具:

  • 概化;
  • 简化建筑物;
  • 规则化建筑物覆盖区;
  • ……

 

这些处理可能会消耗大量计算机性能,毕竟图斑、折点数都很多。

 

(听到笔记本的风扇声,我怕它起飞后找不回来了 😢)

 

去除道路

放大看效果,道路也被提取出来,如何剔除道路呢?

 

 

通常道路的值都是统一的,毕竟颜色都一致,或者就两三种颜色,通常与建筑物的颜色不一样,我们可以在确定道路的值后,使用栅格计算器将该值得像元抹去,或者设置为其他会被排除的值。

更具体的处理步骤是在二值化后,整个栅格数据的像元值分布在0-255,确定道路的值后,使用栅格计算器精确打击!

 

完毕!建议收藏!

抛砖引玉,欢迎大家后台提出建议!

目录
相关文章
|
3月前
|
定位技术
GEE(CCDC-3)——根据CCDC segment分割后的影像进行地类变化统计和绘制土地覆被变化地图
GEE(CCDC-3)——根据CCDC segment分割后的影像进行地类变化统计和绘制土地覆被变化地图
47 0
|
7月前
|
定位技术
ArcGIS地形起伏度+地形粗糙度+地表切割深度+高程变异系数提取
ArcGIS地形起伏度+地形粗糙度+地表切割深度+高程变异系数提取
2137 0
|
存储 前端开发 算法
激光SLAM:ALOAM---后端lasermapping地图栅格化处理与提取
不同于前端的scan-to-scan的过程,ALOAM的后端是scan-to-map的算法,具体来说就是把当前帧和地图进行匹配,得到更准确的位姿同时也可以构建更好的地图.由于是scan-to-map的算法,因此计算量会明显高于scan-to-scan的前端,所以后端通常处于一个低频的运行频率,但是由于scan-to-map的精度往往优于scan-to-scan.因此后端也有比前端更高的精度.为了提高后端的处理速度,所以要进行地图的栅格化处理
激光SLAM:ALOAM---后端lasermapping地图栅格化处理与提取
|
7月前
|
数据采集 数据可视化
国内77个城市建筑物轮廓(带高度)数据分享(附百度网盘)
国内77个城市建筑物轮廓(带高度)数据分享(附百度网盘)
550 1
|
5月前
|
定位技术
基于ENVI实现栅格遥感影像按图层行列号与像元数量划定矩形研究区域并裁剪
基于ENVI实现栅格遥感影像按图层行列号与像元数量划定矩形研究区域并裁剪
|
5月前
|
数据采集 定位技术
R语言geodetector包基于栅格图像实现地理探测器操作
R语言geodetector包基于栅格图像实现地理探测器操作
|
5月前
|
编解码 定位技术
ENVI对不含地理参考信息的栅格图像添加地理或投影坐标系信息
ENVI对不含地理参考信息的栅格图像添加地理或投影坐标系信息
|
5月前
|
定位技术
ENVI自动批量生成地面控制点实现栅格遥感影像自动地理配准
ENVI自动批量生成地面控制点实现栅格遥感影像自动地理配准
|
10月前
|
数据可视化 定位技术
GIS空间分析 数字地形分析1 地势图的制作
在本文中,你将学到如何根据DEM数据制作地势图
96 0
|
10月前
|
编解码 定位技术
GIS空间分析 数字地形分析2 基本地形因子的提取
在本文中,你将学习到如何利用DEM数据在ArcGIS中提取出坡度、坡向等地形因子。
183 0