基于复杂阴影测量方法的建筑物高度反演

doi:10.6046/zrzyyg.2023121

基于复杂阴影测量方法的建筑物高度反演

李志新^,¹, 纪松^,¹, 范大昭¹, 高定¹, 李永建¹, 王刃²

1.信息工程大学地理空间信息学院,郑州 450001

2.山东五征集团,日照 276800

Building height inversion based on the complex shadow measurement method

LI Zhixin^,¹, JI Song^,¹, FAN Dazhao¹, GAO Ding¹, LI Yongjian¹, WANG Ren²

1. Institute of Geospatial Information, Information Engineering University, Zhengzhou 450001, China

2. Shandong Wuzheng Group Co.,Ltd., Rizhao 276800, China

通讯作者: 纪松(1983-),男,博士,副教授,博士生导师,主要从事航天摄影测量以及遥感影像处理与分析方面的研究。Email:jisong_chxy@163.com。

责任编辑: 陈昊旻

收稿日期: 2023-05-9 修回日期: 2023-07-20

基金资助:

国家自然科学基金项目“高分辨率卫星影像协同处理与定位能力星间传递技术研究”(41971427)
嵩山实验室项目“空间智能融合感知技术及精密重建系统”(221100211000-5)
国防科工局项目“高分遥感测绘应用示范系统(二期)”(42-Y30B04-9001-19/21)

Received: 2023-05-9 Revised: 2023-07-20

作者简介 About authors

李志新(1997-),男,硕士研究生,主要从事数字摄影测量研究。Email: 1294622314@qq.com。

摘要

阴影法能以较低的成本测量大范围的城市建筑物高度,但其中的阴影长度测量方法在复杂阴影场景中存在测量效率低、精度差、鲁棒性低等问题。对此,该文提出了一种复杂建筑阴影测量方法。该方法首先通过渔网法和多种约束条件相结合实现对阴影的测量与划分; 其次,统计所有划分区域的阴影长度值,利用四分位法和双向逼近策略确定阴影最优值; 然后,综合评估所有区域最优值确定阴影长度; 最后,利用立体交会的建筑高度对应用该阴影测量方法计算的建筑物高度进行验证。结果表明: 该方法计算建筑物高度有90.6%以上绝对误差在0~5 m之间。因此可以得出该文方法计算的阴影长度有着较高的精度,能够满足不同种类的复杂阴影测量,提高了阴影法反演建筑物高度的精度,为城市建筑物高度的反演提供新的技术支撑。

关键词： 渔网法; 分区统计; 四分位法剔除粗差; 双向逼近策略; 最优值评估

Abstract

Building heights are necessary for urban informatics, providing a significant basis for the planning and early warning of risks for urban construction. The shadow method, which can measure the heights of urban buildings on a large scale at a low cost, faces challenges such as low efficiency, accuracy, and robustness in building height inversion in complex shadow scenes. This study proposed a measurement method for these scenes. First, the shadows were measured and delineated using the fishing net method combined with multiple constraints. Second, the shadow lengths of all the zones divided were obtained, and the optimal values were determined using the quartile method and the bidirectional approximation strategy. Third, the shadow lengths were determined through a comprehensive assessment of the optimal values of all zones. The results show that 90.6% of building heights calculated using the new method exhibited absolute errors ranging from 0 to 5 m. Therefore, this method features elevated accuracy of building height inversion for various complex shadow scenes, laying a basis for research into the inversion and expansion of urban building heights.

Keywords： fishing net method; zonal statistics; quartile method to eliminate gross error; bidirectional approximation strategy; optimal value evaluation

PDF (11331KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

李志新, 纪松, 范大昭, 高定, 李永建, 王刃. 基于复杂阴影测量方法的建筑物高度反演[J]. 自然资源遥感, 2024, 36(3): 108-116 doi:10.6046/zrzyyg.2023121

LI Zhixin, JI Song, FAN Dazhao, GAO Ding, LI Yongjian, WANG Ren. Building height inversion based on the complex shadow measurement method[J]. Remote Sensing for Land & Resources, 2024, 36(3): 108-116 doi:10.6046/zrzyyg.2023121

0 引言

近年来我国城市化进程不断推进,不仅体现在城市面积上的增长,也体现在建筑物高度的增长^[1]。建筑物高度增长一方面能尽量克服城市土地资源匮乏的瓶颈,另一方面为优化城市结构及城市功能作出贡献。在城市遥感研究领域,对于城市建筑物高度的提取也成为研究的重点。城市建筑物高度的估计与测量,已成为城市规划与扩张和城市灾害风险预警与评估的重要参数^[2-3],同时也为数字城市三维模型的建立提供了基础测绘资料。因此,需要一种经济有效的方法来确定大规模的建筑高度信息,以支持城市设计和管理。

近些年来,卫星图像、建筑足迹矢量数据和激光雷达点云等众多地理空间数据的开源^[4],为低成本方式生成大型三维城市模型创造机会。影像中的阴影对遥感解译工作带来许多负面影响^[5],但某些阴影可用来计算目标地物的高度^[6]。在一些无人机无法到达或立体像对缺乏的地区,通过建筑物阴影反演高度信息成了常用的技术手段。对于单张图像来说,无法直接获取目标地物的高度,但可通过目标地物阴影与太阳、卫星高度角与方位角之间的特定关系,计算地物的高度信息。自1989年以来,在航空摄影测量中,研究人员一直使用阴影信息来估算建筑物高度^[7]。Shettigara等^[8]利用SPOT全色图像上的阴影信息构建模型,获得建筑物高度; Wang等^[9]利用ZY-3图像建立了阴影长度与建筑物高度之间的几何关系模型,并结合阴影长度计算了建筑物高度,对城市建筑物进行了三维建模; Turker等^[10]利用建筑物阴影计算地震中倒塌的建筑物高度。这些研究虽只在局部简单场景下进行了初探,但为利用阴影信息计算建筑物高度提供了很好的理论支撑。之后,诸多学者就阴影提取与阴影测量工作展开了研究。Shao等^[11]提出了一种结合图像对象空间指数的方法来提高阴影提取的精度,并以IKONOS图像为例,利用阴影长度来估计建筑物高度; Liasis等^[12]利用卫星图像的频谱和空间分析信息实现了一种新的主动轮廓模型,优化了建筑物的阴影分割过程,提高了阴影提取的精度,并通过阴影长度估算了建筑物高度; 张豪^[13]通过改进U-Net提高了建筑物阴影的检测精度,且能较好区分建筑物阴影与水体。

优化阴影提取工作对建筑物高度的估算精度提升有限,获得准确的阴影长度才是建筑物高度反演的关键。目前,阴影测量算法相对较少,主要有面积与周长反算法^[14]、角点最近距离法^[15]和渔网法^[16]3种。面积与周长反算法和角点最近距离法仅局限于规则、无遮挡的阴影区域,对于形状不规则与受遮挡的建筑物阴影并不适用。渔网法结合像元大小、太阳方位线等条件,在阴影区域生成一系列等间隔的渔网线,通过多次测量减少误差,相对精度较高,适用范围较广。Xie等^[16]结合建筑物轮廓与渔网法,解决了建筑密集区域阴影相互粘连问题,实现了建筑物与阴影的一一对应,但未探究复杂场景下阴影边界不平滑与遮挡对测量结果的影响。

针对目前阴影测量算法在复杂阴影场景中存在测量效率低、精度差、鲁棒性低等问题,本文提出了一种复杂建筑物阴影测量方法,旨在通过阴影的最优区域计算阴影长度的最优测量值。测量过程充分考虑了对测量区域、统计数据,以及最优值的择优,能够最大限度地降低阴影边界不平滑与遮挡对测量结果的影响,提高阴影法计算建筑物高度的准确性与鲁棒性,从而实现大规模、大面积的城市建筑物高度测量,为城市建筑物高度的反演提供基础研究数据。

1 原理与方法

复杂建筑物阴影测量算法主要包含分区测量与最优值评估2个部分,如图1所示。首先,采用渔网法^[17]对阴影矢量进行布点测量; 其次,利用最邻近搜索算法确定每个建筑物矢量与之对应的阴影矢量,并依据建筑物四至范围和采样点间距确定分割点; 然后,结合太阳方位分割所对应的阴影矢量形成多个统计区域; 进而,统计每个区域内的阴影长度值,运用四分位法剔除粗差,并通过双向逼近策略确定每个区域的最优值; 最后,评估每个区域的最优值,确定该建筑物的阴影长度。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 建筑物阴影测量技术路线

Fig.1 Building shadow measurement technique route

1.1 分区测量

现有基于阴影信息反演建筑物高度的方法通常是建立在阴影规则化、地形平坦化和分布理想化3个假设之上^[18⇓-20],理想化假设也导致多数算法在实际应用中鲁棒性较低,究其原因在于实际阴影情况复杂,建筑物相连、遮挡严重,使得阴影测量结果不准确,如图2所示。

图2-1

新窗口打开| 下载原图ZIP| 生成PPT

图2-1 建筑物阴影场景

Fig.2-1 Building shadow scene

图2-2

新窗口打开| 下载原图ZIP| 生成PPT

图2-2 建筑物阴影场景

Fig.2-2 Building shadow scene

图2分别展示了真实场景中10种常见的阴影形态: 其中图2(a)—(b)展示了较为理想的阴影分布情况; 图2(c)—(e)展示阴影边界受植被影响,造成阴影内部孔洞、阴影边界不平滑、阴影遮挡与延长等问题; 图2(f)—(j)展示阴影受建筑物不同程度遮挡,以及阴影相互粘连等问题。可以看出,对于受遮挡的建筑物阴影,只有部分区域能够用来计算建筑物高度,故称这部分区域为最优测量区域,用其计算得到的值为最优值。因此,需重点考虑阴影的最优测量区域,以此获得准确而可靠的阴影长度值。

本文以建筑物多边形为附加条件,设计了一种分区测量方式,如图3所示,解决建筑物密集区域阴影横连、纵连和阴影受遮挡导致测量不准确的问题。首先,采用渔网法对整个阴影掩模进行测量,原理为: 以一定的采样间隔(以像素为单位)在阴影范围内布设种子点,如图3(a)所示,并测量阴影边界在太阳方位上的距离; 其次,利用最邻近搜索算法确定每个建筑物与之对应的阴影; 接着,依据建筑物外包矩形的四至范围和采样点间距确定分割比例点。对于分割比例点的确定,以图3(b)1号统计区域为例,若分割比例太大,则统计区域过大,将有更多受遮挡影响的测量值计入集合当中,造成测量结果偏低; 若分割比例太小,将仅有少量或没有测量值计入集合,导致计算结果不可靠。因此,为了保证在水平方向上有一定数量的测量结果计入统计集合中,本文以水平测量间隔加一定的缓冲距离所占建筑物外包矩形宽度的比值作为最小分割比例,公式为:

(1)

s e g = \frac{d \cdot S R + η}{{X^{B}}_{m a x} - {X^{B}}_{m i n}}

式中: seg为分割比例; d为采样间隔; $η$ 为缓冲距离; ${X^{B}}_{m a x}$ 与 ${X^{B}}_{m i n}$ 分别为建筑物外包矩形的最大经度和最小经度; SR为影像空间分辨率。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 分区测量示意图

Fig.3 Partition measurement

为了降低遮挡测量值对最终统计结果的影响,共定义了4个测量区域: 以最小分割比例确定左右测量区域(分割点div1和div3的经纬度分别为 $d i v 1_X$ , $d i v 1_Y$ , $d i v 3_X$ , $d i v 3_Y$ ),来降低测量区域内遮挡测量值的比值,提高遮挡情形的测量能力。将中间测量区域一分为二(分割点div2经纬度为 $d i v 2_X$ 和 $d i v 2_Y$ ),可有效避免因中间测量区域过大或更多受遮挡影响的测量值计入集合当中,造成测量结果偏低。以此为依据划分的测量区域,可有效兼顾不同程度的遮挡。分割点的经纬度坐标计算公式分别为:

(2)

d i v 1_X = {X^{B}}_{m i n} + ({X^{B}}_{m a x} - {X^{B}}_{m i n}) \cdot s e g

(3)

d i v 2_X = ({X^{B}}_{m a x} + {X^{B}}_{m i n}) / 2

(4)

d i v 3_X = {X^{B}}_{m a x} - ({X^{B}}_{m a x} - {X^{B}}_{m i n}) \cdot s e g

(5)

d i v 1_Y = d i v 2_Y = d i v 3_Y = ({Y^{B}}_{m a x} + {Y^{B}}_{m i n}) / 2

式中: ${Y^{B}}_{m a x}$ 和 ${Y^{B}}_{m i n}$ 分别为建筑物外包矩形的最大纬度和最小纬度。

然后,通过建筑物四至坐标、分割比例点坐标确定分割起始点,并依据太阳方位、分割长度(阴影长度加上建筑物长度的一半)形成分割区域,分割区域与阴影相交的部分为测量区域。

定义 $P o l y ⌈·⌉$ 为一个封闭的多边形,4个测量区域计算公式分别为:

(6)

H_{b u f} = ({Y^{S}}_{m a x} - {Y^{S}}_{m i n}) + ({Y^{B}}_{m a x} - {Y^{B}}_{m i n}) / 2

(7)

S 1 = P o l y_{shadow} \cap P o l y ⌈X^{B}_{min}, d i v 1_X, d i v 1_Y, γ, H_{buf}⌉

(8)

S 2 = P o l y_{shadow} \cap P o l y ⌈d i v 1_X, d i v 1_Y, d i v 2_X, d i v 2_Y, γ, H_{buf}⌉

(9)

S 3 = P o l y_{shadow} \cap P o l y ⌈d i v 2_X, d i v 2_Y, d i v 3_X, d i v 3_Y, γ, H_{buf}⌉

(10)

S 4 = P o l y_{shadow} \cap P o l y ⌈X^{B}_{max}, d i v 3_X, d i v 3_Y, γ, H_{buf}⌉

式中: $S i (i = 1,2, 3,4)$ 为测量统计区域; γ为太阳方位角; $H_{b u f}$ 为分割高度; ${Y^{S}}_{m a x}$ 和 ${Y^{S}}_{m i n}$ 分别为阴影外包矩形的最大与最小纬度, $P o l y ⌈X^{B}_{min}, d i v 1_X, d i v 1_Y, γ, H_{buf}⌉$ 为以 ${X^{B}}_{m i n}$ , $d i v 1_X$ , $d i v 1_Y$ ,γ, $H_{b u f}$ 所组成的分割区域,其余同理; $P o l y_{s h a d o w}$ 为阴影多边形。

最后,统计每个区域内的测量值,并构建阴影长度集合 $L_{i} = {l_{1}, l_{2}, \dots, l_{n}} (i = 1,2, 3,4)$ ,如图3(c)所示。

1.2 最优值评估

本文采用统计学中四分位法^[21]剔除L_i中的异常值,四分位法可以分析数据的离散程度,或者两极分化程度,通过对数据分割,剔除极大极小异常值的影响,并通过判断数据密集度,判断数据主要集中的区间。为了进一步得到较为可靠的测量值,采用双向逼近策略,通过设置择优阈值T,选出一组高度相近的测量值,取其平均值作为该统计区域的最优值,如图4所示。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 确定最优值

Fig.4 Determine the optimal value

具体步骤为: 将L_i数值从小到大排列,四分位法剔除粗差; 取L_i的中位数作为参考值l_med,设置择优阈值T; 计算L_i中l_max与l_min之差,若大于T,则进一步剔除,若 $l_{m a x} - l_{m e d} > l_{m e d} - l_{m i n}$ ,则移除 $l_{m a x}$ ; 若 $l_{m a x} - l_{m e d} < l_{m e d} - l_{m i n}$ ,则移除 $l_{m i n}$ 。

迭代步骤,直至L_i中l_max与l_min之差小于T,取其均值作为该区域的最优值 $o p t$ ,公式为:

(11)

o p t = \frac{1}{n} \sum_{i = 1}^{n} l_{i}

式中n为经双向逼近策略筛选后测量值的数量。

通过分区测量的方式,可以有效解决现有阴影测量算法在阴影顶部边界受周边建筑物部分遮挡情形下测量精度差的问题,但是阴影顶部边界受周边建筑物全部遮挡情形,因图像无法呈现完整的阴影长度,故用其测量的建筑物高度值并不可靠。在实际应用中,需将阴影全遮挡的目标建筑物筛选出来,通过人工或结合其他辅助信息进一步分析。

本文以4个测量区域 $o p t$ 值的差异情况来对阴影场景进行分类,从图2的阴影场景中可以直观发现,若使用上述分区测量方法对每个类型的阴影场景进行测量,阴影无遮挡与全遮挡的4个区域边界相对平滑, $o p t$ 较为相近,而阴影受建筑物、树木等其他地物部分遮挡的4个区域的 $o p t$ 会有明显差异,且最大的统计结果最可信。故此,根据不同统计区域 $o p t$ 差异情况,将阴影场景分为2类: 阴影无遮挡与全遮挡为第一类场景; 阴影受部分遮挡为第二类场景。以此背景,本文设计了一种最优值评估策略,其技术路线如图5所示。图5中,opt_max $与$ opt_min分别为4次统计结果的最大值与最小值; T'为差异阈值。本文以建筑物高度允许误差对应的阴影长度作为差异阈值T'。公式为:

(12)

T' = \frac{δ}{t a n β}

式中: δ为建筑高度的允许误差; β为太阳高度角。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 最优值评估技术路线

Fig.5 Optimal value evaluation technique route

针对第二类场景,直接取 $o p t_{m a x}$ 所对应的测量区域作为最优区域, $o p t_{m a x}$ 也将代表该建筑物的阴影长度值。对于第一类场景,以中间的统计区域建立缓冲区,搜索是否有相交建筑物,若无,则取4个统计区的平均值,作为该建筑物的阴影长度值。若有,则进一步向四周扩大搜索范围,以最邻近的无遮挡建筑物对应的 $o p t$ 或第二类建筑物所对应的 $o p t_{m a x}$ 作为该建筑物的阴影长度值,并对其进行标记,为人工校正提供便利。

2 实验与分析

2.1 实验数据

GF-7是我国目前空间分辨率最高的民用立体测绘卫星,可有效获取20 km幅宽,空间分辨率优于0.8 m,本文将以GF-7后视全色影像为例,验证本文方法的可行性。实验数据位于中国河南省郑州市,影像中心经纬度为E113.8°,N34.6°,拍摄时间2021年9月20日上午11: 30,共从中选择了8个研究区,如图6所示,其中简单场景(建筑物阴影规则完整、无遮挡)4个,复杂场景4个,研究区影像尺寸在600~1 000像素之间。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 研究区可视化图

Fig.6 A visualization of the study area

以U-Net提取的阴影掩模作为测量对象,来验证本文方法在不同场景的有效性。U-Net作为最具代表性的全卷积神经网络(fully convolutional network,FCN),被广泛应用于目标分割及目标分类等领域,已有研究证明了U-Net在阴影提取方面具有出色的性能^[13,22-23]。使用1 450张GF-7后视全色图像作为阴影边界提取模型的数据集,手动对阴影进行标记,然后构建用于阴影提取的U-Net模型。在数据增强后,按7∶2∶1的比例被随机分为训练集、测试集和验证集。在训练过程中,采用SGD优化器对损失函数进行优化,训练100轮次,初始学习率为0.001,批大小为16。

2.2 参考数据与评价指标

建筑物阴影长度的参考值不易直接获得,可根据阴影长度、太阳与卫星参数、影像成像时间和经纬度等信息计算建筑物高度,通过分析建筑物高度的精度间接反映阴影测量算法的精度。GF-7后相机倾角为-5°,太阳、传感器方位角在0°~180°以内,传感器、太阳和建筑物之间的几何关系(为了方便计算方便,方位角的基准方向设置为东),如图7所示。图7中β为太阳高度角,γ为太阳方位角,α为传感器高度角,ε为传感器方位角,L为建筑物阴影的实际长度,L₂为在遥感图像上观测到的建筑物阴影长度,d为建筑物阴影受建筑物本体遮挡的长度。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 建筑物与阴影几何关系示意图

Fig.7 Schematic diagram of geometric relationship between building and shadow

太阳与传感器的高度角和方位角可从影像元数据中获得。建筑物高度的计算公式为:

(13)

H = \frac{L_{2} \cdot t a n β \cdot t a n α \cdot s i n γ}{t a n α \cdot s i n γ - t a n β \cdot s i n ε}

。

通过GF-7立体像对以人工刺点的方式交会出每栋建筑物的高度,每栋建筑物保证3次以上测量,取平均值作为参考高度。本文选取绝对误差、平均绝对误差(mean absolute error,MAE)、平均相对误差(mean relative error,MRE)和均方根误差(root mean square error,RMSE)为评价指标,计算公式分别为:

(14)

M A E = \frac{1}{m} \sum_{i = 1}^{m} |H_{i} - H_{i}^{0}|

(15)

M R E = \frac{1}{m} \sum_{i = 1}^{m} \frac{|H_{i} - H_{i}^{0}|}{H_{i}^{0}} \times 100 %

(16)

R M S E = \sqrt{\frac{1}{m} \sum_{i = 1}^{m} (H_{i} - H_{i}^{0})^{2}}

式中 $H_{i}^{0}$ 为参考高度值; m为测量次数。

2.3 建筑物阴影分区测量结果

使用本文的复杂阴影测量方法对图2所示的10个阴影场景进行了测量。在本文实验中,以10个像素作为采样间隔进行布点测量,并结合建筑物四至确定测量区域,不同场景的分区情况如图8所示。统计每个区域内的测量值,设置择优阈值T为3 m,对择优后阴影长度取平均值作为该区域的最优值。设置建筑高度的允许误差δ为5 m,进行场景判断,选出最能代表该阴影长度的最优值。从图8可以看出,本文的阴影测量方法可以自动分成4个不同的测量区域,尽可能地减少了遮挡和不平滑边界对测量结果的影响,较为准确地测量了阴影的长度,能够适应不同类型阴影场景的测量。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 分区测量结果

Fig.8 Schematic diagram of partition measurement results

2.4 多场景建筑物阴影测量精度分析

采用不同的阴影测量算法与本文方法进行比较,文献[9]与文献[12]分别以整个测区的平均值和中值作为建筑物的阴影长度,文献[16]剔除了3倍中误差之外的测量值,并取剩余测量值的平均值作为建筑物的阴影长度,而本文方法相较于对比方法,进行了分区测量与最优值评估。在精度评估方面,对比不同方法计算的建筑物高度与参考值之间的差异,并统计不同场景类型的绝对误差和不同方法的MAE,MRE和RMSE,如表1和图9所示。从表1可知,本文的阴影测量算法精度最高,MAE为1.332 m,相较于文献[9]、文献[12]、文献[16]分别提升了6.162 m,4.343 m和5.804 m,MRE与RMSE也取得了最优的结果。结合图9可以看出,前5个阴影场景,不同测量算法的精度相当。对于阴影边界不平滑、受植被少量遮挡等场景(图8(c)—(e)),文献[9]与文献[16]会有不同程度的低估,文献[12]所使用的中值法受影响较小。对于阴影边界受建筑物部分遮挡等情形(图8(f)—(i)),对比方法因统计了大量不完整的阴影长度而失效,但本文方法仍能保持良好的稳定性与较高的精度。对于阴影边界受建筑物全部遮挡的场景(图8(j)),因无法测得完整的阴影长度,不予比较。

表1 测量误差统计

Tab.1 Statistic of measurement error

测量方法	MAE/m	MRE/%	RMSE/m
文献[9]	7.494	8.8	10.498
文献[12]	5.675	6.7	9.318
文献[16]	7.136	8.3	10.565
本文方法	1.332	1.8	1.657

新窗口打开| 下载CSV

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 不同阴影场景建筑物高度估算误差

Fig.9 Error of building height estimation in different shadow scenes

2.5 建筑物高度估算结果与分析

通过对图6所示的8个研究区进行建筑高度估算,进一步验证本文方法的有效性。8个研究区共139栋建筑物,高度分布均匀。简单场景共87栋,半遮挡情形2栋,全遮挡情形1栋。复杂场景共52栋,受半遮挡情形20栋,全遮挡情形10栋。建筑物高度误差统计结果如表2和图10所示。

表2 建筑物高度误差统计

Tab.2 Building height error statistics

研究区	MAE/m	MRE/%	RMSE/m
简单场景1	3.179	5.3	2.292
简单场景2	2.995	9.8	4.362
简单场景3	1.949	2.1	1.388
简单场景4	1.204	3.5	4.014
复杂场景1	3.455	3.8	4.589
复杂场景2	2.362	5.4	2.745
复杂场景3	4.512	6.3	5.526
复杂场景4	3.733	10.1	4.763

新窗口打开| 下载CSV

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 建筑物高度误差统计

Fig.10 Statistical chart of building height error

从图10和表2可以看出,使用本文阴影测量算法计算的建筑物高度值具有较高的精度,简单场景的MAE在1.204~3.179 m之间,MRE在2.1%~9.8%之间,RMSE在1.388~4.362 m之间,误差曲线相对稳定,其中简单场景1区域最大的误差值为全遮挡阴影场景对应的建筑物高度误差。复杂场景的MAE在2.362~4.512 m之间,MRE在3.8%~10.1%之间,RMSE在2.745~5.526 m之间。复杂场景2区域精度最高,复杂场景3区域精度较差,分析其原因为: 测区内有2栋全遮挡情形的建筑物,因与最邻近建筑物高度差异较大而未得到较好的补偿(复杂场景1区域2栋全遮挡情形的建筑物同理),最左侧3栋半遮挡情形的建筑物,因所在地势有明显的落差,可视阴影被延长,使得建筑物高度值被高估。复杂场景4区域内有6栋全遮挡情形的建筑物,有3栋搜索到与之高度相近的建筑物而得到有效补偿。图10(c)统计了所有研究区内除全遮挡情形外的建筑物高度误差。从结果可以看出,88栋建筑的绝对误差小于3 m,116栋建筑的绝对精度小于5 m。虽然,在建筑物高度反演过程中,包含了太阳、卫星参数等其他误差,但从结果来看,可以间接反映出本文阴影测量算法具有较高的精度。一般情况下,该方法提取的建筑物高度与立体交会的建筑物高度相比,有90.6%的建筑物高度绝对误差在0~5 m之间,可以满足城市规划的要求,具有一定的实用价值。

3 结论

本文提出了一种以建筑物多边形为附加约束条件的阴影分区测量方法,有效提升了遮挡情况下阴影测量鲁棒性。该方法通过渔网法和多种约束条件相结合实现对阴影的最优划分。然后统计所有划分区域的阴影长度值,利用四分位法和双向逼近策略确定阴影最优值。

使用本文阴影测量方法估算的建筑物高度与立体交会的建筑物高度相比,90.6%的建筑物高度绝对误差优于5 m,提高了阴影法反演建筑物高度的精度,可以满足城市规划的要求,具有一定的实用价值。

同时,本文方法也具有一定的局限性: 虽然有效优化了复杂场景中阴影长度提取的精度,但对于阴影全部遮挡情形仍无法有效解决。为了进一步提升阴影长度测量的鲁棒性,可以考虑利用神经网络学习该情形建筑物及其邻域特征来弥补阴影信息严重缺失的问题。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

Huang

, Chen

, Xu

, et al.

Estimating building height in China from ALOS AW3D30

[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2022, 185:146-157.