第十四章 地理信息系统应用实例 导读：本章介绍了地理信息系统在一些具体领域的应用。 由于GIS是用来管理、分析空间数据的信息系统，所以几乎所有的使用空间数据和空间信息的部门都可以应用GIS。由于各个部门的不同，GIS在具体业务系统中所占的比重、应用方式也各异，结合比较紧密的如城市规划、环境领域等，相对松散的有商业、医疗卫生领域。 本章简单介绍了一些GIS应用实例，针对不同实例，介绍方法也不同，有的概述一个具体区域的应用，有的叙述工作流程，有的介绍了工作原理和分析方法，可以对相关领域的GIS建设提供借鉴，也可以作为其它领域建设GIS的参考。 本章也引入了一些新的概念，如AM/FM、GIS/T以及一些具体的分析方法，如动态分段模型，它们处于GIS与具体领域的交叉部分，同时也是GIS理论和技术的有益补充。 1．城市规划、建设管理 城市是人类活动高度集中的区域，同时也是信息、物质高度集中的区域。随着科技的进步和经济的发展，城市系统越来越复杂，数据和信息越来越多，服务要求越来越高。城市管理面临着新的挑战，为了城市的现代化、生态平衡和持续发展，城市需要全面的规划，而地理信息系统给城市的规划和管理带来了新的工具。 在城市管理中，可以应用GIS的方面非常多，如土地、道路、管网、环境、人口等等诸多要素都可以通过地理信息系统进行管理*，并且形成了应用于专门领域的GIS工具，如土地信息系统（LIS-Land Information System）和AM/FM（自动制图/设施管理，Automatic Mapping/Facility Management），前者侧重于管理地块的位置、面积、权属以及地籍的动态变化；后者则主要对电力、通讯等设施的管理和分析、制图。 城市建设规划涉及的因素非常多，开发新城要征用土地，改建旧城要拆迁安置，同时需要基础设施、公共服务设施的配套。在开发建设活动中，如果不注意各工程项目之间的协调，就可能造成混乱，而采用GIS对各种信息进行管理，并基于此进行分析和辅助决策，可以有效地防止这种混乱局面的出现。由于城市在不断地建设发展，所以需要随时更新城市基础数据库，这就要求应用GIS管理日常城市建设活动，以保证信息的时效性。下面介绍GIS在密尔沃基（Milwaukee）市城市建设日常管理方面的应用。 美国威斯康星州的密尔沃基从1976年开始建设城市地理信息系统，并不断发展、完善。由于数据库比较完整，地理信息系统在很多方面发挥了作用，主要可以分为三个方面的活动。 1）地图更新 密尔沃基市常用地图分为三种： 第一种是地籍房产图（当地称为Plot Map），由政府税务部门负责编制。这种地图通常在以下情况下需要修改更新： 地契内容改变，包括地块边界的变化以及地块权属的变化； 在一个大的地块内再划分为若干小块，产生新的土地所有者； 公共道路的拓宽或者连通； 小地块合并成为大地块； 其它变化，如道路名称、地名的改变以及错误的纠正。 第二种是一般的测绘图（当地称为Quarter-Section Map）。这种图由政府市政工程部门负责编制，图面上标有主要的地块边界（不作为法律依据），但是没有普通的建筑物，通常在下列情况下需要更新： 城市建设造成地面重要物体的边界变化； 重新测量，纠正过时的内容； 道路打通、拓宽或封死； 大的地块合并或重新划分； 街道名称改变； 地图上的错误纠正。 第三种是土地使用图（当地称为Land Use Map）。这种地图反映土地的实际用途，也反映了主要的地块边界。土地使用图由城市规划部门负责编制，主要为城市规划服务，通常在以下条件下需要更新： 房屋的拆除和重建； 房屋改变用途； 地块边界改变； 地块重新划分； 地块兼并； 道路变化； 道路名称的修改； 地图错误纠正。 以上三种地图由三个相互独立的部门各自负责，在信息上存在很多重复内容，在地图更新时也需要大量的重复工作。在建立城市地理信息系统之后，三种地图上各类信息的更新工作明确地分配到三个部门，三种数字化地图集中地存放于数据库中，实现数据共享，每个部门对数据的修改内容可以被其它部门得到，减少了数据冗余和重复工作。 2）土地区划管理 在美国，大多数地方政府都用土地区划（Zoning）来指导、限制城市土地利用和城市发展，但是具体做法各地有许多差异，在密尔沃基市，先由城市规划部门编制土地使用规划，经批准后再制定区划，区划的形式为文本和图件，它是城市建设和审批修建申请的依据。 密尔沃基市的土地区划图分为用途图、范围图和高度图，用途图用于限制和规定相冲突的土地利用，如不准贴近主要交通干道修建大型商业设施等；范围图对建筑之间的距离以及建筑后退道路等作出了规定；高度图则规定了建筑本身的高度、相邻建筑在高度上的相互关系的限制。 在采用GIS后，这三种区划图分为三个图层进行管理，并可以叠加在一起显示，避免各个部门数据不一致的情况，此外，可以很方便地检查规划图的错误，避免法律上的纠纷；也便于税收部门查询地产信息，调整对土地所有者的税收。 3）建筑审批处的内部工作管理 每个建筑物在修建之前，必须经过政府主管部门批准，由于申请非常多，并且需要实地勘察，工作繁重，如何分配并确定每个经办人员的工作，以充分调动每个人的积极性成为一个难题。在建立地理信息系统后，修建项目申请卡片被记录到属性数据库中，其中的建筑物地址属性与地块空间数据相关联，这样可以很方便地将一定时期内所有的修建申请的位置分布显示在地图上（图14-1），可以方便的根据申请量的多少、到市中心的距离、前往勘察的边界程度，划分每个工作人员的负责范围，提高了工作效率，改善了建筑审批处的内部管理。  图14-1：利用GIS进行建筑审批工作协调 2．农业气候区划 2．1项目背景与需求 自八十年代以来，中国农业生产环境、气候环境发生了巨大变化，包括： 1）农业科学技术进步，高新技术引进，农业朝高产、优质、高效、低耗方向发展；种植业向粮食—经济—饲料作物三元结构转变；小生产、大市场向规模经营方向发展，逐步实现种养加、产供销、贸工农一体化。农村产业结构调整向农业气候区划提出了新的要求。 2）气候条件与气候资源本身发生了变化。如东北地区平均气温升高，北方地区干旱范围扩大，长江流域洪涝增多，特别是九十年代以来，异常气候事件呈现明显增多的趋势。有必要重新认识气候资源的变化及其合理利用和保护问题。 3）农业气候区划技术条件发生了巨大变化。“3S”（GIS、RS、GPS）与网络技术在农业气候资源动态监测与开发应用中展示出广阔前景。 4）农业生产的发展，对气候资源开发利用和保护提出了更高要求。气象部门必需依靠科技进步，技术创新适应社会的需求。建立一套基于“3S”、网络平台的区划信息系统，有助于广大气象台站在气象服务手段和技术上得到提升。 5）面对21世纪农业新技术革命和可持续农业战略，农业气候资源作为一种重要的投资环境，有必要进行科学、客观的评价。 中国气象局提出建设“第三次农业气候区划”项目的目的是：采用新技术、新方法、新资料，开发“农业气候区划信息系统（Agriculture & Climate Distributed Information System,简称ACDIS）”软件，建立气候资源开发利用和保护监测体系，实行资源平面与立体，时间与空间全方位优化配置；发挥区域气候优势，趋利避害减轻气候灾害损失，提高资源开发的总体效益。为各级政府分类指导农业生产，农村产业结构调整，退耕还林防止水土流失等提供决策依据，为地方政府服务。 2．2农业气候区划信息系统（ACDIS）系统结构及工作流程 2．2．1系统结构 农业气候区划信息系统（ACDIS）是基于GIS工具平台，建立起面对专业技术人员的专用工具，适用于农业气候资源监测评价、气候资源管理与分析，小网格气候资源推算与空间查询、省地县三级区划产品制作等。具体实现了以下功能： 1）地理基础信息管理：管理工作区的基础地理数据，如行政区划，水系，交通数据等； 2）小网格资源信息管理：管理栅格格式的、与农业气候区划有关的信息，包括各种DEM数据； 3）小网格资源推算与区划产品制作； 4）农业气候资源监测与评价； 5）农业气候区划成果演示。 2．2．2系统工作流程 根据“3S”等新技术在区划中应用需求，建立了农业气候区划工作基本流程（如图14-2）： 图14-2：农业气候区划工作基本流程图 2．3“3S”在项目中的应用 1）利用GIS对农业气候区划综合要素空间查询和管理 ACDIS利用GIS工具提供的基本功能，对其它子系统输入、处理和生成的气候资源等数据进行综合查询和管理。生成不同查询条件下的区划产品，作为区划专题内容，进而实现对区划各类产品的矢量图、栅格图、DEM、注记和属性数据进行以地理表达式为条件的逻辑查询以及不同图件和属性数据的综合查询、管理，并把查询结果制图输出到绘图仪或打印机，或保存为其它格式文件。 2）利用GIS对气候资源进行小网格推算模式研究 气候资源小网格推算模式研究是区划的一项基础性工作。在收集山区气候研究成果基础上，辅助以气象哨、水文站雨量资料后，通过GIS平台可快速计算和获取测点地理参数（高程、坡向、坡度）。采用统计分析方法，根据要素的统计特征值和地理特征将全省划分若干气候区，分别建立各区域气候要素推算统计模式，通过验证和残差订正，应用到气候资源小网格推算中。 3）GIS在农业气候资源分析中的应用 应用GIS来定量采集、管理、分析具有空间特性的气候资源，建立GIS分析气候资源的思想、方法、步骤，包括数字高程模型建立、GIS农业气候资源数据建立、空间分析模型建立、气候资源分析计算、气候分区及定量分析等。 4）利用“3S”提取农业背景信息参与区划计算 农业气候区划是根据农作物生长发育过程中对气候条件的要求和气候资源的地理分布特征来进行分区划片的，在某种农作物的气候可种植区内还有不同的地物类型，不同的农作物要求不同的地理环境。为使农业气候区划对农业生产更具有指导作用，将非气象因子引入到农业气候区划中。农业气候区划对象中往往对土壤PH值要求很高，根据土壤类型分布可以得出土壤PH值的分布，将其作为区划的一个关键指标，使得区划更加有实际应用意义。利用GIS将土壤分类图作为一项数据层参与气候资源数据层集运算，得出包含土壤类型信息的区划结果。 江西省在全省优质早稻种植气候区划和万安县脐橙种植综合区划中（图14-3），除了应用1∶25万的地理数据，还结合了TM影像数据，辅助GPS定位抽样，把早稻、脐橙的可能种植区（农田、荒山荒坡）提取出来，排除了山体、水体、居民点、道路等不能种植脐橙和早稻的区域，把可能种植区与农业气候区划图做逻辑交集运算，得到了全省优质早稻和万安县脐橙种植规划图。  图14-3：江西省万安县脐橙种植区划 5）利用GIS建立集区划—资源动态监测—高产栽培技术为一体的信息服务系统 利用GIS工具，开发出基于Web的、对气候条件敏感的两系杂交稻制种气候资源空间配置信息系统，将区划服务于作物生产基地选择、关键生育期气象服务、农业生产技术指导等各个方面。 3．大气污染监测管理 随着经济的发展，环境污染直接影响了人们的生活质量，环境质量问题也得到了越来越多的重视。污染环境包括水污染、大气污染、固体废弃物污染等，其中就大气污染而言，城市区域由于受到工业生产、居民生活的影响，成为大气污染发生的集中区域，历史上几次严重的污染事故，如伦敦烟雾事件（1952）、洛杉矶光化学烟雾事件（1943），都是发生在大城市。近几十年来，研究者对大气污染问题进行了大量研究，并且通过实验或计算来建立适合于特定区域的大气污染物扩散模式以及确定相关参数的计算方法。 城市大气污染的来源主要包括点状污染源和线状污染源，前者主要包括烟囱，后者则指汽车尾气排放。而污染的扩散的影响因子，除了排放量、排放物之外，还受到气象条件、下垫面等因素的影响。通过大气动力学的研究，研究者给出了用来描述污染物在大气中扩散规律的各种解析方程。这些大气扩散公式为空气污染预报提供了一定的基础，可以为现有污染源条件在不利的气象条件下减少有害污染物的排放、在城市规划中合理进行区域规划以及环境影响评价提供了切实有效的信息。 无论是点源污染，还是线源污染，其空间分布以及属性可以通过地理信息系统进行管理，而污染扩散的影响因子的空间分布同样可以作为GIS的空间数据组成部分，所以，基于GIS可以建立大气污染扩散模型，进而，GIS也提供了丰富的功能以表现污染物强度空间分布，可以查询强度分布状况，并可以结合其它社会经济数据，进行更加细致的评价分析。 下面介绍内蒙古自治区包头市在利用GIS进行大气污染扩散模拟中的应用。 包头市是中国最大的稀土工业基地和著名的钢铁机械工业基地之一，是门类齐全，体系较为完善的现代化工业城市。包头市工业能源消耗以燃煤为主，占56%，其次为焦碳和电力，在能源结构中，重中污染能源占70.8%，包头工业排放的SO2、烟粉尘和氟化物都较大，造成包头市煤烟型和氟污染的特点。其中21家重点空气污染源占全市工业排放的95%，85%，96%（分别为SO2、烟粉尘和氟化物）。其中包钢是最大污染厂家，其SO2、粉尘、氟化物的排放量高居榜首。 从自然条件上看，包头市地处内陆，属内陆半干旱、中温带大陆性季风气候，全年干燥，无霜期短。年均气温6.5℃，取暖期6个月。由于受地形影响，冬季包头多北风和西北风，夏季多东南风，各月夜间均盛行NWW风。近年来城市建筑的快速发展及自然植被，人工绿化面积的不断扩大，市区平均风速呈逐年下降趋势。1995年年均风速为2.2m/s。一年中，春季风速最大，秋季风速最小，冬季风速居中。除七月份主导风向与最大风速风向恰好相反外，常年盛行风向与最大风速风向重合。 逆温层影响着污染物的扩散、稀释，包头市冬季逆温频率很高，夏季逆温频率较低，逆温的强度与厚度，冬季均明显大于夏季。 包头是我国大气污染治理的重点城市，包头关于大气污染扩散的研究工作较多。冶金部建筑研究总院1982针对包钢地区的烟气综合治理规划，利用风洞模拟试验、现场实验等提出了“大气输送气候学模式”（ATCM）。1989年包头市环境监测站，针对包头新市区大气扩散模式和SO2容量计算，提出了基于美国EPA 的ISC（工业复合源大气扩散模式）的城市多源高斯模式。这些模式的建立为包头市的大气污染治理和管理提供了可靠的依据。 基于GIS建立大气污染扩散模型，可以 1）模拟污染物的空间分布，评价不同区域的环境质量； 2）将污染物空间分布与人口密度空间进行复合分析，确定受污染影响的人口数目； 3）预测在给定气象条件下污染物的空间分布； 4）确定不同点源对整个研究区污染总量的贡献；进而 5）为污染整治，如降低排放量、甚至关闭某些污染源，提供决策依据； 6）如果要增加污染点源，可以比较不同的方案（如烟囱的位置，高度等），从中选择最优方案； 7）在城市规划时，作为确定不同用地（居住、工业、商业等）的分布。 在包头市的研究工作中，利用1月份平均风速、风向、频率，并将其换算为风频表，对包头市的37个高架点源造成的地面SO2浓度的空间变化进行模拟，结果如图14-4所示。  图14-4：包头市大气二氧化硫分布图 将模型预测结果图与包头市环境监测站绘制的等值线图相比较。模型在工业区的预测值比较切合实际，在昆都仑区，预测值偏小，原因应该是由于包头地区在一月份特定风向条件下，工业污染对该区域的影响比较小，相比之下居民取暖燃煤造成的二氧化硫污染就较突出。总体来看，无论从工业区还是居民区，模型预测的二氧化硫的浓度不存在数量级上的差别。 在大气环境规划中，对规划中欲添加的新污染源需要预测新增污染源造成的污染、对新源的进行合理选址、确定新源的排污量大小以及点源烟囱的高度等。在研究中，拟在青年农场（图上三角点位置）增加一个新的污染源，有三种方案，分别为源强为100000mg/s、烟囱高150米；源强1000000mg/s、烟囱高150米；源强1000000mg/s、烟囱高200米。在一月风向风频条件下进行预测，将预测结果叠加在原来的污染分布状况上，结果如图14-5。  图14-5：增加不同参数污染源的污染分布图 大气污染的受害人群是大气污染管理和控制中最重要的关注对象，结合包头市新市区（青山区和昆都仑区）的人口分布，研究中分析了添加新源前后大气中SO2污染受害人群的变化。研究方法是将模型的预测结果与包含人口数据的行政区划图进行叠加复合分析，获得各个污染区的面积和人口，结果见表14-1。 表14-1:新添加污染源造成的污染人群变化结果表（单位：万人） 浓度单位：mg/m3 0.00-0.05 0.05-0.15 0.15-0.25 0.25-1.6  新源强：100000mg/s 新源高：150米 46.05 15.33 2.10 0.80  新源强：1000000mg/s 新源强：150米 36.20 24.30 2.73 1.04  新源强：1000000mg/s 新源强：200米 41.02 20.02 2.33 0.89  未添加前大气状况 49.83 15.61 2.03 0.80  从图和表可以看出，源强为1000000mg/s的污染源造成的高浓度受害人群数量明显增加，而源强为100000mg/s的新添污染源造成的高浓度污染受害人群数量变化不大。同时可以看出，增高污染源的高度来降低污染的危害是以低浓度受害人群数量的增加为代价的。从新增污染点源的污染分析来看可以考虑在青年农场位置新建立一个源强为100000mg/s、烟囱高为150米的新点源。 4．道路交通管理 近年来，GIS在交通方面的应用得到了广泛的重视，并形成了专门的交通地理信息系统GIS-T，以满足道路交通管理方面的要求。下面分几个方面介绍GIS在公路方面的具体应用。 4．1路廓设计 路廓设计是公路设计中的一个重要环节，是定出公路最终线向的一个步骤。在路廓设计中，要综合分析多种空间数据，包括大比例尺的土地利用图、地形图以及现有的道路网等。 从各方面收集来的资料合并到GIS资料库中，建造约束多边形，确定公路需要避开的区域，通常这些多边形可以分为三级： 1）第一级约束是必须要避开或减至最少； 2）第二级约束是尽量避开或减至最少； 3）第三级约束是可以避开或减至最少。 根据需要可以用权重表示不同的约束等级，并在地图上分别显示，以便设计员在设计路线时尽量绕开这些区域。然后将每条试验性线路输入到GIS中，利用缓冲区分析得到路廓多边形，与约束多边形进行叠加分析，计算路廓与每个多边形相交的总面积，并乘以权重，得到该试验线路的约束影响的加权总面积（Total Impact Weighted Area, TIWA），用于比较这些线路设计。 在计算得到TIWA之后，可以根据DEM进行线向的详细设计，其主要操作是根据DEM得到道路横断面，进而进行填挖方的优化，得到最终的路廓设计结果。 4．2道路管理 4．2．1道路管理涉及的属性 在公路管理中，通常需要根据其类别、车道数目和路面种类进行分类，其管理的属性包括： 1）公路属性：路网标识，道路标识，路段标识，参考点等； 2）几何属性：路面宽度，路向，车道数目，坡度等； 3）设施属性：中央隔离带和隔离物，护轨，桥梁，交通标志，路灯，交通量计算仪； 4）公用设施：收费所，交叉口，涵洞，交通管理站等； 5）建筑材料属性：面层，覆盖层，基层，路基； 6）道路使用统计属性：交通量，每天平均流量，事故，车重统计，车类统计，进路要求； 7）验收记录属性：强度，粗糙度，障碍物，路面损坏程度，限速区段，车速调查； 8）合同属性：预算，费用，图则，日期等； 9）维修工程属性：工程类别，预算，费用，数量，日期。 这些属性是通过线性定位方法与路网建立关联，通常线性定位的方法有里程碑法和控制段 法，前者采用路名和公路上的里程点来确定物体的位置；后者则将公路分成相连的、长度不等的控制段，每段有一个编码，并且其属性一致（表2）。 表14-2：采用控制段方法记录属性 控制段 表层 车道数目 状态 起点 终点  CS1 砾石 2 坏    CS2 沥青 4 良好 30.08 30.16  CS3 沥青 2 尚好 30.16 30.2  CS4 沥青 2 良好    4．2．2动态分段 交通模型的一个特点是，多种线性物体，重叠在同一个路网上，如限速区段，公共交通线路都与路网重叠，这使得GIS在应用于交通时遇到困难，而动态分段（Dynamic Segmentation）方法较好地解决了该问题。 动态分段是在数据库中纪录道路的每种属性的起止点到道路原点的距离，并不是真的将道路切断存储，适合于动态的分析，顾名动态分段。 采用动态分段之后，一个路段（Segment），是路网上两个交点间连线或者弧的一部分，路段的长度用其占连线的比例来表示，具有唯一的标识码。在GIS数据库中，路段是依附于路网数据，本身没有坐标。由于采用动态分段将道路的各种属性以及其分布集中在一个图层中进行管理，采用线性定位方法，因而容易实现各种“点线”以及“线线”的叠加查询分析，形如： 交通流量>800,000 AND 路面状况=“良好” AND 道路分级=“省道” 4．3流量和路径分析 4．3．1道路网络拓扑 一个路网的表现方式可以偏重其“几何性”或者“拓扑性”，在道路设计中，需要用到其几何表现，而网络分析则着重于拓扑关系（图14-6）。  图14-6：路网的两种表现方式 在交通网络中，大多数线与线之间的交点是具有拓扑性的交点，即在道路的交点处车辆可以转向；但是，由于在GIS中，使用的是平面图形（Planar Graph），而有些道路的关系需要在三维空间中才能够正确表现，如立交桥，实际上道路并不直接相交，但是在平面图形上却表现为一个交点，这需要在应用时进行判定，区分拓扑交点和非拓扑交点，以防止出现错误的网络分析结果。 4．3．2用于道路网络分析的数据结构 进行道路网络分析时，可以只关注其拓扑数据表现，其数据结构和相关算法可以采用已经成熟的“图”数据结构，一种最为简单的实现方式就是“连通矩阵”，记录了结点与结点之间的连接关系（图14-7）。  图14-7：网络连通矩阵 在连通矩阵中，单元格数值为1表示两个结点直接相连，反之表示不直接相连，这样的连通矩阵可以用于进行连通性、可达性分析；矩阵中每个单元格的数值还可以是结点之间的距离，不直接相连的结点间距离可以用无穷大表示，利用该矩阵可以进行最短路径搜寻。在不考虑连线的方向性，即图是无向的，连通矩阵是对称阵，而在城市道路交通中，常常会出现“单行线”的情况，即连线是有方向的，这样矩阵是非对称的。 在城市交通应用中，还会出现“禁止左转”、甚至“禁止右转”的情况，这时，除了连通矩阵之外，还需要记录每个结点——即路口——的连通属性，同样可以用连通矩阵实现（图14-8）。  图14-8：路口的转向和道路连通矩阵，其中第2，3条道路禁止左转 在实际的网络分析中，需要综合利用结点的连通矩阵和结点上道路的连通矩阵。 4．3．3流量和网络分析的具体应用和相关模型 流量和网络分析一般较多地应用于城市交通中，为城市交通管理和道路规划提供科学的决策支持，其中，有如下四类主要的模型： 1）“出行产生（Trip Generation）”模型 根据交通分析区（Traffic Analysis Zone,TAZ）的土地利用形态以及其它社会经济数据来估计各个区域所产生和吸引的交通量，例如，居住区会产生出行，而商业区则吸引出行。 2）“出行分布（Trip Distribution）”模型 根据出行产生分析的结果，确定各个区域之间的交通量。 3）“交通工具选择（Mode Choice）”模型 将交通量计算分配到交通工具上，计算车流量。 4）“出行分派（Trip Assignment）”模型 根据交通线路的容量和速度，将车流量分配到各条道路上。 上述的四种模型，包含多种传统GIS工具所不支持的统计分析运算以及具体的经济地理模型，在具体实现时，可以在GIS中运算，将结果输出到其它系统中，进行专业模型计算。 5．地震灾害和损失估计 对地震灾害以及地震次生灾害的评估对于一个区域的降低危险，资源分配以及紧急响应规划具有重要的意义，而通过存储和分析地质构造信息，利用GIS可以预测地震发生的“场景”并估计该区域由于地震引发的潜在损失。此外，GIS也提供了有力的工具使得在地震实际发生时，分析灾害严重程度的空间分布，帮助政府分配紧急响应资源。 进行地震灾害评估时要综合考虑地质构造等各种信息的空间分布，通常包括以下几个步骤（图14-9）： 1）估计地表震动灾害 需要识别地震源点，然后建立在该点发生地震以及地震波传播的模型，最后根据地表的土壤条件得到最终的震动强度。 2）估计次生的地震灾害 次生的地震灾害包括液化（Liquefaction），滑坡，断裂等，评估这些灾害需要收集相应区域的地质构造信息，计算地表运动的强度和持续时间以及在以前的地震发生过程中这些灾害发生的情况。 3：估计对于建筑物的损害 需要收集地震区域内建筑和生命线的分布状况，然后对每种建筑的建立损害模型，该模型是一个函数，与地表震动强度以及潜在的次要灾害有关。 4：估计可以用金钱衡量和不可以用金钱衡量的损失 可以用金钱衡量的损失包括受损建筑的修复和重建，而不可以用金钱衡量的损失包括人员伤亡。估算这些损失需要相应的社会经济信息，此外，清除垃圾和重新安置费用、失业、精神影响以及其它的长期或短期的影响，需要建立不同的模型，分别加以确定。  图14-9：地震损失评估过程 如上图所示，在地震损失评估中，用到了多种空间信息，如地质构造，建筑等等，因此GIS是非常理想的进行地震损失评估的工具，图14-10，13-11描述了采用GIS进行评估的过程。  图14-10：基于GIS的地震损失评估过程（1）  图14-11：基于GIS的地震损失评估过程（2） 通常，地表震动强度可以根据震源位置以及地震波传播公式计算；而次生灾害以及建筑物的损害要根据相关的图件进行计算，并基于上述计算的结果来评估金钱损失和非金钱损失。在分析过程中，由于地震强度以及破坏程度随着到震源的距离增大而衰减，所以要采用缓冲区计算模型；而在计算金钱损失以及非金钱损失时，因为要综合考虑多个因素，要使用叠加复合模型。 6．地貌 地貌学理论发展和生产实践需要加强计量地貌研究。然而，由于地貌现象的复杂性、地貌数据的庞大等多方面的原因，需要在地貌研究中采用GIS工具，使其成为地貌定量研究的一个有效途径。 在地貌研究中采用GIS技术具体实现步骤和内容归纳为： 1）研究地貌信息内涵和地貌系统的特点。 2）根据地貌信息及其分析方法的特点，结合GIS工具建立地貌信息专题分析系统。 3）在地貌信息系统的支持下，建立综合的定量地貌分析模型，利用尽可能多的信息源，分析地貌信息流及其相关物质流、能量流规律，研究地貌形态、物质组成、成因机制、分布特征、发生发展规律及其对人类生产、生活环境的影响作用。 4）开发机制地貌制图系统作为地貌信息系统的输出系统，将以上的分析结果编制成各种平面或立体的地貌分析图件，实现地貌制图自动化，将在资源开发、环境整治、生产建设的实际目标中发挥更大的作用。 6．1地貌信息的内涵、地貌系统的特点与分析方法 地貌是地壳表层内外营力共同作用到某一发展阶段所形成的三维地表形态，是内、外营力，介质性质和时间的函数，可表述如下：  其中M表示地貌形态；F表示内外营力对地表的作用；x，y，z为空间坐标；Fi表示内营力；Fe表示外营力；m表示构成地貌的介质性质，包括岩性与构造两个方面；t表示作用时间。 内外营力的类型、强度、作用方式、介质性质与作用时间的不同组合，直接影响着地表形态特征。也就是说，在地貌信息中蕴含着归属于不同内外营力和不同发展阶段的特征地貌信息，以此为信息源，提取所需各种特征参数进行复合分析，可以得到不同专题的综合评价结果。 地貌系统是一个开放性的动态系统，是由各种地貌形态类型和要素构成的复杂综合体，它不仅在内部各子系统之间进行着物质、能量和信息流的迁移和转化，而且在系统与外部环境之间也进行着物质、能量和信息流的交换，从而形成一个复杂的物质、能量和信息的传递网络。地貌系统总处于输入——转换——输出的动态过程中，即经常有能量自源区经系统流向耗散区，使系统出现导致有序化的熵减过程，从而具有一定的结构性和相应的熵位，存贮容量和自调节能力。地貌的发展、发生和演化，实质上就是在地貌系统外部环境能量流内、外营力作用下，构成地貌要素的地表物质经历的变形变位（抬升、断错等）分离组合（侵蚀、搬运、堆积等）的复杂物质流传输过程。由于地貌物质传输过程与热力学传导有相似性，许多地貌演化模型，如斜坡演化模型都是热传导方程的形式：  其中H=H(x,y,t)表示地表高程。 由于实际区域地貌系统的过于复杂与庞大，计量研究中往往不能取得数学意义上的精确解释。如以上的演化模型往往是针对某一外力过程，限于极为简单、理想的初始条件并且不考虑介质空间差异的情形；而实际区域内外营力多样化，介质不均一，甚至连初始高程H0都无法表达成确定的函数式。因此只能采用适当的方式离散化，化无限为有限，求得尽量逼近实况的近似解。 具体地说，就是引入有限元素法与有限差分法的思想，将整个区域地貌体分解成许多“微小”的单元（姑且称之为地貌体，其顶部地表面称之为地貌面），如图14-12所示。在小的区域空间内以平代曲，由单个到总体，最终建立整个区域上的种种普适性更强的分析模型以求算各种空间地貌参数和研究复杂地貌过程等问题，可以说离散化是区域地貌信息系统定量分析的前提与特色。  图14-12：区域地貌离散分析 地貌是构造、过程和时间的函数，现存地貌形态的定量指标及其物质组成等特点集中反映了地貌系统以往内外营力共同作用的总和。地貌信息系统分析方法就是提取地貌系统各不同侧面的特征参数，采用系统综合分析方法取得对区域地貌发展规律的总体认识。地理信息系统方法的发展，使得建立地貌信息系统，对庞大复杂的地貌空间数据进行定量分析处理和信息复合与综合评价成为可能。 6．2地貌信息系统 地貌信息系统基于GIS工具开发，集成了GIS的通用功能，并包括如下地貌分析模块： 1）地貌形态分析 用于提取由地表形态反映出来的地貌特征参数，包括DEM以及基于DEM的各种地貌参数：坡度、坡向、起伏度、河网密度、沟谷参数等等。 2）地貌信息复合 包括地貌参数复合模型和地貌信息/专题信息复合模型，前者将各个地貌参数复合，发现其间的规律性，得到综合的分析结果；后者考虑到地貌格局还受到外动力地质作用的影响，将地貌信息与地质、水文、气象等信息复合，发现其联系。 3）地貌区域评价 根据单个或者多个地貌特征的分析，评价地貌环境对农业、建筑、工程、旅游、居住条件的影响以及评价地貌条件对于各种自然灾害发生、发展以及防治的影响。 4）地貌成因分析 包括： 内动力成因分析模型：内动力是地壳内部能量产生的营力，决定地貌的格局与骨架，可以从几何效应和力学机制两个方面进行分析，有限元方法是研究地貌内动力成因机制的有效方法。 外动力成因分析模型：外动力是地球表面受太阳能和重力以及生物活动的营力，主要包括重力作用、风力作用、水利作用、温差作用和生物作用等，影响着地貌的具体形态。 5）地貌预测 可以采用两种模型进行地貌预测，即理论地貌预测模型和统计地貌预测模型，前者以成因分析为基础，采用数理模型进行宏观和微观的模拟，预测地貌的发展；后者在地貌发展历史数据的基础上，采用“黑箱”或“灰箱”的方法进行统计分析，预测地貌发展趋势。 7．医疗卫生 无疑，健康问题是人类生活中最为关注的问题之一，考虑到以下因素，在医疗卫生领域应用GIS，将有助于解决健康问题以及制订相应的政策。 1）某些疾病只是限于一些地区（地方病），探讨其局部的自然和社会条件，对于研究疾病的原因以及治疗手段有着重要的意义。 2）医疗设施（医院、诊所、急救中心等）具有空间分布特性，在一个区域内合理分布医疗设施有利于资源的有效利用。 7．1 GIS与流行病研究 由于流行病是用于描述和解释某种疾病的发病率，从空间的角度来看，流行病学需要很好地描述流行病发病率空间分布特征的手段，进而可以研究发病率模型，以发现流行病和周围环境的关系。通常，GIS在流行病研究中主要提供了如下三个方面的功能： 7．1．1流行病数据的可视化 假定要使用一组点来表示某个人群的发病率，而通常无法直接获得点数据，而只能获得整个区域——如人口普查区或其它行政区划单元——的数据，这些数据包括了疾病的发病率及各疾病的年龄结构。 由于各个空间区域的大小各异，使得表现发病率的空间分布时不直观，一个解决的办法是采用比较统计地图（Cartogram），在比较统计地图上，各个单元的大小与具有发病危险的人口数目成比例，然后将发病率绘制在该图上，这样可以利用GIS，清楚地将流行病数据的空间分布可视化（图14-13）。  图14-13：宾西法尼亚州1981-1988年爱滋病的扩散情况， 图上灰色的区域成几何级数增长生动地表现了爱滋病的扩散 7．1．2空间数据分析 实际上，空间数据可视化与空间数据分析的界限是十分模糊的，如制作比较统计地图需要对地理空间进行变换。由于各个区域人口的分布是不均匀的，为了分析流行病病例的空间分布，通常采用的一种方法称为密度估算（Density Estimation）。该方法是采用一个移动窗口（Moving Window）覆盖于栅格化的位置点上，计算每个窗口内的密度，采用这种方法，关键是确定窗口的大小。此外，计算发病率密度与人口密度的关系；以及对病例空间分布进行聚类分析（Cluster Analysis），并探讨每一特定聚集区域具体特征，都可以应用于流行病数据的空间分析。 7．1．3流行病模型 在识别出每个流行病的空间分布聚集区之后，一个重点的研究是确定流行病和其周围环境的关系，这将有助于验证对流行病发病原因的假定。做法之一是将空气、水质的空间分布与流行病的空间分布进行对比研究，例如在非洲疟疾传染的研究中，使用利用遥感图像计算的NDVI指数来进行建模。此外为了得到流行病的源地，采用距离分析，通常是根据预先假定的病源地进行缓冲区计算，然后得到每个距离范围内的病例数目，进而确定真正的病源地。 利用空间分析手段确定疾病原因，最著名的案例是英国医生琼·斯诺利用地图发现霍乱病病源： 1854年8月到9月英国伦敦霍乱病流行时，当局始终找不到发病原因，后来医生琼·斯诺博士在绘有霍乱流行地区所有道路、房屋、饮用水机井等内容的1:6500城区地图上（图14-14），标出了每个霍乱病死者的住家位置，得到了霍乱病死者居住位置分布图，他分析了这张图，马上明白了霍乱病源之所在——死者住家都集中于饮用“布洛多斯托”井水的地区及周围。根据斯诺博士的分析和请求，当局于9月8日摘下了这个水井的水泵，再往后就没有出现新的霍乱病人了。  图14-14：1954年英国霍乱病流行地区霍乱病死者居住位置分布图 7．2 GIS与医疗设施分布 7．2．1医疗设施规划 不论是在发达国家还是发展中国家，随着“规划需从当地的实际情况出发”的观点的不断被认可，一种以医疗中心为主的观点也不断被采纳。这意味这不再为了所谓声望而投资大型医疗基地，而代之以建造更多面向小社区的诊所，以方便大众。这时，显然需要一种办法来提供人口统计和发病率的详细资料医疗设施需求。在进行这种医疗设施定位时，GIS可以发挥的重要作用。图14-15显示了在英国west Sussex的研究中，通过确定病人到医生的流量，从而对该地区医疗设施规划的合理性做出判断。  图14-15：病人到全科医师的主要流动方向 反应了每个全科医师的“吸引范围”（west Sussex） 在定义医疗需求区域后，为了评估当地居民的需求，通常需要将已获得的调查统计数据进行社会经济分类，以便更好地描述那些小区域。这样，拥有不同生活方式的不同阶层，或不同生活方式与发病率和死亡率之间联系的数据就可获得。然后将标准化发病率通过这一阶层的情况来推及个人，从而产生综合疾病指数应用于实践。这个综合疾病指数可用于预测病人药品的开销等方面。 7．2．2可达性，可用性及结果 GIS在理解医疗结果与医疗设施可达性的关系方面很有帮助。对于一些疾病，如交通事故，迅速送往医院对于病人的抢救至关重要，这需要医院具有良好的可达性，对于可达性与医疗关系的探讨，同样有助于规划医疗设施。 7．3联系流行病学与医疗实施规划—空间决策支持系统 在医学地理学中，地理流行病学和医疗设施规划是两个最主要的领域。建立两者间的桥梁，需要包括空间决策支持系统，如果分析显示在特殊的局部区域存在健康问题，此时，空间决策支持系统将提供工具及策略来处理这一问题。在研究中，发现以下问题是值得注意的： 1）地图的比例尺，小比例尺的地图会掩盖许多细部的重要问题，如小区域间发病率的变化。大比例尺地图则可使有关当局集中力量于最需要他们的区域。 2）疾病的传播是无国界的，随着人类活动的扩展，疾病也随之传播，通过对人口流制图，可以研究疾病的传播规律。 8．军事 军事是以准备和实施战争为中心的社会活动。一切军事行动都是在一定的地理环境中进行的，地理环境对军事行动有着极其重要的影响与作用。随着人类社会向信息化迅速发展，未来高技术战争中信息对抗的含量将越来越高，特别是高技术条件下的局部战争，由于战争爆发突然，战争进程加快、战机稍纵即逝等特点，对作战指挥的时效性有了更高的要求。指挥决策智能化、作战指挥自动化、武器装备信息化成为未来战争取胜的关键。在这种需求下，出现了数字化战场，数字化的地理环境信息已成为指挥决策的必要条件之一。因此，作为空间军事信息保障的军事地理信息系统已成为现代化军事斗争的一项重要内容。 军事地理信息系统MGIS（Military Geographic Information System）是地理信息系统技术在军事方面的应用，是指在计算机软硬件的支持下，对军事地形、资源与环境等空间信息进行采集、存储、检索、分析、显示和输出的技术系统。它在军事地理信息保障和指挥决策中起着重要的作用。 军事地理信息系统和遥感、全球定位系统关系密切，同时和指挥自动化系统 C3I（Command，指挥；Control，控制；Communication，通信；Information，情报）紧密地联系在一起，形成一个多功能的统一系统。它一般由六个子系统组成：信息收集子系统、信息传递子系统、信息处理子系统、信息显示子系统、决策监控子系统和执行子系统。其中，情报是军事决策的基础；信息收集、处理和显示是系统的核心；通信和控制是信息传输和决策过程的保证；指挥使军事决策具体执行。地理信息系统技术在情报的收集、处理、显示和指挥决策方面发挥着重要的作用。 另外，由军事技术革命引发的数字化战场建设已成为未来战场发展的主流，建设数字化战场和数字化部队已成为21世纪军队发展的大趋势，引起了各国的普遍关注。美国著名未来学家托夫勒指出，建设数字化战场是一项比研制原子弹的“曼哈顿工程”更具挑战性的系统工程，“数字化战场是打赢信息战的关键”。战场数字化就其内容来讲，主要是战场地理环境的数字化、作战部队的数字化、各种武器的数字化和士兵装备的数字化。从某种意义上来讲，战场地理环境的数字化是其他数字化的基础，它为作战部队和各种武器装备的数字化提供了必需的战场背景环境和空间定位基础。 8．1国外军事地理信息系统（MGIS）现状 军事地理信息系统在海湾战争及以后的战争中发挥了重要的作用，受到了各国军方的普遍重视。世界上大部分国家都建立了用途不同、规模大小不等的军事地理信息系统。报道较多的是美国、俄罗斯、英国、澳大利亚等国的MGIS。其应用领域包括： 1）基础地理信息 包括地形图、DEM、DTM等。 2）航海、航空管理 航海图、制定计划航线、障碍物、禁区、助航设施、导航管理、空中交通控制等。 3）地形分析 包括战场模拟、行军路线、应急线路分析、越野机动、涉水分析、通视点分析、距离量测、面积量测、武器打击轨迹分析等。 4）任务规划（战略层次） 包括军事基地规划、军事基础设施管理、打击效果评估、巡航导弹支持、战区规划、入侵应急规划、目标分折、轨道建模等。 5）战争管理（战术层次） 包括战场监测、战场管理、小战区规划、登陆计划、战术模拟、后勤保障规划、交通规划等。 6）基础作业支持 包括拦截应用、环境应用、军事设施分类规划等。 7）边界控制 包括边界巡逻和交叉分析、毒品禁运、移民控制等。 8）情报 包括反毒品活动、反恐怖主义活动、武器监视与跟踪、情报收集等。 8．2国外军事应用系统示例 美国国防制图局（DMA，Defense Mapping Agency）是美国军事GIS的管理部门，1989年以来投入大量的人力和物力来开发军事GIS系统。例如： 1）数字航海图 数字化海图（Digital Navigation Chart，DNC）和GPS结合用于海军的导航，来取代以前使用的4500张纸质海图。这些数字地图可以支持广泛的空间分析和地理查询。数字化海图一般包含有12个特征数据库：地名注记、水文要素、界限、港口设施、水下地形、航道、航海标志、救援标志、海洋环境、陆地地形地物、航海障碍物、数据的优先级。 2）数字地图 VMAP（Vector smart MAP）数据分两个层次，分别对应于1：25万比例尺和1：5万比例尺的地形图。该系统包括有10个数据层：交通、网络、界限、水文、城市、公用设施、植被、等高线、工业设施。 3）水文资源评估系统 用于水文和测深数据的收集。 实际上，由于军事保密的特殊性，我们对国外许多功能强大的MGIS了解得都很少。作为军事情报的一个设施，美国正在建立和完善军事应用的全球空间信息和服务（Global Geospatial Information and Services，GGI&S）系统。其中的数据包含两个层次，国家性的和全球性的，符合NOTA DIGEST、ISO TC/211以及OGIS标准。 另外，英国国家遥感中心（British National Remote Sensing Center， NRSC）应用法国的SPOT图像与MGIS结合，模拟敌方的三维地形，对军方飞行员进行模拟训练，取得了较好的效果。 海湾战争的实践证明，利用MGIS技术和遥感技术相结合，可为战时提供实时的地理信息保障。美国国防制图局（DMA）在海湾战争的战场上的MGIS实时服务，主要包括利用自动影像匹配和自动目标识别技术处理卫星和高低空侦察机实时获得的数字影像，及时地为军事决策提供24小时的实时服务。对于冲突地区，最基本的地理信息保障是为作战部队提供作战地图。与传统的地图制图技术相比，MGIS的制图功能无疑是出色的，它既可以提供数字地图，也可以输出精美的印刷地图，其生产速度比传统的制图方法要快得多，特别是在遥感技术支持下，可实时地获得制图数据。在海湾战争期间，美国国防制图局进入一天24小时生产状态，共开发了12000套新的地图产品，其中600套数字地图，印刷了100万幅战地地图，这些地图覆盖科威特、沙特阿拉伯、伊拉克、叙利亚等国和地区，比例尺为1：5万不等。由于对这些战地地图的迫切需要，多国部队在这次战争中使用C130远程运输机运送军事地图到沙特阿拉伯，优先于药品的运输，仅次于爱国者导弹发射架部件的运输，由此可见，战时的实时地理信息保障是何等重要。 除上述地图产品外，美国国防制图局和工程地形实验室还使用了数字地形高程数据、栅格图形数据、地名字典等。在没有现成地图资料的地区，利用MGIS处理SPOT和LANDSAT影像及机载雷达影像，得到可直接应用的影像地图。因此，MGIS和RS相结合，可较好地完成战时的实时地理信息保障任务。 8．3 MGIS的发展方向 概括而言，军事地理信息系统有以下的发展趋势，以便更能够满足现代战争的要求： 1）基础数据、基础设施大型化。 2）各领域应用专业化。 3）实用化，简单化，易学易用。 4）小型化。 5）标准化。