1.一种活动断裂解译方法,其特征在于,所述方法包括:获取工作区域的热红外遥感数据、光学遥感数据、地形资料和地质资料;基于遥感图像处理软件,对所述热红外遥感数据进行辐射标定;基于大气校正法进行地表温度反演,并基于山区的坡向、坡度和太阳辐射,对所述地表温度进行修正,获取反演修正的地表温度;基于数值仿真,计算地面净辐射,获取数值仿真的地表温度;基于所述反演修正的地表温度、所述数值仿真的地表温度,获取所述工作区域的地表热异常;基于所述地形资料、所述地质资料和所述光学遥感数据,提取活动断裂光学遥感解译成果;基于GIS软件,叠置分析所述活动断裂光学遥感解译成果和所述地表热异常,获取活动断裂影响范围的面状信息。
2.根据权利要求1所述的活动断裂解译方法,其特征在于,所述热红外遥感数据包括ASTER和/或Landsat TM和/或Landsat ETM+和/或LandsatTIRS;所述光学遥感数据包括分辨率为2.5m的SPOT和/或分辨率为2m的高分一号和/或分辨率为1m的高分二号和/或分辨率更高的卫星遥感数据和/或航空遥感数据;所述地形资料包括数字高程模型DEM和/或数字表面模型DSM和/或地形图;所述地质资料包括区域地质资料和/或水文地质资料和/或地震资料和/或地灾资料和/或矿产资料。
3.根据权利要求2所述的活动断裂解译方法,其特征在于,所述进行地表温度反演的步骤,具体为,计算植被覆盖度Pv,Pv=(NDVI-NDVI Soil )/(NDVI Veg -NDVI Soil ),其中,NDVI为归一化植被指数,NDVI Soil 为完全是裸土或无植被覆盖区域的NDVI值,NDVI Veg 为完全被植被所覆盖的像元的NDVI值;基于所述植被覆盖度Pv计算地表比辐射率ε,ε=0.004Pv+0.986;获取大气剖面参数,所述大气剖面参数包括大气在热红外波段的透过率τ、大气向上辐射亮度L↑、大气向下辐射亮辐射亮度L↓;基于所述大气剖面参数,计算卫星传感器接收到的热红外辐射亮度值L λ ,L λ =(εB(T S )+(1-ε)L↓)τ+L↑,其中,ε为地表比辐射率,T S 为地表真实温度(K),B(T S )为黑体热辐射亮度;计算温度为T的黑体在热红外波段的辐射亮度B(T S ),B(T S )=(L λ -L↑-τ(1-ε)L↓)/τε;计算地表温度T rs ,T rs =K 2 /ln(K 1 /B(T S )+1),其中,对于Landsat TM,K 1 =607.76W/(m 2 ·μm·sr),K 2 =1260.56K;对于Landsat ETM+,K 1 =666.09W/(m 2 ·μm·sr),K 2 =1282.71K;对于Landsat TIRS Band10,K 1 =774.89W/(m 2 ·μm·sr),K 2 =1321.08K。
4.根据权利要求3所述的活动断裂解译方法,其特征在于,设置NDVI Veg =0.70、NDVI Soil =0.05,当某个像元的NDVI大于0.70时,Pv取值为1;当某个像元的NDVI小于0.05,Pv取值为0。
5.根据权利要求4所述的活动断裂解译方法,其特征在于,所述对所述地表温度进行修正的步骤,具体为,利用坡向、坡度和太阳辐射的函数关系对山区地表温度差进行修正,设置所述反演修正的地表温度为T s ,T s =T rs +ΔT,其中,ΔT=f(ΔQ)=2aRcot(h)sinα|cos(β-η)|+b,ΔT为地表斜坡温差修正,ΔQ为太阳辐射差异,R为平面单位太阳辐射吸收率,h为太阳高度角,α为坡度,β为太阳方位角,η为坡向,a、b为回归常数。
6.根据权利要求5所述的活动断裂解译方法,其特征在于,针对林地、灌丛两种类型,统计地表温度与归一化植被指数NDVI的线性关系进行修正。
7.根据权利要求6所述的活动断裂解译方法,其特征在于,所述获取数值仿真的地表温度的步骤,具体为,基于地面净辐射Q和stefan-boltzmann方程,其中,Q=(Qa+Qb)-(Qr+Qc+Qd+Qw),计算出所述数值仿真的地表温度 其中,ε为黑体辐射系数,σ为斯特藩常量,取5.67×10 -8 W·m -2 ·K -4 ,a、b为回归常数,Qa为太阳直射辐射,Qb为太阳散射辐射,Qr为自身辐射,Qc为空气对流换热,Qd为热传导,Qw为潜热交换。
8.根据权利要求7所述的活动断裂解译方法,其特征在于,在S时刻所述太阳直射辐射的计算方法如下:其中,Qa(i,j)表示坐标为(i,j)的点太阳直射辐射,纬度为 坡度为α,坡向为β,I a 为太阳常数,δ为太阳倾角,w sr ,w ss 为日升、日落的时角,u、v、w分别定义如下:w=sinαcsinβ所述太阳散射辐射的计算公式为Qb=c 1 (sinδ) c2 ,式中,c 1 、c 2 分别为经验参数,δ为太阳倾角;所述自身辐射Qr=εσT 4 ,ε为散热率,T为地表温度,σ为斯特藩常量;所述空气对流换热Qc=H(T air -T),H为对流换热系数,T air 为低空大气温度;所述热传导Qd=∈(T-T inr )/d,∈为导热率,T in 为物体内部温度,d为物体厚度;所述潜热交换Qw=(ρC p /γ)·(e a (T)-e s (T))/R,式中,ρ为材料密度,C p 为定压比热容,γ为比湿常数,e a (T)、e s (T)分别为气温为T时的实际气压和饱和气压,R为不同类型地物的空气动力学总阻抗。
9.根据权利要求8所述的活动断裂解译方法,其特征在于,所述获取所述工作区域的地表热异常的步骤,具体为,基于地表异常热的计算公式ΔT=(T s -T e )/T s ,分别计算昼间地表异常热ΔT昼间、夜间地表热异常ΔT夜间,并将所述昼间地表异常热ΔT昼间和所述夜间地表热异常ΔT夜间进行叠置分析,获取所述工作区域的综合地表热异常ΔT,ΔT=(ΔT昼间+ΔT夜间)。
10.根据权利要求9所述的活动断裂解译方法,其特征在于,所述活动断裂光学遥感解译成果,具体包括,光学遥感影像上活动断裂的线状影像和/或色调差异和/或线状微地貌特征。
11.一种活动断裂解译装置,其特征在于,包括:数据资料获取模块,用于获取工作区域的热红外遥感数据、光学遥感数据、地形资料和地质资料;辐射标定模块,用于基于遥感图像处理软件,对所述热红外遥感数据进行辐射标定;地表温度反演修正模块,用于基于大气校正法进行地表温度反演,并基于山区的坡向、坡度和太阳辐射,对所述地表温度进行修正,获取反演修正的地表温度;地表温度数值仿真模块,用于基于数值仿真,计算地面净辐射,获取数值仿真的地表温度;地表热异常获取模块,用于基于所述反演修正的地表温度、所述数值仿真的地表温度,获取所述工作区域的地表热异常;光学遥感解译成果获取模块,用于基于所述地形资料、所述地质资料和所述光学遥感数据,提取活动断裂光学遥感解译成果;面状信息获取模块,用于基于GIS软件,叠置分析所述活动断裂光学遥感解译成果和所述地表热异常,获取活动断裂影响范围的面状信息。