chapter7.htm

arcinfo workstation资料

<p><font size="2">4 　　　　　　32 　　-74.0 　　　41.0500 </font></p>
<p><font size="2">5 　　　　　　33 　　-73.9875 　41.0500 </font></p>
<p><font size="2">6 　　　　　　34　　 -73.9750　 41.0500</font></p>
<p><font size="2"> ENTER COMMAND><b> Q STOP</b></font></p>
<p><font size="2"> <b>第三步: 定义COVERAGE的坐标系统 </b></font></p>
<p><font size="2">ARC: <b>USAGE</b>　<b> PROJECTDEFINE</b></font></p>
<p><font size="2"> USAGE: PROJECTDEFINE &lt;COVER | GRID | FILE | TIN> &lt;TARGET></font></p>
<p><font size="2"> <cover | GRID FILE TIN> <target> ARC: <b>PROJECTDEFINE </b>　<b>COVER</b>　<b> 
  GEOREF01</b> </font></p>
<p><font size="2">DEFINE PROJECTION</font></p>
<p><font size="2"> PROJECT: <b>PROJECTION</b>　<b> GEOGRAPHIC </b></font></p>
<p><font size="2">PROJECT:<b> UNITS</b>　<b> DD</b></font></p>
<p><font size="2"> PROJECT: <b>PARAMETERS</b></font></p>
<p><font size="2"> ARC:<b>DESCRIBE</b>　<b> GEOREF01</b></font></p>
<p><font size="2"> ..........</font></p>
<p><font size="2"> coordinate system description </font></p>
<p><font size="2">projection　　</font><font size="2">GEOGRAPHIC </font></p>
<p><font size="2">UNITS　　　　　 DD　　　　　　 SPHEROID 　　　CLAK186 </font></p>
<p><font size="2">PARAMETERS: </font></p>
<p>&nbsp;</p>
<p><font size="2"> 当前TIC点的坐标是以参考坐标(经纬度)来表示的, 我们需要把它们投影到地图的坐标系统中, 才能进行下一步的坐标转换. 
  </font></p>
<p><font size="2">我们采用的地图的坐标投影方式是UTM系统(universal Transverse Mecator). 因此, 我们将要把TIC点坐标由经纬度投影转换为UTM的米单位坐标. 
  投影转换可以从键盘输入命令方式进行, 也可以采用读入文件方式进行. </font></p>
<p><font size="2">ARC: <b>USAGE</b>　<b> PROJECT </b></font></p>
<p><font size="2">USAGE: PROJECT <cover | GRID FILE> &lt;COVER | GRID | FILE> 
  &lt; <output> INPUT> &lt;OUTPUT&gt; {PROJECTION_FILE> {NEAREST | BILINEAR | 
  CUBIC} {OUT_CELLSIZE} </font></p>
<p><font size="2">ARC: <b>PROJECT</b>　<b> COVER </b>　<b>GEOREF01</b>　<b> GEOUTM02 
  </b></font></p>
<p><font size="2">**************************************************************** 
  </font></p>
<p><font size="2">* THE INPUT PROJECTION HAS BEEN DEFINED *</font></p>
<p><font size="2"> ****************************************************************</font></p>
<p><font size="2"> USE OUTPUT TO DEFINE THE OUTPUT PROJECTION AND END TO FINISH</font></p>
<p><font size="2"> PROJECT: <b>OUTPUT</b></font></p>
<p><font size="2"> PROJECT:<b> PROJECTION </b>　<b>UTM</b></font></p>
<p><font size="2"> PROJECT:<b> UNITS</b>　<b> METERS</b> </font></p>
<p><font size="2">PROJECT: <b>ZONE 18 </b></font></p>
<p><font size="2">PROJECT: <b>XSHIFT </b>　<b>-580000</b></font></p>
<p><font size="2"> PROJECT:<b> YSHIFT</b>　<b> -4540000</b></font></p>
<p><font size="2"> PROJECT:<b> PARAMETERS</b></font></p>
<p><font size="2"> PROJECT:<b> END </b></font></p>
<p><font size="2">ARC:</font></p>
<p><font size="2"> 现在再列出投影后的TIC点坐标:</font></p>
<p><font size="2"> ARC: <b>LIST</b>　<b> GEOUTM02.TIC</b></font></p>
<p><font size="2"> RECORD　 IDTIC 　　XTIC 　　　YTIC</font></p>
<p><font size="2"> 　1　　　　　 22　　 4024.895 　　5964.872</font></p>
<p><font size="2"> 　2 　　　　　23 　　5074.965　　 5976.987 </font></p>
<p><font size="2">　3 　　　　　24　　 6125.677　　 5989.261</font></p>
<p><font size="2"> 　4 　　　　　32 　　4040.805 　　4577.147 </font></p>
<p><font size="2">　5 　　　　　33 　　5091.074 　　4589.262 </font></p>
<p><font size="2"> 　6　　　　　34　 　6141.984 　　4601.535</font></p>
<p><font size="2"> 同时, 我们可以看到与COVERAGE一起存储的坐标信息也自动地更新, 正确地反映投影后的坐标系统. </font></p>
<p><font size="2">坐标系统信息也可以用命令PROJECTCOPY来从一个数据集拷贝到另外一个数据集. 显然, 当有很多相同坐标系统的数据集, 
  拷贝比一个一个地定义快的多.. 定义坐标信息有利于一些需要坐标比较的分析操作命令如CLIP, INTERSECT的执行. </font></p>
<p><font size="2"><b>第五步 将COVERAGE从数字化坐标转换到真实世界坐标</b> </font></p>
<p><font size="2">TRANSFORM命令可以进行投影转换和仿射转换. 仿射转换是根据三个或三个以上TIC点来计算出输出COVERAGE在比例尺, 
  X偏移, Y偏移, 以及任何可能的旋转变换. </font></p>
<p><font size="2">在进行转换之前, 我们需要先生成一个只包含投影后的TIC点的空COVERAGE. </font></p>
<p><font size="2">ARC: <b>USAGE </b>　<b>CREATE</b> </font></p>
<p><font size="2">　USAGE: CREATE &lt;OUT_COVER&gt; {TIC_BND_COVER}<out_cover> 
  <tic_bnd_cover> </font></p>
<p><font size="2">ARC: <b>CRETE </b>　<b>LANDTR06</b>　<b> GEOUTM02 </b></font></p>
<p><font size="2">这样, 我们将以LANDAT05作为输入COVERAGE, 其单位是数字化仪单位, 输出COVERAGE为LANDTR06, 
  其单位是真实世界坐标(UTM投影). </font></p>
<p><font size="2">ARC:<b> USAGE</b>　<b> TRANSFORM </b></font></p>
<p><font size="2">USAGE: TRANSFORM <in_cover> <out_cover> { &lt;IN_COVER&gt; &lt;OUT_COVER&gt; 
  AFFINE | PROJECTIVE | SIMILARITY}</font></p>
<p><font size="2"> ARC: <b>TRANSFORM </b>　<b>LANDAT05 </b>　<b>LANDTR06 </b></font></p>
<p><font size="2">TRANSFORMING COORDINATES FOR COVERAGE LANDAT05... </font></p>
<p><font size="2">TRANSFORM生成并显示结果报告, 包括输出COVERAGE的TIC, 转换参数, 两个COVERAGE匹配的误差. 
  </font></p>
<p><font size="2">SCALE (X,Y) = (246.145, 245.874) 　SKEW (DEGREES = (-0.037)</font></p>
<p><font size="2"> ROTATION (DEGREES) = (0.347) 　TRANSLATION = (2903.706, 3714.036)</font></p>
<p><font size="2"> RMS ERROR (INPUT, OUTPUT) = (0.006, 1.522) </font></p>
<p><font size="2">AFFINE　　X = AX + BY + C</font></p>
<p><font size="2"> 　　　　　　Y = DX + EY + F </font></p>
<p><font size="2">A = 246.141 　B = -1.648　 C = 2903.706　</font></p>
<p><font size="2"> D = 1.489　 E = 245.868　 F = 3714.036</font></p>
<p><font size="2"> TICID　　 INPUT X 　　　INPUT Y</font></p>
<p><font size="2">　　　　　 OUTPUT X 　 OUTPUT Y 　　X ERROR 　　Y ERROR </font></p>
<p><font size="2">22 　　　4.614 　　　　　　　9.125 </font></p>
<p><font size="2">　　　　4024.895　　　　 5964.872　　　　 -0.531　　　 -0.419</font></p>
<p><font size="2">......................... </font></p>
<p><font size="2">如果RMS误差比较高(如大于0.005英寸), 说明: </font></p>
<ul>
  <li><font size="2">原始的COVERAGE没有准确地数字化</font></li>
  <li><font size="2">地图有较大的变形</font></li>
  <li><font size="2">真实世界坐标没有准确地记录, 例如, 某个点的X,Y坐标记录反了.</font></li>
</ul>
<p><font size="2"> ARC: <b>DESCRIBE</b>　<b> LANDTR06</b></font></p>
<p><font size="2"> 如果你遭遇到大的RMS误差, 需要删除生成的COVERAGE,仔细查出原因. </font></p>
<p><font size="2">ARC:<b> KILL </b>　<b>GEOUTM02 </b>//一旦删除就没有办法恢复 </font></p>
<p><font size="2"><b><font size="3">合并相邻地图</font></b></font></p>
<p><font size="2"> 通常, 每个数据层的地图都是由很多幅相邻的地图组成的. 在开始建立COVERAGE时候, 我们一般都采用边界来代表地图的边界. 
  但是, 完成前面的工作以后, 相邻的COVERAGE应该被合并称为一个完整的COVERGE, 起先的地图边界应该被消去. </font></p>
<p><font size="2">合并地图的第一步工作通常是要检查相邻的COVERGE特征在边界处是否很好地吻合. 这可以通过绘制出同比例尺的地图来与相邻地图比较完成, 
  也可以在ARCPLOT下生成不同的地图来比较相邻的COVERAGE. </font></p>
<p><font size="2">现在我们将合并原有的土地利用图LUAT10和新生成的LANDTR COVERAGE. </font></p>
<p><font size="2"><b>第一步: 比较检查相邻的COVERAGE</b> </font></p>
<p><font size="2">在合并地图之前, 需要确保所有在相邻边界上的特征很好地吻合, 属性项名字和定义应该完全相同, 相连通的特征的属性值也应该相同. 
  </font></p>
<ul>
  <li><font size="2">边界上所有特征是否很好地吻合?? </font></li>
  <li><font size="2">所有的项定义是否相同? </font></li>
  <li><font size="2">边界两边特征属性是否连贯? </font></li>
</ul>
<p><font size="2">ARC:<b> ARCPLOT</b></font></p>
<p><font size="2"> ARCPLOT:<b> MAPEXTENT</b>　<b> LANDTR </b></font></p>
<p><font size="2">ARCPLOT: <b>ARCS</b>　<b> LUAT10</b> </font></p>
<p><font size="2">ARCPLOT: <b>LINESYMBOL</b>　<b> 2 </b></font></p>
<p><font size="2">ARCPLOT: <b>ARCS 　LANDTR</b></font></p>
<p><font size="2"> 可以看到LANDTR和LUAT10部分地重叠. 这是因为我们数字化时候向外延伸了约3毫米. </font></p>
<p><font size="2">ARCPLOT: <b>LINESYMBOL</b>　<b> 3</b></font></p>
<p><font size="2"> ARCPLOT: <b>ARCS </b>　<b>CLIPCOV </b>//CLIPCOV定义了LANDTR的外边界. 
  它应该能够准确地与原有底图重叠, 通常采用拷贝生成或用真实坐标GENERATE生成)</font></p>
<p><font size="2"> ARCPLOT: <b>QUIT </b></font></p>
<p><font size="2">ARC:</font></p>
<p><font size="2"> <b>第二步: 接边COVERAGE准备</b> LANDTR与LUAT10重叠的部分需要删除掉, 我们使用CLIP命令来完成. 
  </font></p>
<p><font size="2">ARC: <b>USAGE</b>　<b> CLIP </b></font></p>
<p><font size="2">USAGE: CLIP &lt;<in_cover>IN_COVER&gt; &lt;CLIP_COVER&gt; &lt;OUT_COVER&gt; 
  <clip_cover> <out_cover> {POLY | LINE | POINT | NET | LINK | RAW} {FUZZY_TOLERANCE} 
  </font></p>
<p><font size="2">ARC: <b>CLIP</b>　<b>LANDTR</b>　<b>CLIPCOV</b>　<b>LANDCP07</b>　<b>POLY</b></font></p>
<p><font size="2">CLIPPING LANDTR WITH CLIPCOV TO CREATE LANDCP07. </font></p>
<p><font size="2">SORTING... </font></p>
<p><font size="2">INTERSECTING... </font></p>
<p><font size="2">ASSEMBLING POLYGONS... </font></p>
<p><font size="2">CREATE NEW LABELS...</font></p>
<p><font size="2"> CREATING DRIVE>:\YOURNAME\LAND\LAND07.PAT... </font></p>
<p><font size="2">ARC: <b>COPY</b>　<b>LANDCP07</b>　<b>LANDEM08</b> //使用拷贝来进行下一步工作EDGEMATCH 
  </font></p>
<p><font size="2">下面使用ARCEDIT来看看CLIP的结果: </font></p>
<p><font size="2">ARC:<b> ARCEDIT </b></font></p>
<p><font size="2">ARCEDIT:<b> EDIT</b>　<b>LANDEM08</b></font></p>
<p><font size="2"> ARCEDIT: <b>DRAWENVIRONMENT</b>　<b>ARC</b>　<b>NODE</b></font></p>
<p><font size="2"> ARCEDIT: <b>DRAW </b></font></p>
<p><font size="2">ARCEDIT: <b>BACKCOVERAGE　LUAT10</b>　<b>4</b> //采用LUAT10作为背景来显示</font></p>
<p><font size="2"> ARCEDIT: <b>BACKENVIRONMENT</b>　<b>ARC</b>　<b>NODE</b></font></p>
<p><font size="2"> ARCEDIT: <b>DRAW </b></font></p>