GIS学习笔记（四）
//先一个一个的来
Arc: &run srtmgrid.aml n39e114.hgt n39e114
Copyright (C) 1982-2006 Environmental Systems Research Institute,
Inc.
All rights reserved.
GRID 9.2 (Sun Sep 17 16:05:34 PDT 2006)
Reading n39e114.hgt as 3-arc-second DEM ...
Copying image file ...
Image parameters:
BYTEORDER M
LAYOUT BIL
NROWS 1201
NCOLS 1201
NBITS 16
ULXMAP 114
ULYMAP 40
XDIM 0.0008333333333
YDIM 0.0008333333333
Converting Image to Grid ...
Running... 100%
Define Projection
Local:
ST
0
Local:
CH
1
Description of Grid N39E114
Cell Size
=
0.000833
Data
Type:
Integer
Number of
Rows
=
1201
Number of Values
=
2422
Number of Columns
=
1201
Attribute Data (bytes)
=
8
BOUNDARY STATISTICS
Xmin
=
113.999583
Minimum Value
=
220.000000
Xmax
=
115.000417
Maximum Value
=
2729.000000
Ymin
=
38.999583
Mean
=
1163.222439
Ymax
=
40.000417
Standard Deviation
=
360.877811
COORDINATE SYSTEM DESCRIPTION
Projection
GEOGRAPHIC
Datum
WGS84
Zunits
METERS
Units
DD
Spheroid
WGS84
Parameters:
Killed XX00000G with the ALL option
Processing completed.
Grid: quit
Leaving GRID...
Arc: lg
Workspace:
D:\GISECO2010\LECTURE2\EXCISE02
Available GRIDs
-------------------
LAODEM15
LAODEM15W
N39E114
Arc: describe n39e119
Geo-dataset n39e119 does not exist or is not readable
Arc: describe n39e114
Description of Grid n39e114
Cell Size
=
0.001
Data
Type:
Integer
Number of
Rows
=
1201
Number of Values
=
2422
Number of Columns
=
1201
Attribute Data (bytes)
=
8
BOUNDARY STATISTICS
Xmin
=
114.000
Minimum Value
=
220.000
Xmax
=
115.000
Maximum Value
=
2729.000
Ymin
=
39.000
Mean
=
1163.222
Ymax
=
40.000
Standard Deviation
=
360.878
COORDINATE SYSTEM DESCRIPTION
Projection
GEOGRAPHIC
Datum
WGS84
Zunits
METERS
Units
DD
Spheroid
WGS84
Parameters:
Arc: &run srtmgrid.aml n39e115.hgt n39e115
Copyright (C) 1982-2006 Environmental Systems Research Institute,
Inc.
All rights reserved.
GRID 9.2 (Sun Sep 17 16:05:34 PDT 2006)
Reading n39e115.hgt as 3-arc-second DEM ...
Copying image file ...
Image parameters:
BYTEORDER M
LAYOUT BIL
NROWS 1201
NCOLS 1201
NBITS 16
ULXMAP 115
ULYMAP 40
XDIM 0.0008333333333
YDIM 0.0008333333333
Converting Image to Grid ...
Running... 100%
Define Projection
Local:
ST
0
Local:
CH
1
Description of Grid N39E115
Cell Size
=
0.000833
Data
Type:
Integer
Number of
Rows
=
1201
Number of Values
=
2789
Number of Columns
=
1201
Attribute Data (bytes)
=
8
BOUNDARY STATISTICS
Xmin
=
114.999583
Minimum Value
=
4.000000
Xmax
=
116.000417
Maximum Value
=
2836.000000
Ymin
=
38.999583
Mean
=
466.910157
Ymax
=
40.000417
Standard Deviation
=
467.708517
COORDINATE SYSTEM DESCRIPTION
Projection
GEOGRAPHIC
Datum
WGS84
Zunits
METERS
Units
DD
Spheroid
WGS84
Parameters:
Killed XX00000G with the ALL option
Processing completed.
Grid: q
Leaving GRID...
Arc: lg
Workspace:
D:\GISECO2010\LECTURE2\EXCISE02
Available GRIDs
-------------------
LAODEM15
LAODEM15W
N39E114
N39E115
Arc: grid
Copyright (C) 1982-2006 Environmental Systems Research Institute,
Inc.
All rights reserved.
GRID 9.2 (Sun Sep 17 16:05:34 PDT 2006)
Grid: bjdem =
mosaic(n39e114,n39115)
Unable to open input grid(s)
Grid: bjdem =
mosaic(n39e114,n39e115)
//把多个图像，grid文件拼在一起
Mosaicing grids ... 100%
Grid: quit
Leaving GRID...
Arc: lg
Workspace:
D:\GISECO2010\LECTURE2\EXCISE02
Available GRIDs
-------------------
BJDEM
LAODEM15
LAODEM15W
N39E114
N39E115
Arc: kill n39e114 all
Killed n39e114 with the ALL option
Arc: kill n39e115 all
Killed n39e115 with the ALL option
//写个有两个run
srtmgrid.aml命令的bat.aml，可以两个一起来了，可以理解，像编程的思想
Arc:
&run bat.aml
Copyright (C) 1982-2006 Environmental Systems Research Institute,
Inc.
All rights reserved.
GRID 9.2 (Sun Sep 17 16:05:34 PDT 2006)
Reading n39e114.hgt as 3-arc-second DEM ...
Copying image file ...
Image parameters:
BYTEORDER M
LAYOUT BIL
NROWS 1201
NCOLS 1201
NBITS 16
ULXMAP 114
ULYMAP 40
XDIM 0.0008333333333
YDIM 0.0008333333333
Converting Image to Grid ...
Running... 100%
Define Projection
Local:
ST
0
Local:
CH
1
Description of Grid N39E114
Cell Size
=
0.000833
Data
Type:
Integer
Number of
Rows
=
1201
Number of Values
=
2422
Number of Columns
=
1201
Attribute Data (bytes)
=
8
BOUNDARY STATISTICS
Xmin
=
113.999583
Minimum Value
=
220.000000
Xmax
=
115.000417
Maximum Value
=
2729.000000
Ymin
=
38.999583
Mean
=
1163.222439
Ymax
=
40.000417
Standard Deviation
=
360.877811
COORDINATE SYSTEM DESCRIPTION
Projection
GEOGRAPHIC
Datum
WGS84
Zunits
METERS
Units
DD
Spheroid
WGS84
Parameters:
Killed XX00000G with the ALL option
Processing completed.
Reading n39e115.hgt as 3-arc-second DEM ...
Copying image file ...
Image parameters:
BYTEORDER M
LAYOUT BIL
NROWS 1201
NCOLS 1201
NBITS 16
ULXMAP 115
ULYMAP 40
XDIM 0.0008333333333
YDIM 0.0008333333333
Converting Image to Grid ...
Running... 100%
Define Projection
Local:
ST
0
Local:
CH
1
Description of Grid N39E115
Cell Size
=
0.000833
Data
Type:
Integer
Number of
Rows
=
1201
Number of Values
=
2789
Number of Columns
=
1201
Attribute Data (bytes)
=
8
BOUNDARY STATISTICS
Xmin
=
114.999583
Minimum Value
=
4.000000
Xmax
=
116.000417
Maximum Value
=
2836.000000
Ymin
=
38.999583
Mean
=
466.910157
Ymax
=
40.000417
Standard Deviation
=
467.708517
COORDINATE SYSTEM DESCRIPTION
Projection
GEOGRAPHIC
Datum
WGS84
Zunits
METERS
Units
DD
Spheroid
WGS84
Parameters:
Killed XX00000G with the ALL option
Processing completed.
Grid: lg
Workspace:
D:\GISECO2010\LECTURE2\EXCISE02
Available GRIDs
-------------------
BJDEM
LAODEM15
LAODEM15W
N39E114
N39E115
Grid: q
Leaving GRID...
Arc: grid
Copyright (C) 1982-2006 Environmental Systems Research Institute,
Inc.
All rights reserved.
GRID 9.2 (Sun Sep 17 16:05:34 PDT 2006)
Grid: basin
Usage: (I) BASIN
(<dir_grid>)
//下面的命令都在subwatershed.aml中，是提取小流域命令集
Grid: fill
Usage: FILL
<in_grid>
<out_grid> {SINK | PEAK} {z_limit}
{out_dir_grid}
Grid: fill laodem15 dem
Filling...
Number of sink(s): 333
Number of sink(s): 52
Number of sink(s): 4
Number of sink(s): 1
Number of sink(s): -1
Grid: dir = flowdirection(dem)
Computing flow direction...
Grid: basn = basin(dir)
Delineating drainage basin...
Grid: accum = FLOWACCUMULATION(dir)
Computing flow accumulation...
Grid: stmgrd = con(accum >
4000,1)
//生成流域
Running... 100%
Grid: stmlnk =
streamlink(stmgrd,dir)
//河流标记
Labeling stream links...
Grid: stmlnk_acumax = zonalmax(stmlnk,accum)
Getting data ..
Computing Statistics ...
Percentage of Cells Processed: 100%
Writing Output ...
Percentage of Cells Processed: 100%
Grid: outlet = con(stmlnk_acumax == accum,stmlnk)//输出出水口的位置
Running... 100%
Grid: subwatershed =
watershed(dir,outlet)
Delineating watershed...
Grid: gridline
Usage: (C) GRIDLINE
(<in_grid>, {POSITIVE | DATA}, {THIN
| NOTHIN},
{NOFILTER | FILTER}, {ROUND | SHARP}, {item},
{thickness}, {dangle}, {weed})
Grid: stmcov =
gridline(stmlnk) //grid-->line
参数（dem数据）
Converting GRID STMLNK to Coverage
STMCOV //每一个河段对应一个小流域
Thinning...
Vectorizing...
Generalized coverage
D:\GISECO2010\LECTURE2\EXCISE02\XX03664000
From: 9058 Vertices and 159 Arcs
To : 2390 Vertices and 159 Arcs
Removing dangling arcs...
151 unique nodes built for
D:\GISECO2010\LECTURE2\EXCISE02\XX03664000
Eliminating pseudo nodes...
Number of
Arcs
(Input,Output) = 149 141
Grid: lg
Workspace:
D:\GISECO2010\LECTURE2\EXCISE02
Available GRIDs
-------------------
ACCUM
BASN
BJDEM
DEM
DIR
LAODEM15
LAODEM15W
N39E114
N39E115
OUTLET
STMGRD
STMLNK
STMLNK_ACUMAX
SUBWATERSHED
Grid: streamorder
Usage: (I) STREAMORDER
(<net_grid>,
<dir_grid>, {STRAHLER |
SHREVE})
Grid: stmord =
streamorder(stmgrd,dir)
//使每一个河段都有一个编码
Computing stream order...
Grid: streamline
Usage: (C) STREAMLINE
(<net_grid>,
<dir_grid>, {out_item}, {weed})
Grid: stmcovln =
streamline(stmord,dir)
//让编码后的图层转化为矢量的
Converting GRID STMORD to Coverage STMCOVLN
Grid: q
Leaving GRID...
Arc: q
Leaving ARC ...
D:\GISECO2010\lecture2\excise02>
