基于TauDEM进行水文响应单元划分的方法研究-QGIS集成SWAT插件

QGIS中swat插件，执行命令摘录。

用于制作一个完整的水文响应单元。

pitremove -z srtm_30m.tif -fel srtm_30mfel.tif
d8flowdir -fel srtm_30mfel.tif -sd8 srtm_30msd8.tif -p srtm_30mp.tif
dinfflowdir -fel srtm_30mfel.tif -slp srtm_30mslp.tif -ang srtm_30mang.tif
aread8 -p srtm_30mp.tif -nc -ad8 srtm_30mad8.tif
areadinf -ang srtm_30mang.tif -nc -sca srtm_30msca.tif
gridnet -p srtm_30mp.tif -plen srtm_30mplen.tif -tlen srtm_30mtlen.tif -gord srtm_30mgord.tif
threshold -ssa srtm_30mad8.tif -thresh 1000 -src srtm_30msrcStream.tif
threshold -ssa srtm_30mad8.tif -thresh 100 -src srtm_30msrcChannel.tif
streamnet -fel srtm_30mfel.tif -p srtm_30mp.tif -ad8 srtm_30mad8.tif
streamnet -fel srtm_30mfel.tif -p srtm_30mp.tif -ad8 srtm_30mad8.tif -src srtm_30msrcChannel.tif -ord srtm_30mordChannel.tif -tree srtm_30mtreeChannel.dat -coord srtm_30mcoordChannel.dat -net srtm_30mchannel.shp -w srtm_30mwChannel.tif
moveOutletstostreams -p srtm_30mp.tif -src srtm_30msrcStream.tif -o outlet.shp -om drawoutlets_sel_snap.shp
aread8 -p srtm_30mp.tif -o drawoutlets_sel_snap.shp -nc -ad8 srtm_30mad8.tif
gridnet -p srtm_30mp.tif -o drawoutlets_sel_snap.shp -plen srtm_30mplen.tif -tlen srtm_30mtlen.tif -gord srtm_30mgord.tif
threshold -ssa srtm_30mad8.tif -thresh 1000 -src srtm_30msrcStream.tif
threshold -ssa srtm_30mad8.tif -thresh 100 -src srtm_30msrcChannel.tif
streamnet -fel srtm_30mfel.tif -p srtm_30mp.tif -ad8 srtm_30mad8.tif -src srtm_30msrcStream.tif -o drawoutlets_sel_snap.shp -ord srtm_30mordStream.tif -tree srtm_30mtreeStream.dat -coord srtm_30mcoordStream.dat -net srtm_30mstream.shp -w srtm_30mwStream.tif
streamnet -fel srtm_30mfel.tif -p srtm_30mp.tif -ad8 srtm_30mad8.tif -src srtm_30msrcChannel.tif -o drawoutlets_sel_snap.shp -ord srtm_30mordChannel.tif -tree srtm_30mtreeChannel.dat -coord srtm_30mcoordChannel.dat -net srtm_30mchannel.shp -w srtm_30mwChannel.tif

　　

其中，在动态划分水文响应单元时，可以先整体走一遍，把流向分析、累积量分析什么的，都走完，后续只执行包含shp点的步骤，即可加速水文单元制作。

以下是非公用部分：

moveOutletstostreams -p srtm_30mp.tif -src srtm_30msrcStream.tif -o outlet.shp -om drawoutlets_sel_snap.shp
aread8 -p srtm_30mp.tif -o drawoutlets_sel_snap.shp -nc -ad8 srtm_30mad8.tif
gridnet -p srtm_30mp.tif -o drawoutlets_sel_snap.shp -plen srtm_30mplen.tif -tlen srtm_30mtlen.tif -gord srtm_30mgord.tif
threshold -ssa srtm_30mad8.tif -thresh 1000 -src srtm_30msrcStream.tif
threshold -ssa srtm_30mad8.tif -thresh 100 -src srtm_30msrcChannel.tif
streamnet -fel srtm_30mfel.tif -p srtm_30mp.tif -ad8 srtm_30mad8.tif -src srtm_30msrcStream.tif -o drawoutlets_sel_snap.shp -ord srtm_30mordStream.tif -tree srtm_30mtreeStream.dat -coord srtm_30mcoordStream.dat -net srtm_30mstream.shp -w srtm_30mwStream.tif
streamnet -fel srtm_30mfel.tif -p srtm_30mp.tif -ad8 srtm_30mad8.tif -src srtm_30msrcChannel.tif -o drawoutlets_sel_snap.shp -ord srtm_30mordChannel.tif -tree srtm_30mtreeChannel.dat -coord srtm_30mcoordChannel.dat -net srtm_30mchannel.shp -w srtm_30mwChannel.tif