文章信息
- 段祝庚, 肖化顺
- Duan Zhugeng, Xiao Huashun
- 利用单幅影像及DEM数据进行森林火点定位及精度分析
- Forest Fire-Point Location and Precision Analysis Based on Single Image and DEM
- 林业科学, 2013, 49(4): 77-84
- Scientia Silvae Sinicae, 2013, 49(4): 77-84.
- DOI: 10.11707/j.1001-7488.20130411
-
文章历史
- 收稿日期:2012-06-05
- 修回日期:2012-10-30
-
作者相关文章
2. 中南大学地球科学与信息物理学院 长沙 410083
2. School of Geo-Sciences and Info-Physics, Central South University Changsha 410083
森林火灾是一种最具破坏性的森林灾害,直接危及立木、土壤、微生物和野生动物,降低森林的更新能力,引起土壤贫瘠,破坏森林涵养水源的功能,导致生态环境恶化,造成大气污染等(关百钧等,1992; 舒立福等,1997)。林火监测是森林火灾预防和灭火的重要环节和手段。随着现代科技的发展,视频监测已经成为一种重要的森林火灾监测技术和方法。火点位置是林火视频监测与森林火灾扑救必须获取的重要信息,它直接影响火灾报警的准确程度,林火扑救方案选择与决策,林火扑救行为的实施,扑火人员、扑火器材、林火扑救路线的选择等,将最终影响林火的燃烧面积和危害程度(林其钊等,2003)。
基于视频监测技术的火点位置确定是林火视频监测的关键技术。通过图像或视频监测技术进行森林火点定位的本质是从二维平面图像恢复真实世界的三维空间信息,主要的方法有基于双像或多像的立体重构、基于图像序列的三维重构(孙敏,2007)和基于单幅影像(刘亚文等,2005)的火点定位。前2种方法虽然精度高,但由于通视条件的影响,在森林地区难以获取相同目标的2幅影像和图像序列,且其算法较为复杂。而基于单幅影像的火点定位在地形复杂的森林地区易于实现,由于森林火点定位精度要求不高(本文森林火点定位精度要求为100 m以内),利用单幅影像进行火点定位不失为一种简单而快捷的方法。本文针对利用单幅影像和DEM数据进行森林火点定位的原理、数据模型、定位方法及其精度进行研究与分析,为利用视频监测技术进行森林火点监测火点位置确定提供一种新的方法和途径。
1 单幅影像火点定位原理由摄影测量解析理论可知,利用单幅影像不能由二维平面影像恢复真实世界的三维空间信息,由像点、物点及摄影中心所确定的共线方程只能确定像点、物点及摄影中心所决定的方向线,无法确定物点位置; 但将共线方程与摄影区域的数字高程模型DEM(digital elevation model)结合,即可确定物点位置,其本质是由共线方程所决定的像点、物点及摄影中心所在的直线与DEM所确定的曲面的交点。如图 1所示,S为安装在地势较高的林火监测云台上的视频设备,可以对较低区域的森林进行视频监测,获取视频影像数据。同时,利用固定安装摄像机上的三维电子罗盘获取摄影瞬间摄像机的三维姿态角方位角(azimuth angle)、倾角(pitch)、侧滚角(roll),并将获取的视频影像数据和摄像机的三维姿态角通过无线传输系统传回到森林防火监测指挥中心,则可由指挥中心林火监测系统的3D空间分析功能,对DEM数据作摄影视线方向的剖面图得到剖面线。这样火点定位问题可转化为摄影视线与山体表面的交线,即摄影视线方向与该方向DEM剖面线的交点。
为了建立火点定位的数学模型,建立如图 2所示坐标系统。
1)像平面坐标系统o-xy是以影像中心o为原点,x轴平行于影像像元水平排列方向,y轴平行于影像像元垂直排列方向。定位时需将图像像素坐标转换为像平面坐标o-xy。
2)像空间辅助直角坐标S-XYZ是以摄影中心S为原点、向东方向为X轴、垂直向上方向为Y轴、向南方向为Z轴的三维空间坐标系。
3)摄影测量坐标系O-XYZ是以森林区域西南角某点为原点、实地东方向为X轴、垂直向上方向为Y轴、实地南方向为Z轴的右手三维空间坐标系。摄影中心S在摄影测量坐标系O-XYX中的坐标为(XS,YS,ZS)。
2.2 数学模型视频或摄影属于中心投影,像点、物点及摄影中心满足共线方程(李德仁等,2001),如式(1)所示:
$\begin{array}{l} X = {X_S} + \left({Y - {Y_S}} \right)\frac{{{a_1}\left({x - {x_0}} \right)+ {a_2}\left({y - {y_0}} \right)- {a_3}f}}{{{b_1}\left({x - {x_0}} \right)+ {b_2}\left({y - {y_0}} \right)- {b_3}f}};\\ Z = {Z_S} + \left({Y - {Y_S}} \right)\frac{{{c_1}\left({x - {x_0}} \right)+ {c_2}\left({y - {y_0}} \right)- {c_3}f}}{{{d_1}\left({x - {x_0}} \right)+ {d_2}\left({y - {y_0}} \right)- {d_3}f}}。 \end{array}$ | (1) |
式中: x,y为像点在像平面坐标系o-xy中坐标; x0,y0,f为相机内方位元素; XS,YS,ZS为摄影中心在摄影测量坐标系O-XYZ中坐标; X,Y,Z为像点对应物点在摄影测量坐标系O-XYZ中坐标; ai,bi,ci(i=1,2,3)为以Y-X-Z为主轴的旋转角构成的方向余弦。在瞭望塔上安装数字云台,数字云台携带了三维电子罗盘,可利用固定安装在摄影机上的三维电子罗盘实时测定拍摄像片时的三维姿态角———方位角(φ)、倾角(ω)、侧滚角(κ)。
由式(1)可知,若知道Y 值,也就是物点高程,即可通过共线方程求解物点平面位置(X,Z)。但此时Y 值未知,不能直接利用共线方程求解物点平面位置,需要利用DEM数据作为辅助数据求解物点平面位置(X,Z)。由式(1)可导出式(2),即像点、物点位于过摄影中心(XS,YS,ZS)的铅垂面内,且X,Z满足式(2)。沿X轴方向每隔一定间隔ΔX=Xi+1-Xi得到若干Xi(i=1,2,…,n),利用式(2)可计算出对应的Zi(i=1,2,…,n),从而得到一组坐标(Xi,Zi),这组坐标位于过视线的铅垂面内。将Zi代入式(1)即可求得视线方向Yi,则点(Xi,Yi,Zi)位于像点、物点及摄影中心所决定的视线上,Yi即为视线上对应点高程,即H(Li)。
$Z = {Z_S} + \left({X - {X_S}} \right)\frac{{{c_1}\left({x - {x_0}} \right)+ {c_2}\left({y - {y_0}} \right)- {c_3}f}}{{{a_1}\left({x - {x_0}} \right)+ {a_2}\left({y - {y_0}} \right)- {a_3}f}}$ | (2) |
依据计算出的坐标(Xi,Zi),利用空间数据查询功能在DEM数据中查询对应的高程值H(Pi),若不在DEM格网上时需采用双线性法等内插方法内插得到。当ΔH=H(Li)-H(Pi)小于一定的阈值时,该点(Xp,Yp,Zp)即为摄影视线方向与该方向DEM剖面线的交点,也就火点所在位置。阈值的选择可根据DEM数据高程精度确定。
2.3 基于单幅影像及DEM火点定位流程基于单幅影像及DEM数据火点定位按以下步骤实现火点定位:
1)在地势高处S安装视频监测设备,获取影像数据及摄影时三维姿态角(φ,ω,κ),利用手持GPS测量摄影中心在地面测量坐标系中坐标,并转换为摄影测量坐标系O-XYZ中坐标(XS,YS,ZS);
2)通过3D空间分析功能对监测区域视线方向进行空间搜索,视线方向高程H(Li)与DEM数据高程H(Pi)之差ΔH=H(Li)-H(Pi)小于一定的阈值时,即认为视线方向与沿视线方向的剖面线相交,交点就是定位点。将Y=H(Pi)代入式(1)中,计算得到待定点P在摄影测量坐标系O-XYZ中坐标(Xp,Yp,Zp),并转换为地面测量坐标系坐标。
3 精度分析 3.1 误差源分析由单幅影像火点定位数学模型可知,火点定位的精度主要取决于2方面的精度: 一是视线方向的精度; 二是DEM数据的高程精度。具体而言,火点位置(X,Z)定位精度与像点坐标(x,y)、相机内方位元素(x0,y0,f)、摄影中心坐标即像片外方位线元素(XS,YS,ZS)、像片外方位角元素(φ,ω,κ)及DEM高程精度即Y坐标的精度等因素有关。
由于普通摄像机或相机属于非量测相机,存在相机畸变且相机内方位元素(x0,y0,f)未知,需利用纠正校正方法(Heikkil et al.,1996; Tsai,1987)如Heikkila的四步校正法(Heikkil et al.,1997;Faugeras et al.,1987; Shih et al.,1993; Wei et al.,1993; Weng et al.,1992)对相机进行校正获取内方位元素(x0,y0,f)及相机变形参数。经纠正后像元坐标(x,y,f)误差很小。拍摄像片时采用手动调焦模式,相机的主距f 与相机纠正时主距一致。
摄影中心坐标(XS,YS,ZS)可由带广域差分功能的手持GPS接收机测定,其精度可达0.5 m左右。摄影中心坐标(XS,ZS)的误差只对火点位置产生平移影响,而摄影中心高程YS的影响与DEM高程精度综合一起影响火点定位的精度。目前固定安装在摄影机上的三维电子罗盘实时测定拍摄像片时的三维姿态角方位角(φ)精度可达0.3°左右,倾角(ω)、侧滚角(κ)的精度可达0.2 °左右。
DEM可由1: 10 000地形图经矢量化后建立,DEM高程精度取决于1: 10 000地形图的高程精度。对于山区1: 10 000地形图基本等高距为10 m,设DEM数据高程精度为2倍基本等高距,即DEM数据的高程误差为±20.0 m。
3.2 仿真分析为了定量分析各参数对火点定位精度的影响,对式(1)两边分别求全微分可得:
$\begin{array}{l} {\rm{d}}X = {\rm{d}}{X_S} + {a_{11}}{\rm{d}}Y - {a_{11}}{\rm{d}}{Y_S} + {a_{12}}{\rm{d}}x - {a_{12}}{\rm{d}}{x_0} + \\ \;\;\;\;\;\;\;\;{a_{13}}{\rm{d}}y - {a_{13}}{y_0} + {a_{14}}{\rm{d}}f + {a_{15}}{\rm{d}}\varphi + {a_{16}}{\rm{d}}\omega + \\ \;\;\;\;\;\;\;\;{a_{17}}{\rm{d}}\kappa ;\\ {\rm{d}}Z = {\rm{d}}{Z_S} + {a_{21}}{\rm{d}}Y - {a_{21}}{\rm{d}}{Y_S} + {a_{22}}{\rm{d}}x - {a_{22}}{x_0} + \\ \;\;\;\;\;\;\;{a_{23}}{\rm{d}}y - {a_{23}}{y_0} + {a_{24}}{\rm{d}}f + {a_{25}}{\rm{d}}\varphi + {a_{26}}{\rm{d}}\omega + \\ \;\;\;\;\;\;\;{a_{27}}{\rm{d}}\kappa . \end{array}$
式中: a1i,a2i(i=1~7)分别为X,Z对各参数的偏导数。利用式(3),通过Matlab程序语言编程对各参数对火点位置定位精度的影响进行仿真分析。
3.2.1 DEM高程精度对火点定位精度的影响假设像点坐标(x,y)量测误差为1个像元,相机内方位元素(x0,y0,f)已经过纠正,其误差忽略不计。摄影中心坐标(XS,YS,ZS)误差为0.5 m,像片外方位角元素φ的误差为0.3°,ω,κ的误差为0.2°,且φ=κ=0; 并设摄影机安装于地势较高的S 处,摄影中心与火点位置高差为+150 m,拍摄影像时主距为100 mm,火点离摄影中心距离约为1 000,2 000 m时,此时像片倾角(ω)分别为-4° 17' 21″,-8°31'51″。设DEM数据的高程误差为±20.0 m,像点位置为(1 000,1 000)处火点定位X及Z方向精度如图 3所示。计算结果显示,当火点离摄影中心距离约为1 000 m时,即像片倾角(ω)为-8'31'51″时,DEM 高程误差对火点位置X方向及Z方向的影响分别为-19.2~21.3,2.5~-4.5 m。
当火点离摄影中心距离约为2 000 m时,即像片倾角(ω)为-4°17'21″时,DEM高程误差对火点位置X方向及Z方向的影响分别为-19.0~21.1 m,0.9~-3.0 m(图 4)。由此可以看出,DEM高程误差对火点位置X方向及Z方向的影响呈线性关系,且随着摄影距离的增加,DEM高程误差对火点位置X方向及Z方向的影响约为减弱。
为了分析像平面内火点定位精度分布,假设像点坐标(x,y)量测误差为1个像元,摄影中心坐标(XS,YS,ZS)误差为0.5 m,像片外方位角元素φ的误差为0.3°,ω,κ的误差为0.2°,且φ=κ=0; 并设摄影机安装于地势较高的S处,摄影中心与火点位置高差为+150 m,拍摄影像时主距为100 mm,火点离摄影中心距离约为1 000,2 000 m时,此时像片倾角(ω)分别为-4°17'21″,-8°31'51″; DEM数据的高程精度约为± 20 m。按式(3)计算各项误差综合影响X,Z坐标方向在像平面内的分布。如图 5所示分别为在x,y≤[-500,500],[-200,200],[-10,10]像平面内dX,dZ 变化情况,采样间隔分别为50,20,1像元,由图可以看出:
1)整体上dX,dZ值在像平面内绝大部分区域小于100 m;
2)相对而言,各精度影响因子在Z方向产生的定位误差dZ明显大于X方向定位的误差dX,即沿摄影主光轴方向的定位精度要明显低于垂直于主光轴方向的定位精度;
3)在x,y≤[-500,500]像平面内,当像平面坐标y为[0,300],而x取较大值时,dX,dZ值较大,达到±100 m以上,如图 6。此区域本文称为粗差分布区域。当y取较小值(0<y<300)时,随着x值的增加而定位误差dX,dZ迅速增大,以致无法确定物点位置。在粗差分布区域内利用本文方法无法进行准确定位,而像点离开粗差分布区域以外一定像素位置处,dX,dZ快速收敛,降到森林火点定位精度要求100 m以内。在x,y≤[-500,500]像平面内,dX大于100 m粗差分布区域约占整个像平面的9.9%,dZ大于100 m粗差分布区域约占整个像平面的29.7%。粗差分布区域占整个像平面的比率随着像平面像素的增加而减少。粗差分布区域主要位于像平面的上半部分,而被摄物体大多成像于像平面的下半部分,因此,粗差分布区域占火点定位总区域比率很少。
仿真分析表明,利用摄像机拍摄单幅影像,同时利用三维电子罗盘仪测定摄影时三维姿态角,结合DEM高程数据确定森林火点位置的方法,在粗差分布区域以外的像平面内量测像点坐标,依据共线方程及DEM高程数据能快速准确地确定森林火点位置,且精度可保证在100 m以内,但在粗差分布区域无法通过本方法进行准确位置,需要采取其他措施才能准确定位。
4 试验验证 4.1 试验地点试验地点为大围山国家森林公园,位于湖南省浏阳市东北部,地处湘东幕阜山—罗霄山接壤地带的大围山支脉,地理位置为114°2'—114°12'E,28°21'—28°26'N,属于中亚热带常绿阔叶林地带。最高海拔七星岭1 607.9 m; 最低海拔花门电站230 m,相对高差1 377.9 m。年平均气温11~16 ℃,1月平均温度-4~2.5 ℃,7月平均温度20~28 ℃。年平均降雨量1 200~2 000 mm。面积约3 700 hm2,园内森林茂密,森林覆盖率高达99%,近年旅游资源开发迅速,旅游观光人员多,森林火险高。摄站点设在细棚子附近的山头上。
4.2 试验仪器设备及试验方法1)本次试验影像采集采用佳能EOS 550D单反相机,配EF50(F1.8)镜头。其主要技术参数为:分辨率为5 184×3 456,芯片大小22.300 mm×14.900 mm,像素大小为4.3 μm。摄影时采用手动调焦,焦距为相机纠正时的焦距f(1 279个像素单位)。被摄区域主要为森林,影像中应有部分明显地物,以便精度验证。
2)将佳能EOS 550D单反相机和三维电子罗盘仪用胶带固定,三维电子罗盘仪方位向与像平面的法线指向物镜方向一致,三维电子罗盘仪倾角为零时,像平面的法线位于水平面内。此时,三维电子罗盘仪测定的三维姿态角方位角、倾角、侧滚角对应于摄影测量坐标系统中像片外方位角元素(φ,ω,κ)。三维电子罗盘仪测量摄影瞬间的三位姿态角方位角(φ)为186.7°,倾角(ω)为-5.7°、侧滚角(κ)为1.3°。三维电子罗盘仪采用瑞芬DCM600,三维姿态角方位角(φ)、倾角(ω)、侧滚角(κ)的标称精度分别为0.3°,0.2°,0.2°。
3)摄站点位置采用能接收广域差分信号的Trimble手持GPS接收机测定,信号解算状态为固定解(fixed),并将摄影点WGS84坐标转换为北京54坐标系统下的相机中心坐标(XS,YS,ZS),坐标为(3 144 397.1,510 332.9,935.6) m。
4)对1: 10 000纸质地形图沿等高线采集高程点,建立TIN并重采样为规则格网DEM,格网间距为5 m。
5)按2.3节所述流程量测像点坐标,并通过共线方程及DEM数据求得物点位置。
4.3 试验结果假定被摄区域内明显地物点如道路拐点、山顶点为森林火点。利用本文方法计算被摄区域内10个明显地物点坐标(图 7,8),并利用高精度手持GPS测量进行精度验证,其坐标差如表 1所示。试验结果显示,利用本文方法计算出坐标和手持GPS实测坐标差均在±50 m,能满足火点定位要求,证实本方法可行。
本文基于单幅影像摄影测量共线方程,利用数码相机拍摄森林区域的影像和三维电子罗盘测定拍摄影像时三维姿态角,结合被摄区域的DEM数据,确定森林火点位置,通过理论仿真分析和现场试验证明,除粗差分布区域外该方法能准确确定森林火点位置,火点位置精确可控制在100 m以内,能满足森林防火对火点定位精度的要求。火点定位的精度与三维电子罗盘测定三维姿态角的精度、相机的内方位元素校正精度及DEM数据的高程精度等有关,且沿摄影主光轴方向的定位精度要明显低于垂直于主光轴方向的定位精度。DEM高程误差与火点定位误差呈线性关系,提高DEM的精度能有效地提高火点定位精度。
该方法存在定位误差较大的区域,即粗差分布区域。粗差分布区域位于像平面坐标y为[0,300]而x取较大值时区域。在x,y≤[-500,500]像平面内,dX大于100 m粗差分布区域约占整个像平面的9.9%,dZ大于100 m粗差分布区域约占整个像平面的29.7%。粗差分布区域占整个像平面的比率随着像平面像素的增加而减少。在粗差分布区域内,利用本文方法可能无法获得准确的定位结果,需要采取其他措施才能准确定位。粗差分布区域位置分布还与相机焦距、三维罗盘测定三维姿态角的精度等有关。
另外,当多个山峰连线与搜索视线方向平行或接近平行时,在火点定位搜索时可能会出现多值现象,即有多个点为待定点火点位置,需要利用其他信息对待定点进行筛选,最终确定正确的火点位置。
[1] | 关百钧, 魏宝麟. 1992. 世界林业发展概论. 北京: 中国林业出版社.(1) |
[2] | 李德仁, 周月琴, 金为铣. 2001. 摄影测量与遥感概论. 北京: 测绘出版社, 28.(1) |
[3] | 林其钊, 舒立福. 2003. 林火概论. 合肥: 中国科学技术大学出版社.(1) |
[4] | 刘亚文, 张祖勋, 张剑清, 等. 2005. 利用地图与单影像进行建筑物三维重建的新途径. 武汉大学学报: 信息科学版, 30 (2): 146-149.(1) |
[5] | 舒立福, 田晓瑞. 1997. 国外森林防火工作现状及展望. 世界林业研究, 10(2): 28-35.(1) |
[6] | 孙敏. 2007. 多视几何与传统摄影测量理论. 北京大学学报: 自然科学版, 43(4): 453-459.(1) |
[7] | Faugeras O D, Toscani G. 1987. Camera calibration for 3D computer vision. Proc International Workshop on Industrial Applications of Machine Vision and Machine Intelligence, Silken, Japan, 240-247.(1) |
[8] | Heikkil J, Silven O. 1996. Calibration procedure for short focal length off-the-shelf CCD cameras. Proc 13 th International Conference on Pattern Recognition. Vienna, Austria, 166-170.(1) |
[9] | Heikkil J, Silven O. 1997. A four-step camera calibration procedure with implicit image correction. IEEE Computer Society Conference on Computer Vision and Pattern Recognition. San Juan: IEEE Comput Soc Press, 1106-1112.(1) |
[10] | Shih S W, Hung Y P, Lin W S. 1993. Accurate linear technique for camera calibration considering lens distortion by solving an eigenvalue problem. Optical Engineering , 32(1): 138-149.(1) |
[11] | Tsai R Y. 1987. A versatile camera calibration technique forhighaccuracy 3D machine vision metrology using off-the-shelf TV cameras and lenses. IEEE Journal of Robotics and Automation, 3(4): 323-344.(1) |
[12] | Wei G Q, Ma S D. 1993. A complete two-plane camera calibration method and experimental comparisons. Proc 4 th International Conference on Computer Vision, Berlin, Germany, 439-446.(1) |
[13] | Weng J, Cohen P, Herniou M. 1992. Camera calibration with distortion models and accuracy evaluation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 14(10): 965-980.(1) |