空间数据的转换与处理ppt课件

1、第四章第四章 空间数据的转换与处理空间数据的转换与处理l投影变换l数据格式转换l数据处理lArcToolbox1l4.1定义投影 在GIS中，空间数据是一个重要的部分。整个GIS都是围绕空间数据的采集、加工、存储、分析和表现展开的。而原始空间数据本身通常在数据结构、数据组织、数据表达上和用户自己的信息系统不一致，就需要对原始数据进行转换与处理，如投影变换，不同数据格式之间的相互转换，以及数据的裁切、拼接等处理。以上所述的各种数据转换与处理均可以利用ArcToolbox中的工具实现。 2l投影变换（projection transformation）是将一种地图投影点的坐标变换为另一种地图投影点

2、的坐标的过程,是研究投影点坐标变换的理论和方法。 3 GIS中的坐标系定义是GIS系统的基础，正确定义GIS系统的坐标系非常重要。GIS中的坐标系定义由基准面和地图投影两组参数确定，而基准面的定义则由特定椭球体及其对应的转换参数确定，因此欲正确定义GIS系统坐标系，首先必须弄清地球椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念及它们之间的关系。 45 大地水准面是由静止海水面并向大陆延伸所形成的不规则的封闭曲面。它是重力等位面，即物体沿该面运动时，重力不做功（如水在这个面上是不会流动的） 67 地球椭球体地球椭球体只不过是一个具有长半轴，

3、短半轴和变率的椭球体，可以任意放置的，它没有为我们规定度量的起点，所以就有基准面的产生，而基准面就是规定了度量标准。 8 基准面是利用特定椭球体对特定地区地球表面的逼近，因此每个国家或地区均有各自的基准面，我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面.椭球体与基准面之间的关系是一对多的关系，也就是基准面是在椭球体基础上建立的，但椭球体不能代表基准面，同样的椭球体能定义不同的基准面，一般意义上基准面与参考椭球体是同一个概念。9 对于地理坐标，只需要确定两个参数，即椭球体和大地基准面。10l地球上同一位置的坐标在不同的基准面上是不一样的，而基准面是构成坐标系的一个部

4、分，因为基准面在定位的时候牵扯到了相对地心的平移或旋转等，所以对于地理坐标转换无法直接进行，需要一个转换参数，而这些参数也是基于不同的模型的，常用的有3参数和7参数，3参数是比较简单的也是比较容易理解的，3参数是在两个基准面之间进行了X，Y，Z轴的平移113参数之间两个基准面的关系12如果知道了这三个平移的参数外加个基准面上的点，那么另外一点的坐标就是13 7参数的模型比较复杂，这种复杂的同时让精度大为提高，7参数不仅仅考虑了两个基准面之间的平移，还考虑了旋转外加一个比例因子(椭球体的大小可能不一样)14知道平移三参数, 旋转三参数以及比例因子S，外加一个基准上的坐标就可按照下 面的公式求出另

5、外一个基准面上的坐标：15l4.1.1 定义投影 定义投影（Define Projection），指按照地图信息源原有的投影方式，为数据添加投影信息。在ArcGIS中利用【数据管理工具】工具箱, 【投影变换】工具集中的【定义投影】命令，能够为数据定义投影。16定义投影有三种方法定义投影有三种方法l选择固定的投影l导入投影（当已知原始数据与某一数据的投影相同时）l自定义投影17 1 . 栅格（Raster）数据的投影变换：利用【数据管理工具】工具箱,【投影和变换】中【栅格】工具集,【投影栅格】命令对栅格数据进行投影变换。输入的栅格必须已具有投影信息，若没有则在【输入坐标系】中指定数据的原始投影信

6、息【地理坐标变换】：当输入和输出坐标系的基准面相同时，地理（坐标）变换为可选参数。如果输入和输出基准面不同，则必须指定地理（坐标）变换。4.1.2投影变换182 . 要素类（Feature）数据的投影变换 ：利用【数据管理工具】工具箱, ,【投影和变换】中的【要素】工具集,【投影】命令，对矢量数据进行投影变换。【批量投影】可以进行多个要素类的投影变换。19 20一次处理多一次处理多个投影（不个投影（不经常用到）经常用到）21选择Geographic Transformation时，可以选择一种。因为没有识别出坐标系，所以提示错误。22选择投影坐选择投影坐标系标系2324使用已有使用已有的投影坐

7、的投影坐标系投影标系投影25以以t25_Project投影图像投影图像的坐标系作的坐标系作为投影坐标为投影坐标系系26新建一新建一个投影个投影坐标系坐标系2728l4.1.3 数据变换 数据变换是指对数据进行诸如放大、缩小、翻转、移动、扭曲等几何位置、形状和方位的改变等的操作。对矢量数据的相应操作在ArcMap中【编辑器】工具条的若干工具实现（详见第三章）。而栅格数据的相应操作则集中于ArcMap中的地理配准工具条和ArcToolbox中的投影和变换工具集来实现，以下分别就栅格数据的翻转（Flip）、镜像（Mirror）、重设比例尺（Rescale）、旋转（Rotate）、移动（Shift）和

8、扭曲（Warp）等分别介绍。 291.翻转（Flip）：是指将栅格数据沿着通过数据中心点的水平轴线，将数据进行上下翻转。利用Flip命令，对数据进行翻转。FlipFlip对对话框话框翻转翻转（翻转（FlipFlip）的图解表达）的图解表达302.镜像（Mirror）：是指将栅格数据沿着通过数据中心点的垂直轴线，将数据进行左右翻转。利用Mirror命令，对数据进行镜像。镜像（镜像（Mirror）的图解表达）的图解表达镜像Mirror对对话框话框313.重设比例尺（Rescale）：是指将栅格数据按照指定比例分别沿X轴和Y轴放大或缩小。Rescale命令，对数据重设比例尺。重设比例尺重设比例尺（重

9、设比例尺（Rescale）的图解表达）的图解表达Rescale对对话框话框324.旋转（Rotate）：是指将栅格数据沿着指定的中心点旋转指定的角度。利用Rotate命令，对数据进行旋转。 图图4.13 Rotate对对话框话框 旋转旋转（旋转（Rotate）的图解表达）的图解表达335.移动（Shift）：是指将栅格数据分别沿X轴和Y轴移动指定的距离。利用Shift命令，对数据进行移动。Shift对对话框话框移动移动（移动（Shift）的图解表达）的图解表达346.扭曲（Warp）：是指将栅格数据通过输入的控制点进行多项式变换。利用Warp命令，对数据进行扭曲变换。Warp对对话框话框 原数

10、据一次多项式二次多项式三次多项式扭曲（扭曲（Warp）的图解表达）的图解表达354.1.3数据变换数据变换l一、矢量数据的空间校正【空间校正】工具用于矢量数据的空间位置匹配，它提供用于对齐和整合数据的交互式方法，看执行的一些任务包括：将数据从一个坐标系转换到另一坐标系中、纠正几何变形、将沿着某一图层的边的要素与邻接图层的要素对齐及在图层之间的复制属性。由于空间校正在编辑会话中执行，因此可使用现有的编辑功能来增强校正效果。36 当要素的坐标超出要素类的属性域范围时会发生此错误。这种情况会在创建新要素或编辑现有要素的坐标时发生。 坐标必须位于要素类 x,y 属性域以及 z 或 m 属性域（如果要素

11、类可存储这些值）的范围内。属性域在创建要素类时设置并且无法在以后进行编辑。如果需要在该 x,y 位置或者使用这些 z 或 m 值创建新要素，可将要素类导出为新要素类并增加属性域值。37二、地理配准二、地理配准l栅格数据通常是通过扫描纸质地图或采集航空及卫星照片获得 。 l通过扫描获取的影像不包含定义其地理空间位置所需的信息。而航空及卫星照片所使用的坐标系统是相对于通用GIS平台软件所使用的坐标系统是独立的。l为了能够将这些影像数据与其它的数据集成，以便进行分析， 就必须对其进行处理：用户需要事先将这些数据校准（配准）到一个指定的地图坐标系38扫描地图的误差来源扫描地图的误差来源 GIS中数据的

12、来源主要是对地图图纸的数字化,扫描数字化过程中引起的误差主要决定于l要素对象l软件处理技术l扫描仪39l扫描要素对象扫描要素对象要素本身的宽度、复杂程度、粘连以及图面的整洁和清晰程度都对扫描数字化误差有一定的影响。例如,线条的粘连,结合处易出现较大的误差;线条发虚,会得到多个实体;图面不整洁引入了噪声,易引起软件误判;线条不光滑,易出现毛刺等。40l软件处理技术软件处理技术在扫描数字化过程中,三个主要参数:分辨率、门槛值(灰度值或对比度值)和滤波值的确定将对扫描图的质量产生重大的影响;而图像处理、几何校正和矢量化等后处理技术,其功能的强弱、模型的优化将直接影响扫描数字化的精度。如采用合适的校正

13、模型(仿射变换、双线形变换、多项式变换等)、定向点的自动对中、采样点的自动对中等,将有效的提高扫描数字化的精度。41l扫描仪扫描仪在扫描过程中,由于使用CCD扫描仪,会引入一些误差。主要包含有:扫描仪的分辨率;光学误差;电信号传输过程中造成的辐射误差;沿导轨扫描过程中,由于机械运动、速度不均或其它原因所造成的直线性误差;线阵方向与扫描方向不垂直所引起的CCD线阵的直线性误差;外界因素影响产生的误差;随机性误差等等。42所有以上这些扫描的误差引起的几何变形,可看成平移、旋转、缩放、仿射、弯曲以及各种更高变形的综合作用结果。在实际操作过程中,很难对这些误差一一进行变形改正,只能综合考虑它们的影响,

14、综合校正。 43l输入到计算机中的图形，实际上都是通过其位置坐标(x,y)来表示，因此校正过程实质上是找一种数学关系(或函数关系)，描述变换前图形坐标(x,y)与变换后图形坐标(x,y)之间的换算，其数学关系一般描述为 x=f1(x,y) y=f2(x,y)44多项式变换多项式变换(Polynomial)l这个数学关系常表示为二元多项式一次、二次或三次及更高次表达式. 其中A、B代表二次以上高次项之和。上式是高次变换方程，符合上式的变换称为高次变换。在进行高次变换时，需要有6对以上控制点的坐标和理论值，才能求出待定系数。 45l当不考虑高次变换方程中的A和B时，则变成二次变换方程，称为二次变换

15、。二次变换适用于原图有非线性变形的情况，至少需要5对控制点的坐标及其理论值，才能求出待定系数。 46仿射变换（仿射变换（1次多项式）次多项式）l仿射变换是使用最多的一种几何纠正方式，只考虑到x和y方向上的变形，仿射变换的特性是：直线变换后仍为直线；平行线变换后仍为平行线；不同方向上的长度比发生变化。47l对于仿射变换，只需知道不在同一直线上的三对控制点的坐标及其理论值，就可求得待定系数。但在实际使用时，往往利用4个以上的点进行纠正，利用最小二乘法处理，以提高变换的精度。 48l地理配准：是为了使得影像数据可以和GIS矢量数据集成在一起，而为影像数据指定一个参考坐标系的过程。49l将扫描地图配准

16、到坐标系下50影像配准的步骤影像配准的步骤（RegisterRectify）校准栅格数据 （选择控制点）坐标变换 （求解二元多项式n次方程）检查均方差（计算控制点误差）重采样矫正（Rectify）：生成新的影像文件 （三种重采样算法）5152校准栅格数据校准栅格数据l通常,你会将栅格数据校准到已经存在具有坐标信息的空间数据 （矢量数据） 。首先假定矢量化数据中的一些空间要素 (目标数据)也同时存在于要进行配准的栅格图像上比如： 街道、建筑物、河流. l地理配准的基本过程是在栅格图像中选取一定数据的控制点，将它们的坐标指定为矢量数据中对应点的坐标（在空间数据中，这些点的坐标是已知的，坐标系统为地

17、图坐标系）53控制点控制点l在配准中我们需要知道一些特殊点的坐标，即控制点。l控制点可以是经纬 线网格的交点、公里网格的交点或者一些典型地物的坐标。l我们可以从 图中均匀的取几个点。如果我们知道这些点在我们矢量坐标系内坐标， 则直接输入控制点的坐标值，如果不知道它们的坐标，则可以采用间接方法获取从矢量数据中选取。 54选取控制点选取控制点l控制点的数目取决于 你打算使用哪一种数学方法来实现坐标转换. 但是，过多的控制点并不一定能够保证高精度的配准 。要尽可能使控制点均匀分布于整个格格图像，而不是只在图像的某个较小区域 选择控制点。l通常,先在图像的四个角选择4个控制点，然后在中间的位置有规律地

18、选择一些控制点能得到较好的效果55l地形图中，读取控制点的坐标，图中红色控制点的坐标为（564000,2776000），单位：米56l输入控制点坐标57坐标变换l一旦你选取了足够的控制点, 你就可以将栅格数据变换（或转换）到地图坐标系统下. 转换(Transformation) 运用一种数学变换方法来重新确定栅格数据中每个像元的灰度值。l一次多项式：仿射（affine）变换可以将栅格数据平移, 缩放, 及旋转. 栅格图像中的一条直线变换后然后为直线。 矩形和正方形变换后为平行四边形5859检查均方差(RMS)l坐标转换的准确程度可以通过比较某一点在地图中的实际坐标与根据变换公式得到的坐标来判断

19、。l这两个点之间的距离之差称为残差(residual error). 通过计算均方差(RMS)获取控制点总误差。l均方差(RMS) 的大小描述了变换公式在不同控制点间的一致性。l可以将残差特别大控制点删除，然后添加新的控制点。60lRMS 比较小时，说明控制点的选取是比较准确的，但也要注意有可能存在残差非常大的控制点，但由于其它控制点是很精确的，所以总的RMS误差比较小的情况。61矫正栅格数据-重采样l你可能会认为一旦实现了栅格数据到地图坐标变换之后，每一个像元都被转换到了新的地图坐标。事实并非如此。在地理配准过程中，将基于地图坐标生成一个“空的” 矩阵, 矩阵中每个元素的值（表示颜色）将通过

20、重采样重新计算。62有三种通用的图像重采样技术 l最小紧邻l双线性内插l立方卷积Nearest neighbor assignment takes the value from the cell closest to the transformed cell as the new value. Its the fastest resampling technique and is appropriate for categorical, or thematic, data. Bilinear interpolation and cubic convolution techniques comb

21、ine a greater number of nearby cells (4 and 16, respectively) to compute the value for the transformed cell. These two techniques use a weighted averaging method to compute the output transformed cell value and thus are only appropriate for continuous data such as elevation, slope, and other continu

22、ous surfaces.63Rectify：矫正栅格数据lRectify 生成一个新的已经过地理配准的的格文件。可以保存为 ESRI GRID、 TIFF或者 ERDAS IMAGINE的格式.64要素矢量化要素编辑的一般步骤创建新图层（要素类）：点、线、多边形编辑环境及工具点要素编辑线要素编辑多边形要素编辑65要素编辑的一般步骤要素编辑的一般步骤1.在ArcCatlog中创建新的要素类2.在ArcMap中加载新建的要素类，如有必要，加载经过配准的扫描地图3.在ArcMap中打开“编辑器”工具栏，执行其中的 “开始编辑”命令，进入编辑状态。4.利用“编辑器”的功能完成地图要素的分层提取，打开

23、图层的属性表，输入要素的相关属性5.在编辑过程中，点击“编辑器”中的“保存编辑”可以随时保存修改的结果，点“停止编辑”完成编辑。66创建新图层（要素类）创建新图层（要素类）lArcCatalog中可以创建矢量数据集(shape文件和地理数据库geodatabase = “要素类”) 点：Point线：Line多边形：Polygonl在ArcMap中提供了全面的工具，可以完成矢量数据的编辑空间数据属性数据67l新建图层要素类 指定要素类型：点？线？多边形？68l在创建图层时指定坐标系统 在ArcMap中指定数据框的坐标系统69l在创建图层（要素类）的时候除指定坐标系外，还需要指定要素类的坐标范围

24、l可以从已有的数据导入或手动指定这个范围70l添加新建的图层，添加已配准的影像图层，在编辑器中，点击“开始编辑”，并在目标图层中选中要编辑的图层。71屏幕跟踪数字化屏幕跟踪数字化l从“编辑器”工具栏中选中草图工具，根据扫描地图上的内容分层提取地图要素。72编辑图层中要素的属性编辑图层中要素的属性l属性字段可以在创建图层是定义，也可以在ArcMap中添加Fields are added; definel字段名称：field namel数据类型：data typel字段宽度：widthl精度：decimal precision73l打开图层（要素类）的属性表，添加新字段添加属性字段添加属性字段74

25、编辑要素属性信息编辑要素属性信息l在属性字段添加后，就可以输入或修改要素的属性75编辑环境及工具编辑环境及工具l主要编辑环境“编辑器”工具栏76“编辑器编辑器”工具栏工具栏77捕捉环境捕捉环境Snappingl捕捉或悬挂操作:how features align during creation/editing线段之间的连接 (在两条线段连接的地方加一个结点)completion of polygonsavoid overshoots/undershootsavoid slivers or gaps78捕捉线段捕捉线段79点要素的编辑点要素的编辑l在图层中新建、删除、移动一个点要素80线要素的编

26、辑线要素的编辑l悬挂功能l修改线段的形状（编辑顶点）l分割线要素、公共顶点的编辑l合并线要素l封闭线段81打断线要素打断线要素打断工具82封闭线段封闭线段83多边形要素的编辑多边形要素的编辑l多边形顶点编辑（改变形状）l分割多边形l公共边的编辑l要素合并（Union）：求并集，比如：由多个岛屿组成的行政区，飞地。l要素联合（Combine): 中心有岛屿的湖泊l要素相减（Subtract）：去除重叠部分（landuse)l要素相交（Intersection): 获取多个多边形的公共部分844.2 数据格式转换数据格式转换 空间数据的来源有很多，如地图、工程图、规划图、照片、航空与遥感影像等，因

27、此空间数据也有多种格式。根据应用需要，对数据的格式要进行转换。转换是数据结构之间的转换，而数据结构之间的转化又包括同一数据结构不同组织形式间的转换和不同数据结构间的转换。其中，不同数据结构间的转换主要包括矢量到栅格数据的转换和栅格到矢量数据的转换。数据格式转换工具数据格式转换工具85数据格式转换数据格式转换 ArcGIS ArcGIS中的空间数据主要有两种类型：一是基中的空间数据主要有两种类型：一是基于文件的空间数据于文件的空间数据, , 如如coverageShapefileGridImageTIN coverageShapefileGridImageTIN 。GeodabaseGeodab

28、ase数据模型也可以在数据库中管理同数据模型也可以在数据库中管理同样类型的空间数据；样类型的空间数据；二是基于数据库的空间数据二是基于数据库的空间数据86基于文件的空间数据基于数据库的空间数据CoveragesOracleShapfilesOracle with SpatialGridsDB2 with its Spatial TypeTINsInformix with its Spatial TypeImages(各种格式)SQL ServerVector Product Format(VPF)filesPersonal Geodatabase(微软的Access)CAD文件表（各种格式）A

29、rcGIS中的数据类型中的数据类型87l4.2.1 数据结构转换 地理信息系统的空间数据结构主要有栅格结构和矢量结构，这两种数据结构是模拟地理信息的两种不同的方法。在地理信息系统中栅格数据与矢量数据各具特点与适用性，为了在一个系统中可以兼容这两种数据，以便有利于进一步的分析处理，常常需要实现两种结构的转换。 881、栅格数据向矢量数据的转换、栅格数据向矢量数据的转换 目的：是为了将栅格数据分析的结果，通过目的：是为了将栅格数据分析的结果，通过矢量绘图装置输出，或者为了数据压缩的需要，矢量绘图装置输出，或者为了数据压缩的需要，将大量的面积栅格数据转换成由少量数据表示将大量的面积栅格数据转换成由少

30、量数据表示的多边形边界。的多边形边界。 根据描述对象的不同，栅格数据可以转换根据描述对象的不同，栅格数据可以转换为点状、线状和面状的矢量数据。为点状、线状和面状的矢量数据。89Simplify polygons 是默认的，可以简化面状数据的边界形是默认的，可以简化面状数据的边界形状状栅格图像90参数为VALUE参数为COUNT912、矢量数据向栅格数据的转换、矢量数据向栅格数据的转换 原因：诸如行政边界、交通干线、土地利用类型、土壤类型等数据都是用矢量数字化的方法输入计算机或者以矢量数据的方式存在计算机中，表现为点、线、多边形数据。但是，矢量数据直接用于多种数据的复合分析等处理比较复杂，特别是

31、不同数据要在位置上一一配准，寻找交点并进行分析。92 相比之下，栅格数据模式处理要容易的很多。 注意：矢量数据的基本坐标是直角坐标X、Y,其坐标原点在图的左下角，网格数据的基本坐标是行和列（i,j），其坐标原点在图的左上角，变换时要让直角坐标的X、Y分别与行和列平行。93Field窗口选择数据转换时所依据的属性值，Output cell size栅格的大小，也可以通过浏览某一个栅格数据，来定义栅格大小949596974.2.2数据格式转换数据格式转换1、CAD数据的转换数据的转换 ArcGIS中的要素类、Shapfile数据可以转换成CAD数据，CAD数据也可以转换成要素类和地理数据库。可以将

32、CAD格式的工程图和规划图等放到ArcGIS软件中处理。98数据输出为数据输出为CAD格式格式,Conversion Tools Output Type窗口选择输出的CAD文件的版本 Ignore Paths in Table （可选），在选择的状态下输出单一的CAD文件 Appeed to Existing Files （可选），选择的状态下将输出地数据添加到已有的CAD文件中Seed File在此步为选中状态时浏览确定所需的已有CAD文件99CAD的输入转换的输入转换 Spatial Reference是可选项，用于设置输出地理数据库的空间属性100数据格式在ArcCatalog的转换 通

33、过ArcToolBox/Conversion tools 输入的非ArcGIS格式数据，统一输入到Geodatabase(to Geodatabase) ESRI目前将Shapefile作为交换格式。 在ArcCatalog下面提供了其他的工具：工具条-右键-ArcView8x tools1011024.3数据处理数据处理 在实际应用研究中，根据研究区域的特点，首先需要对空间数据进行处理，以便获取需要的数据。主要包括：裁切、拼接等，借助于ArcToolbox中的工具进行多种空间数据处理操作。1034.3.1数据的裁切数据的裁切1、矢量数据的裁剪、矢量数据的裁剪从大范围里面裁剪出从大范围里面裁剪

34、出 一部分区域一部分区域Tolerance是可选是可选项，用于确定容项，用于确定容差的大小差的大小104 2、栅格数据的裁切，在ArcToolBox中提供了很多剪切方式，其中最常用的是利用现有的栅格或矢量数据裁切栅格数据。105以以Extract by Rectangle为例为例Extraction area是可选项，定义裁切矩形内部还是外部的数据（默认是内部（INSIDE）1064.3.2数据拼接数据拼接 数据拼接是指将相邻的数据拼接为一个完整的目标数据。 原因：研究区域是个非常大的范围，跨越若干相邻数据，而空间数据是分幅存储的，因此要对这些相邻的数据进行拼接。 前提：矢量数据经过了严格的接

35、边。 方法：利用Spatial Adjustment工具可完成数据接边处理。107矢量数据的拼接矢量数据的拼接 Arctoolbox【数据管理工具】 【常规】 【合并】1082、栅格数据的拼接、栅格数据的拼接1094.3.3数据提取数据提取 数据提取是从已有的数据中，根据属性表内容选择符合条件的数据，构成新的数据层。可以通过设置SQL表达式进行条件选择。 包括矢量数据的提取和栅格数据的提取。110矢量数据的提取矢量数据的提取1111122、栅格数据的提取、栅格数据的提取ArctoolboxSpatial AnalysisExtractionExtract by Attributes.11311

36、4ArcToolbox ArcToolBox包含了ArcGIS地理处理的大部分分析工具和数据管理工具。ArcToolBox底下根据不同的功能，分了许多的工具箱。 在地理处理框架中，工具箱是根据核心和扩展功能来划分的。只有激活了License以后，才能使用这些工具箱。 115 可以单独使用这些工具，可以单独使用这些工具，也可以把这些工具连接也可以把这些工具连接起来建空间分析工作流起来建空间分析工作流模型。模型。 可以创建自己的工具和可以创建自己的工具和工具箱。地理处理框架工具箱。地理处理框架还为工作环境提供组织还为工作环境提供组织和管理功能，允许执行和管理功能，允许执行简单和复杂的分析，也简单和复杂的分析，也可以让其他人共享自己可以让其他人共享自己定制的工具定制的工具 116