ArcGIS-在DEM的基础上划分小流域

1、ArcGIS10.0在DEM勺基础上划分小流域根据deMI取河流网络，计算流水累积量、流向、根据指定的流域面积大小自动划分流域1.加载DEM打开ArcMap,加载数据DE嗷据文件如果DE汕有很多异常值，可能刚加载进来的时候是全灰或者全白，这个时候需要右键t图层属性,进行调整。2.FillDEMg面存在着一些凹陷的区域（DE冲身是插值计算的，很难与现实情况完全符合）由于这些区域异常低值存在，使得该区域在进行水流流向计算时得到不合理的水流方向。因此，应该首先对原始DE嗷据进行洼地填充，得到无洼地的DEIM打开ArcToolbox中的水文分析模块。使用SpatialAnalystToolsTHydr

2、ologytFill如果没有水文分析模块，请用CustomizerExtensionsTSpatialAnalyst添加扩展模块。?&&CjonversionTools?&dDataInteroperabilityTookDataManagementToolsEditingToolsGeocodingToolsGeortatistical:2尺aAnalystToolsLinearReferencingToolsMultidimensionToolsNetworkAnalystToolsParcelFabric?筋?TookSchematicsToolsServerTo

3、ols?SpatialAnalystTools霉Conditional甩Density&Distance屉Extraction3Generalization2?d7Groundwater裁Hydrology卜-%Basin?：0卜FlowAccumulation理卜FlowDirectionh衣*FlowLength卜气Sink卜? ? 5-.? ? E.Stream Order-%SnapPourPoint卜.$StreamLink卜旅、StreamtoFeature?EE-SBFillinputsurfaceraster|DEM.mgOutputsurfacerasterF:Hes

4、eadikTEMPzhangyaafa-igningxiaDEMFIZlimitOptionalCancat-rorm&PShawHeA>>3.FlowDirectionArcGIS中的水流方向利用D8算法计算。针对每一个栅格，将其高程与周围八个栅格进行比较，得到水流方向。最终得到的结果是一个栅格图，每个象元上的值的含义如下图所示。如果某象元算得的水流方向是正东方向，则值为1?如果西北，则为128,以此类推。注意，这里要用的DE汜经不是原始的二步生成的结果DEM（文件命名不能超过使用SpatialAnalystToolsfHydrologyDEM了，需要用经过Fill后的D

5、EM也就是第13个字符？好像是）fFlowDirection啊ffiPIGeocodingToolsj?GeostatisticalAnalystTools±令LinearReferencingTools田A1MultidimensionTools±?NetworkAnalystToolsS)商IParcelFabricTools习野SchematicsTools由令ServerTools令SpatialAnalystToolsl±氐Conditional庄aaDensity由邑？Distancel±Extraction?霉Generalization色

6、甩Groundwater曰Hydrologyl气BasinAllFlowAccumulationFlowDirection.?_%FlowLength一冬SinkSnapPourPoint?-气StreamLink?-气StreamOrderStreamtoFeature*.)*?飞Watershed>Interpolation)-.LocalS*:?M叩Algebra4.FlowAccumulationFlowAccumulation由FlowDirection数据计算而来。每一i个栅格FlowAccumulation的值代表着其上游有多少个栅格的FlowDirection最终汇流经过

7、该栅格。一般而言，计算出来的Accumulation的数值数值越大，代表越有可能是河谷。使用SpatialAnalystToolsfHydrologyfFlowAccumulationg-QlGeostatisticalAnalystToolsj令LinearReferencingTools±j令MultidimensionToolsTNetworkAnalystTools囱令ParcelFabricTools2)ISchematicsTools士野ServerToolsj野SpatialAnalystTools砂Conditional?(2-电DensityDistancei自Ex

8、tractionl±GeneralizationiI田电GroundwateryHydrology?iABasin-气FillFlowAccumulation?一氏？FlowDirection?"飞FlowLength一气Sink?-%SnapPourPoint%StreamLink?;?、，StreamOrder'、StreamtoFeatureWatershed由Interpolation也电Localipr屉MapAlgebrag)乘MathFlowAccumulationoErvrMFFtEShowHgip>>5.河网计算也就是说，如果在某个栅格

9、的上游有 定义为河流。由于本文用的栅格是当FlowAccumulation达到一定值的时候，就会产生地表水流,格那么所有那些Accumulaten大于阈值的栅格就是潜在的水流路径，由这些水流与定义的汇路径构成的网络，就是河网。这个阈值是需要人为设定的，10000个栅格的水流会流经这个栅格，则将这个栅30米X30米，也就是说，如果某个栅格点，其上游的汇水面积超过9平方公里，则认为这里是河流使用SpatialAnalystToolsfMapAlgebrafRasterCalculatorT NetworkAnalyst ToolsGeostatisticalAnalystToolsj文寸Linea

10、rReferencingTools±令MultidimensionToolsAnalystTools囱令ParcelFabricTools2)商ISchematicsTools±0ServerTools自野Spatial阿Conditional?电DensityDistancei自Extraction&?电GeneralizationiI田？，Groundwater商HydrologyE%?Interpolation由屉Local?心“叩Algebra(3?甩Math±屉MultivariateIT辛Neighborhood审心OverlayIF?电Ras

11、terCreation国？甩Reclass庄电？SolarRadiation审？心Surface曲区罔ZonalM令SpatialStatisticsTools0令TrackingAnalystToolsRasterCalculatorMapAlgebraexpressionLayersandvariables?HowAccum?flowdir<>DEMFillConditionalConPickConCFIowAccunT10000fl)SetNulMathinnrnMnAbsExpOutputrasterF:ResearchTEMPzhangyaofangningxjaxsSt

12、reamnetCancel|Envronments.ShowHelp?QUri山d*ArMi?-占Urtttled-ArcMap-ArcJnfoFileEditViewBookmArksInsertSelectionGeoprocessingCustomiieWindowsHelp+1:250.mrtMjfiEditorTTableOf匚onttntslayersSSF:RetearchTEr/lPzh;jngyaoaSitreamnet1Floi/vAccumValueHigh:4.18仔7E+007Low0口flowdir1124QJ2_430DEMFiHValueHrah:3537Htr

13、r1B3685.8543836041.J03Mete6.StreamlinkStreamlink记录着河网中节点之间的连接信息。Streamlink的每条弧段要么连接着两个作为出水点或汇合点的结点，要么连接着作为出水点的结点和河网起始点。河网节点是计算流域的输入参数。使用SpatialAnalystToolsfHydrologyfStreamLinkt)-QlSchematicsToolsj)?ServerToolsT射ISpatialAnalystTools(5*;ConditionalE甩Density国屉Di&tance中苍ExtractionGeneralization阿Gro

14、undwaterQ甩Hydrology!?一叫aBasinj%FillFlowAccumulation卜-冬FlowDirection卜%FlowLength|ASink$?%SnapPourPointStreamLink;StreamOrder!jr严、StreamtoFeature!?-气Watershed(TInterpolation由$Local阿？MapAlgebra审电Math国Multivariate由心Neighborhood由屉Overlay手屉RasterCreation色乐ReclassStreamLinkinputstreamTBaier|strearrmeiInput

15、flowdirectionrasterflowdirOutputrasterF:ResearGhTEMPjhangydofangYiinAciaftreamlinkCancelEnvrDnrnentsShowHet>>7.Watershed先确定一个出水点，然后结合水流方向数据，分析搜索出该出水点上游所有流过该出水口的栅格,化为这个出水点之上的流域。这里，我们将第6步算出来的河流节点作为这一步的出水点数据。当然，这样的话计算出的流域将会很破碎，而且得不到流域之间的隶属关系如果需要得到更为精确的流域，那么就需要人工去判读出水点。使用SpatialAnalystToolsfHydrol

16、ogyfWatershed田翻ISchematicsToolsj?ServerToolsT?SpatialAnalystTools田ConditionalgDensity匡屉DistancegmExtraction由？电Generalisation焉电Groundwater"?E氐HydrologyIBasinj|?%Fill"？气FlowAccumulation卜-气FlowDirectionFlowLengthSinkSnapPourPointStreamW-iWatershe:LCLinkStreamOrderStreamtoFeature/L气L，*SInterpo

17、lation?心Local!?MapAlgebra每甩Math庄岭？MultivariateffiNeighborhood时屉Overlayi4霉RasterCreationE令ReclassI CD IIS I'%WatershedInputflowdirectionrasterpowdirInputrasterorfeaturepourpointdata|streamlinkPourpointfield(optional)VALUEOutputrasterF:flesearchTEMPvihangydofang'ningxjaA3tershedOKCancai|Envron

18、ments.ShowHelpFlc*fi?AccuiT14?eaS42Edrtcn口 - ArcMip Ardmf QSeltcticnGecpAcefEingiCutTomizeWindowH*lpFileEditVlcvvBsckrnihcEImeit：Zt"dbleOfCiantenti:J县|I.耳Lav*rtB因FARjSB0rrhTEMPzhengQr?gdLow:1回：VAIU?High;?27tLciw:口1TtrcotrlinkValutKgr:3271rkAnaHetT'1Higr;A16117A017917?4.BA+01A£|3.59&

19、;为了显示的需要，我们还需要将河网加进来到这一步，就把所有的集水区划分出来了。以显示每个小流域隶属哪一段河道。8.将河网转成矢量图第五步生成的河网图是栅格的，这里需要将其转化为矢量的，以作为图层查看流域。使用SpatialAnalystToolsfHydrologyfStreamtoFeature田命SchematicsToolsj令ServerToolsT?SpatialAnalystTools田;ConditionalgDensity匡甩Distanceg野ExtractionLt)*%；Generalisation焉电Groundwater?E电JHydrology!Basinj卜Fil

20、lPFlowAccumulation卜-气FlowDirection卜-%FlowLength|卜气Sink$?气SnapPourPointStreamLinkStreamOrderL-气Watershed(T.辛Interpolation每屉Local焉StreamtoFeature电MapAlgebra14Math£EMultivariateit心Neighborhood庄屉Overlay?E霉RasterCreation住乐Reclass到这里为止，就基本符合要求了。下面第9步之后的东西，是为了以后其他计算的需要，例如给流域添加属性，计算流域的面积等用的。9.将流域的栅格转矢量

21、ConversionToolpFromRasterRastertoPolygonArcToolbox?YDAralystTools瞬AnalyticToolriCartographyTcok?ConversionToolsL±?崎FrarrisKMLi耐FrcirnRasterAx:RastertoASCEfRastertoFloatL-xRastertoPointji-ARastertoPolygon*RastertoP&lylins"RasterToVideoFrom'MF5MetadataE±-ToCAD3 ToColada4 ToCover

22、age3-TodBASE$ToGeodatatodseToKMLToRaster5 -ToShapefile6&DataIntenoperataiIrtyTooIs轲DaitaManagementToolsEditingTOols0GuQAdingTq-qIS尊GeostatiSticaIAnalystToolsLinEdrReferEncingTOCe4MuHidmensionT001s现在就可以做专题图来显示不同的流域了Gfe1nmAaA.miiEl-ALrcMAp-Arc|nFa-I'内FiltEdit甘由斡BoAknitfrlUkittrLSwpiP盂血in号CuaAD

23、mAtiWindow*HtlpU已拙3找；c*耳劭0:!*事TA:baC*CcAtcn?'：便砂mSILayers¥>|命*|l:2JOOjDOO呼QUl圉二；tJetworkm*Mfft亘AiiialysFrirtrr*二F口Etreeimret_GuanEtre=imre£E口UjHtershedSuxuanI*口GuymnR曲attuqhdnn12g111nloEjl*EI"*IJ2仙5皆5:q.fiJ丁卡m即带曲5Wh10.通过裁剪获得指定范围的栅格图给定数据，1）栅格图。2）给定范围的面状图层加载数据QlArcWc-AtcMap-FibEditViewBAukmaritiiIrbeHSdudUnGe&pTutcAaingiCusftxjirriiiEWindMdpi巾-71E葩血.:"YliNetworkAniaAEtT帕OflCeOtenos-J"|&国I言T：