干度指数NDBSI和湿度WET

按照徐涵秋老师的文章，干度指标的建筑指标选择的是IBI建筑指标，徐涵秋老师2008年提出的，和土壤指数SI二者合成。
NDBSI=(IBI+SI)/2
以Landsat7为例，以下表格可知Landsat4/5 与之相同
IBI={2b5/(b5+b4)- [b4/(b4+b3) + b2/(b2+b5)]} / {2b5/(b5+b4) + [b4/(b4+b3) + b2/(b2+b5)]}
SI=[ (b5+b3) - (b4+b1) ] / [ (b5+b3) + (b4+b1) ]

Wet=0.2626 b1+0.2141 b2+0.0926 b3+0.0656 b4-0.7629 b5-0.5388 b7
引文：Huang C Derivation of a tasselled cap transformation based on Landsat 7 at satellite reflectance

landsat4-5	landsatETM7	landsat8
-	-	Coastal（海岸波段） b1
Blue b1	Blue b1	Blue b2
Green b2	Green b2	Green b3
Red b3	Red b3	Red b4
NIR近红外 b4	NIR b4	NIR b5
SWIR短波红外b5	SWIR短波红外 b5	SWIR1短波红外1 b6
LWIR热红外b6	LWIR热红外 b6	TIRS 1 热红外1 b10、TIRS 2 热红外1 b11
SWIR短波红外b7	SWIR短波红外 b7	SWIR2短波红外2 b7
-	Pan全色波段 b8	Pan全色波段 b8
-	-	Cirrus圈云波段短波红外1 b9
-	-

对于Landsat 8 而言
IBI={2b6/(b6+b5)- [b5/(b5+b4) + b3/(b3+b6)]} / {2b6/(b6+b5) + [b5/(b5+b4) + b3/(b3+b6)]}
SI=[ (b6+b4) - (b5+b2) ] / [ (b6+b4) + (b5+b2) ]
WET=0.2626 b2+0.2141 b3+0.0926 b4+0.0656 b5-0.7629 b6-0.5388 b7