《地理信息系统基础》实验指导书

《地理信息系统基础》实验指导书实验一GIS演示与操作实验目的通过MapInfo实例的演示与操作，理解GIS的基本概念和普通GIS平台应具备的功能。实验内容1.地理数据可视化2.操作电子地图3.地理查询与分析实验步骤1.进入MapInfo，打开World.tab。2.浏览人口数据。3.显示人口结构分布或老年人口比例专题地图。4.控制图层在缩放范围内显示。5.使用不同的投影方法显示地图。6.使用Info工具查询地理对象的信息。7.建立统计图窗口，作各国的人口增长率、工业增长率及其它统计图。8.调入世界河流图层(Rivers.tab)，查询几条著名河流(如多瑙河，Danube)流经的国家。方法：打开SQL对话框，在表文框中键入World，在条件文本框中键入表达式：objIntersects(selectobjfromRiverswhereR=”Danube”)9.调入各国首都图层(Worldcap.tab)，选择四个城市作半径为1000km的缓冲区，查询这些缓冲区覆盖的人口。实验报告格式1.实验名称、时间、地点1.实验目的2.实验内容3.操作步骤4.实验结果5.实验心得、体会

实验二MapInfo的数据导入实验目的掌握用MIF(MapInfoInterchangeFormat)文件将自己的数据导入MapInfo系统的方法。实验内容1.查看MIF文件格式2.将例题的图形进行手工矢量化3.将手工矢量化的结果构造成MIF文件，导入MapInfoMIF和MID文件简介MapInfoInterchangeFormat(MIF)是一种能完整描述MapInfo数据库的ASCII文件格式。从MapInfo导出数据时，图形数据放在.MIF文件中，表格数据放在.MID(MapInfoInterchangeData)文件中。(一)MIF文件1.MIF文件的文件头MIF文件有一个文件头，头三行是版本信息、字符集、分界符；如果对数据库作了索引，第四行是索引信息，接下来为坐标投影信息(包括边界)，然后是表格文件(数据库)的字段(列)数、各字段的名字及类型，最后一行以Data结束，表明下面各行为数据。Version450Charset"WindowsSimpChinese"Delimiter","Index1,3CoordSysEarthProjection1,0Columns7NameChar(10)NoIntegerCode1SmallintXFloatYDecimal(10,1)DMYDateYNLogicalData上面的文件头中，第四行表明对表格文件的第1、3字段作了索引。第五行是地理坐标投影信息：CoordSysEarthProjection1,0说明采用的是Longitude/Latitude~Longitude/Latitude投影，如是其它投影，则投影代码及表示形式不同，还可包含边界信息，如：CoordSysEarthProjection9,104,"m",110,25,15,40,3800000,2500000Bounds(-16525939.2309,-5282952.71727)(24125939.2309,35368925.7445)对于非地理坐标信息，可使用如下的投影信息(包括边界)：CoordSysNonEarthUnits"m"Bounds(0,0)(500,500)第六行说明表格文件有7个字段，随后是每个字段的名字和类型。MapInfo的字段类型有7种：Char(n)字符型，括号中为字符长度，在MID文件中用“”界定。Integer长整型Smallint短整型Float实数型Decimal(n,m)十进位数型，括号中为字长和小数点后的位数。Date日期型，如2003年3月1日，属性表显示和输入形式是2003-03-01,在MID文件中为20030301。Logical逻辑型，T或F。文件头后即是数据体。2.点对象(符号)对于每个点对象，有2行，第1行以Point开始，后为一对坐标数据，随后一行以Symbol开始，指明点对象(符号)的类型、颜色和大小。如无Symbol，则采用默认类型、颜色和大小。Point3799847.114256235.84Symbol(52,255,9)Point2.09009-0.530879Symbol(34,16719904,12)除用MIF文件外，有地理坐标的点数据还可直接用<表(图表文件)>菜单中的“创建点(制造图点)”选项导入；对于没有地理坐标的点数据，可通过地理编码赋予这些点地理坐标，成为点对象。3.线对象线对象包括直线、折线和圆弧。直线用Line标识，后跟直线的起点(x1,y1)和终点(x2,y2),如：Line-1.1905540.301614-0.791748-0.00516Pen(2,7,16776960)其中Pen指明了线的宽度(代码)、类型和颜色。如无Pen，则采用默认宽度、类型和颜色：Pen(1,2,0)宽度代码为1(像素为1)，类型为2(实线)，颜色为黑色。折线包括两种类型的折线，一种为无分叉的折线，用PLINE标识，后该线的点数，随后是每个点的坐标(每个点1行)；另一种为有分叉的折线，如有支流的河流，这些支流可组合在一起作为一个线对象。后一类型以PLINEMULTIPLE为标识，后跟折线数，如PLINEMULTIPLE5表示该线对象有5条折线组成，随后为每条折线的点数和点坐标。也可用Pen指明线的宽度(代码)、类型和颜色。PLINE3124.00678253173856.4725761413574123.25955963134856.6561851501465123.84651947021557.5234107971191PLINEMULTIPLE45-61.2772979736328-3.60705971717834-60.6846580505371-3.53694128990173-60.4775352478027-3.28227639198303-59.9512214660645-3.26529860496521-59.6999473571777-3.085333585739143-58.5025100708008-3.27193522453308-59.4667854309082-3.006098985672-59.6999473571777-3.085333585739144-67.7270431518555-3.05758810043335-67.7542037963867-3.28567099571228-69.1090316772461-3.43507599830627-69.295783996582-4.036087989807132-64.5731887817383-3.14008784294128-63.9511756896973-3.63880848884583圆弧圆弧用Arc标识，后面的数据依次是圆弧椭圆最右(或最左)点的X坐标，圆弧椭圆最上(或最下)点的Y坐标，圆弧椭圆最左(或最右)点的X坐标、圆弧椭圆最下(或最上)点的Y坐标，圆弧在其椭圆中的起始角、终止角。如：Arc0.636941-0.3495130.8648330.1763832700Pen(1,2,0)Pen的意义同上。4.面对象包括多边形(区域)、椭圆、矩形、圆角矩形。多边形多边形对象以Region标识，后跟该多边形对象的多边形个数，然后是每个多边形的点数和各点的坐标。如:Version300Charset"WindowsSimpChinese"Delimiter","CoordSysEarthProjection12,62,"m",0Bounds(0,0)(500,500)Columns2PolygonChar(10)HighDecimal(5,0)DataRegion152506436012525023316212525064Pen(1,2,0)Brush(16,65280,16777215)Center250130Region2525099297125250172213125250995250117265125250140238125250117Pen(1,2,0)Brush(16,65280,16777215)Center250130每个多边形对象后的Pen指明多边形外框线的特征，意义同上；Center指明该多边形对象中心的X、Y坐标；Brush中的数据依次为填充的图案、前景色、背景色。如没有Pen、Center、Brush数据，均采用默认值。椭圆椭圆以Ellipse标识，后跟的数据意义同圆弧中头四个数据，Pen、Brush意义同上。Ellipse-0.6578070.0697670.750830.707641Pen(1,2,0)Brush(2,16777215,16777215)矩形椭圆以Rect标识，后跟两个角点的坐标(X1,Y1,X2,Y2)，Pen、Brush意义同上。Rect-1.089701-0.50038-0.24832-0.189369Pen(1,2,0)Brush(2,16777215,16777215)圆角矩形圆角矩形以Roundrect标识，头四个数据亦为两个角点的坐标(X1,Y1,X2,Y2)，最后一个数据为圆角半径，Pen、Brush意义同上。Roundrect-0.870633-1.375010.224989-0.9279990.26295Pen(1,2,0)Brush(2,16777215,16777215)5.文本对象文本对象以Text标识，下跟3-6行内容，第1行是文本内容(用“”界定，如换行，加\n)，第2行为文本对象未旋转前两个角点的坐标(X1,Y1,X2,Y2)，第3行Font指明了字符集、字体效果(粗体、斜体、加底线等)、样式、字体颜色、文本背景样式与颜色(也可没有该项)。如文本旋转了角度，则用一行指明Angle；当文本不是左对齐时，用对齐方式Justify说明；还可用Spacing、LabelLine指明行距、标注线特征等。例："Xyz\n内蒙古"1.3732-1.1939091.954369-0.761465Font("幼圆",289,0,16719904,65535)Angle45JustifyCenterSpacing1.5LabelLineSimple1.740617-0.699837(二)MID文件MID文件存放表格数据，项目数与类型要与MIF文件中的定义的字段对应。如上例，在MIF文件中有：Columns2PolygonChar(10)HighDecimal(5,0)则MID文件内容可为："多边形1",450"多边形2",350"多边形34",250实验步骤1.进入MapInfo，打开World.tab、Rivers.tab、Prc.tab、Prc_hiwy.tab、Prc_air.tab。2.将上述文件中的数据分别转出到MIF文件(在图表文件菜单中选择“导出图表文件”)。3.通过这些MIF文件观测点、线、多边形等对象数据的格式。4.将下图中的线和多边形在屏幕采样(手工矢量化)，另加2个点数据POINT285400和POINT250125和Text对象。5.根据采样数据构造MIF文件。6.将构造的MIF文件调入MapInfo显示出来。实验三空间数据和地理编码实验目的理解GIS空间数据编码和地理编码的原理，掌握建立链状双重独立式数据文件的方法和用MapInfo进行地理编码的方法；掌握将点对象直接导入MapInfo的方法。实验内容1.对实验二中的多边形数据建立链状双重独立式数据结构中的四个数据文件，即多边形文件、弧段文件、弧段坐标文件、节点文件。2.利用实验二中的数据进行地理编码。3.将深海钻孔数据(allcore.txt)用点对象直接导入法导入MapInfo.三、地理编码方法简介1.何谓地理编码地理编码是向文件或数据库增加地理信息从而使对象可被显示在地图上的过程。该文件必须包含具有地理性质的文本数据（如州、县、邮政编码或街道地址）。在地理编码中，MapInfo取得这些文本数据并将它们与(已地图化的表的)地理信息关联，从而使文件中的对象可被显示在地图上。2.对表地理编码(1) 至少打开一张已地图化的表。该表在地理编码操作中用作地理信息源表（如用地址作索引的街道地图）。(2) 将源表中要用于地理编码的字段作索引。有关索引，详见帮助主题“创建表或修改表”。在“文件”＞“新建表”和“表”＞“维护”＞“表结构”中都有选项对表建索引和使其可地图化。(3)打开要地理编码的表。该表必须与源表不同。(4)选择“表”＞“地理编码”，出现“地理编码”对话框。(5)指定合适的地理编码选项。3.关于表MapInfo数据库是以表来组织的。表由行和列组成。每行包含关于特定地理特征、事件等的信息。每列包含关于表中的数据项的特定类型的信息。MapInfo表可以包含表示地理对象的图形对象。这样的表可被显示为地图。参见“源表”和“查询表”。四、实验步骤1.设在实验二中利用MIF文件已经将数据导入MapInfo，建立了test1.tab等文件，先将test1.mid文件完善。如"多边形1",450"多边形2",350"多边形3-4",250"线1",1"线2",2"线3",3"线4",4"线5",5"线6",6"线7",7"线8",8"线9",9"线10",10"线11",11"线12",12"线13",13"线14",14"线15",15"线16",16"线17",17"点1",1"点2",22.另建立一个文本文件(如test11.txt)，用于地理编码，如："对象","编号""多边形1","450""多边形2","350""多边形3-4","250-150""线1","L1""线2","L2""线3","L3""线4","L4""线5","L5""线6","L6""线7","L7""线8","L8""线9","L9""线10","L10""线11","L11""线12","L12""线13","L13""线14","L14""线15","L15""线16","L16""线17","L17""点1","1001""点2","1002"3.进入MapInfo系统，打开test1.tab，对该表的第1列作索引,即以该表作为源表。4.打开test11.txt(从数据库中调入)，对该表作地理编码(在“图表文件”菜单中)。地理编码完成后，就生成test11.tab。5.在图层控制中调入test11.tab，将该表的第2列作为标注显示在图上。6.打开world.tab、grid15.tab，选择mollweide(equal-area)投影。7.调入allcore.txt文件，在“图表文件”菜单中用“制造图点”项生成allcore.tab，保存工作区文件(worldcore.wor)。8.比较、总结三种将点对象导入MapInfo方法(通过MIF文件、通过地理编码、直接导入法)的特点。实验四利用MapInfo进行栅格图像配准和建立GIS数据库实验目的1.掌握栅格图像配准方法2.掌握利用MapInfo建立GIS数据库的基本方法3.熟悉MapInfo的空间SQL查询方法二、实验内容1.将图像文件SG1.JPG在MapInfo下配准2.对SG1.JPG中的高速公路层跟踪矢量化3.对高速公路层建立属性库，具体属性项如下道路编码(ID)、道路名称、道路类型、道路长度4.利用SQL对话框进行SQL查询三、图像配准简介栅格图象是由多行微小的点（象素）组成的一种计算机化的图象。如果手头有扫描仪及扫描软件，可以通过扫描一幅纸张地图来创建栅格图象。完成地图扫描并将其保存于文件中后，即可在MapInfo中显示该文件。有多种不同的栅格图象文件格式。MapInfo能够处理以下格式的栅格图象文件：JPEG、GIF、TIFF、PCX、BMP、TGA（Targa）和BIL（SPOT卫星图片）。配准一幅栅格图象时，要输入地图坐标（如经度／纬度），并指定栅格图象上与该坐标对应的点(至少3个点)。要在MapInfo中显示栅格图象前必须进行配准，以使MapInfo在显示图象时能够完成地理计算，如计算距离和面积等。在MapInfo中首次打开一幅栅格图象时，MapInfo显示“配准栅格图象”对话框。填写该对话框以告知MapInfo如何配准图象。MapInfo将栅格图象配准信息保存在表文件中以供以后使用。下一次打开该栅格图象表时就不必再进行配准了。这样，只须对栅格图象进行一次配准。四、实验步骤1.打开Prc.tab、Prc_ctya.tab、Prc_ctyb.tab、Prc_hiwy.tab。(采用Longitude/Latitude投影)2.调入SG1.JPG进行配准，参考配准点坐标：北京116.90,39.90天津117.19,39.13唐山118.18,39.62沧州116.85,38.30廊坊116.67,39.51任丘116.08,38.703.新建道路表，对高速公路层进行屏幕矢量化(构造MIF文件)4.完善道路表的属性数据库(构造MID文件)5.进行SQL查询

实验五空间分析综合应用实例一、实验目的了解并初步掌握GIS空间分析原理，学会运用MapInfo软件进行空间分析操作并解决综合性的实际问题。二、实验内容利用MapInfo和Access数据库软件进行空间分析、导出必要的数据、计算分析问题所需的证据，然后分析计算结果，得出自己的结论或认识。三、问题及分析1.问题Atlanta市的公园分布与居民收入的关系(是否偏向富人居住区?)。2.资料分析现有主要资料：Atlanta市的普查地段(Tract)图、公园分布图。普查地段图表示了Atlanta市482个普查地段的位置与范围，共有47个属性项，其中与本问题有关的主要属性有:No_tractDecimal(5,0);普查地段编号，为普查地段唯一标识符POP_93Decimal(6,0);93年人口数POP_98Decimal(6,0);98年人口数AVG_INCDecimal(10,4);93年普查地段户均收入MED_INCDecimal(10,4);93年普查地段户收入中值PER_CAPINCDecimal(10,4);93年普查地段人均收入根据这些属性可以得到全市93年人均收入，据此确定每个小区的相对贫富位置。公园分布图表示了Atlanta市内及周边47个公园的位置与范围，主要属性是：No_parkDecimal(5,0);公园编号，为公园唯一标识符缺少资料：每个地段的几何中心，每个公园的几何中心。可以用它们代表普查地段和公园的位置及计算地段与公园的距离，并据此将地段分为离公园远的和近的。这些资料可通过MapInfo导出得到。3.解答方案根据现有资料，采用1993年的人口及收入数据。先将各地段根据人均收入分为收入高(4/3全市人均收入)的和中低的(<4/3全市人均收入)两类，它们的总数分别用TH和TML标识，再根据各地段离各公园的平均距离或中值距离分为远和近的两类，其总数分别用TF和TN标识。需计算的证据数据为：收入中低且离公园远的地段占收入中低地段总数的百分比(P(TMLTF))，收入高且离公园远的地段占收入高地段总数的百分比(P(THTF))；收入中低且离公园远的人口占收入中低人口总数的百分比(P(POPMLPOPF))，收入高且离公园远的人口占收入高人口总数的百分比(P(POPHPOPF))(注：由于缺乏详细的人口收入数据，这里的收入高低人口并不是真正意义上的各类收入人口统计数，而是根据地段类型将各地段的所有人口都划归某一收入类型);离公园远和近的居民的人均收入(和)。4.结果分析参考P154的思路分析计算结果，得出自己的结论或认识。四、实验步骤(一)导入Atlanta市的基础图件1.通过At_tracts.mif和At_tracts.mid生成Atlanta市的普查地段(小区)图，查看其属性数据，记录数据缺失点。2.导入其它图件At_parks、At_streets、At_high。(二)导出Atlanta市的普查地段和公园的几何中心数据1.导出Atlanta市的普查地段几何中心数据和其它相关字段(用下面的SQL查询语句)Fields{No_tract,Val(Format$(CentroidX(obj),"###.######"))"经度",Val(Format$(CentroidY(obj),"###.######"))"纬度",area(obj,"sqkm")"面积",Perimeter(obj,"km")"周长",POP_93,POP_98,AVG_INC,MED_INC,PER_CAPINC}Tables{At_tracts}Where{}Group{}Order{}Into{selection}Browse生成查询文件后导出到Access数据库(表名At_tracts2)。2.导出Atlanta市公园的几何中心数据(用下面的SQL查询语句)Fields{No_park,Val(Format$(CentroidX(obj),"###.######"))"经度",Val(Format$(CentroidY(obj),"###.######"))"纬度",area(obj,"sqkm")"面积",Perimeter(obj,"km")"周长"}Tables{At_parks}Where{}Group{}Order{}Into{selection}Browse生成查询文件后导出到Access数据库(表名At_parks2)。(三)在Access中整理数据并导出txt文件并在MapInfo导入之1.删除At_tracts2中多余的点，导出文本文件At_tracts2.txt。2.导出文本文件At_park2.txt。3.在MapInfo中将该两文件导入(用制造图点法)(四)将几何中心数据用长度单位(m)表示并导出到Access由于原几何中心是用经纬度表示的，不利于计算距离，故要将At_tracts2、At_park2用Non_Earth/Non_Earth(meters)投影，再用SQL语句Fields{No_tract,Val(Format$(CentroidX(obj),"###.######"))"X",Val(Format$(CentroidY(obj),"###.######"))"Y",面积,周长}Tables{At_traks2}Wher