ArcGIS專題操作之-空間分析建模

實(shí)驗(yàn)七、3D數(shù)據(jù)集成分析專業(yè)年級：地信071姓名：王媛媛學(xué)號：06407024一、實(shí)驗(yàn)?zāi)康呐c要求：（1）掌握DOM(數(shù)字正射影像)、DEM與DLG(數(shù)字線劃圖形)三種數(shù)據(jù)的格式與組合方式。（2）掌握利用GIS空間分析方法解決實(shí)際問題。（3）要求利用現(xiàn)有的數(shù)據(jù)進(jìn)行深層次的數(shù)據(jù)挖掘與分析，得出有價(jià)值的信息。（1）綜合利用幾種不同的數(shù)據(jù)源，求取實(shí)驗(yàn)區(qū)域中坡度大于25度的耕地面積以及繪制相關(guān)的圖件，作出必要分析。（2）結(jié)合DOM數(shù)據(jù)進(jìn)行坡度大于25度的耕地區(qū)域的三維顯示，并制作三維飛行動(dòng)畫。二、實(shí)驗(yàn)原理實(shí)驗(yàn)技術(shù)路線如下：計(jì)算25度以上的耕地面積：首先統(tǒng)計(jì)出整個(gè)區(qū)域中滿足條件的柵格的個(gè)數(shù)，柵格分辨率是5米。由幾何關(guān)系可以知道，一個(gè)傾斜平面的表面積等于它的水平投影面積除以它和水平面的夾角的余弦值，即“柵格平面面積÷cos（slope）”，由此計(jì)算出25注意：1.土地利用類型編碼為141、143、144、146的屬于耕地類型對DEM數(shù)據(jù)進(jìn)行“fill”處理以減少數(shù)據(jù)誤差，得到的坡度圖層的四個(gè)邊界的數(shù)據(jù)是不確定的，在分析中要避免使用。三、實(shí)驗(yàn)數(shù)據(jù)本次實(shí)驗(yàn)選選取我國陜北地區(qū)的3D數(shù)據(jù)，分別為DEM（柵袼分辨率5米）、土地利用線劃圖(tdly.shp)，數(shù)字正射影像(35c.img)四、實(shí)驗(yàn)內(nèi)容及步驟首先激活各個(gè)擴(kuò)展功能模塊1、原始dem添洼點(diǎn)擊ArcToobox下的SpatialAnalystTools/Hydrology/fill，將原始dem添洼，注意設(shè)置Z值為1，不設(shè)Z值可能不能執(zhí)行。2、提取坡度打開空間分析工具條，點(diǎn)擊SurfaceAnalysis/Slope，設(shè)置空間分辨率即cell為53、利用柵格計(jì)算器計(jì)算坡度大于25°以上的柵格，打開柵格計(jì)算器，輸入以下公式：con([SlopeofFill_DEM]>=25,1)，然后執(zhí)行。此方法只對坡度大于等于25的柵格賦1，其他柵格無值4、將土地利用矢量圖“以土地利用編號”為關(guān)鍵字檢索耕地土地類型打開屬性表，點(diǎn)擊右下方的option,選擇selectbyAttribute，在出現(xiàn)的對話框中，輸入圖上顯示：5、將耕地區(qū)域轉(zhuǎn)化為柵格打開空間分析工具條，點(diǎn)擊SpatialAnalyst/Covert/FeaturetoRaster，注意選擇屬性段為LU-CODE6、利用柵格計(jì)算器計(jì)算25°以上的耕地分布區(qū)域打開柵格計(jì)算器，將剛生成的耕地柵格圖乘以坡度大于25的結(jié)果圖，即可輸入公式：[polygon2]*[Calculation2]7、計(jì)算25°以上的耕地面積公式：柵格平面面積÷cos（slope）slope可能為平均坡度值644081/cos(/180*π)=7324358、三維顯示并錄制動(dòng)畫打開3DAnalyst工具條，點(diǎn)擊Arcscence，加載25°以上的耕地分布區(qū)域圖、填充的dem、35c.img打開25°以上的耕地分布區(qū)域圖及35c.img的LayerProperties/BaseHeights，選擇obtainheightforlayerfromsurface,選擇基準(zhǔn)為dem在空白處右擊選擇Animation,打開錄制對話框點(diǎn)擊錄制點(diǎn)擊錄制開始錄制開始錄制速