資源描述:
《arcgis教程:坡度》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、地理國情監(jiān)測云平臺(tái)ArcGIS教程:坡度一、了解坡度 坡度可表明表面上某個(gè)位置的最陡下坡傾斜程度??舍槍IN中的每個(gè)三角形和柵格中的每個(gè)像元計(jì)算坡度。對于不規(guī)則三角網(wǎng)(TIN),坡度為各三角形中的最大高程變化率。對于柵格,坡度為每個(gè)像元及與其相鄰的八個(gè)像元中的最大高程變化率?! ∑露让羁商崛≥斎氡砻鏂鸥瘢⒂?jì)算出包含各個(gè)像元坡度的輸出柵格。坡度值越小,地勢越平坦;坡度值越大,地勢越陡峭??墒褂冒俜直葐挝挥?jì)算輸出坡度柵格,也可以以度為單位進(jìn)行計(jì)算?! ∑露冉堑扔?5度時(shí),高程增量等于水平增量。如果以百分比形式表示,此角的坡度為100%。如果坡度接近垂直(90度),則百分比坡度接近無窮大?!?/p>
2、 坡度函數(shù)最常用于處理高程格網(wǎng),如下圖所示。較為陡峭的坡在輸出坡度地圖中顯示為紅色。二、坡度北京數(shù)字空間科技有限公司地理國情監(jiān)測云平臺(tái) 1、摘要 判斷柵格表面的各像元中的坡度(梯度或z值的最大變化率)?! ?、插圖 3、用法 ·坡度是指各像元中z值的最大變化率?! ぎ?dāng)表面z單位用其他單位而非地面x,y單位表示時(shí),需要使用Z因子對坡度計(jì)算進(jìn)行校正?! ぽ敵鲋档姆秶Q于測量單位的類型?! 烊绻麨槎?,坡度值的范圍為0至90。 §如果為高程增量,范圍為0至無窮大。平坦表面為0%,45度表面為100%,隨著表面變得越來越接近垂直,高程增量百分比將變得越來越大?! と绻苯余徲?3x3
3、窗口)中的中心像元為NoData,則輸出將為NoData。 ·如果有任何鄰域像元為NoData,則會(huì)向這些像元分配中心像元的值,然后再計(jì)算坡度?! ?、語法 Slope(in_raster,{output_measurement},{z_factor})參數(shù)說明數(shù)據(jù)類型in_raster輸入表面柵格。RasterLayeroutput_measurement(可選)確定輸出坡度數(shù)據(jù)的測量單位(度或百分比)?!?????????DEGREE?—?坡度傾角將以度為單位進(jìn)行計(jì)算?!?????????PERCENT_RISE?—?輸出高程增量百分比的關(guān)鍵字,也稱為百分比坡度。Stringz_fac
4、tor(可選)一個(gè)表面z單位中地面x,y單位的數(shù)量。z?單位與輸入表面的x,y單位不同時(shí),可使用z因子調(diào)整z單位的測量單位。計(jì)算最終輸出表面時(shí),將用z因子乘以輸入表面的z值。如果x,y單位和z單位采用相同的測量單位,則z因子為1。這是默認(rèn)值。如果x,y單位和z單位采用不同的測量單位,則必須將z因子設(shè)置為適當(dāng)?shù)囊蜃?,否則會(huì)得到錯(cuò)誤的結(jié)果。例如,如果z單位是英尺而x,y單位是米,則應(yīng)使用?z因子0.3048將z單位從英尺轉(zhuǎn)換為米(1英尺=0.3048米)。Double 5、返回值北京數(shù)字空間科技有限公司地理國情監(jiān)測云平臺(tái)名稱說明數(shù)據(jù)類型out_raster輸出坡度柵格。Raster 6、代碼實(shí)
5、例 坡度示例1(Python窗口) 在本例中將確定輸入表面柵格的坡度值。import?arcpyfrom?arcpy?import?envfrom?arcpy.sa?import?*env.workspace=?"C:/sapyexamples/data"outSlope=Slope("elevation",?"DEGREE",?0.3043)outSlope.save("C:/sapyexamples/output/outslope01") 坡度示例2(獨(dú)立腳本) 在本例中將確定輸入表面柵格的坡度值。#Name:_Ex_02.py#Description:Identifiesther
6、ateofmaximumchange#?????????????inz-valuefromeachcell.#Requirements:SpatialAnalystExtension#Importsystemmodulesimport?arcpyfrom?arcpy?import?envfrom?arcpy.sa?import?*#Setenvironmentsettingsenv.workspace=?"C:/sapyexamples/data"#SetlocalvariablesinRaster=?"elevation"outMeasurement=?"DEGREE"zFactor=?0.
7、3043#CheckouttheArcGISSpatialAnalystextensionlicensearcpy.CheckOutExtension("Spatial")#ExecuteSlopeoutSlope=Slope(inRaster,outMeasurement,zFactor)#SavetheoutputoutSlope.save("C:/sapyexamples/output/ou