1 引 言
大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] 。快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定。震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性。震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] 。作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] 。
21世纪以来,遥感卫星的数量和时—空—谱分辨率不断提升,为地震应急提供了丰富的数据支撑。应用于地震灾情调查的遥感数据可分为3类,即光学影像、SAR和LiDAR[5 , 6 ] ,它们有着各自的特性和优缺点,可以从不同的角度反映灾情。其中光学影像具有分辨率高、直观、可靠等优点,便于目视解译,在震后灾情调查这样形势紧急、不容有误的情况下更能快速有效地获取初步且珍贵的灾情信息,因此光学影像(包括与其他数据的融合)一直是震后灾情调查中主要使用的数据[6 -8 ] 。如王晓青等[3 ] 在汶川地震发生后,利用高分辨率光学遥感影像,对地震灾害及其损失进行了快速评估,弥补了地面调查的不足,为地震应急指挥和救援决策提供了重要依据。
依据是否采用震前的遥感影像,地震灾害信息可通过仅采用震后影像进行解译和同时采用震前和震后影像进行变化检测这两种方法进行提取[6 ,9 ] 。震后地面的基础设施和土地覆盖相对于震前发生变化的部分就可能是地震造成损害的区域[7 ] ,因此通过对震前、震后的两幅影像的变化检测就可以提取出震害信息。由于震前影像可获取地物在震前的状态,因此变化检测的方法能够更准确地获取灾情信息,受到了更为广泛的应用[6 ] 。但该方法在实际应用中可能会受到局限:首先震前影像可能无法获取或震前影像太过久远,尤其是在欠发达国家和地区;其次,震前和震后影像来自不同的时相、不同的传感器,光照条件、传感器特性等的差异会增大数据处理的难度[10 ] 。所幸的是,随着卫星遥感数据的日益丰富和遥感图像处理技术的进步,这些局限都已得到适当缓解。
因此,基于地震灾情调查中的遥感数据,以及灾情信息提取的方法,对地震应急中使用的光学遥感数据及其产品进行归纳总结,进而对地震灾情提取中的变化检测算法进行综述,详细阐述了各种方法的原理、优势和实际应用情况。
2 高分辨率光学遥感卫星数据
1972年,世界上第一颗遥感卫星Landsat 1发射成功,开启了遥感卫星发展的序幕;20世纪80年代末呈现出了百家争鸣的景象,法国、日本等其他国家相继发射遥感卫星,光学卫星的全色分辨率达到2 m;21世纪以来,一系列高分辨率的商业遥感卫星发射,并广泛应用于对地观测中,时—空—谱分辨率全面提升,且卫星的载荷形式、运行模式、数据处理技术以及应用场景愈加多样和成熟,为地震灾情调查提供了有力的数据支撑。
表1 按全色分辨率由高到低的顺序列出了数据影像已广泛用于地震灾情调查的高分辨率光学遥感卫星[6 ,11 -12 ] 。可以看出,空间分辨率优于1 m的卫星多达十数颗,其中WorldView系列卫星分辨率可达0.31 m。单星的时间分辨率已经可以达到1 d,而且目前的商业卫星均在大力构建卫星星座,获取数据的能力和数据质量显著提升,如美国的WorldView-2–4与GeoEye-1组成的4星星座可达到每日4.5次的重访频率,法国的Pléiades-1A/1B和SPOT-6/7组成的4星星座可实现每日2次的重访(SPOT-6/7提供大幅宽的普查图像,Pléiades-1A/B提供0.5 m分辨率的详查图像),中国的北京2号三星星座可实现每日2次的重访。此外,震后进行多星联合调度观测,可充分利用卫星资源、最大化灾情观测效益,极大地提高了震后应急响应的速度。在传感器配置上,这些卫星多为1个全色波段加4个或8个多光谱波段,波长范围覆盖了可见光到近红外波段(多为0.4~0.9 μm),像元的动态范围可达11 bit。
1999年法国航天局国家空间研究中心(CNES)和欧洲航天局(ESA)倡议发起了“空间与重大灾害国际宪章”(The International Charter Space and Major Disasters),加入该宪章的成员可以无偿获取其他成员的卫星数据用于应急救灾,解决了卫星资源的共享问题[13 ] ,对遥感技术在地震应急中的应用有着深远影响。
使用卫星影像进行灾情调查,需首先进行预处理,包括辐射校正、大气校正、几何校正、图像融合等。表2 列出了美国DigitalGlobe、欧洲Astrium和中国二十一世纪空间技术应用股份有限公司(21at)旗下的遥感影像产品[14 -16 ] 。可以看出用户拿到的数据已经过不同程度的预处理,因此用户可以根据实际需求选择合适等级的数据产品,无需考虑数据预处理的问题。此外,数据服务商可为用户提供加急处理选项(Rush option)以确保震后数据的时效性。
3 震害提取的变化检测算法
地震发生后,基于获取的遥感影像,需要快速查明3类震害信息:倒塌的建筑物、损毁的生命线工程和次生地质灾害等[9 ] 。其中建筑物震害是造成人员伤亡的主因,也是灾情调查的重点和难点;生命线工程包括道路、电线、水管、通讯等设施,关系到受灾人员的生命维持和现场第一手资料的获取;滑坡、泥石流、堰塞湖、土壤液化等次生地质灾害严重威胁灾区居民的安全和灾区的生态系统。这3类震害在发生的区域、尺度大小、形态特征和表现形式等方面存在着明显的差别,需要采用不同的数据和方法进行提取。本研究将不拘于这3种震害类别,而是着眼于震害信息提取的变化检测算法,对地震灾情调查的现状进行综述,以期为未来地震应急中的灾情提取工作提供参考和裨益。
3.1 变化检测算法概述
变化检测的基本流程主要包括影像预处理、变化检测和精度评价[17 -18 ] 。
3.1.1 震前和震后影像预处理
在地震应急情况下,震前、震后的影像可能来自不同的传感器,获取影像的条件(如太阳高度角、光照条件和大气状况等)往往存在着差异,需要消除这些因素造成的虚假变化。几何配准和相对辐射校正是影响变化检测精度的关键步骤[18 ] 。
几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数。几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] 。Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测。面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性。Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”。
相对辐射校正是指通过比较和修改震前、震后影像的灰度统计量,使其具有趋于一致的光谱特征,以减少影像由于接收条件的不同而造成的灰度与对比度的差异[28 , 29 ] 。地震灾情调查中常用的相对辐射校正方法有直方图匹配法[29 -32 ] 、暗集—亮集法[33 ] 、线性回归法[34 ] 等。
3.1.2 变化检测
完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情。对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类。最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述。
3.1.3 精度评价
在地震应急情况下,不易及时获取实地调查的真实数据,因此对变化检测结果的精度评价通常需要借助高分辨率的震后影像,以目视解译的结果作为真值,构造混淆矩阵,计算精度评价指标,如总体精度、生产者精度、使用者精度、κ 系数等。
3.2 基于像元的变化检测算法
基于像元的变化检测方法是以像元作为基本处理单元,优点是简单、直接,能精确定位变化的位置,但对震前、震后影像的配准精度,以及相对辐射校正和变化阈值的选取有着很高的要求。依据提取变化信息的原理的不同,基于像元的算法又可分为直接比较法、特征提取法和分类后比较法。
3.2.1 直接比较法
(1)代数法
代数法是通过震前和震后影像对应波段间的代数运算来提取变化信息,常用的代数运算方法包括:
差值: I d x , y = I p o s t x , y - I p r e x , y (1)
比值: I r x , y = I p o s t x , y I p r e ( x , y ) (2)
归一化差值: I n d x , y = I p o s t x , y - I p r e ( x , y ) I p o s t x , y + I p r e ( x , y ) (3)
变化向量: I ⃗ d x , y = I ⃗ p o s t x , y - I ⃗ p r e x , y = R p o s t ( 1 ) - R p r e 1 , ⋯ , R p o s t ( n ) - R p r e n T (4)
其中:I p r e ( x , y ) 和I p o s t x , y 分别表示震前和震后影像在坐标x , y 处的像素值(如DN值、辐射亮度或反射率等);R p r e i 和R p o s t i 的上标i 表示第i 个波段。
差值法的思想来自于对“变化”二字的直观理解,是变化检测算法的基础,其他的方法都或多或少地用到了差值法。由于差值法太过简单,仅采用差值法提取震害信息的文献很少,如Yusuf等[38 ] 在2001年印度古吉拉特邦地震后,采用差值法计算出了震前、震后两幅Landsat-7 ETM+影像的差值影像,设定阈值提取了建筑物倒塌的区域和土壤液化区域;Ishii等[39 ] 通过震前和震后航空影像的色调、亮度的差值提取出了建筑物的损毁区域。归一化差值是对影像的差值进行归一化,凸显了相对变化量,在震害提取中应用较为广泛,如Romaniello等[23 ,40 -41 ] 人的研究。变化向量法是在差值法的基础上扩展而来的,通过震前、震后影像的所有波段的差值构建变化向量,能够避免单一波段比较带来的信息不完整,如Zhao等[34 ] 在利用Landsat 8影像提取尼泊尔震后的滑坡时,采用7个波段的反射率和3种植被指数(NDVI、EVI、SAVI)构建了两组变化向量。
代数法优点在于简单、易操作,便于解释变化检测结果,但在实际应用中存在着较多不足:①代数法对震前、震后影像的几何配准和相对辐射校正有着很高的要求,用于高分辨率影像时会产生明显的“椒盐效应”;②对于地震应急中常用的多光谱数据,简单的代数运算没有考虑波段之间的相关性,存在着信息重叠和冗余,且波段间的运算使得不同的数值可能会产生同样大小的差值或比值,忽略了不同地物在不同敏感波段存在的差异,丢失潜在的可利用信息[42 ] ;③由于存在“同质异谱”和“异质同谱”的现象,代数法提取的震害信息可能存在误检或遗漏。
(2)回归法
回归法考虑了地物自身随时间的变化,震后的理论影像I ̂ p o s t 可以通过震前一系列影像的最小二乘回归得到,获取的真实影像I p o s t 与理论影像之差即为地震灾害引起的变化:
I d x , y = I p o s t x , y - I ̂ p o s t ( x , y ) (5)
地震应急情况下,通常不容易获取震前多个时相的数据,因此几乎没有回归法应用于地震灾情的案例。Kohiyama和Yamazaki[22 ] 在提取2003年巴姆地震的建筑物震害时,提出了一种类似于回归法的方法——图像波动模型法(Image Fluctuation Model,IMF),认为不同时刻同一像元的DN值由于传感器、大气、光照等差异而产生微小波动,利用震前的17幅影像获得其波动模型,计算震后影像的置信水平,从而得到建筑物发生损毁的可能性。回归法的优点在于采用了大量的震前影像对变化的趋势进行拟合或假设,其变化检测的结果较为可靠,但对数据要求过高,且每个时相的影像均要进行配准,计算量大,因此在地震应急中并不常用。
(3)相关系数法
相关系数法是通过计算震前、震后两幅影像的相关系数来判断变化的程度。假设选定的移动窗口大小为m ×n ,则两幅影像对应窗口的相关系数为:
γ = ∑ i , j I p r e i , j - I ¯ p r e I p o s t i , j - I ¯ p o s t ∑ i , j I p r e i , j - I ¯ p r e 2 ∑ i , j I p o s t i , j - I ¯ p o s t 2 (6)
其中:i = 1,2 , ⋯ , m ; j = 1,2 , ⋯ , n ;I ¯ p r e 和I ¯ p o s t 分别表示两个时相窗口内像元的均值。相关系数越大,表示震前震后影像的变化越小,即损毁程度越小。相关系数不随影像灰度的线性变化而变化,可以消除地震前后由于影像接收条件的不同所造成的灰度与对比度的差异[29 ] 。
Zhang等[43 ] 对1999年台湾南投地震的震前震后影像,通过设定相关系数的阈值范围,将建筑物的损毁程度划分成了4个等级;Rathje等[44 ] 采用了相关系数法提取巴姆城区震后的建筑物震害,并发现灰度共生矩阵(Gray Level of Co-occurrence Matrix,GLCM)的相关系数比DN值的相关系数能更好地区分损毁和完好的房屋;Chesnel等[26 ] 以屋顶范围作为窗口计算相关系数,实现了震前震后屋顶矢量的精确配准;张伟等[24 ] 为了减小配准误差对相关系数的影响,通过在对应点周围搜索相关系数最大的潜在同名点,使得倒塌和未倒塌建筑的分类精度(κ 系数)由0.60提高到了0.85。可见,相关系数不仅可以作为变化检测的特征变量,也可以作为几何配准的相似测度。
3.2.2 特征提取法
(1)植被指数法
为了凸显植被信息,通过强化不同波段的反射和吸收特性,构造了各种植被指数,其中最常用的是归一化的植被指数(NDVI):
N D V I = R N I R - R r R N I R + R r (7)
其中:R N I R 、R r 分别表示近红、红波段的反射率。
在地震应急中,通常并不关注植被的变化情况,甚至需要剔除植被的影响,如Yamazaki等[22 ,44 -46 ] 在提取城市的建筑物震害信息时,对植被进行了掩膜处理,以消除虚假的变化和干扰;Chini等[41 ] 发现若不对植被进行甄别,20%~50%未受损的建筑物会被误判为“倒塌”。但在某些情况下可通过植被的变化提取需要的灾情信息,如热带地区或山区植被的明显减少表明可能发生了山体滑坡[47 ] 。Yang等[33 ] 使用震前的Landsat 5 TM和震后的Terra ASTER影像通过波段运算得到NDVI,通过对比试验表明了相对辐射校正的必要性和植被指数法对滑坡检测的有效性。
(2)主成分变换法
主成分变换是将具有相关性的多波段影像压缩到完全独立的较少的几个主成分上,使新影像更便于解译。主成分变换法进行变化检测原理是,未变化的区域具有较高的相关性,而发生变化的区域相关性相对较低[48 ] 。在提取变化信息时,主成分变换法常与差值法结合使用,依据差值和主成分变换的先后可分为差值主成分法、主成分差值法和多波段主成分变换法[48 -49 ] 。主成分变换的缺点是变换后各分量所突出的信息有所不同且物理意义并不明确,具体哪个分量可用来表征变化信息需要一定经验和分析[17 ] 。
Tomowski等[31 ] 以2010年海地地震后太子港总统府区域为研究对象,从GLCM能量值第一主成分的差值中提取了损毁的建筑、植被和新增的救灾帐篷等变化信息;李成龙等[50 ] 采用主成分差值法对2003年阿尔及利亚地震的建筑物损毁情况进行检测,然而检测结果中存在很多零散杂乱的图斑,尤其受建筑物阴影变化的影响很大,需要较多的人工干预;Gong等[51 ] 对比了主成分分析法与差值法、面向对象法的震害提取效果,发现主成分分析法并无优势,提取的变化多为建筑物阴影的变化而非建筑物的损毁。此外,穗帽变换作为一种特殊的主成分变换方法,可以得到具有物理意义的主成分,也在震害提取中得到应用,如Oommen等[52 ] 基于2001印度古吉拉特邦地震震前、震后的Landsat 7 ETM+影像,利用穗帽变换得到的湿度分量,探测出了灾区土壤液化的分布情况。
3.2.3 分类后比较法
分类后比较法是对震前和震后影像中的像元进行独立分类,然后比较分类结果,如果对应的类别相同,则判定没有发生变化,否则认为发生了变化。分类后比较法在一定程度上借鉴了面向对象变化检测方法的思想,只对感兴趣的对象进行比较,可以直接获得变化的位置、类型和数量等信息,而且通过对震前和震后影像的单独分类,克服了影像由于获取条件、分辨率等因素造成的差异[18 ] 。
Xu等[53 ] 采用分类后比较法探测出了2008年汶川地震中的唐家山堰塞湖,采用决策树的方法分别对配准好的震前、震后影像进行分类,将震前影像提取出的水体和震后影像提取出的水体、滑坡进行叠加,发现滑坡阻塞河道导致河道变宽形成堰塞湖,进而定量估算了河道增加的宽度和被淹区域的面积。
然而,分类后比较法的缺点限制了其在震害信息提取中的应用:①对每幅影像单独分类,增加了计算量和复杂度,影响时效性;②两次分类的误差会传递到变化检测的结果中而被放大[18 , 42 ] 。此外,震害信息提取中,变化检测的关注点主要在于“是否”发生了变化或发生变化的程度,至于具体类型的变化(from-to)并非必须获取的信息。
3.3 面向对象的变化检测算法
基于像元的变化检测算法在应用于高分辨率影像时常常会受到局限。首先,分辨率越高,像元越小,空间不一致或者几何配准误差等造成的影响越加明显;其次,高分辨率影像在承载更丰富的地表信息的同时,也导致了同一地物内部像元灰度起伏明显,在使用基于像元的方法时,过多无关的变化会被检测出来,产生“椒盐效应”[18 , 54 ] 。
面向对象的变化检测方法能很好地克服基于像元的方法在配准和噪声等方面的局限[55 ] ,通过影像中地物对象的提取,直接比较这些地物的光谱特征、纹理特征、形状特征和上下文关系等的变化,更加高效可靠[56 ] ,此外还能减少光照条件、阴影等差异造成的影响[57 ] 。面向对象方法的关键在于对象的提取(包括分割和分类),而对变化信息的提取基本上都是借鉴了基于像元的方法,所不同的只是将像元间的比较变成了对象间的比较[57 ] 。Bitelli等[25 ,55 ,58 ] 对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害。
基于Hussain等[18 ] 的分类方法,结合地震灾情调查中的实际需求和应用情况,本研究将地震中面向对象的变化检测方法分为了3类:目视解译法、直接比较法、多时相叠加法。
3.3.1 目视解译法
震后的灾损评估可分为两大类,即基于算法的定量分析和基于人工解译的定性分析,应当根据具体的任务目标选择合适、有效的分析方法[59 ] 。
人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译。人工解译虽然费事费力,但精度和可靠性相对更高[58 ] 。对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] 。
Saito等[59 ] 采用震后的IKONOS全色影像,并结合现场实地调查的照片,对目视解译的效果进行评价,表明了目视解译的快速性和对部分建筑损毁解译的局限性;Yamazaki等[62 ] 采用巴姆地震震前和震后的QuickBird影像,对12 063个建筑物的损毁程度逐个进行了目视解译;Ehrlich等[7 ] 通过震后映秀城区的WorldView影像、震前和震后都江堰城区的SPOT-5影像识别了滑坡灾害和建筑物的损毁情况;吴玮莹等[63 ] 利用震后获取的高分二号、北京二号影像和震前的Google Earth影像,通过人工目视解译的方法,以道路整齐边界遭受的破坏和山体植被覆盖的变化作为解译标志,初步建立了九寨沟地震的同震滑坡分布图。为便于目视解译,Zhang等[43 ] 在提取1999年南头地震的灾情信息时,将震前和震后共2幅灰度影像进行假彩色合成显示,在新合成的彩色图像中,未发生变化的区域为白色—灰色,发生变化区域表现为彩色,这些色彩就包含了震害信息,据此可判断震害分布及其性质。
此外,许多分析算法也是半自动的(semi-automated algorithm)或人机交互的,特别是阈值的设定、震害提取效果的评定等关键问题,都或多或少的需要人工目视解译的参与。然而不可否认的是,如今计算机技术和图像处理算法的愈加成熟,对于大范围的灾情调查,自动化的变化检测算法更快速,是现今震害信息检测的发展趋势。
3.3.2 多时相叠加法
多时相叠加法是指将震前和震后影像叠加在一起,只对震前的影像进行对象的分割和分类,将分类后对象的掩膜直接应用到所有的时相中,仅对掩膜所包含的区域进行比较。多时相叠加法的优点在于对多个时相的影像只需进行一次对象分割和分类,且同一对象具有相同的几何特征[18 ] ,但该方法需要进行影像的几何配准。
Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%。Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类。Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图。Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势。Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物。Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物。Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化。Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰。
3.3.3 直接比较法
直接比较法是指将震前和震后影像中分别分割出来的对象进行直接比较,以提取变化信息。直接比较法在提取变化信息时借鉴了基于像元的方法,直接比较震前、震后同一地物对象在光谱、纹理、形状等特征上的变化,因此可以使用前文所述的基于像元的算法。相比于多时相叠加法,直接比较法对配准的要求相对较低,只需要保证震前震后影像中同一对象间的匹配,但直接比较法通常需要进行两次分割和分类。
Matsumoto等[55 ,67 ] 对震前、震后的影像进行分类,提取出两个时相的建筑物,通过建筑物的减少来判别损毁。类似地,Gong等[51 ] 根据房屋的长宽比、纹理、尺度等特征建立分类的规则集提取出了震前和震后影像中的房屋对象,认为震后相比于震前减少的房屋对象即为被损毁的房屋。Tiede等[27 ] 在未对震前、震后影像配准的情况下,通过对比对应建筑物阴影的面积和形状的变化来判断建筑物是否损毁;Zhu等[68 ] 对震后影像进行分割,通过决策树分类提取出人造地表(如道路、停车场)和裸土,采用多波段主成分变换法判断这些地表是否发生液化。Anniballe等[32 ] 借助于拉奎拉市已有的GIS底图,利用建筑物高度和传感器姿态等数据精确计算了震前和震后影像中建筑物的矢量多边形,测试了4类共13种影像特征(光谱特征、信息论特征、纹理特征、HSV颜色特征)和2种分类方法(最大后验概率和支持向量机)对建筑物震害的提取效果。
如表3 所示,总结和对比了前文所述的地震灾情调查中的变化检测算法,从中可以看出各个算法之间的差异和优缺点。
4 讨 论
4.1 灾区的数据获取问题
在利用变化检测方法进行震害信息提取时,用到的数据主要包括GIS底图、震前和震后影像、验证数据等。
GIS底图可以极大地便利震害信息提取,降低问题的复杂度,减少计算量和工作量,提高效率和精度,如Samadzadegan和Rastiveisi等[30 ,32 ,64 -66 ] 直接采用了GIS底图提供的建筑物矢量,省去了对象分割和建筑物分类等工作。因此,在地震危险度高的城市应大力开展城市GIS底图的绘制和更新工作。
进行变化检测时,震前和震后的影像是必须的,且时效性尤为关键。研究统计了相关文献所采用的数据的获取时间,同时搜集了CHARTER机制响应地震灾害而生成的光学影像产品,统计了这些数据的获取时间,如图1 所示,图中横坐标表示地震发生的日期,纵坐标表示震前或震后影像获取时间相对于地震发生时间所间隔的天数。可以看出,绝大多数震后影像都能在1~3 d内获取,基本上可以满足地震应急的时效性需求;震前数据则相对久远,多为一年以前的,期间地表覆盖可能发生了变化(如植被的生长更替、建筑物的拆迁和新建等),不利于震后的变化检测。因此,应对地震危险度高的城市定期观测和更新本地的遥感影像。总体上,震后数据获取的快慢随发震日期的推移并无明显变化,可能与具体的区域、地震的影响程度等因素有关。
图1
图1
各地震前后的光学卫星影像获取时间
Fig.1
Acquisition time of pre- and post-earthquake images
验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] 。但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差。Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”。 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] 。因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性。因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性。
4.2 变化检测的方法
为了系统地综述变化检测算法在地震灾情调查中的应用,本研究对这些算法的类别进行了较为细致的划分,但这种类别的划分并非唯一的[18 -19 ,35 ,73 ] ,不同的类别的界限不是严格分明的,可能存在着交叉、借鉴或套用。在实际研究和应用这些方法时,也不会拘泥于算法的类别,如Janalipour等[74 ] 综合采用了差值法、变化向量法、植被指数法、主成分分析法等多种方法,构建了多准则判据来判断建筑物是否损毁,使检测结果更加稳定、精确。
从实际应用的角度,基于像元的方法简单、直观、快速,在震后的灾情调查尤其是早期的快速普查中应用广泛。本研究对基于像元的方法分类均是基于其提取变化信息的原理,从类别的名称基本上可以知晓其进行灾情提取的操作方法。
从学术研究的角度,面向对象的方法能够克服基于像元的方法在配准和噪声等方面局限,且涉及到地物分类、变化提取等诸多问题,可研究的内容更丰富、更具挑战性,因而更受研究学者的青睐,从而发展出了丰富繁杂的面向对象变化检测算法。正因如此,对面向对象的方法的进一步分类是非常困难的。通过删繁就简,结合地震应急中的应用情况,从两个时相影像中对象分割角度对这些方法分门别类,模糊了这些算法提取变化信息的原理。
此外,就地震灾情调查而言,无论是损毁的建筑物和道路,还是新出现的滑坡、堰塞湖等,在震后影像中都相对于震前发生了变化,因此理论上变化检测算法都是可以检测出这些震害类型的。然而不同的震害类型在尺度、光谱、纹理等特征上存在着各自的特点,因此在实际应用中应针对性地选择合适的数据和方法。
5 结 论
大地震后的灾情调查是变化检测算法的一个极为特殊的应用场景,对变化检测精度、地表类型变化等方面的要求相对较低,对方法的实用性和结果的可靠性、时效性有着非常高的要求。因此目前在实际应用中对复杂算法的引入较为谨慎,当震害信息提取的精度达到可接受的程度时,多数研究并不执着于算法精度的提升和横向对比,也很少对已有的方法进行深度挖掘和优化,因此已有的研究虽提出了各种各样的震害提取方法,但系统性和条理性相对较弱,对这些方法的分类较为困难。
在采用高分辨率影像进行震后的灾情提取时,基于像元的方法简单直接,但需要精确配准,且未考虑像元的邻域信息,存在椒盐效应。面向对象的变化检测方法可以很好地克服上述局限,具有较高的检测精度和稳健性,代表了变化检测的发展方向。尤其是随着数字城市的建设,包括GIS底图在内的基础数据越加完备,将显著提高面向对象方法提取震害信息的效率和精度。
在后续的研究和应用中,需要进一步开展和加强的工作包括:①提高卫星遥感数据的时效性,除了卫星星座组网,还可以加强卫星调度技术、去云算法等的研究;②进行GIS底图的绘制和更新,尤其是地震危险度较高的城市;③当前震害信息提取中的关键问题(包括阈值的选取、配准和精度评价等)都倾向于选择简单直观或人机结合的办法,需要引入更多自动化的算法解决这些问题;④基于遥感影像的建筑物损毁等级存在“低估”,需要建立遥感影像与实地调查间的建筑物损毁程度的对应关系,使不同的研究成果之间具有可比性。
参考文献
View Option
[1]
Ma Yuhong , Xie Lili . A study on factors influencing earthquake casualties
[J]. Journal of Nature Disasters , 2000 , 9 (3 ): 84 -90 .
[本文引用: 1]
马玉宏 , 谢礼立 . 关于地震人员伤亡因素的探讨
[J]. 自然灾害学报 , 2000 , 9 (3 ): 84 -90 .
[本文引用: 1]
[2]
Liu Jiahang , Yang Jianfeng , Wei Chengjie , et al . Acquisition of earthquake damage information based on remote sensing technology: history, current situation and trend
[J]. Journal of Natural Disasters , 2004 , 13 (6 ): 48 -54 .
[本文引用: 2]
柳稼航 , 杨建峰 , 魏成阶 , 等 . 震害信息遥感获取技术历史、现状和趋势
[J]. 自然灾害学报 , 2004 , 13 (6 ): 48 -54 .
[本文引用: 2]
[3]
Wang Xiaoqing , Wang Long , Wang Yan , et al . Study on the damage assessment of Wenchuan Ms8.0 Earthquake based on remote sensing imagery
[J]. Technology for Earthquake Disaster Prevention , 2008 , 3 (3 ): 251 -8 .
[本文引用: 2]
王晓青 , 王龙 , 王岩 , 等 . 汶川8.0级大地震应急遥感震害评估研究
[J]. 震灾防御技术 , 2008 , 3 (3 ): 251 -258 .
[本文引用: 2]
[4]
Liu Zhi , Jiang Hui , Chen Kunhao . Research on quasi real-time comprehensive estimation system for earthquake disaster
[J]. Technology for Earthquake Disaster Prevention , 2017 , 12 (4 ): 834 -44 .
[本文引用: 1]
刘智 , 姜慧 , 陈焜浩 . 准实时地震灾情综合评估系统的研发
[J]. 震灾防御技术 ,2017 ,12 (4 ):834 -844 .
[本文引用: 1]
[5]
Rathje E M , Adams B . The role of remote sensing in earthquake science and engineering: Opportunities and challenges
[J]. Earthquake Spectra , 2008 , 24 (2 ): 471 -92 .
[本文引用: 2]
[6]
Dong L , Shan J . A comprehensive review of earthquake-induced building damage detection with remote sensing techniques
[J]. ISPRS Journal of Photogrammetry Remote Sensing , 2013 , 84 : 85 -99 .
[本文引用: 5]
[7]
Ehrlich D , Guo H , Molch K , et al . Identifying damage caused by the 2008 Wenchuan Earthquake from VHR remote sensing data
[J]. International Journal of Digital Earth , 2009 , 2 (4 ): 309 -326 .
[本文引用: 3]
[8]
Matin S S , Pradhan B . Challenges and limitations of Earthquake-induced Building damage mapping techniques using remote sensing images-a systematic review
[J]. Geocarto International ,2021 :1 -27 . DOI:10.1080/10106049.2021.1933213 .
[本文引用: 1]
[9]
Dell'acqua F , Gamba P . Remote sensing and earthquake damage assessment: Experiences, limits, and perspectives
[J]. Proceedings of the IEEE , 2012 , 100 (10 ): 2876 -2890 .
[本文引用: 3]
[10]
Kaya G , Musaoglu N , Ersoy O K . Damage assessment of 2010 Haiti Earthquake with post-earthquake satellite image by support vector selection and adaptation
[J]. Photogrammetric Engineering Remote Sensing , 2011 , 77 (10 ): 1025 -1035 .
[本文引用: 1]
[11]
Joyce K E , Belliss S E , Samsonov S V , et al . A review of the status of satellite remote sensing and image processing techniques for mapping natural hazards and disasters
[J]. Progress in Physical Geography , 2009 , 33 (2 ): 183 -207 .
[本文引用: 1]
[12]
Geiß C , Taubenböck H . Remote sensing contributing to assess earthquake risk: from a literature review towards a roadmap
[J]. Natural Hazards , 2013 , 68 (1 ): 7 -48 .
[本文引用: 1]
[13]
Bessis J L , Béquignon J , Mahmood A . The international charter “Space and Major Disasters” initiative
[J]. Acta Astronautica , 2004 , 54 (3 ): 183 -190 . DOI: 10.1016/S0094-5765(02)00297-7 .
[本文引用: 1]
[14]
Curiel A D S , Gomes L , Purl D , et al . Small satellites for earth oservation [M]. Dordrecht : Springer Netherlands , 2008 .
[本文引用: 1]
[15]
Astrium ,EADS .Pléiades imagery user guide
[EB/OL]. , 2019-06-17 . Accessed 2019-09-07 .
URL
[16]
DigitalGlobe . Digitalglobe core imagery products guide
[EB/OL],, 2013 . Accessed 2019-09-07 .
URL
[本文引用: 1]
[17]
Lu D , Mausel P , Brondízio E , et al . Change detection techniques
[J]. International Journal of Remote Sensing , 2004 , 25 (12 ): 2365 -2401 . DOI: 10.1080/0143116031000139863 .
[本文引用: 3]
[18]
Hussain M , Chen D , Cheng A , et al . Change detection from remotely sensed images: from pixel-based to object-based approaches
[J]. ISPRS Journal of Photogrammetry and Remote Sensing ,2013 ,80 :91 -106 . DOI:10.1016/j.isprsjprs. 2013.03.006 .
[本文引用: 9]
[19]
Gong J , Sui H , Ma G , et al . A review of multi-temporal remote sensing data change detection algorithms
[J]. The International Archives of the Photogrammetry , Remote Sensing Spatial Information Sciences, 2008 , 37 (B7 ): 757 -762 .
[本文引用: 3]
[20]
Gusella L , Adams B J , Bitelli G , et al . Object-oriented image understanding and post-earthquake damage assessment for the 2003 Bam, Iran, Earthquake
[J]. Earthquake Spectra , 2005 , 21 (): 225 -238 . DOI: 10.1193/1.2098629 .
[本文引用: 3]
[21]
Huyck C K , Adams B J , Cho S , et al . Towards rapid citywide damage mapping using neighborhood edge dissimilarities in very high-resolution optical satellite imagery—application to the 2003 Bam, Iran, Earthquake
[J]. Earthquake Spectra , 2005 , 21 (): 255 -266 . DOI: 10.1193/1.2101907 .
[22]
Kohiyama M , Yamazaki F . Damage detection for 2003 Bam, Iran, Earthquake using Terra-Aster satellite imagery
[J]. Earthquake Spectra ,2005 ,21 ():267 -274 . DOI:10.1193/1.2098947 .
[本文引用: 3]
[23]
Romaniello V , Piscini A , Bignami C , et al . Earthquake damage mapping by using remotely sensed data: The Haiti case study
[J]. Journal of Applied Remote Sensing , 2017 , 11 (1 ): 016042 .
[本文引用: 2]
[24]
Zhang Wei , Wang Xiaoqing , Sun Guoqing , et al . An improved change detection algorithm for building damage extraction based on high resolution RS images
[J]. Earthquake , 2013 , 33 (2 ): 87 -95 .
[本文引用: 2]
张伟 , 王晓青 , 孙国清 , 等 . 高分遥感影像震害提取的变化检测算法改进
[J]. 地震 2013 , 33 (2 ): 87 -95 .
[本文引用: 2]
[25]
Bitelli G , Camassi R , Gusella L , et al . Image change detection on urban area: The earthquake case
[C]∥ XXth ISPRS Congress . Istanbul, Turkey , 2004 .
[本文引用: 2]
[26]
Chesnel A , Binet R , Wald L . Object oriented assessment of damage due to natural disaster using very high resolution images
[C]∥ 2007 IEEE International Geoscience and Remote Sensing Symposium . Barcelona, Spain , 2007 .
[本文引用: 3]
[27]
Tiede D , Lang S , Füreder P , et al . Automated damage indication for rapid geospatial reporting
[J]. Photogrammetric Engineering and Remote Sensing , 2011 , 77 (9 ): 933 -942 .
[本文引用: 4]
[28]
Yang X , Lo C . Relative radiometric normalization performance for change detection from multi-date satellite images
[J]. Photogrammetric engineering and remote sensing , 2000 , 66 (8 ): 967 -980 .
[本文引用: 1]
[29]
Zhang Jingfa , Xie Lili , Tao Xiaxin . Change detection of remote sensing image for earthquake-damaged buildings and its application in seismic disaster assessment
[J]. Journal of Natural Disasters , 2002 , 11 (2 ): 59 -64 .
[本文引用: 3]
张景发 , 谢礼立 , 陶夏新 . 建筑物震害遥感图像的变化检测与震害评估
[J]. 自然灾害学报 , 2002 , 11 (2 ): 59 -64 ].
[本文引用: 3]
[30]
Samadzadegan F , Rastiveisi H . Automatic detection and classification of damaged buildings, using high resolution satellite imagery and vector data
[J]. The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences , 2008 , 37 (Part B8 ): 415 -420 .
[本文引用: 3]
[31]
Tomowski D , Ehlers M , Klonus S . Colour and texture based change detection for urban disaster analysis
[C]∥ 2011 Joint Urban Remote Sensing Event . Munich, Germany , 2011 .
[本文引用: 2]
[32]
Anniballe R , Noto F , Scalia T , et al . Earthquake damage mapping: An overall assessment of ground surveys and vhr image change detection after L'aquila 2009 earthquake
[J]. Remote Sensing of Environment , 2018 , 210 : 166 -178 . DOI: 10.1016/j.rse.2018.03.004 .
[本文引用: 6]
[33]
Yang X , Chen L . Using multi-temporal remote sensor imagery to detect earthquake-triggered landslides
[J]. International Journal of Applied Earth Observation and Geoinformation , 2010 , 12 (6 ): 487 -495 . DOI: 10.1016/j.jag.2010.05.006 .
[本文引用: 4]
[34]
Zhao W , Li A , Nan X , et al . Postearthquake landslides mapping from Landsat-8 data for the 2015 Nepal earthquake using a pixel-based change eetection method
[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing ,2017 ,10 (5 ):1758 -1768 . DOI:10.1109/JSTARS. 2017.2661802 .
[本文引用: 3]
[35]
Tewkesbury A P , Comber A J , Tate N J , et al . A critical synthesis of remotely sensed optical image change detection techniques
[J]. Remote Sensing of Environment , 2015 , 160 : 1 -14 . DOI: 10.1016/j.rse.2015.01.006 .
[本文引用: 2]
[36]
Zhang Liangpei , Wu Chen . Advance and future development of change detection for multi-temporal remote sensing imagery
[J]. Acta Geodaetica et Cartographica Sinica , 2017 , 46 (10 ): 1447 -1459 .
[本文引用: 1]
张良培 , 武辰 . 多时相遥感影像变化检测的现状与展望
[J]. 测绘学报 , 2017 , 46 (10 ): 1447 -59 ].
[本文引用: 1]
[37]
Si Salah H , Ait-Aoudia S , Rezgui A , et al . Change detection in urban areas from remote sensing data: A multidimensional classification scheme
[J]. International Journal of Remote Sensing ,2019 ,40 (17 ):6635 -6679 . DOI:10.1080/01431161. 2019.1583394 .
[本文引用: 1]
[38]
Yusuf Y , Matsuoka M , Yamazaki F . Damage assessment after 2001 gujarat earthquake using Landsat 7 satellite images
[J]. Journal of the Indian Society of Remote Sensing , 2001 , 29 (1-2 ): 17 -22 .
[本文引用: 2]
[39]
Ishii M , Goto T , Sugiyama T , et al . Detection of earthquake damaged areas from aerial photographs by using color and edge information
[C]∥ The 5th Asian Conference on Computer Vision . Melbourne, Australia , 2002 .
[本文引用: 1]
[40]
Stramondo S , Bignami C , Chini M , et al . Satellite radar and optical remote sensing for earthquake damage detection: results from different case studies
[J]. International Journal of Remote Sensing , 2006 , 27 (20 ): 4433 -4447 . DOI: 10.1080/01431160600675895 .
[本文引用: 2]
[41]
Chini M , Pierdicca N , Emery W J . Exploiting SAR and VHR optical images to quantify damage caused by the 2003 Bam Earthquake
[J]. IEEE Transactions on Geoscience and Remote Sensing ,2009 ,47 (1 ):145 -152 . DOI:10.1109/TGRS. 2008.2002695 .
[本文引用: 6]
[42]
Singh A . Review article digital change detection techniques using remotely-sensed data
[J]. International Journal of Remote Sensing ,1989 ,10 (6 ):989 -1003 . DOI:10.1080/01431168 908903939 .
[本文引用: 2]
[43]
Zhang J , Xie L , Tao X . Change detection of earthquake-damaged buildings on remote sensing image and its application in seismic disaster assessment
[C]∥ IGARSS 2003 2003 IEEE International Geoscience and Remote Sensing Symposium Proceedings (IEEE Cat No03CH37477) . Toulouse, France : IEEE , 2003 .
[本文引用: 3]
[44]
Rathje E M , Kyu-Seok W , Crawford M , et al . Earthquake damage identification using multi-temporal high-resolution optical satellite imagery
[C]∥ Proceedings 2005 IEEE International Geoscience and Remote Sensing Symposium , 2005 IGARSS ' 05 . Seoul, KoreaSouth, 2005 .
[本文引用: 4]
[45]
Wang X , Li P . Extraction of earthquake-induced collapsed buildings using very high-resolution imagery and airborne Lidar data
[J]. International Journal of Remote Sensing , 2015 , 36 (8 ): 2163 -2183 . DOI: 10.1080/01431161.2015.1034890 .
[46]
Liu J , Li P . Extraction of earthquake-induced collapsed buildings from Bi-Temporal VHR images using object-level Homogeneity Index and histogram
[J]. IEEE Journal of Selected To-pics in Applied Earth Observations and Remote Sensing , 2019 ,12 (8 ):2755 -2770 .DOI:10.1109/JSTARS.2019. 2904670 .
[本文引用: 2]
[47]
Guzzetti F , Mondini A C , Cardinali M , et al . Landslide inventory maps: New tools for an old problem
[J]. Earth-Science Reviews , 2012 , 112 (1 ): 42 -66 . DOI: 10.1016/j.earscirev.2012.02.001 .
[本文引用: 1]
[48]
Deng J S , Wang K , Deng Y H , et al . Pca‐based Land‐use change detection and analysis using multitemporal and Multisensor Satellite Data
[J]. International Journal of Remote Sensing ,2008 ,29 (16 ):4823 -4838 . DOI:10.1080/0143116080 1950162 .
[本文引用: 2]
[49]
Byrne G F , Crapper P F , Mayo K K . Monitoring land-cover change by principal component analysis of multitemporal landsat data
[J]. Remote Sensing of Environment , 1980 , 10 (3 ): 175 -184 . DOI: 10.1016/0034-4257(80)90021-8 .
[本文引用: 1]
[50]
Li Chenglong , Zhang Jingfa . Method and application of earthquake damage change detection from remote sensing images based on principal component analysis
[J]. Earthquake , 2013 , 33 (2 ): 105 -110 .
[本文引用: 1]
李成龙 , 张景发 . 基于主成分分析的遥感震害变化检测方法与应用
[J]. 地震 , 2013 , 33 (2 ): 105 -110 ].
[本文引用: 1]
[51]
Gong L , Li Q , Zhang J . Earthquake building damage detection with object-oriented change detection
[C]∥ 2013 IEEE International Geoscience and Remote Sensing Symposium-IGARSS: IEEE , 2013 .
[本文引用: 3]
[52]
Oommen T , Baise L G , Gens R , et al . Documenting earthquake-induced liquefaction using satellite remote sensing image transformations
[J]. Environmental and Engineering Geoscience , 2013 , 19 (4 ): 303 -318 . DOI: 10.2113/gseegeosci.19.4.303 .
[本文引用: 2]
[53]
Xu M , Cao C , Zhang H , et al . Change detection of an earthquake-induced Barrier Lake based on remote sensing image classification
[J].International Journal of Remote Sensing , 2010 ,31 (13 ):3521 -3534 . DOI:10.1080/01431161003727689 .
[本文引用: 2]
[54]
Niemeyer I , Marpu P R , Nussbaum S . Change detection using object features
[M]∥ Object-Based Image Analysis: Spatial Concepts for Knowledge-Driven Remote Sensing Applications . Berlin : Springer Berlin Heidelberg . 2008 : 185 -201 .
[本文引用: 1]
[55]
Matsumoto K , Vu T T , Yamazaki F . Extraction of damaged buildings using high resolution satellite images in the 2006 central java earthquake
[C]∥ Proc 27th Asian Conference on Remote Sensing . New Delhi, India , 2006 .
[本文引用: 3]
[56]
Johansen K , Arroyo L A , Phinn S , et al . Comparison of geo-object based and pixel-based change detection of riparian environments using high spatial resolution multi-spectral imagery
[J]. Photogrammetric Engineering Remote Sensing of Environment ,2010 ,76 (2 ): 123 -136 .
[本文引用: 1]
[58]
Li P , Xu H , Guo J . Urban building damage detection from very high resolution imagery using ocsvm and spatial features
[J]. International Journal of Remote Sensing ,2010 ,31 (13 ): 3393 -3409 . DOI: 10.1080/01431161003727705 .
[本文引用: 3]
[59]
Saito K , Spence R J S , Going C , et al . Using high-resolution satellite images for post-earthquake building damage assessment: A study following the 26 January 2001 Gujarat earthquake
[J]. Earthquake Spectra , 2004 , 20 (1 ): 145 -169 . DOI: 10.1193/1.1650865 .
[本文引用: 2]
[60]
Chini M , Cinti F R , Stramondo S . Co-seismic surface effects from very high resolution panchromatic images: The case of the 2005 Kashmir (Pakistan) earthquake
[J]. Natural Hazards and Earth System Sciences , 2011 , 11 (3 ): 931 -943 .
[本文引用: 1]
[61]
Bhattacharjee G , Soden R , Barns K , et al . Factors affecting earthquake responders’ building damage information needs and use
[J]. Earthquake Spectra , 2021 , 38 (1 ): 56 -80 . DOI: 10.1177/87552930211030297 .
[本文引用: 1]
[62]
Yamazaki F , Yano Y , Matsuoka M . Visual damage interpretation of buildings in bam city using QuickBird images following the 2003 bam,Iran, Earthquake
[J]. Earthquake Spectra , 2005 , 21 (): 329 -336 . DOI: 10.1193/1.2101807 .
[本文引用: 3]
[63]
Wu Xiaoying , Wang Xiaoqing , Deng Fei . Compilation and spatial analysis of co-seismic landslide inventory by using high-resolution remote sensing images in earthquake emergency response: An example of the Jiuzhaigou Ms
7 .0 Earthquake on August 8, 2017[J]. Technology for Earthquake Disaster Prevention , 2017, 12 (4 ): 815 -825 .
[本文引用: 1]
吴玮莹 , 王晓青 , 邓飞 . 基于高分卫星遥感影像的地震应急滑坡编目与分布特征探讨——以2017年8月8日九寨沟7.0级地震为例
[J]. 震灾防御技术 , 2017 , 12 (4 ): 815 -825 .
[本文引用: 1]
[64]
Tong X , Lin X , Feng T , et al . Use of shadows for detection of earthquake-induced collapsed buildings in high-resolution satellite imagery
[J]. ISPRS Journal of Photogrammetry and Remote Sensing , 2013 , 79 : 53 -67 . DOI: 10.1016/j.isprsjprs.2013.01.012 .
[本文引用: 3]
[65]
Kahraman F , Imamoglu M , Ates H F . Disaster damage assessment of buildings using adaptive self-similarity descriptor
[J]. IEEE Geoscience and Remote Sensing Letters , 2016 , 13 (8 ): 1188 -1192 . DOI: 10.1109/LGRS.2016.2574960 .
[本文引用: 1]
[66]
Ranjbar H R , Dehghani H , Ardalan A R A , et al . A GIS based approach for earthquake loss estimation based on the immediate extraction of damaged buildings
[J]. Geomatics, Natural Hazards and Risk , 2017 , 8 (2 ): 772 -791 . DOI: 10.1080/19475705.2016.1265013 .
[本文引用: 2]
[67]
Yamazaki F , Suzuki D , Maruyama Y . Detection of damages due to earthquakes using digital aerial images
[C]∥ 6th International Workshop on Remote Sensing for Disaster Applications: Citeseer , 2008 .
[本文引用: 1]
[68]
Zhu J , Baise L G , Koch M . Mapping earthquake induced liquefaction surface effects from the 2011 Tohoku Earthquake using satellite imagery
[C]∥ 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS) , 2016 .
[本文引用: 1]
[69]
Yang W , Wang M , Shi P . Using MODIS NDVI time series to identify geographic patterns of landslides in vegetated regions
[J]. IEEE Geoscience and Remote Sensing Letters , 2013 , 10 (4 ): 707 -710 . DOI: 10.1109/LGRS.2012.2219576 .
[本文引用: 1]
[70]
Zhang W , Lin J , Peng J , et al . Estimating Wenchuan Earthquake induced landslides based on remote sensing
[J]. International Journal of Remote Sensing , 2010 , 31 (13 ): 3495 -3508 . DOI: 10.1080/01431161003727630 .
[本文引用: 1]
[71]
Pierdicca N , Anniballe R , Noto F , et al . Triple collocation to assess classification accuracy without a ground truth in case of earthquake damage assessment
[J]. IEEE Transactions on Geoscience and Remote Sensing , 2018 , 56 (1 ): 485 -496 . DOI: 10.1109/TGRS.2017.2750770 .
[本文引用: 1]
[72]
Zheng Z , Zhong Y , Wang J , et al . Building damage assessment for rapid disaster response with a deep object-based semantic change detection framework: From natural disasters to man-made disasters
[J]. Remote Sensing of Environment , 2021 , 265 : 112636 . DOI: 10.1016/j.rse.2021.112636 .
[本文引用: 1]
[73]
Sun Xiaoxia , Zhang Jixian , Gao Yan , et al . A summary on current techniques and prospects of remote sensing change detection
[J]. Remote Sensing Information , 2011 ,26 (1 ): 119 -123 .
[本文引用: 1]
孙晓霞 , 张继贤 , 燕琴 , 等 . 遥感影像变化检测方法综述及展望
[J]. 遥感信息 , 2011 ,26 (1 ): 119 -123 .
[本文引用: 1]
[74]
Janalipour M , Taleai M . Building change detection after earthquake using multi-criteria decision analysis based on extracted information from high spatial resolution satellite images
[J]. International Journal of Remote Sensing , 2017 , 38 (1 ): 82 -99 . DOI: 10.1080/01431161.2016.1259673 .
[本文引用: 1]
关于地震人员伤亡因素的探讨
1
2000
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
关于地震人员伤亡因素的探讨
1
2000
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
震害信息遥感获取技术历史、现状和趋势
2
2004
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
... [2 ,4 ]. ...
震害信息遥感获取技术历史、现状和趋势
2
2004
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
... [2 ,4 ]. ...
汶川8.0级大地震应急遥感震害评估研究
2
2008
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
... 21世纪以来,遥感卫星的数量和时—空—谱分辨率不断提升,为地震应急提供了丰富的数据支撑.应用于地震灾情调查的遥感数据可分为3类,即光学影像、SAR和LiDAR[5 , 6 ] ,它们有着各自的特性和优缺点,可以从不同的角度反映灾情.其中光学影像具有分辨率高、直观、可靠等优点,便于目视解译,在震后灾情调查这样形势紧急、不容有误的情况下更能快速有效地获取初步且珍贵的灾情信息,因此光学影像(包括与其他数据的融合)一直是震后灾情调查中主要使用的数据[6 -8 ] .如王晓青等[3 ] 在汶川地震发生后,利用高分辨率光学遥感影像,对地震灾害及其损失进行了快速评估,弥补了地面调查的不足,为地震应急指挥和救援决策提供了重要依据. ...
汶川8.0级大地震应急遥感震害评估研究
2
2008
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
... 21世纪以来,遥感卫星的数量和时—空—谱分辨率不断提升,为地震应急提供了丰富的数据支撑.应用于地震灾情调查的遥感数据可分为3类,即光学影像、SAR和LiDAR[5 , 6 ] ,它们有着各自的特性和优缺点,可以从不同的角度反映灾情.其中光学影像具有分辨率高、直观、可靠等优点,便于目视解译,在震后灾情调查这样形势紧急、不容有误的情况下更能快速有效地获取初步且珍贵的灾情信息,因此光学影像(包括与其他数据的融合)一直是震后灾情调查中主要使用的数据[6 -8 ] .如王晓青等[3 ] 在汶川地震发生后,利用高分辨率光学遥感影像,对地震灾害及其损失进行了快速评估,弥补了地面调查的不足,为地震应急指挥和救援决策提供了重要依据. ...
准实时地震灾情综合评估系统的研发
1
2017
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
准实时地震灾情综合评估系统的研发
1
2017
... 大型的地震灾害因其瞬时突发性、强破坏性和难以预报性,常常造成大面积的人员伤亡和巨额的财物损失[1 ] .快速高效地实施地震应急,可以最大程度减少人员伤亡和经济损失,维护社会稳定.震后灾情信息的快速获取是应急救援的前提,否则灾情无从评估,决策指挥缺乏针对性和适用性.震后灾区的交通、通讯和电力等设施会遭到不同程度的损毁,传统的人工实地勘测虽然精确可靠,但存在着工作量大、工作范围小、时效性差、效率低、信息不直观等问题,难以满足应急救援的需求[2 -3 ] .作为一种先进、成熟的对地观测手段,卫星遥感技术可快速、动态地获取大范围的灾区影像,已成为震后灾情快速调查的主要技术手段之一[2 ,4 ] . ...
The role of remote sensing in earthquake science and engineering: Opportunities and challenges
2
2008
... 21世纪以来,遥感卫星的数量和时—空—谱分辨率不断提升,为地震应急提供了丰富的数据支撑.应用于地震灾情调查的遥感数据可分为3类,即光学影像、SAR和LiDAR[5 , 6 ] ,它们有着各自的特性和优缺点,可以从不同的角度反映灾情.其中光学影像具有分辨率高、直观、可靠等优点,便于目视解译,在震后灾情调查这样形势紧急、不容有误的情况下更能快速有效地获取初步且珍贵的灾情信息,因此光学影像(包括与其他数据的融合)一直是震后灾情调查中主要使用的数据[6 -8 ] .如王晓青等[3 ] 在汶川地震发生后,利用高分辨率光学遥感影像,对地震灾害及其损失进行了快速评估,弥补了地面调查的不足,为地震应急指挥和救援决策提供了重要依据. ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
A comprehensive review of earthquake-induced building damage detection with remote sensing techniques
5
2013
... 21世纪以来,遥感卫星的数量和时—空—谱分辨率不断提升,为地震应急提供了丰富的数据支撑.应用于地震灾情调查的遥感数据可分为3类,即光学影像、SAR和LiDAR[5 , 6 ] ,它们有着各自的特性和优缺点,可以从不同的角度反映灾情.其中光学影像具有分辨率高、直观、可靠等优点,便于目视解译,在震后灾情调查这样形势紧急、不容有误的情况下更能快速有效地获取初步且珍贵的灾情信息,因此光学影像(包括与其他数据的融合)一直是震后灾情调查中主要使用的数据[6 -8 ] .如王晓青等[3 ] 在汶川地震发生后,利用高分辨率光学遥感影像,对地震灾害及其损失进行了快速评估,弥补了地面调查的不足,为地震应急指挥和救援决策提供了重要依据. ...
... [6 -8 ].如王晓青等[3 ] 在汶川地震发生后,利用高分辨率光学遥感影像,对地震灾害及其损失进行了快速评估,弥补了地面调查的不足,为地震应急指挥和救援决策提供了重要依据. ...
... 依据是否采用震前的遥感影像,地震灾害信息可通过仅采用震后影像进行解译和同时采用震前和震后影像进行变化检测这两种方法进行提取[6 ,9 ] .震后地面的基础设施和土地覆盖相对于震前发生变化的部分就可能是地震造成损害的区域[7 ] ,因此通过对震前、震后的两幅影像的变化检测就可以提取出震害信息.由于震前影像可获取地物在震前的状态,因此变化检测的方法能够更准确地获取灾情信息,受到了更为广泛的应用[6 ] .但该方法在实际应用中可能会受到局限:首先震前影像可能无法获取或震前影像太过久远,尤其是在欠发达国家和地区;其次,震前和震后影像来自不同的时相、不同的传感器,光照条件、传感器特性等的差异会增大数据处理的难度[10 ] .所幸的是,随着卫星遥感数据的日益丰富和遥感图像处理技术的进步,这些局限都已得到适当缓解. ...
... [6 ].但该方法在实际应用中可能会受到局限:首先震前影像可能无法获取或震前影像太过久远,尤其是在欠发达国家和地区;其次,震前和震后影像来自不同的时相、不同的传感器,光照条件、传感器特性等的差异会增大数据处理的难度[10 ] .所幸的是,随着卫星遥感数据的日益丰富和遥感图像处理技术的进步,这些局限都已得到适当缓解. ...
... 表1 按全色分辨率由高到低的顺序列出了数据影像已广泛用于地震灾情调查的高分辨率光学遥感卫星[6 ,11 -12 ] .可以看出,空间分辨率优于1 m的卫星多达十数颗,其中WorldView系列卫星分辨率可达0.31 m.单星的时间分辨率已经可以达到1 d,而且目前的商业卫星均在大力构建卫星星座,获取数据的能力和数据质量显著提升,如美国的WorldView-2–4与GeoEye-1组成的4星星座可达到每日4.5次的重访频率,法国的Pléiades-1A/1B和SPOT-6/7组成的4星星座可实现每日2次的重访(SPOT-6/7提供大幅宽的普查图像,Pléiades-1A/B提供0.5 m分辨率的详查图像),中国的北京2号三星星座可实现每日2次的重访.此外,震后进行多星联合调度观测,可充分利用卫星资源、最大化灾情观测效益,极大地提高了震后应急响应的速度.在传感器配置上,这些卫星多为1个全色波段加4个或8个多光谱波段,波长范围覆盖了可见光到近红外波段(多为0.4~0.9 μm),像元的动态范围可达11 bit. ...
Identifying damage caused by the 2008 Wenchuan Earthquake from VHR remote sensing data
3
2009
... 依据是否采用震前的遥感影像,地震灾害信息可通过仅采用震后影像进行解译和同时采用震前和震后影像进行变化检测这两种方法进行提取[6 ,9 ] .震后地面的基础设施和土地覆盖相对于震前发生变化的部分就可能是地震造成损害的区域[7 ] ,因此通过对震前、震后的两幅影像的变化检测就可以提取出震害信息.由于震前影像可获取地物在震前的状态,因此变化检测的方法能够更准确地获取灾情信息,受到了更为广泛的应用[6 ] .但该方法在实际应用中可能会受到局限:首先震前影像可能无法获取或震前影像太过久远,尤其是在欠发达国家和地区;其次,震前和震后影像来自不同的时相、不同的传感器,光照条件、传感器特性等的差异会增大数据处理的难度[10 ] .所幸的是,随着卫星遥感数据的日益丰富和遥感图像处理技术的进步,这些局限都已得到适当缓解. ...
... 人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译.人工解译虽然费事费力,但精度和可靠性相对更高[58 ] .对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] . ...
... Saito等[59 ] 采用震后的IKONOS全色影像,并结合现场实地调查的照片,对目视解译的效果进行评价,表明了目视解译的快速性和对部分建筑损毁解译的局限性;Yamazaki等[62 ] 采用巴姆地震震前和震后的QuickBird影像,对12 063个建筑物的损毁程度逐个进行了目视解译;Ehrlich等[7 ] 通过震后映秀城区的WorldView影像、震前和震后都江堰城区的SPOT-5影像识别了滑坡灾害和建筑物的损毁情况;吴玮莹等[63 ] 利用震后获取的高分二号、北京二号影像和震前的Google Earth影像,通过人工目视解译的方法,以道路整齐边界遭受的破坏和山体植被覆盖的变化作为解译标志,初步建立了九寨沟地震的同震滑坡分布图.为便于目视解译,Zhang等[43 ] 在提取1999年南头地震的灾情信息时,将震前和震后共2幅灰度影像进行假彩色合成显示,在新合成的彩色图像中,未发生变化的区域为白色—灰色,发生变化区域表现为彩色,这些色彩就包含了震害信息,据此可判断震害分布及其性质. ...
Challenges and limitations of Earthquake-induced Building damage mapping techniques using remote sensing images-a systematic review
1
2021
... 21世纪以来,遥感卫星的数量和时—空—谱分辨率不断提升,为地震应急提供了丰富的数据支撑.应用于地震灾情调查的遥感数据可分为3类,即光学影像、SAR和LiDAR[5 , 6 ] ,它们有着各自的特性和优缺点,可以从不同的角度反映灾情.其中光学影像具有分辨率高、直观、可靠等优点,便于目视解译,在震后灾情调查这样形势紧急、不容有误的情况下更能快速有效地获取初步且珍贵的灾情信息,因此光学影像(包括与其他数据的融合)一直是震后灾情调查中主要使用的数据[6 -8 ] .如王晓青等[3 ] 在汶川地震发生后,利用高分辨率光学遥感影像,对地震灾害及其损失进行了快速评估,弥补了地面调查的不足,为地震应急指挥和救援决策提供了重要依据. ...
Remote sensing and earthquake damage assessment: Experiences, limits, and perspectives
3
2012
... 依据是否采用震前的遥感影像,地震灾害信息可通过仅采用震后影像进行解译和同时采用震前和震后影像进行变化检测这两种方法进行提取[6 ,9 ] .震后地面的基础设施和土地覆盖相对于震前发生变化的部分就可能是地震造成损害的区域[7 ] ,因此通过对震前、震后的两幅影像的变化检测就可以提取出震害信息.由于震前影像可获取地物在震前的状态,因此变化检测的方法能够更准确地获取灾情信息,受到了更为广泛的应用[6 ] .但该方法在实际应用中可能会受到局限:首先震前影像可能无法获取或震前影像太过久远,尤其是在欠发达国家和地区;其次,震前和震后影像来自不同的时相、不同的传感器,光照条件、传感器特性等的差异会增大数据处理的难度[10 ] .所幸的是,随着卫星遥感数据的日益丰富和遥感图像处理技术的进步,这些局限都已得到适当缓解. ...
... 地震发生后,基于获取的遥感影像,需要快速查明3类震害信息:倒塌的建筑物、损毁的生命线工程和次生地质灾害等[9 ] .其中建筑物震害是造成人员伤亡的主因,也是灾情调查的重点和难点;生命线工程包括道路、电线、水管、通讯等设施,关系到受灾人员的生命维持和现场第一手资料的获取;滑坡、泥石流、堰塞湖、土壤液化等次生地质灾害严重威胁灾区居民的安全和灾区的生态系统.这3类震害在发生的区域、尺度大小、形态特征和表现形式等方面存在着明显的差别,需要采用不同的数据和方法进行提取.本研究将不拘于这3种震害类别,而是着眼于震害信息提取的变化检测算法,对地震灾情调查的现状进行综述,以期为未来地震应急中的灾情提取工作提供参考和裨益. ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Damage assessment of 2010 Haiti Earthquake with post-earthquake satellite image by support vector selection and adaptation
1
2011
... 依据是否采用震前的遥感影像,地震灾害信息可通过仅采用震后影像进行解译和同时采用震前和震后影像进行变化检测这两种方法进行提取[6 ,9 ] .震后地面的基础设施和土地覆盖相对于震前发生变化的部分就可能是地震造成损害的区域[7 ] ,因此通过对震前、震后的两幅影像的变化检测就可以提取出震害信息.由于震前影像可获取地物在震前的状态,因此变化检测的方法能够更准确地获取灾情信息,受到了更为广泛的应用[6 ] .但该方法在实际应用中可能会受到局限:首先震前影像可能无法获取或震前影像太过久远,尤其是在欠发达国家和地区;其次,震前和震后影像来自不同的时相、不同的传感器,光照条件、传感器特性等的差异会增大数据处理的难度[10 ] .所幸的是,随着卫星遥感数据的日益丰富和遥感图像处理技术的进步,这些局限都已得到适当缓解. ...
A review of the status of satellite remote sensing and image processing techniques for mapping natural hazards and disasters
1
2009
... 表1 按全色分辨率由高到低的顺序列出了数据影像已广泛用于地震灾情调查的高分辨率光学遥感卫星[6 ,11 -12 ] .可以看出,空间分辨率优于1 m的卫星多达十数颗,其中WorldView系列卫星分辨率可达0.31 m.单星的时间分辨率已经可以达到1 d,而且目前的商业卫星均在大力构建卫星星座,获取数据的能力和数据质量显著提升,如美国的WorldView-2–4与GeoEye-1组成的4星星座可达到每日4.5次的重访频率,法国的Pléiades-1A/1B和SPOT-6/7组成的4星星座可实现每日2次的重访(SPOT-6/7提供大幅宽的普查图像,Pléiades-1A/B提供0.5 m分辨率的详查图像),中国的北京2号三星星座可实现每日2次的重访.此外,震后进行多星联合调度观测,可充分利用卫星资源、最大化灾情观测效益,极大地提高了震后应急响应的速度.在传感器配置上,这些卫星多为1个全色波段加4个或8个多光谱波段,波长范围覆盖了可见光到近红外波段(多为0.4~0.9 μm),像元的动态范围可达11 bit. ...
Remote sensing contributing to assess earthquake risk: from a literature review towards a roadmap
1
2013
... 表1 按全色分辨率由高到低的顺序列出了数据影像已广泛用于地震灾情调查的高分辨率光学遥感卫星[6 ,11 -12 ] .可以看出,空间分辨率优于1 m的卫星多达十数颗,其中WorldView系列卫星分辨率可达0.31 m.单星的时间分辨率已经可以达到1 d,而且目前的商业卫星均在大力构建卫星星座,获取数据的能力和数据质量显著提升,如美国的WorldView-2–4与GeoEye-1组成的4星星座可达到每日4.5次的重访频率,法国的Pléiades-1A/1B和SPOT-6/7组成的4星星座可实现每日2次的重访(SPOT-6/7提供大幅宽的普查图像,Pléiades-1A/B提供0.5 m分辨率的详查图像),中国的北京2号三星星座可实现每日2次的重访.此外,震后进行多星联合调度观测,可充分利用卫星资源、最大化灾情观测效益,极大地提高了震后应急响应的速度.在传感器配置上,这些卫星多为1个全色波段加4个或8个多光谱波段,波长范围覆盖了可见光到近红外波段(多为0.4~0.9 μm),像元的动态范围可达11 bit. ...
The international charter “Space and Major Disasters” initiative
1
2004
... 1999年法国航天局国家空间研究中心(CNES)和欧洲航天局(ESA)倡议发起了“空间与重大灾害国际宪章”(The International Charter Space and Major Disasters),加入该宪章的成员可以无偿获取其他成员的卫星数据用于应急救灾,解决了卫星资源的共享问题[13 ] ,对遥感技术在地震应急中的应用有着深远影响. ...
1
2008
... 使用卫星影像进行灾情调查,需首先进行预处理,包括辐射校正、大气校正、几何校正、图像融合等.表2 列出了美国DigitalGlobe、欧洲Astrium和中国二十一世纪空间技术应用股份有限公司(21at)旗下的遥感影像产品[14 -16 ] .可以看出用户拿到的数据已经过不同程度的预处理,因此用户可以根据实际需求选择合适等级的数据产品,无需考虑数据预处理的问题.此外,数据服务商可为用户提供加急处理选项(Rush option)以确保震后数据的时效性. ...
Pléiades imagery user guide
0
Digitalglobe core imagery products guide
1
2013
... 使用卫星影像进行灾情调查,需首先进行预处理,包括辐射校正、大气校正、几何校正、图像融合等.表2 列出了美国DigitalGlobe、欧洲Astrium和中国二十一世纪空间技术应用股份有限公司(21at)旗下的遥感影像产品[14 -16 ] .可以看出用户拿到的数据已经过不同程度的预处理,因此用户可以根据实际需求选择合适等级的数据产品,无需考虑数据预处理的问题.此外,数据服务商可为用户提供加急处理选项(Rush option)以确保震后数据的时效性. ...
Change detection techniques
3
2004
... 变化检测的基本流程主要包括影像预处理、变化检测和精度评价[17 -18 ] . ...
... 完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情.对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类.最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述. ...
... 主成分变换是将具有相关性的多波段影像压缩到完全独立的较少的几个主成分上,使新影像更便于解译.主成分变换法进行变化检测原理是,未变化的区域具有较高的相关性,而发生变化的区域相关性相对较低[48 ] .在提取变化信息时,主成分变换法常与差值法结合使用,依据差值和主成分变换的先后可分为差值主成分法、主成分差值法和多波段主成分变换法[48 -49 ] .主成分变换的缺点是变换后各分量所突出的信息有所不同且物理意义并不明确,具体哪个分量可用来表征变化信息需要一定经验和分析[17 ] . ...
Change detection from remotely sensed images: from pixel-based to object-based approaches
9
2013
... 变化检测的基本流程主要包括影像预处理、变化检测和精度评价[17 -18 ] . ...
... 在地震应急情况下,震前、震后的影像可能来自不同的传感器,获取影像的条件(如太阳高度角、光照条件和大气状况等)往往存在着差异,需要消除这些因素造成的虚假变化.几何配准和相对辐射校正是影响变化检测精度的关键步骤[18 ] . ...
... 完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情.对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类.最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述. ...
... 分类后比较法是对震前和震后影像中的像元进行独立分类,然后比较分类结果,如果对应的类别相同,则判定没有发生变化,否则认为发生了变化.分类后比较法在一定程度上借鉴了面向对象变化检测方法的思想,只对感兴趣的对象进行比较,可以直接获得变化的位置、类型和数量等信息,而且通过对震前和震后影像的单独分类,克服了影像由于获取条件、分辨率等因素造成的差异[18 ] . ...
... 然而,分类后比较法的缺点限制了其在震害信息提取中的应用:①对每幅影像单独分类,增加了计算量和复杂度,影响时效性;②两次分类的误差会传递到变化检测的结果中而被放大[18 , 42 ] .此外,震害信息提取中,变化检测的关注点主要在于“是否”发生了变化或发生变化的程度,至于具体类型的变化(from-to)并非必须获取的信息. ...
... 基于像元的变化检测算法在应用于高分辨率影像时常常会受到局限.首先,分辨率越高,像元越小,空间不一致或者几何配准误差等造成的影响越加明显;其次,高分辨率影像在承载更丰富的地表信息的同时,也导致了同一地物内部像元灰度起伏明显,在使用基于像元的方法时,过多无关的变化会被检测出来,产生“椒盐效应”[18 , 54 ] . ...
... 基于Hussain等[18 ] 的分类方法,结合地震灾情调查中的实际需求和应用情况,本研究将地震中面向对象的变化检测方法分为了3类:目视解译法、直接比较法、多时相叠加法. ...
... 多时相叠加法是指将震前和震后影像叠加在一起,只对震前的影像进行对象的分割和分类,将分类后对象的掩膜直接应用到所有的时相中,仅对掩膜所包含的区域进行比较.多时相叠加法的优点在于对多个时相的影像只需进行一次对象分割和分类,且同一对象具有相同的几何特征[18 ] ,但该方法需要进行影像的几何配准. ...
... 为了系统地综述变化检测算法在地震灾情调查中的应用,本研究对这些算法的类别进行了较为细致的划分,但这种类别的划分并非唯一的[18 -19 ,35 ,73 ] ,不同的类别的界限不是严格分明的,可能存在着交叉、借鉴或套用.在实际研究和应用这些方法时,也不会拘泥于算法的类别,如Janalipour等[74 ] 综合采用了差值法、变化向量法、植被指数法、主成分分析法等多种方法,构建了多准则判据来判断建筑物是否损毁,使检测结果更加稳定、精确. ...
A review of multi-temporal remote sensing data change detection algorithms
3
2008
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... 完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情.对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类.最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述. ...
... 为了系统地综述变化检测算法在地震灾情调查中的应用,本研究对这些算法的类别进行了较为细致的划分,但这种类别的划分并非唯一的[18 -19 ,35 ,73 ] ,不同的类别的界限不是严格分明的,可能存在着交叉、借鉴或套用.在实际研究和应用这些方法时,也不会拘泥于算法的类别,如Janalipour等[74 ] 综合采用了差值法、变化向量法、植被指数法、主成分分析法等多种方法,构建了多准则判据来判断建筑物是否损毁,使检测结果更加稳定、精确. ...
Object-oriented image understanding and post-earthquake damage assessment for the 2003 Bam, Iran, Earthquake
3
2005
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
... Gusella等[20 ] ...
Towards rapid citywide damage mapping using neighborhood edge dissimilarities in very high-resolution optical satellite imagery—application to the 2003 Bam, Iran, Earthquake
0
2005
Damage detection for 2003 Bam, Iran, Earthquake using Terra-Aster satellite imagery
3
2005
... 地震应急情况下,通常不容易获取震前多个时相的数据,因此几乎没有回归法应用于地震灾情的案例.Kohiyama和Yamazaki[22 ] 在提取2003年巴姆地震的建筑物震害时,提出了一种类似于回归法的方法——图像波动模型法(Image Fluctuation Model,IMF),认为不同时刻同一像元的DN值由于传感器、大气、光照等差异而产生微小波动,利用震前的17幅影像获得其波动模型,计算震后影像的置信水平,从而得到建筑物发生损毁的可能性.回归法的优点在于采用了大量的震前影像对变化的趋势进行拟合或假设,其变化检测的结果较为可靠,但对数据要求过高,且每个时相的影像均要进行配准,计算量大,因此在地震应急中并不常用. ...
... 在地震应急中,通常并不关注植被的变化情况,甚至需要剔除植被的影响,如Yamazaki等[22 ,44 -46 ] 在提取城市的建筑物震害信息时,对植被进行了掩膜处理,以消除虚假的变化和干扰;Chini等[41 ] 发现若不对植被进行甄别,20%~50%未受损的建筑物会被误判为“倒塌”.但在某些情况下可通过植被的变化提取需要的灾情信息,如热带地区或山区植被的明显减少表明可能发生了山体滑坡[47 ] .Yang等[33 ] 使用震前的Landsat 5 TM和震后的Terra ASTER影像通过波段运算得到NDVI,通过对比试验表明了相对辐射校正的必要性和植被指数法对滑坡检测的有效性. ...
... Zhao等
[34 ] 回归法 稳定可靠 需要多个时相的震前影像,需要分别配准,通常难以实现 Kohiyama等[22 ] 相关系数法 可以消除地震前后影像的灰度和对比度的差异;可用于配准 需确定窗口大小和阈值;配准误差会降低相关性 Rathje等[44 ] ...
Earthquake damage mapping by using remotely sensed data: The Haiti case study
2
2017
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... 差值法的思想来自于对“变化”二字的直观理解,是变化检测算法的基础,其他的方法都或多或少地用到了差值法.由于差值法太过简单,仅采用差值法提取震害信息的文献很少,如Yusuf等[38 ] 在2001年印度古吉拉特邦地震后,采用差值法计算出了震前、震后两幅Landsat-7 ETM+影像的差值影像,设定阈值提取了建筑物倒塌的区域和土壤液化区域;Ishii等[39 ] 通过震前和震后航空影像的色调、亮度的差值提取出了建筑物的损毁区域.归一化差值是对影像的差值进行归一化,凸显了相对变化量,在震害提取中应用较为广泛,如Romaniello等[23 ,40 -41 ] 人的研究.变化向量法是在差值法的基础上扩展而来的,通过震前、震后影像的所有波段的差值构建变化向量,能够避免单一波段比较带来的信息不完整,如Zhao等[34 ] 在利用Landsat 8影像提取尼泊尔震后的滑坡时,采用7个波段的反射率和3种植被指数(NDVI、EVI、SAVI)构建了两组变化向量. ...
高分遥感影像震害提取的变化检测算法改进
2
2013
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... Zhang等[43 ] 对1999年台湾南投地震的震前震后影像,通过设定相关系数的阈值范围,将建筑物的损毁程度划分成了4个等级;Rathje等[44 ] 采用了相关系数法提取巴姆城区震后的建筑物震害,并发现灰度共生矩阵(Gray Level of Co-occurrence Matrix,GLCM)的相关系数比DN值的相关系数能更好地区分损毁和完好的房屋;Chesnel等[26 ] 以屋顶范围作为窗口计算相关系数,实现了震前震后屋顶矢量的精确配准;张伟等[24 ] 为了减小配准误差对相关系数的影响,通过在对应点周围搜索相关系数最大的潜在同名点,使得倒塌和未倒塌建筑的分类精度(κ 系数)由0.60提高到了0.85.可见,相关系数不仅可以作为变化检测的特征变量,也可以作为几何配准的相似测度. ...
高分遥感影像震害提取的变化检测算法改进
2
2013
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... Zhang等[43 ] 对1999年台湾南投地震的震前震后影像,通过设定相关系数的阈值范围,将建筑物的损毁程度划分成了4个等级;Rathje等[44 ] 采用了相关系数法提取巴姆城区震后的建筑物震害,并发现灰度共生矩阵(Gray Level of Co-occurrence Matrix,GLCM)的相关系数比DN值的相关系数能更好地区分损毁和完好的房屋;Chesnel等[26 ] 以屋顶范围作为窗口计算相关系数,实现了震前震后屋顶矢量的精确配准;张伟等[24 ] 为了减小配准误差对相关系数的影响,通过在对应点周围搜索相关系数最大的潜在同名点,使得倒塌和未倒塌建筑的分类精度(κ 系数)由0.60提高到了0.85.可见,相关系数不仅可以作为变化检测的特征变量,也可以作为几何配准的相似测度. ...
Image change detection on urban area: The earthquake case
2
2004
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... 面向对象的变化检测方法能很好地克服基于像元的方法在配准和噪声等方面的局限[55 ] ,通过影像中地物对象的提取,直接比较这些地物的光谱特征、纹理特征、形状特征和上下文关系等的变化,更加高效可靠[56 ] ,此外还能减少光照条件、阴影等差异造成的影响[57 ] .面向对象方法的关键在于对象的提取(包括分割和分类),而对变化信息的提取基本上都是借鉴了基于像元的方法,所不同的只是将像元间的比较变成了对象间的比较[57 ] .Bitelli等[25 ,55 ,58 ] 对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害. ...
Object oriented assessment of damage due to natural disaster using very high resolution images
3
2007
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... Zhang等[43 ] 对1999年台湾南投地震的震前震后影像,通过设定相关系数的阈值范围,将建筑物的损毁程度划分成了4个等级;Rathje等[44 ] 采用了相关系数法提取巴姆城区震后的建筑物震害,并发现灰度共生矩阵(Gray Level of Co-occurrence Matrix,GLCM)的相关系数比DN值的相关系数能更好地区分损毁和完好的房屋;Chesnel等[26 ] 以屋顶范围作为窗口计算相关系数,实现了震前震后屋顶矢量的精确配准;张伟等[24 ] 为了减小配准误差对相关系数的影响,通过在对应点周围搜索相关系数最大的潜在同名点,使得倒塌和未倒塌建筑的分类精度(κ 系数)由0.60提高到了0.85.可见,相关系数不仅可以作为变化检测的特征变量,也可以作为几何配准的相似测度. ...
... Chesnel等[26 ] ...
Automated damage indication for rapid geospatial reporting
4
2011
... 几何配准的关键在于建立两幅影像对应像元坐标转换的映射函数.几何配准误差是变化检测主要的误差来源之一,因此基于像元的变化检测方法需达到亚像元级的配准精度[19 ] .Gusella等[20 -23 ] 均是手动选取控制点实现了震前、震后影像的配准,这也是最常用的配准方法;张伟等[24 ] 采用相关系数法同时完成了影像的配准和变化检测.面向对象的变化检测方法对配准精度的要求较为宽松,只需要保证震前和震后影像中同一对象的对应性.Bitelli等[25 ] 对比了基于像元和面向对象的分类方法在震后建筑物震害识别中的效果,指出在同样的配准条件下面向对象方法精度更高、更稳定;Chesnel等[26 ] 通过计算相关系数实现了屋顶矢量的精确配准,显著提高了建筑物震害的提取精度;Tiede等[27 ] 通过在设定的虚拟的缓冲区域内搜索震前和震后影像中形状、大小相近的建筑物阴影实现了同一对象的“链接”. ...
... 人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译.人工解译虽然费事费力,但精度和可靠性相对更高[58 ] .对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] . ...
... Matsumoto等[55 ,67 ] 对震前、震后的影像进行分类,提取出两个时相的建筑物,通过建筑物的减少来判别损毁.类似地,Gong等[51 ] 根据房屋的长宽比、纹理、尺度等特征建立分类的规则集提取出了震前和震后影像中的房屋对象,认为震后相比于震前减少的房屋对象即为被损毁的房屋.Tiede等[27 ] 在未对震前、震后影像配准的情况下,通过对比对应建筑物阴影的面积和形状的变化来判断建筑物是否损毁;Zhu等[68 ] 对震后影像进行分割,通过决策树分类提取出人造地表(如道路、停车场)和裸土,采用多波段主成分变换法判断这些地表是否发生液化.Anniballe等[32 ] 借助于拉奎拉市已有的GIS底图,利用建筑物高度和传感器姿态等数据精确计算了震前和震后影像中建筑物的矢量多边形,测试了4类共13种影像特征(光谱特征、信息论特征、纹理特征、HSV颜色特征)和2种分类方法(最大后验概率和支持向量机)对建筑物震害的提取效果. ...
... Tiede等[27 ] ...
Relative radiometric normalization performance for change detection from multi-date satellite images
1
2000
... 相对辐射校正是指通过比较和修改震前、震后影像的灰度统计量,使其具有趋于一致的光谱特征,以减少影像由于接收条件的不同而造成的灰度与对比度的差异[28 , 29 ] .地震灾情调查中常用的相对辐射校正方法有直方图匹配法[29 -32 ] 、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
建筑物震害遥感图像的变化检测与震害评估
3
2002
... 相对辐射校正是指通过比较和修改震前、震后影像的灰度统计量,使其具有趋于一致的光谱特征,以减少影像由于接收条件的不同而造成的灰度与对比度的差异[28 , 29 ] .地震灾情调查中常用的相对辐射校正方法有直方图匹配法[29 -32 ] 、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
... [29 -32 ]、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
... 其中:i = 1,2 , ⋯ , m ; j = 1,2 , ⋯ , n ;I ¯ p r e 和I ¯ p o s t 分别表示两个时相窗口内像元的均值.相关系数越大,表示震前震后影像的变化越小,即损毁程度越小.相关系数不随影像灰度的线性变化而变化,可以消除地震前后由于影像接收条件的不同所造成的灰度与对比度的差异[29 ] . ...
建筑物震害遥感图像的变化检测与震害评估
3
2002
... 相对辐射校正是指通过比较和修改震前、震后影像的灰度统计量,使其具有趋于一致的光谱特征,以减少影像由于接收条件的不同而造成的灰度与对比度的差异[28 , 29 ] .地震灾情调查中常用的相对辐射校正方法有直方图匹配法[29 -32 ] 、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
... [29 -32 ]、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
... 其中:i = 1,2 , ⋯ , m ; j = 1,2 , ⋯ , n ;I ¯ p r e 和I ¯ p o s t 分别表示两个时相窗口内像元的均值.相关系数越大,表示震前震后影像的变化越小,即损毁程度越小.相关系数不随影像灰度的线性变化而变化,可以消除地震前后由于影像接收条件的不同所造成的灰度与对比度的差异[29 ] . ...
Automatic detection and classification of damaged buildings, using high resolution satellite imagery and vector data
3
2008
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
... Samadzadegan等[30 ] ...
... GIS底图可以极大地便利震害信息提取,降低问题的复杂度,减少计算量和工作量,提高效率和精度,如Samadzadegan和Rastiveisi等[30 ,32 ,64 -66 ] 直接采用了GIS底图提供的建筑物矢量,省去了对象分割和建筑物分类等工作.因此,在地震危险度高的城市应大力开展城市GIS底图的绘制和更新工作. ...
Colour and texture based change detection for urban disaster analysis
2
2011
... Tomowski等[31 ] 以2010年海地地震后太子港总统府区域为研究对象,从GLCM能量值第一主成分的差值中提取了损毁的建筑、植被和新增的救灾帐篷等变化信息;李成龙等[50 ] 采用主成分差值法对2003年阿尔及利亚地震的建筑物损毁情况进行检测,然而检测结果中存在很多零散杂乱的图斑,尤其受建筑物阴影变化的影响很大,需要较多的人工干预;Gong等[51 ] 对比了主成分分析法与差值法、面向对象法的震害提取效果,发现主成分分析法并无优势,提取的变化多为建筑物阴影的变化而非建筑物的损毁.此外,穗帽变换作为一种特殊的主成分变换方法,可以得到具有物理意义的主成分,也在震害提取中得到应用,如Oommen等[52 ] 基于2001印度古吉拉特邦地震震前、震后的Landsat 7 ETM+影像,利用穗帽变换得到的湿度分量,探测出了灾区土壤液化的分布情况. ...
... Tomowski等[31 ] ...
Earthquake damage mapping: An overall assessment of ground surveys and vhr image change detection after L'aquila 2009 earthquake
6
2018
... 相对辐射校正是指通过比较和修改震前、震后影像的灰度统计量,使其具有趋于一致的光谱特征,以减少影像由于接收条件的不同而造成的灰度与对比度的差异[28 , 29 ] .地震灾情调查中常用的相对辐射校正方法有直方图匹配法[29 -32 ] 、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
... Matsumoto等[55 ,67 ] 对震前、震后的影像进行分类,提取出两个时相的建筑物,通过建筑物的减少来判别损毁.类似地,Gong等[51 ] 根据房屋的长宽比、纹理、尺度等特征建立分类的规则集提取出了震前和震后影像中的房屋对象,认为震后相比于震前减少的房屋对象即为被损毁的房屋.Tiede等[27 ] 在未对震前、震后影像配准的情况下,通过对比对应建筑物阴影的面积和形状的变化来判断建筑物是否损毁;Zhu等[68 ] 对震后影像进行分割,通过决策树分类提取出人造地表(如道路、停车场)和裸土,采用多波段主成分变换法判断这些地表是否发生液化.Anniballe等[32 ] 借助于拉奎拉市已有的GIS底图,利用建筑物高度和传感器姿态等数据精确计算了震前和震后影像中建筑物的矢量多边形,测试了4类共13种影像特征(光谱特征、信息论特征、纹理特征、HSV颜色特征)和2种分类方法(最大后验概率和支持向量机)对建筑物震害的提取效果. ...
... Anniballe等[32 ] ...
... GIS底图可以极大地便利震害信息提取,降低问题的复杂度,减少计算量和工作量,提高效率和精度,如Samadzadegan和Rastiveisi等[30 ,32 ,64 -66 ] 直接采用了GIS底图提供的建筑物矢量,省去了对象分割和建筑物分类等工作.因此,在地震危险度高的城市应大力开展城市GIS底图的绘制和更新工作. ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
... [32 ]对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Using multi-temporal remote sensor imagery to detect earthquake-triggered landslides
4
2010
... 相对辐射校正是指通过比较和修改震前、震后影像的灰度统计量,使其具有趋于一致的光谱特征,以减少影像由于接收条件的不同而造成的灰度与对比度的差异[28 , 29 ] .地震灾情调查中常用的相对辐射校正方法有直方图匹配法[29 -32 ] 、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
... 在地震应急中,通常并不关注植被的变化情况,甚至需要剔除植被的影响,如Yamazaki等[22 ,44 -46 ] 在提取城市的建筑物震害信息时,对植被进行了掩膜处理,以消除虚假的变化和干扰;Chini等[41 ] 发现若不对植被进行甄别,20%~50%未受损的建筑物会被误判为“倒塌”.但在某些情况下可通过植被的变化提取需要的灾情信息,如热带地区或山区植被的明显减少表明可能发生了山体滑坡[47 ] .Yang等[33 ] 使用震前的Landsat 5 TM和震后的Terra ASTER影像通过波段运算得到NDVI,通过对比试验表明了相对辐射校正的必要性和植被指数法对滑坡检测的有效性. ...
... Yang等[33 ] ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Postearthquake landslides mapping from Landsat-8 data for the 2015 Nepal earthquake using a pixel-based change eetection method
3
2017
... 相对辐射校正是指通过比较和修改震前、震后影像的灰度统计量,使其具有趋于一致的光谱特征,以减少影像由于接收条件的不同而造成的灰度与对比度的差异[28 , 29 ] .地震灾情调查中常用的相对辐射校正方法有直方图匹配法[29 -32 ] 、暗集—亮集法[33 ] 、线性回归法[34 ] 等. ...
... 差值法的思想来自于对“变化”二字的直观理解,是变化检测算法的基础,其他的方法都或多或少地用到了差值法.由于差值法太过简单,仅采用差值法提取震害信息的文献很少,如Yusuf等[38 ] 在2001年印度古吉拉特邦地震后,采用差值法计算出了震前、震后两幅Landsat-7 ETM+影像的差值影像,设定阈值提取了建筑物倒塌的区域和土壤液化区域;Ishii等[39 ] 通过震前和震后航空影像的色调、亮度的差值提取出了建筑物的损毁区域.归一化差值是对影像的差值进行归一化,凸显了相对变化量,在震害提取中应用较为广泛,如Romaniello等[23 ,40 -41 ] 人的研究.变化向量法是在差值法的基础上扩展而来的,通过震前、震后影像的所有波段的差值构建变化向量,能够避免单一波段比较带来的信息不完整,如Zhao等[34 ] 在利用Landsat 8影像提取尼泊尔震后的滑坡时,采用7个波段的反射率和3种植被指数(NDVI、EVI、SAVI)构建了两组变化向量. ...
... Zhao等[34 ] ...
A critical synthesis of remotely sensed optical image change detection techniques
2
2015
... 完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情.对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类.最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述. ...
... 为了系统地综述变化检测算法在地震灾情调查中的应用,本研究对这些算法的类别进行了较为细致的划分,但这种类别的划分并非唯一的[18 -19 ,35 ,73 ] ,不同的类别的界限不是严格分明的,可能存在着交叉、借鉴或套用.在实际研究和应用这些方法时,也不会拘泥于算法的类别,如Janalipour等[74 ] 综合采用了差值法、变化向量法、植被指数法、主成分分析法等多种方法,构建了多准则判据来判断建筑物是否损毁,使检测结果更加稳定、精确. ...
多时相遥感影像变化检测的现状与展望
1
2017
... 完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情.对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类.最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述. ...
多时相遥感影像变化检测的现状与展望
1
2017
... 完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情.对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类.最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述. ...
Change detection in urban areas from remote sensing data: A multidimensional classification scheme
1
2019
... 完成了数据预处理工作后,就可以选择合适的变化检测方法对震前和震后的影像进行比较,以确定发生变化的区域和变化量,进而评估灾情.对于变化检测的方法,不同的学者有着不同的分类方式,如Lu等[17 ] 将变化检测方法分为7类:代数法、变换法、分类法、高级模型法、GIS法、目视解译法和其他方法等;Gong等[19 ] 分成了直接比较法、分类后比较法、统一模型法和时间序列分析法4大类;Tewkesbury等[35 ] 基于分析的单元对象的不同,将变化检测方法分成了7类;张良培和武辰[36 ] 根据采用的数据将变化检测分为中低分辨率影像变化检测、高光谱影像变化检测和高分辨率影像变化检测3大类;Si Salah等[37 ] 从变化检测的目的、输入数据、时间分辨率、分析单元、输出单元、建筑物特征、处理方法、变化类别和精度评估等多个维度上将建筑物的变化检测方法进行了复杂而完备的分类.最常用的分类是根据基本处理单元的不同将变化检测方法分为两大类,即基于像元的变化检测方法和面向对象的变化检测方法[18 ] ,然后再依据算法原理进行更细类别的划分,本研究将采用这种分类方式对地震灾情调查中的变化检测方法进行综述. ...
Damage assessment after 2001 gujarat earthquake using Landsat 7 satellite images
2
2001
... 差值法的思想来自于对“变化”二字的直观理解,是变化检测算法的基础,其他的方法都或多或少地用到了差值法.由于差值法太过简单,仅采用差值法提取震害信息的文献很少,如Yusuf等[38 ] 在2001年印度古吉拉特邦地震后,采用差值法计算出了震前、震后两幅Landsat-7 ETM+影像的差值影像,设定阈值提取了建筑物倒塌的区域和土壤液化区域;Ishii等[39 ] 通过震前和震后航空影像的色调、亮度的差值提取出了建筑物的损毁区域.归一化差值是对影像的差值进行归一化,凸显了相对变化量,在震害提取中应用较为广泛,如Romaniello等[23 ,40 -41 ] 人的研究.变化向量法是在差值法的基础上扩展而来的,通过震前、震后影像的所有波段的差值构建变化向量,能够避免单一波段比较带来的信息不完整,如Zhao等[34 ] 在利用Landsat 8影像提取尼泊尔震后的滑坡时,采用7个波段的反射率和3种植被指数(NDVI、EVI、SAVI)构建了两组变化向量. ...
... Yusuf等[38 ] ...
Detection of earthquake damaged areas from aerial photographs by using color and edge information
1
2002
... 差值法的思想来自于对“变化”二字的直观理解,是变化检测算法的基础,其他的方法都或多或少地用到了差值法.由于差值法太过简单,仅采用差值法提取震害信息的文献很少,如Yusuf等[38 ] 在2001年印度古吉拉特邦地震后,采用差值法计算出了震前、震后两幅Landsat-7 ETM+影像的差值影像,设定阈值提取了建筑物倒塌的区域和土壤液化区域;Ishii等[39 ] 通过震前和震后航空影像的色调、亮度的差值提取出了建筑物的损毁区域.归一化差值是对影像的差值进行归一化,凸显了相对变化量,在震害提取中应用较为广泛,如Romaniello等[23 ,40 -41 ] 人的研究.变化向量法是在差值法的基础上扩展而来的,通过震前、震后影像的所有波段的差值构建变化向量,能够避免单一波段比较带来的信息不完整,如Zhao等[34 ] 在利用Landsat 8影像提取尼泊尔震后的滑坡时,采用7个波段的反射率和3种植被指数(NDVI、EVI、SAVI)构建了两组变化向量. ...
Satellite radar and optical remote sensing for earthquake damage detection: results from different case studies
2
2006
... 差值法的思想来自于对“变化”二字的直观理解,是变化检测算法的基础,其他的方法都或多或少地用到了差值法.由于差值法太过简单,仅采用差值法提取震害信息的文献很少,如Yusuf等[38 ] 在2001年印度古吉拉特邦地震后,采用差值法计算出了震前、震后两幅Landsat-7 ETM+影像的差值影像,设定阈值提取了建筑物倒塌的区域和土壤液化区域;Ishii等[39 ] 通过震前和震后航空影像的色调、亮度的差值提取出了建筑物的损毁区域.归一化差值是对影像的差值进行归一化,凸显了相对变化量,在震害提取中应用较为广泛,如Romaniello等[23 ,40 -41 ] 人的研究.变化向量法是在差值法的基础上扩展而来的,通过震前、震后影像的所有波段的差值构建变化向量,能够避免单一波段比较带来的信息不完整,如Zhao等[34 ] 在利用Landsat 8影像提取尼泊尔震后的滑坡时,采用7个波段的反射率和3种植被指数(NDVI、EVI、SAVI)构建了两组变化向量. ...
... 人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译.人工解译虽然费事费力,但精度和可靠性相对更高[58 ] .对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] . ...
Exploiting SAR and VHR optical images to quantify damage caused by the 2003 Bam Earthquake
6
2009
... 差值法的思想来自于对“变化”二字的直观理解,是变化检测算法的基础,其他的方法都或多或少地用到了差值法.由于差值法太过简单,仅采用差值法提取震害信息的文献很少,如Yusuf等[38 ] 在2001年印度古吉拉特邦地震后,采用差值法计算出了震前、震后两幅Landsat-7 ETM+影像的差值影像,设定阈值提取了建筑物倒塌的区域和土壤液化区域;Ishii等[39 ] 通过震前和震后航空影像的色调、亮度的差值提取出了建筑物的损毁区域.归一化差值是对影像的差值进行归一化,凸显了相对变化量,在震害提取中应用较为广泛,如Romaniello等[23 ,40 -41 ] 人的研究.变化向量法是在差值法的基础上扩展而来的,通过震前、震后影像的所有波段的差值构建变化向量,能够避免单一波段比较带来的信息不完整,如Zhao等[34 ] 在利用Landsat 8影像提取尼泊尔震后的滑坡时,采用7个波段的反射率和3种植被指数(NDVI、EVI、SAVI)构建了两组变化向量. ...
... 在地震应急中,通常并不关注植被的变化情况,甚至需要剔除植被的影响,如Yamazaki等[22 ,44 -46 ] 在提取城市的建筑物震害信息时,对植被进行了掩膜处理,以消除虚假的变化和干扰;Chini等[41 ] 发现若不对植被进行甄别,20%~50%未受损的建筑物会被误判为“倒塌”.但在某些情况下可通过植被的变化提取需要的灾情信息,如热带地区或山区植被的明显减少表明可能发生了山体滑坡[47 ] .Yang等[33 ] 使用震前的Landsat 5 TM和震后的Terra ASTER影像通过波段运算得到NDVI,通过对比试验表明了相对辐射校正的必要性和植被指数法对滑坡检测的有效性. ...
... 人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译.人工解译虽然费事费力,但精度和可靠性相对更高[58 ] .对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] . ...
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
... Chini等[41 ] ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Review article digital change detection techniques using remotely-sensed data
2
1989
... 代数法优点在于简单、易操作,便于解释变化检测结果,但在实际应用中存在着较多不足:①代数法对震前、震后影像的几何配准和相对辐射校正有着很高的要求,用于高分辨率影像时会产生明显的“椒盐效应”;②对于地震应急中常用的多光谱数据,简单的代数运算没有考虑波段之间的相关性,存在着信息重叠和冗余,且波段间的运算使得不同的数值可能会产生同样大小的差值或比值,忽略了不同地物在不同敏感波段存在的差异,丢失潜在的可利用信息[42 ] ;③由于存在“同质异谱”和“异质同谱”的现象,代数法提取的震害信息可能存在误检或遗漏. ...
... 然而,分类后比较法的缺点限制了其在震害信息提取中的应用:①对每幅影像单独分类,增加了计算量和复杂度,影响时效性;②两次分类的误差会传递到变化检测的结果中而被放大[18 , 42 ] .此外,震害信息提取中,变化检测的关注点主要在于“是否”发生了变化或发生变化的程度,至于具体类型的变化(from-to)并非必须获取的信息. ...
Change detection of earthquake-damaged buildings on remote sensing image and its application in seismic disaster assessment
3
2003
... Zhang等[43 ] 对1999年台湾南投地震的震前震后影像,通过设定相关系数的阈值范围,将建筑物的损毁程度划分成了4个等级;Rathje等[44 ] 采用了相关系数法提取巴姆城区震后的建筑物震害,并发现灰度共生矩阵(Gray Level of Co-occurrence Matrix,GLCM)的相关系数比DN值的相关系数能更好地区分损毁和完好的房屋;Chesnel等[26 ] 以屋顶范围作为窗口计算相关系数,实现了震前震后屋顶矢量的精确配准;张伟等[24 ] 为了减小配准误差对相关系数的影响,通过在对应点周围搜索相关系数最大的潜在同名点,使得倒塌和未倒塌建筑的分类精度(κ 系数)由0.60提高到了0.85.可见,相关系数不仅可以作为变化检测的特征变量,也可以作为几何配准的相似测度. ...
... Saito等[59 ] 采用震后的IKONOS全色影像,并结合现场实地调查的照片,对目视解译的效果进行评价,表明了目视解译的快速性和对部分建筑损毁解译的局限性;Yamazaki等[62 ] 采用巴姆地震震前和震后的QuickBird影像,对12 063个建筑物的损毁程度逐个进行了目视解译;Ehrlich等[7 ] 通过震后映秀城区的WorldView影像、震前和震后都江堰城区的SPOT-5影像识别了滑坡灾害和建筑物的损毁情况;吴玮莹等[63 ] 利用震后获取的高分二号、北京二号影像和震前的Google Earth影像,通过人工目视解译的方法,以道路整齐边界遭受的破坏和山体植被覆盖的变化作为解译标志,初步建立了九寨沟地震的同震滑坡分布图.为便于目视解译,Zhang等[43 ] 在提取1999年南头地震的灾情信息时,将震前和震后共2幅灰度影像进行假彩色合成显示,在新合成的彩色图像中,未发生变化的区域为白色—灰色,发生变化区域表现为彩色,这些色彩就包含了震害信息,据此可判断震害分布及其性质. ...
... Zhang等[43 ] ...
Earthquake damage identification using multi-temporal high-resolution optical satellite imagery
4
2005
... Zhang等[43 ] 对1999年台湾南投地震的震前震后影像,通过设定相关系数的阈值范围,将建筑物的损毁程度划分成了4个等级;Rathje等[44 ] 采用了相关系数法提取巴姆城区震后的建筑物震害,并发现灰度共生矩阵(Gray Level of Co-occurrence Matrix,GLCM)的相关系数比DN值的相关系数能更好地区分损毁和完好的房屋;Chesnel等[26 ] 以屋顶范围作为窗口计算相关系数,实现了震前震后屋顶矢量的精确配准;张伟等[24 ] 为了减小配准误差对相关系数的影响,通过在对应点周围搜索相关系数最大的潜在同名点,使得倒塌和未倒塌建筑的分类精度(κ 系数)由0.60提高到了0.85.可见,相关系数不仅可以作为变化检测的特征变量,也可以作为几何配准的相似测度. ...
... 在地震应急中,通常并不关注植被的变化情况,甚至需要剔除植被的影响,如Yamazaki等[22 ,44 -46 ] 在提取城市的建筑物震害信息时,对植被进行了掩膜处理,以消除虚假的变化和干扰;Chini等[41 ] 发现若不对植被进行甄别,20%~50%未受损的建筑物会被误判为“倒塌”.但在某些情况下可通过植被的变化提取需要的灾情信息,如热带地区或山区植被的明显减少表明可能发生了山体滑坡[47 ] .Yang等[33 ] 使用震前的Landsat 5 TM和震后的Terra ASTER影像通过波段运算得到NDVI,通过对比试验表明了相对辐射校正的必要性和植被指数法对滑坡检测的有效性. ...
... Rathje等[44 ] ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Extraction of earthquake-induced collapsed buildings using very high-resolution imagery and airborne Lidar data
0
2015
Extraction of earthquake-induced collapsed buildings from Bi-Temporal VHR images using object-level Homogeneity Index and histogram
2
2019
... 在地震应急中,通常并不关注植被的变化情况,甚至需要剔除植被的影响,如Yamazaki等[22 ,44 -46 ] 在提取城市的建筑物震害信息时,对植被进行了掩膜处理,以消除虚假的变化和干扰;Chini等[41 ] 发现若不对植被进行甄别,20%~50%未受损的建筑物会被误判为“倒塌”.但在某些情况下可通过植被的变化提取需要的灾情信息,如热带地区或山区植被的明显减少表明可能发生了山体滑坡[47 ] .Yang等[33 ] 使用震前的Landsat 5 TM和震后的Terra ASTER影像通过波段运算得到NDVI,通过对比试验表明了相对辐射校正的必要性和植被指数法对滑坡检测的有效性. ...
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
Landslide inventory maps: New tools for an old problem
1
2012
... 在地震应急中,通常并不关注植被的变化情况,甚至需要剔除植被的影响,如Yamazaki等[22 ,44 -46 ] 在提取城市的建筑物震害信息时,对植被进行了掩膜处理,以消除虚假的变化和干扰;Chini等[41 ] 发现若不对植被进行甄别,20%~50%未受损的建筑物会被误判为“倒塌”.但在某些情况下可通过植被的变化提取需要的灾情信息,如热带地区或山区植被的明显减少表明可能发生了山体滑坡[47 ] .Yang等[33 ] 使用震前的Landsat 5 TM和震后的Terra ASTER影像通过波段运算得到NDVI,通过对比试验表明了相对辐射校正的必要性和植被指数法对滑坡检测的有效性. ...
Pca‐based Land‐use change detection and analysis using multitemporal and Multisensor Satellite Data
2
2008
... 主成分变换是将具有相关性的多波段影像压缩到完全独立的较少的几个主成分上,使新影像更便于解译.主成分变换法进行变化检测原理是,未变化的区域具有较高的相关性,而发生变化的区域相关性相对较低[48 ] .在提取变化信息时,主成分变换法常与差值法结合使用,依据差值和主成分变换的先后可分为差值主成分法、主成分差值法和多波段主成分变换法[48 -49 ] .主成分变换的缺点是变换后各分量所突出的信息有所不同且物理意义并不明确,具体哪个分量可用来表征变化信息需要一定经验和分析[17 ] . ...
... [48 -49 ].主成分变换的缺点是变换后各分量所突出的信息有所不同且物理意义并不明确,具体哪个分量可用来表征变化信息需要一定经验和分析[17 ] . ...
Monitoring land-cover change by principal component analysis of multitemporal landsat data
1
1980
... 主成分变换是将具有相关性的多波段影像压缩到完全独立的较少的几个主成分上,使新影像更便于解译.主成分变换法进行变化检测原理是,未变化的区域具有较高的相关性,而发生变化的区域相关性相对较低[48 ] .在提取变化信息时,主成分变换法常与差值法结合使用,依据差值和主成分变换的先后可分为差值主成分法、主成分差值法和多波段主成分变换法[48 -49 ] .主成分变换的缺点是变换后各分量所突出的信息有所不同且物理意义并不明确,具体哪个分量可用来表征变化信息需要一定经验和分析[17 ] . ...
基于主成分分析的遥感震害变化检测方法与应用
1
2013
... Tomowski等[31 ] 以2010年海地地震后太子港总统府区域为研究对象,从GLCM能量值第一主成分的差值中提取了损毁的建筑、植被和新增的救灾帐篷等变化信息;李成龙等[50 ] 采用主成分差值法对2003年阿尔及利亚地震的建筑物损毁情况进行检测,然而检测结果中存在很多零散杂乱的图斑,尤其受建筑物阴影变化的影响很大,需要较多的人工干预;Gong等[51 ] 对比了主成分分析法与差值法、面向对象法的震害提取效果,发现主成分分析法并无优势,提取的变化多为建筑物阴影的变化而非建筑物的损毁.此外,穗帽变换作为一种特殊的主成分变换方法,可以得到具有物理意义的主成分,也在震害提取中得到应用,如Oommen等[52 ] 基于2001印度古吉拉特邦地震震前、震后的Landsat 7 ETM+影像,利用穗帽变换得到的湿度分量,探测出了灾区土壤液化的分布情况. ...
基于主成分分析的遥感震害变化检测方法与应用
1
2013
... Tomowski等[31 ] 以2010年海地地震后太子港总统府区域为研究对象,从GLCM能量值第一主成分的差值中提取了损毁的建筑、植被和新增的救灾帐篷等变化信息;李成龙等[50 ] 采用主成分差值法对2003年阿尔及利亚地震的建筑物损毁情况进行检测,然而检测结果中存在很多零散杂乱的图斑,尤其受建筑物阴影变化的影响很大,需要较多的人工干预;Gong等[51 ] 对比了主成分分析法与差值法、面向对象法的震害提取效果,发现主成分分析法并无优势,提取的变化多为建筑物阴影的变化而非建筑物的损毁.此外,穗帽变换作为一种特殊的主成分变换方法,可以得到具有物理意义的主成分,也在震害提取中得到应用,如Oommen等[52 ] 基于2001印度古吉拉特邦地震震前、震后的Landsat 7 ETM+影像,利用穗帽变换得到的湿度分量,探测出了灾区土壤液化的分布情况. ...
Earthquake building damage detection with object-oriented change detection
3
2013
... Tomowski等[31 ] 以2010年海地地震后太子港总统府区域为研究对象,从GLCM能量值第一主成分的差值中提取了损毁的建筑、植被和新增的救灾帐篷等变化信息;李成龙等[50 ] 采用主成分差值法对2003年阿尔及利亚地震的建筑物损毁情况进行检测,然而检测结果中存在很多零散杂乱的图斑,尤其受建筑物阴影变化的影响很大,需要较多的人工干预;Gong等[51 ] 对比了主成分分析法与差值法、面向对象法的震害提取效果,发现主成分分析法并无优势,提取的变化多为建筑物阴影的变化而非建筑物的损毁.此外,穗帽变换作为一种特殊的主成分变换方法,可以得到具有物理意义的主成分,也在震害提取中得到应用,如Oommen等[52 ] 基于2001印度古吉拉特邦地震震前、震后的Landsat 7 ETM+影像,利用穗帽变换得到的湿度分量,探测出了灾区土壤液化的分布情况. ...
... Matsumoto等[55 ,67 ] 对震前、震后的影像进行分类,提取出两个时相的建筑物,通过建筑物的减少来判别损毁.类似地,Gong等[51 ] 根据房屋的长宽比、纹理、尺度等特征建立分类的规则集提取出了震前和震后影像中的房屋对象,认为震后相比于震前减少的房屋对象即为被损毁的房屋.Tiede等[27 ] 在未对震前、震后影像配准的情况下,通过对比对应建筑物阴影的面积和形状的变化来判断建筑物是否损毁;Zhu等[68 ] 对震后影像进行分割,通过决策树分类提取出人造地表(如道路、停车场)和裸土,采用多波段主成分变换法判断这些地表是否发生液化.Anniballe等[32 ] 借助于拉奎拉市已有的GIS底图,利用建筑物高度和传感器姿态等数据精确计算了震前和震后影像中建筑物的矢量多边形,测试了4类共13种影像特征(光谱特征、信息论特征、纹理特征、HSV颜色特征)和2种分类方法(最大后验概率和支持向量机)对建筑物震害的提取效果. ...
... Gong等[51 ] ...
Documenting earthquake-induced liquefaction using satellite remote sensing image transformations
2
2013
... Tomowski等[31 ] 以2010年海地地震后太子港总统府区域为研究对象,从GLCM能量值第一主成分的差值中提取了损毁的建筑、植被和新增的救灾帐篷等变化信息;李成龙等[50 ] 采用主成分差值法对2003年阿尔及利亚地震的建筑物损毁情况进行检测,然而检测结果中存在很多零散杂乱的图斑,尤其受建筑物阴影变化的影响很大,需要较多的人工干预;Gong等[51 ] 对比了主成分分析法与差值法、面向对象法的震害提取效果,发现主成分分析法并无优势,提取的变化多为建筑物阴影的变化而非建筑物的损毁.此外,穗帽变换作为一种特殊的主成分变换方法,可以得到具有物理意义的主成分,也在震害提取中得到应用,如Oommen等[52 ] 基于2001印度古吉拉特邦地震震前、震后的Landsat 7 ETM+影像,利用穗帽变换得到的湿度分量,探测出了灾区土壤液化的分布情况. ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Change detection of an earthquake-induced Barrier Lake based on remote sensing image classification
2
2010
... Xu等[53 ] 采用分类后比较法探测出了2008年汶川地震中的唐家山堰塞湖,采用决策树的方法分别对配准好的震前、震后影像进行分类,将震前影像提取出的水体和震后影像提取出的水体、滑坡进行叠加,发现滑坡阻塞河道导致河道变宽形成堰塞湖,进而定量估算了河道增加的宽度和被淹区域的面积. ...
... Gong等[51 ]
分类后比较法 克服了影像因获取条件引起差异;可以提供变化的类型 需要训练样本;需要多次分类;分类误差会传递 Xu等[53 ] 面向对象的方法 目视解译 简单实用,精度和可靠性高 需投入大量的人力资源 Yamazaki等[62 ] ...
Change detection using object features
1
2008
... 基于像元的变化检测算法在应用于高分辨率影像时常常会受到局限.首先,分辨率越高,像元越小,空间不一致或者几何配准误差等造成的影响越加明显;其次,高分辨率影像在承载更丰富的地表信息的同时,也导致了同一地物内部像元灰度起伏明显,在使用基于像元的方法时,过多无关的变化会被检测出来,产生“椒盐效应”[18 , 54 ] . ...
Extraction of damaged buildings using high resolution satellite images in the 2006 central java earthquake
3
2006
... 面向对象的变化检测方法能很好地克服基于像元的方法在配准和噪声等方面的局限[55 ] ,通过影像中地物对象的提取,直接比较这些地物的光谱特征、纹理特征、形状特征和上下文关系等的变化,更加高效可靠[56 ] ,此外还能减少光照条件、阴影等差异造成的影响[57 ] .面向对象方法的关键在于对象的提取(包括分割和分类),而对变化信息的提取基本上都是借鉴了基于像元的方法,所不同的只是将像元间的比较变成了对象间的比较[57 ] .Bitelli等[25 ,55 ,58 ] 对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害. ...
... ,55 ,58 ]对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害. ...
... Matsumoto等[55 ,67 ] 对震前、震后的影像进行分类,提取出两个时相的建筑物,通过建筑物的减少来判别损毁.类似地,Gong等[51 ] 根据房屋的长宽比、纹理、尺度等特征建立分类的规则集提取出了震前和震后影像中的房屋对象,认为震后相比于震前减少的房屋对象即为被损毁的房屋.Tiede等[27 ] 在未对震前、震后影像配准的情况下,通过对比对应建筑物阴影的面积和形状的变化来判断建筑物是否损毁;Zhu等[68 ] 对震后影像进行分割,通过决策树分类提取出人造地表(如道路、停车场)和裸土,采用多波段主成分变换法判断这些地表是否发生液化.Anniballe等[32 ] 借助于拉奎拉市已有的GIS底图,利用建筑物高度和传感器姿态等数据精确计算了震前和震后影像中建筑物的矢量多边形,测试了4类共13种影像特征(光谱特征、信息论特征、纹理特征、HSV颜色特征)和2种分类方法(最大后验概率和支持向量机)对建筑物震害的提取效果. ...
Comparison of geo-object based and pixel-based change detection of riparian environments using high spatial resolution multi-spectral imagery
1
2010
... 面向对象的变化检测方法能很好地克服基于像元的方法在配准和噪声等方面的局限[55 ] ,通过影像中地物对象的提取,直接比较这些地物的光谱特征、纹理特征、形状特征和上下文关系等的变化,更加高效可靠[56 ] ,此外还能减少光照条件、阴影等差异造成的影响[57 ] .面向对象方法的关键在于对象的提取(包括分割和分类),而对变化信息的提取基本上都是借鉴了基于像元的方法,所不同的只是将像元间的比较变成了对象间的比较[57 ] .Bitelli等[25 ,55 ,58 ] 对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害. ...
Object-based change detection
2
2012
... 面向对象的变化检测方法能很好地克服基于像元的方法在配准和噪声等方面的局限[55 ] ,通过影像中地物对象的提取,直接比较这些地物的光谱特征、纹理特征、形状特征和上下文关系等的变化,更加高效可靠[56 ] ,此外还能减少光照条件、阴影等差异造成的影响[57 ] .面向对象方法的关键在于对象的提取(包括分割和分类),而对变化信息的提取基本上都是借鉴了基于像元的方法,所不同的只是将像元间的比较变成了对象间的比较[57 ] .Bitelli等[25 ,55 ,58 ] 对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害. ...
... [57 ].Bitelli等[25 ,55 ,58 ] 对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害. ...
Urban building damage detection from very high resolution imagery using ocsvm and spatial features
3
2010
... 面向对象的变化检测方法能很好地克服基于像元的方法在配准和噪声等方面的局限[55 ] ,通过影像中地物对象的提取,直接比较这些地物的光谱特征、纹理特征、形状特征和上下文关系等的变化,更加高效可靠[56 ] ,此外还能减少光照条件、阴影等差异造成的影响[57 ] .面向对象方法的关键在于对象的提取(包括分割和分类),而对变化信息的提取基本上都是借鉴了基于像元的方法,所不同的只是将像元间的比较变成了对象间的比较[57 ] .Bitelli等[25 ,55 ,58 ] 对比了基于像元和面向对象的方法,结果表明面向对象方法能够克服配准和椒盐噪声的影响,能更有效地提取震害. ...
... 人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译.人工解译虽然费事费力,但精度和可靠性相对更高[58 ] .对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] . ...
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
Using high-resolution satellite images for post-earthquake building damage assessment: A study following the 26 January 2001 Gujarat earthquake
2
2004
... 震后的灾损评估可分为两大类,即基于算法的定量分析和基于人工解译的定性分析,应当根据具体的任务目标选择合适、有效的分析方法[59 ] . ...
... Saito等[59 ] 采用震后的IKONOS全色影像,并结合现场实地调查的照片,对目视解译的效果进行评价,表明了目视解译的快速性和对部分建筑损毁解译的局限性;Yamazaki等[62 ] 采用巴姆地震震前和震后的QuickBird影像,对12 063个建筑物的损毁程度逐个进行了目视解译;Ehrlich等[7 ] 通过震后映秀城区的WorldView影像、震前和震后都江堰城区的SPOT-5影像识别了滑坡灾害和建筑物的损毁情况;吴玮莹等[63 ] 利用震后获取的高分二号、北京二号影像和震前的Google Earth影像,通过人工目视解译的方法,以道路整齐边界遭受的破坏和山体植被覆盖的变化作为解译标志,初步建立了九寨沟地震的同震滑坡分布图.为便于目视解译,Zhang等[43 ] 在提取1999年南头地震的灾情信息时,将震前和震后共2幅灰度影像进行假彩色合成显示,在新合成的彩色图像中,未发生变化的区域为白色—灰色,发生变化区域表现为彩色,这些色彩就包含了震害信息,据此可判断震害分布及其性质. ...
Co-seismic surface effects from very high resolution panchromatic images: The case of the 2005 Kashmir (Pakistan) earthquake
1
2011
... 人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译.人工解译虽然费事费力,但精度和可靠性相对更高[58 ] .对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] . ...
Factors affecting earthquake responders’ building damage information needs and use
1
2021
... 人脑是人工智能的模仿对象,可根据自身的认知和经验自动地实现目标对象的提取及其邻域信息的解译.人工解译虽然费事费力,但精度和可靠性相对更高[58 ] .对于地震灾情调查这种情况紧急的实际问题,尤其是在人口密度大的城市范围内提取建筑物震害时,人工目视解译仍然是最常用、最可靠的方法[7 ,27 ,40 -41 ,60 -61 ] . ...
Visual damage interpretation of buildings in bam city using QuickBird images following the 2003 bam,Iran, Earthquake
3
2005
... Saito等[59 ] 采用震后的IKONOS全色影像,并结合现场实地调查的照片,对目视解译的效果进行评价,表明了目视解译的快速性和对部分建筑损毁解译的局限性;Yamazaki等[62 ] 采用巴姆地震震前和震后的QuickBird影像,对12 063个建筑物的损毁程度逐个进行了目视解译;Ehrlich等[7 ] 通过震后映秀城区的WorldView影像、震前和震后都江堰城区的SPOT-5影像识别了滑坡灾害和建筑物的损毁情况;吴玮莹等[63 ] 利用震后获取的高分二号、北京二号影像和震前的Google Earth影像,通过人工目视解译的方法,以道路整齐边界遭受的破坏和山体植被覆盖的变化作为解译标志,初步建立了九寨沟地震的同震滑坡分布图.为便于目视解译,Zhang等[43 ] 在提取1999年南头地震的灾情信息时,将震前和震后共2幅灰度影像进行假彩色合成显示,在新合成的彩色图像中,未发生变化的区域为白色—灰色,发生变化区域表现为彩色,这些色彩就包含了震害信息,据此可判断震害分布及其性质. ...
... Yamazaki等[62 ] ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
基于高分卫星遥感影像的地震应急滑坡编目与分布特征探讨——以2017年8月8日九寨沟7.0级地震为例
1
2017
... Saito等[59 ] 采用震后的IKONOS全色影像,并结合现场实地调查的照片,对目视解译的效果进行评价,表明了目视解译的快速性和对部分建筑损毁解译的局限性;Yamazaki等[62 ] 采用巴姆地震震前和震后的QuickBird影像,对12 063个建筑物的损毁程度逐个进行了目视解译;Ehrlich等[7 ] 通过震后映秀城区的WorldView影像、震前和震后都江堰城区的SPOT-5影像识别了滑坡灾害和建筑物的损毁情况;吴玮莹等[63 ] 利用震后获取的高分二号、北京二号影像和震前的Google Earth影像,通过人工目视解译的方法,以道路整齐边界遭受的破坏和山体植被覆盖的变化作为解译标志,初步建立了九寨沟地震的同震滑坡分布图.为便于目视解译,Zhang等[43 ] 在提取1999年南头地震的灾情信息时,将震前和震后共2幅灰度影像进行假彩色合成显示,在新合成的彩色图像中,未发生变化的区域为白色—灰色,发生变化区域表现为彩色,这些色彩就包含了震害信息,据此可判断震害分布及其性质. ...
基于高分卫星遥感影像的地震应急滑坡编目与分布特征探讨——以2017年8月8日九寨沟7.0级地震为例
1
2017
... Saito等[59 ] 采用震后的IKONOS全色影像,并结合现场实地调查的照片,对目视解译的效果进行评价,表明了目视解译的快速性和对部分建筑损毁解译的局限性;Yamazaki等[62 ] 采用巴姆地震震前和震后的QuickBird影像,对12 063个建筑物的损毁程度逐个进行了目视解译;Ehrlich等[7 ] 通过震后映秀城区的WorldView影像、震前和震后都江堰城区的SPOT-5影像识别了滑坡灾害和建筑物的损毁情况;吴玮莹等[63 ] 利用震后获取的高分二号、北京二号影像和震前的Google Earth影像,通过人工目视解译的方法,以道路整齐边界遭受的破坏和山体植被覆盖的变化作为解译标志,初步建立了九寨沟地震的同震滑坡分布图.为便于目视解译,Zhang等[43 ] 在提取1999年南头地震的灾情信息时,将震前和震后共2幅灰度影像进行假彩色合成显示,在新合成的彩色图像中,未发生变化的区域为白色—灰色,发生变化区域表现为彩色,这些色彩就包含了震害信息,据此可判断震害分布及其性质. ...
Use of shadows for detection of earthquake-induced collapsed buildings in high-resolution satellite imagery
3
2013
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
... GIS底图可以极大地便利震害信息提取,降低问题的复杂度,减少计算量和工作量,提高效率和精度,如Samadzadegan和Rastiveisi等[30 ,32 ,64 -66 ] 直接采用了GIS底图提供的建筑物矢量,省去了对象分割和建筑物分类等工作.因此,在地震危险度高的城市应大力开展城市GIS底图的绘制和更新工作. ...
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Disaster damage assessment of buildings using adaptive self-similarity descriptor
1
2016
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
A GIS based approach for earthquake loss estimation based on the immediate extraction of damaged buildings
2
2017
... Gusella等[20 ] 展示了典型的多时相叠加法在建筑物震害信息提取中的应用,作者首先基于均质性对震前影像进行分割,然后采用最近邻法进行分类,将提取的建筑物掩膜叠加到震后影像中,最后再对震后的建筑物进行分类,得到“受损”和“完好”的两类建筑物,总体精度可达70.5%.Samadzadegan等[30 ] 利用已有的建筑物矢量,比较了震前和震后影像建筑物范围内的8种纹理特征的变化,以此作为输入建立模糊推理系统,将1 133个建筑物识别为未损坏、中度损坏、严重损毁和完全损毁4类.Chini等[41 ] 在影像分割的基础上,利用多尺度的形态学开运算和闭运算解决了分类中的“同谱异物”问题,将提取的建筑物掩膜同时叠加到配准好的震前和震后影像,通过建筑物对象的归一化差值得到了损毁建筑的分布图.Li等[58 ] 采用多时相叠加法提取出了建筑物对象,基于震前、震后影像在光谱和伪交叉变异函数纹理特征上的差异,聚焦于“完全倒塌”的建筑物这一单一类别,引入了单类支持向量机来识别倒塌的建筑物,并对比了基于像元和面向对象方法的分类效果,显示了后者的优势.Tong等[64 ] 基于建筑物的矢量底图计算其理论投影范围内的DN值均值在震前、震后影像中的变化,区分出了完全倒塌、部分倒塌和未倒塌的建筑物.Liu等[46 ] 在对植被和均质区域进行掩膜处理后,通过分析倒塌和完好的建筑物在对象直方图特征的差异,在剩余的区域内提取出了不同损毁程度的建筑物.Kahraman等[65 ] 基于已有的建筑物矢量,引入了自相似描述子(self-similarity sescriptor,SSD)来比较震前、震后影像中房屋对象在形状上的变化.Ranjbar等[66 ] 提取了震前的建筑物对象并叠加到震后影像中,优选了8种特征(包括灰度均值、纹理特征、半方差函数、Gabor特征等),利用模糊推理系统建立了这8特征与建筑物损毁程度之间的关系(if-then),避免了阈值选取问题的困扰. ...
... GIS底图可以极大地便利震害信息提取,降低问题的复杂度,减少计算量和工作量,提高效率和精度,如Samadzadegan和Rastiveisi等[30 ,32 ,64 -66 ] 直接采用了GIS底图提供的建筑物矢量,省去了对象分割和建筑物分类等工作.因此,在地震危险度高的城市应大力开展城市GIS底图的绘制和更新工作. ...
Detection of damages due to earthquakes using digital aerial images
1
2008
... Matsumoto等[55 ,67 ] 对震前、震后的影像进行分类,提取出两个时相的建筑物,通过建筑物的减少来判别损毁.类似地,Gong等[51 ] 根据房屋的长宽比、纹理、尺度等特征建立分类的规则集提取出了震前和震后影像中的房屋对象,认为震后相比于震前减少的房屋对象即为被损毁的房屋.Tiede等[27 ] 在未对震前、震后影像配准的情况下,通过对比对应建筑物阴影的面积和形状的变化来判断建筑物是否损毁;Zhu等[68 ] 对震后影像进行分割,通过决策树分类提取出人造地表(如道路、停车场)和裸土,采用多波段主成分变换法判断这些地表是否发生液化.Anniballe等[32 ] 借助于拉奎拉市已有的GIS底图,利用建筑物高度和传感器姿态等数据精确计算了震前和震后影像中建筑物的矢量多边形,测试了4类共13种影像特征(光谱特征、信息论特征、纹理特征、HSV颜色特征)和2种分类方法(最大后验概率和支持向量机)对建筑物震害的提取效果. ...
Mapping earthquake induced liquefaction surface effects from the 2011 Tohoku Earthquake using satellite imagery
1
2016
... Matsumoto等[55 ,67 ] 对震前、震后的影像进行分类,提取出两个时相的建筑物,通过建筑物的减少来判别损毁.类似地,Gong等[51 ] 根据房屋的长宽比、纹理、尺度等特征建立分类的规则集提取出了震前和震后影像中的房屋对象,认为震后相比于震前减少的房屋对象即为被损毁的房屋.Tiede等[27 ] 在未对震前、震后影像配准的情况下,通过对比对应建筑物阴影的面积和形状的变化来判断建筑物是否损毁;Zhu等[68 ] 对震后影像进行分割,通过决策树分类提取出人造地表(如道路、停车场)和裸土,采用多波段主成分变换法判断这些地表是否发生液化.Anniballe等[32 ] 借助于拉奎拉市已有的GIS底图,利用建筑物高度和传感器姿态等数据精确计算了震前和震后影像中建筑物的矢量多边形,测试了4类共13种影像特征(光谱特征、信息论特征、纹理特征、HSV颜色特征)和2种分类方法(最大后验概率和支持向量机)对建筑物震害的提取效果. ...
Using MODIS NDVI time series to identify geographic patterns of landslides in vegetated regions
1
2013
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Estimating Wenchuan Earthquake induced landslides based on remote sensing
1
2010
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Triple collocation to assess classification accuracy without a ground truth in case of earthquake damage assessment
1
2018
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
Building damage assessment for rapid disaster response with a deep object-based semantic change detection framework: From natural disasters to man-made disasters
1
2021
... 验证数据通常有两个来源,一是高分辨率遥感影像的目视解译结果[33 ,69 -70 ] ;二是地面实地调查的结果[32 ,41 ,44 ,52 ] .但地面评估与遥感评估的视角和依据的信息不同,建筑物震害的等级划分标准也不同,二者之间往往存在着较大的偏差.Anniballe等[32 ] 对比了利用QuickBird遥感影像目视解译和现场实地调查的建筑物损毁程度,结果表明遥感影像在识别建筑物震害时存在着严重的“低估”和“漏判”. 一方面是因为高分辨率遥感影像仍不足以识别更细节的信息[62 ] ,如房屋的裂痕、散落的瓦砾等;另一方面,遥感影像多为天顶或近似天顶观测,对建筑物侧面墙壁的损毁信息获取不足,而实地调查评估中则会较多地利用这些信息,从而导致遥感数据与实地调查更大的偏差[5 , 9 ] ,如低矮楼房的平坍塌或低楼层的垮台(pancake collapse)这类严重的建筑物震害在遥感影像中容易被遗漏[71 ] .因而当前较多研究采用自己的准则定义建筑物损毁程度[64 , 72 ] ,导致与实地调查和其他研究成果之间缺乏对应性和可比性.因此在震害信息提取中,还需要建立起遥感影像与实地调查结果之间的对应关系,以解决遥感影像对震害的“低估”和“漏判”问题,也使不同的研究成果之间具有可比性. ...
遥感影像变化检测方法综述及展望
1
2011
... 为了系统地综述变化检测算法在地震灾情调查中的应用,本研究对这些算法的类别进行了较为细致的划分,但这种类别的划分并非唯一的[18 -19 ,35 ,73 ] ,不同的类别的界限不是严格分明的,可能存在着交叉、借鉴或套用.在实际研究和应用这些方法时,也不会拘泥于算法的类别,如Janalipour等[74 ] 综合采用了差值法、变化向量法、植被指数法、主成分分析法等多种方法,构建了多准则判据来判断建筑物是否损毁,使检测结果更加稳定、精确. ...
遥感影像变化检测方法综述及展望
1
2011
... 为了系统地综述变化检测算法在地震灾情调查中的应用,本研究对这些算法的类别进行了较为细致的划分,但这种类别的划分并非唯一的[18 -19 ,35 ,73 ] ,不同的类别的界限不是严格分明的,可能存在着交叉、借鉴或套用.在实际研究和应用这些方法时,也不会拘泥于算法的类别,如Janalipour等[74 ] 综合采用了差值法、变化向量法、植被指数法、主成分分析法等多种方法,构建了多准则判据来判断建筑物是否损毁,使检测结果更加稳定、精确. ...
Building change detection after earthquake using multi-criteria decision analysis based on extracted information from high spatial resolution satellite images
1
2017
... 为了系统地综述变化检测算法在地震灾情调查中的应用,本研究对这些算法的类别进行了较为细致的划分,但这种类别的划分并非唯一的[18 -19 ,35 ,73 ] ,不同的类别的界限不是严格分明的,可能存在着交叉、借鉴或套用.在实际研究和应用这些方法时,也不会拘泥于算法的类别,如Janalipour等[74 ] 综合采用了差值法、变化向量法、植被指数法、主成分分析法等多种方法,构建了多准则判据来判断建筑物是否损毁,使检测结果更加稳定、精确. ...
/
〈
〉