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

ArcGIS10.0在DEM勺基础上划分小流域根据deMI取河流网络，计算流水累积量、流向、根据指定的流域面积大小自动划分流域1.加载DEM打开ArcMap,加载数据DE嗷据文件如果DE汕有很多异常值，可能刚加载进来的时候是全灰或者全白，这个时候需要右键t图层属性,进行调整。2.FillDEMg面存在着一些凹陷的区域（DE冲身是插值计算的，很难与现实情况完全符合）由于这些区域异常低值存在，使得该区域在进行水流流向计算时得到不合理的水流方向。因此，应该首先对原始DE嗷据进行洼地填充，得到无洼地的DEIM打开ArcToolbox中的水文分析模块。使用SpatialAnalystToolsTHydrologytFill如果没有水文分析模块，请用CustomizerExtensionsTSpatialAnalyst添加扩展模块。?&CjonversionTools?&dDataInteroperabilityTookDataManagementToolsEditingToolsGeocodingToolsGeortatistical:2尺aAnalystToolsLinearReferencingToolsMultidimensionToolsNetworkAnalystToolsParcelFabric?筋?TookSchematicsToolsServerTools?SpatialAnalystTools霉Conditional甩Density&Distance屉Extraction3Generalization2?d7Groundwater裁Hydrology卜-%Basin?：0卜FlowAccumulation理卜FlowDirectionh衣*FlowLength卜气Sink卜? ? 5-.? ? E.Stream Order-%SnapPourPoint卜.$StreamLink卜旅、StreamtoFeature?EE-SBFillinputsurfaceraster|DEM.mgOutputsurfacerasterF:HeseadikTEMPzhangyaafa-igningxiaDEMFIZlimitOptionalCancat-rorm&PShawHeA3.FlowDirectionArcGIS中的水流方向利用D8算法计算。针对每一个栅格，将其高程与周围八个栅格进行比较，得到水流方向。最终得到的结果是一个栅格图，每个象元上的值的含义如下图所示。如果某象元算得的水流方向是正东方向，则值为1?如果西北，则为128,以此类推。注意，这里要用的DE汜经不是原始的二步生成的结果DEM（文件命名不能超过使用SpatialAnalystToolsfHydrologyDEM了，需要用经过Fill后的DEM也就是第13个字符？好像是）fFlowDirection啊ffiPIGeocodingToolsj?GeostatisticalAnalystTools令LinearReferencingTools田A1MultidimensionTools?NetworkAnalystToolsS)商IParcelFabricTools习野SchematicsTools由令ServerTools令SpatialAnalystToolsl氐Conditional庄aaDensity由邑？DistancelExtraction?霉Generalization色甩Groundwater曰Hydrologyl气BasinAllFlowAccumulationFlowDirection.?_%FlowLength一冬SinkSnapPourPoint?-气StreamLink?-气StreamOrderStreamtoFeature*.)*?飞WatershedInterpolation)-.LocalS*:?M叩Algebra4.FlowAccumulationFlowAccumulation由FlowDirection数据计算而来。每一i个栅格FlowAccumulation的值代表着其上游有多少个栅格的FlowDirection最终汇流经过该栅格。一般而言，计算出来的Accumulation的数值数值越大，代表越有可能是河谷。使用SpatialAnalystToolsfHydrologyfFlowAccumulationg-QlGeostatisticalAnalystToolsj令LinearReferencingToolsj令MultidimensionToolsTNetworkAnalystTools囱令ParcelFabricTools2)ISchematicsTools士野ServerToolsj野SpatialAnalystTools砂Conditional?(2-电DensityDistancei自ExtractionlGeneralizationiI田电GroundwateryHydrology?iABasin-气FillFlowAccumulation?一氏？FlowDirection?飞FlowLength一气Sink?-%SnapPourPoint%StreamLink?;?、，StreamOrder、StreamtoFeatureWatershed由Interpolation也电Localipr屉MapAlgebrag)乘MathFlowAccumulationoErvrMFFtEShowHgip5.河网计算也就是说，如果在某个栅格的上游有 定义为河流。由于本文用的栅格是当FlowAccumulation达到一定值的时候，就会产生地表水流,格那么所有那些Accumulaten大于阈值的栅格就是潜在的水流路径，由这些水流与定义的汇路径构成的网络，就是河网。这个阈值是需要人为设定的，10000个栅格的水流会流经这个栅格，则将这个栅30米X30米，也就是说，如果某个栅格点，其上游的汇水面积超过9平方公里，则认为这里是河流使用SpatialAnalystToolsfMapAlgebrafRasterCalculatorT NetworkAnalyst ToolsGeostatisticalAnalystToolsj文寸LinearReferencingTools令MultidimensionToolsAnalystTools囱令ParcelFabricTools2)商ISchematicsTools0ServerTools自野Spatial阿Conditional?电DensityDistancei自Extraction&?电GeneralizationiI田？，Groundwater商HydrologyE%?Interpolation由屉Local?心“叩Algebra(3?甩Math屉MultivariateIT辛Neighborhood审心OverlayIF?电RasterCreation国？甩Reclass庄电？SolarRadiation审？心Surface曲区罔ZonalM令SpatialStatisticsTools0令TrackingAnalystToolsRasterCalculatorMapAlgebraexpressionLayersandvariables?HowAccum?flowdirDEMFillConditionalConPickConCFIowAccunT10000fl)SetNulMathinnrnMnAbsExpOutputrasterF:ResearchTEMPzhangyaofangningxjaxsStreamnetCancel|Envronments.ShowHelp?QUri山d*ArMi?-占Urtttled-ArcMap-ArcJnfoFileEditViewBookmArksInsertSelectionGeoprocessingCustomiieWindowsHelp+1:250.mrtMjfiEditorTTableOf匚onttntslayersSSF:RetearchTEr/lPzh;jngyaoaSitreamnet1Floi/vAccumValueHigh:4.18仔7E+007Low0口flowdir1124QJ2_430DEMFiHValueHrah:3537Htrr1B3685.8543836041.J03Mete6.StreamlinkStreamlink记录着河网中节点之间的连接信息。Streamlink的每条弧段要么连接着两个作为出水点或汇合点的结点，要么连接着作为出水点的结点和河网起始点。河网节点是计算流域的输入参数。使用SpatialAnalystToolsfHydrologyfStreamLinkt)-QlSchematicsToolsj)?ServerToolsT射ISpatialAnalystTools(5*;ConditionalE甩Density国屉Di&tance中苍ExtractionGeneralization阿GroundwaterQ甩Hydrology!?一叫aBasinj%FillFlowAccumulation卜-冬FlowDirection卜%FlowLength|ASink$?%SnapPourPointStreamLink;StreamOrder!jr严、StreamtoFeature!?-气Watershed(TInterpolation由$Local阿？MapAlgebra审电Math国Multivariate由心Neighborhood由屉Overlay手屉RasterCreation色乐ReclassStreamLinkinputstreamTBaier|strearrmeiInputflowdirectionrasterflowdirOutputrasterF:ResearGhTEMPjhangydofangYiinAciaftreamlinkCancelEnvrDnrnentsShowHet7.Watershed先确定一个出水点，然后结合水流方向数据，分析搜索出该出水点上游所有流过该出水口的栅格,化为这个出水点之上的流域。这里，我们将第6步算出来的河流节点作为这一步的出水点数据。当然，这样的话计算出的流域将会很破碎，而且得不到流域之间的隶属关系如果需要得到更为精确的流域，那么就需要人工去判读出水点。使用SpatialAnalystToolsfHydrologyfWatershed田翻ISchematicsToolsj?ServerToolsT?SpatialAnalystTools田ConditionalgDensity匡屉DistancegmExtraction由？电Generalisation焉电Groundwater?E氐HydrologyIBasinj|?%Fill？气FlowAccumulation卜-气FlowDirectionFlowLengthSinkSnapPourPointStreamW-iWatershe:LCLinkStreamOrderStreamtoFeature/L气L，*SInterpolation?心Local!?MapAlgebra每甩Math庄岭？MultivariateffiNeighborhood时屉Overlayi4霉RasterCreationE令ReclassI CD IIS I%WatershedInputflowdirectionrasterpowdirInputrasterorfeaturepourpointdata|streamlinkPourpointfield(optional)VALUEOutputrasterF:flesearchTEMPvihangydofangningxjaA3tershedOKCancai|Envronments.ShowHelpFlc*fi?AccuiT14?eaS42Edrtcn口 - ArcMip Ardmf QSeltcticnGecpAcefEingiCutTomizeWindowH*lpFileEditVlcvvBsckrnihcEImeit：ZtdbleOfCiantenti:J县|I.耳Lav*rtB因FARjSB0rrhTEMPzhengQr?gdLow:1回：VAIU?High;?27tLciw:口1TtrcotrlinkValutKgr:3271rkAnaHetT1Higr;A16117A017917?4.BA+01A|3.59&为了显示的需要，我们还需要将河网加进来到这一步，就把所有的集水区划分出来了。以显示每个小流域隶属哪一段河道。8.将河网转成矢量图第五步生成的河网图是栅格的，这里需要将其转化为矢量的，以作为图层查看流域。使用SpatialAnalystToolsfHydrologyfStreamtoFeature田命SchematicsToolsj令ServerToolsT?SpatialAnalystTools田;ConditionalgDensity匡甩Distanceg野ExtractionLt)*%；Generalisation焉电Groundwater?E电JHydrology!Basinj卜FillPFlowAccumulation卜-气FlowDirection卜-%FlowLength|卜气Sink$?气SnapPourPointStreamLinkStreamOrderL-气Watershed(T.辛Interpolation每屉Local焉StreamtoFeature电MapAlgebra14MathEMultivariateit心Neighborhood庄屉Overlay?E霉RasterCreation住乐Reclass到这里为止，就基本符合要求了。下面第9步之后的东西，是为了以后其他计算的需要，例如给流域添加属性，计算流域的面积等用的。9.将流域的栅格转矢量ConversionToolpFromRasterRastertoPolygonArcToolbox?YDAralystTools瞬AnalyticToolriCartographyTcok?ConversionToolsL?崎FrarrisKMLi耐FrcirnRasterAx:RastertoASCEfRastertoFloatL-xRastertoPointji-ARastertoPolygon*RastertoP&lylinsRasterToVideoFromMF5MetadataE-ToCAD3 ToColada4 ToCoverage3-TodBASE$ToGeodatatodseToKMLToRaster5 -ToShapefile6&DataIntenoperataiIrtyTooIs轲DaitaManagementToolsEditingTOols0GuQAdingTq-qIS尊GeostatiSticaIAnalystToolsLinEdrReferEncingTOCe4MuHidmensionT001s现在就可以做专题图来显示不同的流域了Gfe1nmAaA.miiEl-ALrcMAp-Arc|nFa-I内FiltEdit甘由斡BoAknitfrlUkittrLSwpiP盂血in号CuaADmAtiWindow*HtlpU已拙3找；c*耳劭0:!*事TA:baC*CcAtcn?：便砂mSILayers|命*|l:2JOOjDOO呼QUl圉二；tJetworkm*Mfft亘AiiialysFrirtrr*二F口Etreeimret_GuanEtre=imreE口UjHtershedSuxuanI*口GuymnR曲attuqhdnn12g111nloEjl*EI*IJ2仙5皆5:q.fiJ丁卡m即带曲5Wh10.通过裁剪获得指定范围的栅格图给定数据，1）栅格图。2）给定范围的面状图层加载数据QlArcWc-AtcMap-FibEditViewBAukmaritiiIrbeHSdudUnGe&pTutcAaingiCusftxjirriiiEWindMdpi巾-71E葩血.:YliNetworkAniaAEtT帕OflCeOtenos-J|&国I言T：f*门册I?益腔汾鲁厂Edrt:-rwL自|Iv?iCjRaarcnME回GuytldirlnjyectengXuIF;RsearcnMEOEhtiryHighs3537LewW533II本例中右一副大范围的DEM,有一个小的面状区域。现在需要获得该小面状范围的DEM打开ArcToolboxArcToolboxi-DataLnteropcrabilrtyTools白，？Dmt孑ManagementToo|&庄DataCompari;or|4$Dvtdb诳iDistributedGeodd（abaseLr?鬻DomainEi匡磬Feature口垂B虫鸟Features|由.Fields自HGeodatabase1+?抵，GeneralTJ?甩GeneralirationGraphIndies应.心Joins-hr?UI*LayersandTableAiewe（t?亲隹）?Package-aProjectionsandTransformationsji-E-Raster一itiMosaicCatdset14ARasterCatalogIRasterDatasetiaj每-Raitei1Proc?singtlntl-j、ComposrteBands:a.CreateOrthoCorrectedRatterOat;*匚r&atePan-sharperedRsster门妣丙:气ExtractSubddhseth,.RatterToDTED使用DataManagementtoolRasterARasterProcessinClip设置相关参数注意，如果使用的面不是规则矩形定要点上“UseInputFeaturesforClippingGeometry(optional)查看结果