一、什么是叠加分析

1-1、基本概念

叠加分析是地理信息系统最常用的提取空间隐含信息的手段之一。该方法源于传统的透明材料叠加，即将来自不同的数据源的图纸绘于透明纸上，在透光桌上将其叠放在一起，然后用笔勾出感兴趣的部分或提取感兴趣的信息。
2、叠加分析 - 图1

1-2、计算公式

2、叠加分析 - 图2

1-3、应用场景

比如，现在城市规划需要建一个新的大型商场，将规划的新图层于旧城区居民的图层叠加分析，就可以得出那些小区要拆迁。
如规划一条铁路，将规划的图层和农田图层叠加分析，决定需要征用田地的补偿预算款型。

二、步骤

  1. - 输入数据
  2. - 叠加数据
  3. - 显示结果

三、示例

02.png

  1. <!DOCTYPE html>
  2. <html lang="en">
  4. <head>
  5.     <meta charset="UTF-8">
  6.     <meta name="viewport" content="width=device-width, initial-scale=1.0">
  7.     <title>Document</title>
  8.     <script src="./dist/include-openlayers-local.js"></script>
  9.     <script src="./libs/gaode.js"></script>
  10. </head>
  12. <body>
  13.     <button onclick="overAnalysis()">叠加分析</button>
  14.     <div id="map_container">
  15.     </div>
  16.     <script>
  17.          /* 模拟居民区 */
  18.          const vectorData = 'gdbp://MapGISLocalPlus/gis0517/sfcls/Area'
  19.         const vecLayer = new Zondy.Map.GdbpLayer(
  20.             '',
  21.             [vectorData], {
  22.                 ip: 'localhost',
  23.                 port: '6163'
  24.             }
  25.         )
  26.         var map = new ol.Map({
  27.             target: "map_container",
  28.             layers: [gaode,vecLayer],
  29.             view: new ol.View({
  30.                 projection: 'EPSG:4326',
  31.                 center: [114.30, 30.50],
  32.                 zoom: 4
  33.             })
  34.         })
  36.          /* 叠加分析 */
  37.          let resultname;
  38.         function overAnalysis() {
  39.             //2.1设置要叠加的多边形
  40.             //模拟商圈
  41.             const geoRegion = new Zondy.Object.GRegion(
  42.                 [
  43.                     new Zondy.Object.AnyLine([
  44.                         new Zondy.Object.Arc(
  45.                             [
  46.                                 new Zondy.Object.Point2D(100, 33),
  47.                                 new Zondy.Object.Point2D(100, 22),
  48.                                 new Zondy.Object.Point2D(115, 22),
  49.                                 new Zondy.Object.Point2D(115, 33),
  50.                                 new Zondy.Object.Point2D(190, 33)
  51.                             ],
  52.                             0
  53.                         )
  54.                     ])
  55.                 ],
  56.                 0
  57.             )
  58.             resultname = 'gdbp://MapGISLocalPlus/gis0517/sfcls/Area'+Date.now();
  59.             //2.2 创建叠加分析服务
  60.             var overlayParam = new Zondy.Service.OverlayByPolygon({
  61.                 ip: 'localhost',
  62.                 port: '6163', //访问IGServer的端口号，.net版为6163，Java版为8089,
  63.                 //设置被叠加图层URL
  64.                 srcInfo1: vectorData,
  65.                 //设置结果URL
  66.                 desInfo: resultname,
  67.                 //设置多边形坐标序列化对象
  68.                 strGRegionXML: JSON.stringify(geoRegion),
  69.                 //多边形字符串输入格式
  70.                 inFormat: 'json',
  71.                 //设置结果图层的图形参数信息
  72.                 infoOptType:2,
  73.                 //求交
  74.                 overType: 1,
  75.                 //允许重算面积
  76.                 isReCalculate: true,
  77.                 //容差半径
  78.                 radius: 0.05,
  79.             })
  80.             //2.3 执行叠加分析
  81.             //调用基类的execute方法,执行叠加分析,onSuccess为结果回调函数
  82.             overlayParam.execute(
  83.                 onSuccess,
  84.                 'post',
  85.                 false,
  86.                 'json',
  87.                 onError
  88.             )
  89.         }
  90.         //2-4、分析结果显示到地图容器上
  91.         function onSuccess(result){
  92.             console.log(result)
  93.             if(result.succeed){
  94.                 map.removeLayer(vecLayer)
  95.                 console.log(resultname)
  96.                 const resultLayer = new Zondy.Map.GdbpLayer('',[resultname],{
  97.                     ip:'localhost',
  98.                     port:'6163'
  99.                 })
  100.                 map.addLayer(resultLayer)
  101.             }
  102.         }
  103.         function onError(err){
  104.             console.log(err)
  105.         }
  106.     </script>
  107. </body>
  109. </html>