遥感建筑用地信息的快速提取_徐涵秋

12420108JOURNALOFGEO-INFORMATIONSCIENCEVol.12,No.4Aug.,2010:2009-11-18;:2010-06-02.:(40371107);(JK2009-004)。:(1955-),,,,,。E-mail:[e-mail]徐涵秋,杜丽萍(,,350002):,。,。,,(Index-basedBuilt-upIndex,IBI)。(SAVI、MNDWI、NDBI),。,,。IBI,,,。:(IBI);;1　,。40%[1],。,,,。,101240。2007“18”。20077。,,,18。,。20,。,。,,。:(1)[2-5];(2)[6-8];(3),[9-11]IBIZha[12]NDBI。,,。。,NDBI[12]。IBI,。,。,IBI。2　(IBI),。,,[13-14],[3,12]。,(2005、2007、2008)[9-11],(),,,98%[10]。:SAVI[15],NDBI[12,16]MND-WI[17-18],:SAVI=[(NIR-Red)(1+l)]/(NIR+Red+l)(1)NDBI=(MIR-NIR)/(MIR+NIR)(2)MNDWI=Green-MIR)/(Green+MIR)(3),Green、Red、NIR、MIR、、,LandsatTM2、3、4、5;l,0-1。(2005),,SAVINDVI[9]。NDVI,,30%。SAVI,,15%(Ray,2002)[19]。30%,,SAVI。MNDWI2005[17,18],McFeeters,,。“Photo-grammetricEngineeringandRemoteSensing”,,[20];,,MNDWI(AUC0.90,0.85、0.530.37)[21]。,MNDWI,。(1(a)),,,。(1(b)),NDBISAVIMNDWI,,IBI(Index-basedBuilt-upIndex):IBI=[NDBI-(SAVI+MNDWI)/2]/[NDBI+(SAVI+MNDWI)/2](4)1　(a)NDBI(X)SAVI(Y);(b)Fig.1　Spectralfeaturesofthethreemainland-useclasses.(a)Scatterplot(X-axis:NDBI,Y-axis:SAVI);(b)Spectralsignaturesrepresentedbythemeanvaluesofthethreeclasses　　NDBISAVIMNDWI;SAVIMNDWI(1(b)),,,,,,。2,。3　IBI,5:(1);(2);(3)3;(4)IBI;(5)。,,,,。ERMapper7.1,,ERMapper(3),。21996LandsatTM5754:　　　2　1996LandsatTMFig.2　LandsatTMimageofXiamenIslandin1996anditsenhancedandextractedbuilt-upimage,1,,93.9%(1)。3　ERMapperIBIFig.3　TheIBImoduleintegratedinthegraphicaluserinterfaceoftheERMapperpackage1　Tab.1　Accuracyassessment(%)1781719591.28211311598.26180130310(%)98.8986.92(%)93.9Kappa0.8721,;l(0.5);IBI(2)。4,:#＄RC　　　　　　　#＄numberinput=0.5#SAVI0.5#TheparameteofRC:　　if＄RC==FALSEthengotoCcreatthreeindexband:　　if＄RC!=FALSEthengotoRadiometricCorrection#RadiometricCorrection:　　＄alg=newalgorithm　　set＄algmodergb　　＄count=1increment:　　addredlayer　　＄layer2=currentlayer　　set＄layer2typegreen　　set＄layer2description＄count　　set＄layer2formulatoi1-+＄number＄count+　　set＄layer2datasetto＄RasterDataset　　set＄layer2input1toband＄count　　＄count=＄count+1　　if(＄count=7)gotoincrement　　Deletelayer1#Creatthreeindexbands:　　addredlayer　　＄layer1=firstlayer　　set＄layer1typered　　set＄layer1descriptionNDBI　　set＄layer1formulato(i1-i2)/(i1+i2)　　set＄layer1datasetto＄RasterDataset　　set＄layer1input1toband5　　set＄layer1input2toband4　　addgreenlayer　　＄layer2=currentlayer576　　　　　　　　　　　　　　　　　　　2010　　set＄layer2typegreen　　set＄layer2descriptionSAVI　　set＄layer2formulato　　((i1-i2)＊(1++＄numberinput+))/(i1+i2++＄numberinput+)　　set＄layer2datasetto＄RasterDataset　　set＄layer2input1toband4　　set＄layer2input2toband3　　addbluelayer　　＄layer3=currentlayer　　set＄layer3typeblue　　set＄layer3descriptionMNDWI　　set＄layer3formulato(i1-i2)/(i1+i2)　　set＄layer3datasetto＄RasterDataset　　set＄layer3input1toband2　　set＄layer3input2toband5#InputathresholdvalueforIBI:　　setlayerformulatoif(i1-(i2+i3)/2)/(i1+(i2+i3)/2)+＄number5+theni1elsenull　　setlayerinput1toband1　　setlayerinput2toband2　　setlayerinput3toband34　(:)Fig.4　Flowchartofextractionofbuilt-upland.Note:imagewithsinglequotationmarkrepresentsprocessedcorrespondingimage4　SAVI、MNDWINDBI(IBI)。,,,,。IBI,,,。IBINDBI,,LandsatTM/ETM+、ASTER、SPOT-4,5CBERS-1,2。IKONOS5774:　　　,,。,ERMapper,MATLAB,。:[1]SterlingSandDucharneA.ComprehensiveDataSetofGlobalLandCoverChangeforLandSurfaceModelApplica-tions[J].GlobalBiogeochemicalCycles,2008,22(3):GB3017.[2]GuindonB,ZhangY,DillabaughC.LandsatUrbanMap-pingBasedonaCombinedSpectral-spatialMethodology[J].RemoteSensingofEnvironment,2004,92(2):218-232.[3]MasekJG,LindsayFE,GowardSN.DynamicsofUrbanGrowthinWashingtonDCMetropolitanArea,1973-1996,fromLandsatObservations[J].InternationalJournalofRemoteSensing,2000,21(18):3473-3486.[4]LoCP.LandUseMappingofHongKongfromLandsatIm-ages:AnEvaluation[J].InternationalJournalofRemoteSensing,1981,2:231-252.[5],,,.[J].,2009,11(5):694-690.[6]XianG,CraneM.AssessmentofurbangrowthintheTam-paBaywatershedusingremotesensingdata[J].RemoteSensingofEnvironment,2005,97:203-205.[7].[J].,2002,17(3):86-92.[8],.TM[J].,2000,4(2):146-150.[9].[J].,2005,10(2):223-229.[10]XuHanqiu.Extractionofurbanbuilt-uplandfeaturesfromLandsatimageryusingathematic-orientedindexcombina-tiontechnique[J].PhotogrammetricEngineeringandRe-moteSensing,2007,73(12):1381-1392.[11]XuHanqiu.Anewindexfordelineatingbuilt-uplandfea-turesinsatelliteimagery[J].InternationalJournalofRe-moteSensing,2008,29(14):4269-4276.[12]ZhaY,GaoJ,NiS.Useofnormalizeddifferencebuilt-upindexinautomaticallymappingurbanareasfromTMim-agery[J].InternationalJournalofRemoteSensing,2003,24(3):583-594.[13],,,.[J].,2005,60(1):143-150.[14],,,.[J].,2009,11(5):677-683.[15]HueteAR.ASoil-AdjustedVegetationIndex(SAVI)[J].RemoteSensingofEnvironment,1988,25(3):295-309.[16].-[J].,2000,55(6):671-678.[17].(MNDWI)[J].,2005,9(5):589-595.[18]XuHanqiu.Modificationofthenormaliseddifferencewa-terIndex(NDWI)toenhanceopenwaterfeaturesinre-motelysensedimagery[J].InternationalJournalofRe-moteSensing,2006,27(14):3025-3033.[19]RayTW.VegetationinremotesensingFAQs.ERMap-perApplications[M].Perth:ERMapperLtd.,2002,85-97.[20]JiL,ZhangL,WylieB.Analysisofdynamicthresholdsforthenormalizeddifferencewaterindex[J].Photo