土地利用解译及精度评价操作流程

1、土地利用解译操作流程准备工作：在空间数据云平台（）网站上下载1992年-2012年遥感影像数据（注意：为确保郑州市在图幅范围内，检索时候输入与郑州市相近的经纬度坐标，控制云量大小，确保数据时间在5-10月份之间）考虑到下载的数据质量比较好，数据都自带坐标且云量基本都为0,省去配准和大气校正步骤，直接加载已知 1995年的土地分类矢量图，并对遥感图像进行裁剪。1.由于下载的数据都是分波段的，因此先进行波段合成，具体操作如下：打开，打开 File->Open external file->Landsat->GeoTIFF，选择 7 个 TIF 格式文件打开，选择 RGB显示，勾选

2、 742 波段，Load RGB,打开 Basic Tools->Layer Stacking点击Import File，输入要合成的波段，选择742波段进行合成并输出一个文件（名为1992）2加载1995年郑州市土地利用矢量图，具体操作如下：打开File->Open Vetor Filel浏览到要打开的 SHP格式文件，点击打开，点OK选择刚刚打开的SHP图层并加载到之前打开的RGBS层，如下图选择Load VectDisplay #1:Sal Victor WindowDisplay *1Display #2New Vector Wiixd<owCancel得到如下结果:

3、',1;O #1 Zoom 14. .EE.- 現Q #1 Scroll (0.03307) 口 丨甸丨hi|e Overlay Er nanc?3.根据加载的矢量图对遥感影像进行不规则裁剪，具体操作如下:忘 1 Vector Parameters :匚ureQ.File) Mode Edit Options HelpWindow / Imige V Scroll V Zoom打开图像感兴趣区域（ROI）如下所图示:Hide LayerAnnotation,.,Classification.Contour Line氣“Vecto rsiihDensity Slice.Grid Line

4、s.'f 斑阻;Layw （阳nd 1:斷TTF" Layw：阴nd| Fig fSwrl呂y Enhance Took Windowijhr" «u s #1 ROI ToolFile RO!_Tpe Options HelpW indtjw : 0 TmgeSei* olLOffRDT KameCftlcrFiit&lsFolygonsrBjafi.011 #1lUd00/0oA*EVF： L*r*r：刖VhiU1, 129, 200661/1. 129j£New Iegi. on'jctoDelete TartStitE |

5、 Fiwl | 肌1蝕 ROISslftct AllHid» BOIs Show RDIl点击File->Subset Data via ROIs，对图像进行裁剪RJOIiROIHid* EDIiAllFili ROI Type Options Help0 列 ROI Tool输出文件名为1992裁剪，如图所示:加载裁剪图像，显示裁剪结果:Save ROk.Restore ROIs.Sublet D目饲 via RO鳶.J_"_Output R.Ols to ASCII.Carre elExport ROls to EVF.Export RO1& to Sha

6、pefile.Export ROIs to rr-D Vkualizer.File OptionsEOT Uak (Layer OEDI iMk (Layer CBond 1:517.1BOIOLayer1 :B4.1MI IM 4 sk (Layer 0 如d 】：应.T12B3 乂 1364 CB吐电)EECD1 E pl ay #2 t-La ad R.GE1tET.TIF)11E4.TIF)1：E2.TIF)CBATLdCEiCEsnlAvailable Bdnds LiftD imsGray Scle HGD Color日甸1的3裁剪- 口 KOI Mask (Layer 田and

7、1' IT :“ HOI Mask (Lftyer (Band L: E4 L- ROI Lk訥±1 零 iTifcE-0 1992； 匚砒縉- 匚蝉出-3 Map InUJansf*0.SD0££5.将剪裁图像另存为.img格式文件，命名为1992裁剪.img。QuickMspPreferences.Cancel5吕y Diply Group. Restore Display Group.Create Mew Vector Layer.| File | Overlay Enhance Tools ?*/indowPo5tscript File.Im ag

8、e File?5Save Image 氏5Savt Zoom A5vS 丰 1 RR (ROI Msk (Layer fBand "一.匕,门；.1十："16.将1992裁剪.img和1992 土地利用矢量文件加载到ArcMap中，如下图所示:7.新建Shpfile图层，定义坐标系（与遥感影像坐标系保持一致），打开编辑和捕捉工具, 选择 editor->snapping->snapping toolbar女口下图所示：8.比对1994年与1995年的土地利用变化， 将增加变化类型提取矢量化， 并编辑属性表，打 开新建图层属性表，并添加字段信息，如下图所示1-3代

9、表该增加部分是由类型 1用地转化 成类型3用地的。9.增加部分数字化提取结果如下图所示:（1）在提取增加信息的过程中，由于变化部分图斑面积普遍较小，矢量化过程中记得打 开捕捉工具，避免两个斑块离得太近造成拓扑错误；（2）数字化过程中可以利用查询工具快速查找原本地类编码，便于数字化过程中属性表 变化信息编辑；（3）将1995年土地利用属性表中的地类编码进行合并，如将123、113等全部合并为第1类，将51、52、53等全部合并为第5类，整理属性表信息，删除冗余部分，主要属性信息应包括ID、地类编码、土地面积、变化信息、变化面积等基本信息。说明：对增加部分进行提取时，是以1995年为基准进行推算提

10、取的，对于1995年以前的年份所提取出来的增加信息实际上来讲是返减信息，如1-3类型变化，按照时间顺序来看，应该是3-1的变化，而对于1995年后的数据则无此问题。1994年土地利用图10. 对图层zhengzhou1995和B-1994年变化信息图层做合并操作，得到层打开 ArcToolbox->Analysis Tools->Overlay->Union分析工具歸 ArcToolbosl+i ® 30 Analyst Tools曰 Analysis Took+1 窃弓忍 OverlayIdentity气 Interstct气 Spatial Join豕 Syii

11、rirn'etri cal Differen匚巳 Update_n r同Irmt Fa-H.ti：3r6-=* - Lp aaassUnionCorr pule s a 攣 omein c. union of the Input fest uro e.Ill<LIT TrJ ftH-iQicfi Pcpti onr&l)i|ffeters.w« Hdn Tilt*£riviunaiivB<iL,.QOijlpijl Fftlure Cl&ssC lJikflT kii 3 i tririor 皿血 咄塀土地刊用芸It询)f 氓All fe

12、ihirM gndtii«ir atiributoc will be -*/riner to th& output fea-ur class同如上图所示，将 B-1994年变化图层和图层 zhe ngzhou1995添加到In put Features中，输出结果命名为1994，得到结果如下图所示:Jfa£匚opl7 °总j.裂说明：对于合并后属性信息，删除矢量化过程中造成的明显拓扑错误信息，删除无用字段并添加 1994_class新字段信息来表明1994年新增加的地类信息，保留1995-1994年变化信息和1995_class字段信息，方便变化前信息与变

13、化后信息对比，查找变化信息。（注意：前期可以跳过 11 13步的获取土地转移矩阵部分，直接做14步）11. 分别对1994年、1995年土地利用变化图层做数据融合操作打开 ArcToolbox，选择 Data Man ageme nt Tools ->Ge neralizati on -'Dissolve 工具。In put Feature 选择要融合的图层，Output Feature Class选择输出结果存储的位置及名称，Dissolve Field（s）选择土地利用类型字段（如Type1995），然后勾选Creat multipartfeatures选项，点击 OK完成，

14、如下图所示：重复此过程，对1995年数据也进行融合，杳看屈性表和融合后图层变化如下图:a(=7口0灯&1日吨合：FIDShape *1?94 g Lasfi |0Fn-lrgj-E1占F心塔sn二2Polyga-n3Polre&n4斗512.对融合后数据进行叠置分析在ArcMap中打开两个时相融合后的数据， 在ArcToolbox中选择 Analysis Tools ->Overlay ->1 ntersect 工具，In put Features 选择两个时相的图层， Output Feature Class 选择叠 加结果存储的位置，命名为 A,其余选项可以忽略

15、，单击【 OK完成，如图所示：查看A图层属性表信息，删除无用字段，添加面积字段Area，并在属性表中Area字段上单击右键选择 Calculate Geometry ，在打开的 Calculate Geometry对话框中，Property选择Area , Units选择要使用的面积单位，单击【OKI完成图斑面积计算，属性信息如图所示：Table >目7崔毬1- 2 XFIDStXBL 口 E *19f94 claES |1995 ic loss |Area0F-olF-gon117 17 033MTP-clygon12.631133P-olygan14.S99L1S2P-Dly

16、3;DH15».2S75814P-DLy EC-n1.二总右空日SP-cLfedr223-1. 234731GP-olyE-DE：31丁P-DljiED32.05.1&SGPoly fic-n.2g三P-clys 口 et34.02739310p4Z-ly£DQ411.23092 S11P-Dlyg-DQ42.0006 扫12PdItesel龙a.03501113P-qIiedc町444.66a4414P-olygDci511. £4£93123p-oiygDi33196. 7 3311313.制作转移矩阵保存ARCGIS中所得结果，用Excel

17、打开文件，并另存为.xls格式文件。In.ixj"vI1994_clas=：L995_classArea(kui2)11717.03390678500120.63116252142140. 899U4S7782159. 2875S114&58210016627995662234.2947314926031CL 1349994=943932Q 05188601310337.17454009&75340.0273945436641k 23092828902420. 0006993150生30. 025013871074444.66874920220511.64693054

18、55755196. 795112&8600.在Excel中选中所有数据（不要点左上角，只选择有效数据），点击【插入】选项卡，选择 【数据透视表】|【数据透视表】，点击【确定】。在打开的数据透视表中按图示将字段拖入相应区域。新聲d宇段列表龄睡动区域分祈分享a备粉言工具得到结果如下图所示:1994-1995年丄地村用馬申求和项:Are, Om21993 cla££19加 匚1昭亍12S45总计1717. 03390680. 631182921a suots9. 29"5811<7727, B5L78S320. J166379563C 2M7314934.

19、3113501930. 13499914C, Q518B6Q137.174540La QZ73945U1. 36BS2Q15141.Z399333B9a QQQ旳痔0. QZ5013B7LH., 66S7-19215, 9253303&1.64693051&196. "95112719氐 442032启计720. 0633331SL 78199557.190553971IL叨525進206. 08269 3日1013. 9L'J399I#匚：'E DisJ *zkp时宰芹冋谢输出耍素类襄创建泊料色磬常巨缶慕M 妾素蛙口弓prEEl團疑I十李咸闫坦甩匣幕

20、幻职筲O-工員萃砌I14.遥感解译结果精度验证(1 )对1999_B数据进行融合处理打开 ArcToolbox，选择 Data Man ageme nt Tools ->Ge neralizati on-'Dissolve 工具。In putFeature选择要融合的图层1999_B， Output Feature Class选择输出结果存储的位置及名称，Dissolve Field(s)选择土地利用类型字段(如 B_class )，然后勾选Creat multipart features选项，点击 OK完成，如下图所示：削刍*一 = j 回完成后查看新生成文件 1999_B_D

21、is的属性表，发现全部融合为5类，如下图所示:1FIDShape B clc&s5一1T22Z3'd兗"V5(2) 在ArcGIS中生成随机点->要素类-> 创建随机点如下图所示在工具箱里打开数据管理工具E卑愛詳气创逢渔网弋创韋厢机点气更新莖记要素类%计畀默认XY客走 “计算默认空闫將网盍弓I追尬主记豎素娄如下图填写输出路径、输出名称、约束要素类选择上一步融合处理结果，长整型改成50,下边的线性单位可设为0m。添加XY坐标如下图所示打开工具箱中的数据管理工具-> 要素-> 添加XY坐标ArcToolbox X1*1 0索引Ei抵表0 藍素%XY

22、強*様复几何 %切分、爲就角i 气取酒茨分割J在折磁分割铁 卢注庶处缶割裟 匚复倒费養 r勢部件至单酗半 芦4聂小边界几何 ggjng0 拠才点異幕一 要養包帝柜形轻面要董折君转栏声,委養转点 *.賽養徳 气要薰靖面 工3DZ值 气面显曰0要殊“:创建取属注旺要素类 气创建渔网W占4呼* 一4* »卜.选择上一步生成的图层 1999_B_R,并点击确定。结果如下图所示:& 3 Laym a xmmjlr* 1W_R15岳 a lSs- 13KLJ尿19舒卫£ g甘海观QU 口& HtK£HMfl«klljy!|_gBt*. G ROI h

23、usk 3尹釜岂 & i-UL i'.!a!： |Lzv亡a打开1999_B_R属性表可以发现增加了 X、Y坐标值如下图:表 X国T龟殆野期x（4）结合谷歌地球进行校验并调整地类将上一步得到的格式文件用Excel方式打开，可以得到每种地类 20个随机点的X、Y坐标值。打开谷歌地球，打开工具-> 选项->3D视图，勾选“通用横轴墨卡托投影”g Gcx?号 1匸 ±63，刃亂视團 高通無歹亦菟 导翫馬感-纹理颜邑増强鱼（16位）*頁翔鱼（32位 匡压缩各向貝悄过濾©关囲 a中 e高标签/團标尺寸小倬中C'大OpenGL。DirectX使用安全

24、模式显示纬度/经度 已小数度数度、尔秒 度、一进制分*通用横轴罢卡托投影4i4i 皿£度量单位9系统默认 其尺、英里米、公里宇休述择3D宇体提升高展f会同时缩京3D建筑和树木）；1(0. 01 - 3)使用高清地形團如要快速切换分辨率并呈现地形團，请将苴停用）。 団使用3D團濠（停用则采用旧版3D建航环境使用潼頁的环哒築（实验性功能）总览图地图尺寸；小Q大绻敝关奚：无限 i：iQi：无限|恢复默唾贸匚确定|取匸|匚应用（5 ）对比校验，调整地类点击添加地标按钮分别将不同地类的 X、Y坐标如下图所示输入 Google地球中，并放大每一 个标记后点，查看标记后点所处地类是否与之前随机点地

25、类相一致，如果有变化则在Excel表格中对地类信息进行调整编辑（无变化不用调整，且只调整地类，如1调为2，坐标值无需调整）3820135.27327000000335030259550000003300103.11754000000378960.811530000003827148,539530000003834568.804340000003 323040.S62210000003321796.538610000003857376. 82302000000382727*2971100000038E3Y75. 955540000003542980.896370000003809406.8140

26、30000003843092.7GL090000003856L48. 453570000003853853.71292000000打开感兴趣区,东进：680845. 38 以东 北进：3799222. 92米以.北说明样式/濒色视图海拔高度添加遴接添抑團片. 全部调整完毕，将新变化地类全部进行重新归类，如之前2地类变为1地类，则将其直接移至1地类部分，每个地类单独按如下图格式粘贴到记事本中，存为.TXT格式：1745738.276391000001776313.02375SOOOOO1742167,21186000001742279.27173000001725470. 8149T&000001775694. 2601870000017595%. 6CS1650000017G0851. &3555300000170S390. 051685000001755463* 084082000001703834.91258000001769816.SS1352000001739171.