06高光谱遥感第六讲

Lecture6SpectralFeatureSelectionandExtractionofHyperspectralData™™™™™™™NMNMN/N-MM100/100-33=161,700™10100wain,1978Hughes™™99™9;9.9™™X1X2X3X4X5X2X4™5B1(2.16μm),B2(2.18μm),B3(2.21μm),B4(2.32μm),B5(2.38μm)B1B2B3B4B5nxx1,x2…xnXX=(x1,x2…xn)TX{}∫==dxxXPXEM)(mi:PxXPxiixσ()iiiimxpxdx=∫(){}22iiExMσ=−xixjxixjij2=0nXn()(){}jjiiijmxmxE−−=2σ()(){}TEXMXM∑=−−ijxixjjiijijσσσρ2=HP(i)iK8bitK=25516bitK=65535)(log)(21iPiPHKi∑=−=™™™™™™™ijijkn1,n2……nkRi=Rw+RaiiiRPσ=σiiRwiRai™™™PiFeatureExtractionY=AX™Ym1≤mnn-mmXnm™mn-m(KL)X∑xλIuui()0xIuλ−∑=()nii,,2,1==λnλλλ≥≥≥21[]12,,,Tiiiniuuuu=A=UTKLK-LY……nXYiXiYXn111212122212nnTnnnnuuuuuuyXUXuuu⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎣⎦#x∑KL™™™KLKL™™™98%Maximumnoisefraction(MNF)™MNF,™™Green,A.A.,M.Berm,P.Switzer,andM.D.Craig(1988).ATransformationforOrderingMultispectralDatainTermsofImageQuanlitywithImplicationsforNoiseRemoval.IEEETGRS26,No.1(CA)(CA)CA∑=∑=∑miiwM11iiMEmiimoDTY=DTX()()TAioioEmmmm⎡⎤∑=−−⎣⎦∑==MiiomMm11yTxwwDD∑=∑yTxAADD∑=∑dTd2Txwwddσ=∑2TxAAddσ=∑()0xAwdλ∑−∑=CA∑w∑A,DDTCAnλλλ≥≥≥21™CA™NMCAM-1NM-1M-1CADVegetationInverted-gaussianfunctionfitwithin670-780nmspectralrangeDerivativetodeterminetheslopeofrededge(Srg)andareaofrededge(Arg)D20020()()()exp()2sspRRRRλλλσλλσ−−=−−=+))(max(750680λdrSrg=750680ArgRdλ=∫RsR0σs2ρ2ρρ1s1ρmλ2λ1λms2ρ2ρρ1s1ρmλ2λ1λm™“”λλλλΔρρρ™(AbsorptionPosition,AP)(AbsorptionDepthAD)™(AbsorptionWidthAW)(AbsorptionAsymmetry,AAmmddSAIρρρρρ21)1(/−+==™™NnTnxTnxTTnxnhabaa,...2,1)(11)(01)(00)(=⎪⎩⎪⎨⎧≤≤=™“”:ClarkandRoush1984rriii=0i=011……kk--11wwiii=0i=011……kk--11Step1Step1ii=0=0rriiwwiiStep2Step2i=ki=k--11jj=i+1=i+1StepStep33iijjiijj,,j=kj=k--1,1,,,wwjjrrjj(1)m(1)m=j+1=j+1(2)(2)m=km=k--11jjwwjjrrjji=ji=jStepStep22(3)(3)iijjwwmmr1r1mm..(4)(4)rrmmr1r1mm,,jj,,jj=j+1,=j+1,StepStep33;;rrmmr1r1mm,,i,ji,j,,mm=m+1,=m+1,StepStep22StepStep44wwii,,i=0,1,i=0,1,……kk--11hhiii=0i=011……kk--11hhiirriiStepStep55rrii=r=rii/h/hiii=0i=011……kk--11™r’i=r(i)/h(i)™r(i),01.0Multi-temporalIndexImageCubeModel(MIIC)Tt1t2tn-1tnbTimeReflectanceWavelengthdIndexTimeeXSaYXTcYTLJOHEHASOTIJOHEHASOVI456789101112wheatcornwheat456789101112EarlyriceLatericeVITRTRTIHEHATIHEHASOVI456789101112North:wheatNorth:cornSouth:earlypaddySouth:laterpaddywaterSo-Tr-Tl-He-JO-Ha-HarvestTime2Multi-temporalIndexImageCubeModel(MIIC)™:™——•“”•λλϕλλdf∫=21)(λλϕλλdf∫=21)(AD=1-MIN(RC)⏐λ2λ1∫−=21)1(λλλdRCAA