SAR+图像在震害变化检测中的应用.pdfVIP

SAR图像在震害变化检测中的应用①② 窦爱霞1 张景发： 王晓青- 魏久传3 孟庆强3 (1．中国地震局分析预报中心，北京100036；2．中国地震局地壳应力研究所， 北京100085；3．山东科技大学地球科学信息与工程学院，山东泰安271019) 摘 要：本文主要探讨了利用SAR图像和变化检测技术在震害评估中的作用。提出采用单个 相关系数空间分布检测地震破坏区域和利用多个相关系数检测建筑物破坏程度的变化检测方法， 并以张北地震前后SAR图像进行实验，得到满意效果。 关键词：SAR图像变化检测震害 1．引言 遥感图像震害变化检测技术是利用地震前后的两幅图像确定地震造成的地物是否被破 坏的技术，常用的方法有假彩色合成法、图像差值法、图像比值法、相关系数法等…。以 往使用差值、比值法对中低分辨率的遥感图像效果不是很理想；相关系数法仅仅统计图像 灰度值的相关性，并且只是对统计值进行简单的密度分割，结果误差浮动较大。而随着合 视。同光学图像相比。SAR图像不受天气、云层等因素的影响，可以全天候、全天时获得 遥感数据，是较好的震害检测数据源。 利用SAR图像对强震中建筑物等的受灾情况进行检测，国内外学者进行了大量的研 究。Ronald 12 1。Masashi 和6m分辨率的机载SAR震害影像，获取高烈度区域的房屋毁坏率，并以此为特征参量， 快速确定了极震区14J。 总之，机载和星载SAR图像对检测形变和受灾建筑物都有重要意义。但是，由于SAR 图像在成像机理上与光学图像的差异，常规的变化检测方法直接应用于SAR图像会存在一 些问题。本文研究了基于SAR图像单特征相关空间分布和多特征相关进行震害变化检测的 实验分析，结果证实了特征相关检测震害的有效性。 ①本文是“十五”国家“863”课题(编号2001AAl36040)的研究成果。 ②本文是第一作者硕士论文的部分研究内容。 ．132： 2．SAR图像震害变化检测方法 常用变化检测方法，如比值法、差值法，考虑的仅仅是单个像素包含的信息。然而，相 邻像元的灰度值相互之间有很明显的相关性，因此，考虑相邻区域而不是孤立的像元时，会 更可能检测出图像发生的变化。由此，本文使用了更加准确的基于邻域的变化检测算法——’ 相关系数。 利用地震前后中低分辨率的SAR图像的相关系数，可从整体上检测地震破坏的区域及 其空间分布，也可基于少量先验知识评估地物破坏程度与类型。 2．1相关系数 相关系数反映两个函数在不同的相应位置上互相匹配的程度，它不随灰度的线性变化而 变化，可以消除地震前后由于图像接收条件、接收时间的不同造成的影像之间灰度与对比度 的差异。图像差异越大，相关性越小；’反之越大。变化检测使用的两图像可以是原灰度图 像，也可以是增强后的纹理图像。 以点(1，s)为中心的影像窗口之间地震前后两幅图像的相关系数的算法如下： ∑(札；一菇)(yl厂l，) v=————————j=i————————————————————————————一 。 f m．n ／m．n 、／善(噩。一舶2、／善(yf，s．y)2 式中 m，n——窗口大小； x，y一分别表示主影像和从影像的窗口内灰度平均值。 窗口选择可以是规则窗口，也可以是以村镇、居民区等边界为边的不规则多边形窗口。 相关运算后的图像灰度值小于1，为了便于统计分析，可将结果乘以某常数K。 2．2震害区域的检测 基于相关系数检测地震破坏区域，其方法如下：首先建立居民区边界GIS矢量，计算每 个居民区窗口地震前后SAR图像的相关系数，统计出相关系数文件；最后通过分析相关系 数的空间分布来确定地震破坏情况。流程如图l。 图1．震害区域检测流程 2．3震害程度与类型的检测 使用相关系数获取震害信息，首先建立居民区感兴趣区iROI)，然后选择合适窗口进行 ·133· 特征运算，再对居民区特征图像进行相关运算，训练不同破坏等级下