ENVI遥感影像变化检测

1、t宝此电星!星M看套J1A押.躯亚府01射基章*电里厩增祀公也行更ELL过LL堂型府 剧的本力但翱JUMVI 士化窃割由金屯14股胖依礼!->I 由科暗E&总* 1clim|vH霭+府原堵Mjpd Uue DkfktmZ脂用吗曳器U岳JK./ Wf：g，Zpi3 脩£*但KW1 £ g” q%Lni、w«nE1>-Ml CluW r fl F：色卡区I磊，香也色甘*整X.标MIHMtH用埼丽国 »Mfti. 川四十， 性闫 "A Hi fi-m |. kt ft *4 '? 嫌在 i'，*K 情: I i M

2、 国之.“4 iftn u带白山也.U . 1，f 4 * 瞥4t施 1'M Ii，u.一 以，用，*£岫* IftW *.-( WflLll14L上 博戚牌>iL m+“ F”，。J H 期灯* i >i w F n lb f JR trtra fi TMUftfj f aMi4« r fij/ Meji fi f; p.r-i r utv » 期"吧ii睡it*依r ni<岫.flMtMTr”, .iHMMh MMBWiMiHir MNMM加4.歌0&由 主* tttnAUVH燃靶号修巾 ，i/+M* |i盯，咚.i

3、liIH in启骨" di的曲 MAMI «, JiAilihtmtU 1 t，*r£»4X* W，5 -rHk 珠鞋制基幅-H，1；TFRT1I1M事，中串上媪， Wl； Itt. r!. ULflHIihh隔号叶 u =fhm百金也也 ilt 检2/西诲匕 I li7 hl|i*.'tt*hHB«41帆 比4理.哨看*曙禺事喇/蛇II电I tfetAH rAfUR/ifiHItlti *比总物量.tlLJMKM fl Hitter 此 * >*i»tFJ 5候.t fttfttilr* 打| E Wf，*nW IE宽

4、电的H iff> <.V 4 n ii ith-&W1« i.修田喷油胖IltHT L耳畦计巾忖以 冷亭#|J零嗡岫性侪* J HJtR-N &、*1*借串修闻*曼i Kjr>ML*布懒轿0孽mnd工震小千/明的H U.AHu M曰*乜1整地方基，口得 FMH由"IM 1 "fin AtlilULinh4tl.1 f iLff. *丘也力情，善讣L$ ,*占1； Yi* qrk 也的!1/嘴人群阖，” *%。珥胃五雪中皿*呷曲If T W J K444 mu ；*而 :I W：。值卜席H，市号号时l，ft I南花*结口匕WT科 的

5、*，用时1於才田41«尸岬承网 *用，*号I T-患平易方*HmaMbMdMm碑惠事*打i-«4Wh盘* ml iLM rtiflvti里ml xrirtu*的£皿廿电*北帆、1,川*置他/11是”壶工11 4vl4* nrHMUli.il, iaithrAlUlAA 1 * 将刘乩LI啊。空吃11，鹿育勒11盯工11且.叫口卜染3暨靶111414由*: 口个个乾.直一寸龄*"壁包,ftftW 电宙W曲 一”!»加M；叫 _®_y .i JDft* U lTl TM5 1114J偷 JN，|H4 If 的 TM? 'Wftt,

6、itnutfi 总修 普 协/inf Jin? Il序第，,部ItEL<rLjMH ja0网停# Mrflt 中点依餐一露 塞-|非!！：，曾亶4代：,科卡卜力田XE*"匕£*rih YE稣*P瞿4|，妣，网嘴淡"M争&” 照明，I, %:危G.|tf!糊干, fllt-iLAMii *空M带11*睡1 *qMBMlRr 新厦序41|睢| M lift '"43如业.立雷篁->t 1*1-人*功殳0*域环*. «J|1L- UM4l.丸 曾整也值a魏miiH w,安，廿*.爬，*小ni竞:1号牌宽n to v M V

7、J0.111 *、!¥仁科丝加叱色.魄»1111立，旧叫如UK*的> | > 萨r*川，里用丁明 飞 ，隼/，|t|r|Vb,君才中/中宜胃|度彳¥ IfiTJWr 场勺比寸隽1； 4用值，：，梏由r M-ij« *规Pt3fiH木工小，纯,HidrH除海*/ 方1*1 -七甘口*m 恰用片中0哨立 MH七.*.心制导仲片-瞋*， 1.森林开采监测打开实习数据0-森林开采监测下的实习数据。Compute Difference Map选择 basic tools/change detection/ Compute Difference Maph

8、另选择原始的影像 july_06 与 july_00,在弹出的 Compute Difference Map inputparameters0 口 下，查看 define class thresholds,no chang表示没有变化,change(-1表示减少，change(+1表示增加;其他默认选项不变勾选normalize data range0-1选择输出路径与文件名为 com_diff。选择 classification/post classification/classification to vector在输入图层中选择上 一步生成的结果，弹出窗口中选择全部，保存路径生成结果，转

9、化为矢量。（由于耗时过多，故可以不做? Image Difference打开ENVI Zoom 4.8，将原始的影像导入到其中，在ENVI Zoom窗口下的 toolbox 中选择 image change单出 image change detection勺对话框，将 timeclassification image file选择为 00年影像，点击 OK,time2 classificationimage file中选择06年影像数据，点击OK,选择下一步，保持默认设置，选择下一 步，选择image difference选择下步，选择 difference offwtin IndK* 端域通

10、倏直|脾it)*口* 播 wdlhLerzlcl me i'i 1. eyj I ?li：c 山，j 切，j :!1 1' i 卜力.£曰由神牌9之伴fl.点i&HnMh%*.rr作on up»崎t. 一 小 i诚的此瞿号it察，得山帝t£比怦 ptni!)eniBs/smbQosT-e«t wIhb 'F*鼻 tni<u? vnlue» value field 卜 Tiff出门 Ilim,朝 LiR Jecteafi# %*苴，,；巴.tig itxreasc A r：.'i£. q：h

11、<r 中：f 直打，B"性卜 »rt tw 卜 "|%信-.4:和坏叶心心和三仃泳折，L M地府IfeF至兰曲£¥；刚* SFFAK i Hji> Ihfln'l'i 'n电界Uuhm/XDYI打片，耕30位.2(m 的声”肺惬时第四酸依 聿旺 NDV1 HML :且打理 NDVI UkMriqvts d“dbind iZMW ROTI* 儿 RI 照什融通符型1Ams JQWMfil. 小值医修口八最大笆3用L切RI上间的M曜“分:味一： li him 口* m »«hH doib *

12、ij RDli.” tJ JH . F f 2OU2 H 注棒2042 ROL裨班弊 i m肆惊小摩法与支力范.密林 /*bt 碧EM Ttid IPtttti >4H> ff. rcut F41>中价M"氏两个时旭能否像欧丽鼻,四波片.在>Lwl CUf41 x>»出口it个存itm”一 ：a人儿向晴帝口呼此段.逸推*1曰ii> telrtit白M逸琼，星盘港.而自动d抵ft广电运里,叁口蔺才找n十仲子点站变种干/杯为GB 与用匕也配it肃 曷为1方).凿青小片ul单曲mt.用用控制在tih：剧蚓去 杭二£畅剖点.“寸也 干忙

13、心内，晨大祀，为用基为1. 'H.I ntil也A Owve PdKEivi*u»(w诙枢*在重把常*方式中“RN.应 由nrtt, M"特件*th玷H，rii. DiBTl«y Ktl'lf -j.-DiH 曲i孔 词 HtJYl - ti»e » 2» 如日 Tn * 工F & J：'N1'f pfii Rr<wie* *tii*r*ci iw $u*t thLr*. *J 开翼币国的,方国.电4国国后nA的珈七t事的胃卢. aM*X« 4 D. T拚般上的阜蛤W。纪融打样刈

14、宙横州利，口好他门.血星卜的 ftMHftv -aSI2,落前胃上JK京中国挂 gr】MT|MHn， Cltu ibes. |Et9|r . *，二时至F- '' J, tail xtun il.7 达.身制为一L7W51，-o,iia-aK- f,gC.N ."：' I匕 DEpkr "r rft J,L 'I', flp l*J .小 Gt % CL»X1 f ICBtion Posl CI>99cf icalion- jA II .ti AMdd/l 5ffl “ ilff Kernel 5i:e 陆 5+5r

15、士 用小fi.3.耕地变化监测Change Detection StatisticsThematic Change首先加载两幅影像ag_08_maxlike.img,ag_09_maxlike.img观察影彳t信息，我们发现,其中黄色的是休耕地，绿色的是有农作物覆盖的农田农田，黄褐色的是留荏地，蓝色是水体，灰色是裸岩。操作：BasicTools f Change Detection fChange Detection StatisfCs选择08年数据，然后选择09年数据，因为两幅影像的分类名称一致，系统会自动对应。点击 ok。在 change detection statistics outp

16、u中选中pixels、percent、area选择掩膜输出，选择保存位置。生成统计结果，选才? area选 项，在 option中选择单位变换为 square krnt打开ENVI Zoom,把08和09年的数据加载进去，与前面一样，也可以利用Portal 进行查看两个年份的土地覆盖类型的变化，另外，点击crosshairs移动鼠标可以查看两个年份是否发生了土地类型的变化；点击crosshairs移动鼠标可以查看两个年份是否发生了土地类型的变化。Toolbox-workflows -Thematic Change别输入两个时相的影像数据，接着点击Next,然后进行结果影像进行平滑和聚类处理，平

17、滑处理的参数与上面一样，设置为3;聚类处理的参数设置为5。把结果影像和统计结果输出图2将影像在zoom中生成结果后如上图2所示，影像上红色区域变化加大，表示的 是两幅影像时间差之间被砍伐的森林区域，占大部分的面积，蓝色部分表示该时间段 内植被增长的区域，但是面积相比较被砍伐的区域来说小很多，黑色部分是代表没有 变化的区域。在arcmap中将最后结果输出，如下图3所示，很清楚的可以看清楚红色 表示00年到06年植被大量减少的区域，绿色表示植被增加的区域，浅绿色表示植被 覆盖没有发生变化的区域。图3 Arc map制图输出结果2、林地病虫害遥感动态监测图4 ndvi影像如上图4所示,分别是2002

18、年与2007年经过NDVI计算的影像，从上面的影像中可以清楚的看到 2007年的NDVI影像上有些区域变成了灰色，也就是植被收到病虫害的区域,NDVI影像上亮度很高的植物量表示很大。图5如上图5所示,前者为未经过Density Slice处理的结果，后者为经过Density Slice 处理的结果，在最后生成的结果中(图5中右图绿色表示植物被破坏严重区域，这种区 域面积很大，而红色的斑点区域都是破坏严重程度一般的区域，这种区域占面积不很大，只有不大的区域，蓝色区域表示的是基本没有发生大变化的地方。3、耕地变化监测Area (Square Kmg fields s field B soil wa

19、ter b ground R Total C Total Unclassified 0.00 0.00 0.00 0.00 0.00 0.00 0.00 green fields 676.28 365.15 276.10 0.68 50.40 1368.61 1368.61 S_field 92.23 114.28 73.91 0.01 20.24 300.67 300.67 B_soil 49.57 90.61 491.67 1.30 249.44 882.58 882.58 water 0.32 0.62 0.53 340.42 0.71 342.60 342.60 B_ground 2.

20、99 11.01 31.04 0.72 1559.79 1605.54 1605.54 Class Total 821.38 581.66 873.25 343.11 1880.59 0.00 0.00 C_ Changes 145.10 467.38 381.59 2.69 320.79 0.00 0.00 I_Difference 547.23 -281.00 9.32 -0.51 -275.04 0.00 0.00横向表示的是08年的各土地类型的面积数据，纵向的是09年各土地类型对应 的面积，比较横坐标与纵坐标，发现08年的stubble field ft 09年有90.61Km度成 b

21、are ground/roc% Class Total行表示08年每个类别中的总面积，Class Total列表示 09年每个类别中的总面积，Class Chang续示的是对应类别改变为其它土地类型的 面积;Image Difference表示的是08和09年同种地物面积的差值。对角线上08、09没有发生变化的，非对角线上的发生变化，对列方向最终年份的统计。图6如上图6所示，每种颜色分别代表一种变化类型，如green巾elds变化成stubblefield就有一种颜色表示，而影像中大部分的黑色代表地物类型没有变化的区域。下 图为arcmap出图结果，从图中可以更好的看清楚结果。存在问题与解决办

22、法1、此次实习虽然内容并不复杂，但是中间过程比较多， 需要认真一步接一步进行，需 要搞清楚每一步的目的与原理，否则会将参数设置 错误，导致结果发生错误。2、林地病虫害遥感动态监测最后在进行 Density Slice 处理时，需要设置 一个阈值，是对影像的拉伸窗口的工具中进行的，不知道这个 拉伸世道什么程度就好，是要 将手病虫害的区域拉伸凸显出来为止吗？也不知道 这一步究竟为什么要这样做，做到怎么样 的程度才是好的效果，只是根据自己的理解将手病虫害区域拉伸，取其值-0.512。3、本次操作需要使用到ENVI Zoom 4.8工具，由于以前也没使用该工具，所以第一次使用有些难度，但总体感觉跟ENVI还是有一些相像的，但为什么一定要到 ENVI Zoom 4.8中处理而不是在 ENVI中使用呢？ 4、对影像进行自动配准不是很了解，包括其原理、参数设置， 在什么样的情况下才能 进行自动配准，这个算法是怎么实现的， 那可不可以之前 的影像配准都使用自动配准？是不 是进行自动配准还有什么条件？总结1.根据这次实