操作步骤
 
1. 距离制图－创建缓冲区
 
1.1 点要素图层的缓冲区分析
 
在 ArcMap 中新建地图文档，加载图层：StudyArea ，point 打开 Arctoolbox，执行命令<Spatial Analyst Tools>－<Distance>-<欧几里德距离>，按下图所示设置各参数 在 ArcView GIS 3.x 中 Find Distance 功能与 ArcGIS 中<<欧几里德距离>制图功能相同。
 
 
点击“环境设置”按钮设置“常规选项”中的“输入范围”，使其与StudyArea 相同
 
 
 
1）显示并激活由 point.shp 产生的新栅格主题，eucdist_poin（如上图）。
 
在进行分析时，若选中了 point 图层中的某一个或几个要素，则缓冲区分析只对该要素进行；否则，对整个图层的所有要素进行。
 
1.2 线要素图层的缓冲区分析
 
1）.在 ArcMap 中，新建地图文档，加载 line 图层，点击常用工具栏将地图适当缩小。
 
 
2）.分别选中图层 line 中的两条线，进行缓冲区分析，注意比较线的缓冲区分析与点的缓冲区分析有何不同。 方法：打开 Arctoolbox，执行命令<Spatial Analyst Tools>－<Distance>-<欧几里德距离>；设置“环境设置”中“常规选项”中的“输出范围”为“Same As Display”
 
 
3）取消选定，对整个 line 层面进行缓冲区分析，观察与前两个分析结果的区别（如下图）。
 
 
1.3 多边形图层的缓冲区分析
 
在 ArcMap 中新建地图文档，添加图层：polygon，进行缓冲区分析，观察面的缓冲区分析与点、线的缓冲区分析有何区别。 与创建线的缓冲区相同，先将地图适当缩小，将“环境设置”中“常规选项” 中的“输出范围”为“Same As Display”
 
 
．综合应用实验
 
2.1 水源污染防治
 
Point 图层面表示了水源（如：水井）的位置分布，要求利用缓冲区分析找出水源污染防治的重点区域。步骤：
 
在 ArcMap 中，新建地图文档，添加表示水源分布的点图层数据 point.shp。
打开 ArcToolbox，执行<Spatial Analyst Tools>－<Distance>-<欧几里德距离> 命令。
显示并激活由 point.shp 产生的新栅格图层 eucdist_poi1（如下图）。
右键点击栅格图层 eucdist_poi1，在右键菜单中执行“属性”命令，设置图层显示符号，可调整分类和设置图例显示方式。
 
新的栅格图层显示了区域内每个栅格距最近的水井的距离，其中红色的栅格距各个井的距离最近，对水源的影响最大；蓝色的栅格距各个井的距离最远，影响最小。
 
 
在本例中认为距各个水井 0.1 以内的区域对水质的影响和污染最大，因此，打开“空间分析”工具栏（首先要确保“空间分析”扩展模块已经加载），执行菜单命令<空间分析>-<栅格计算器>，将距离≤0.1 以内的区域提取出作为缓冲区进行专项的污染防治。（如下图）。
 
                  
 
 
2.2 受污染地区的分等定级 
 
假定：point 图层表示的是几个点状污染源，距污染源的远近不同，受污染的状况也不同，距污染源越近，受污染越严重，据此对污染源附近地区进行分等定级。
 
在 ArcMap 中，新建地图文档，添加图层：point.、StudyArea。打开 ArcToolbox，执行<Spatial Analyst Tools>-<Distance>-<欧几里德距离>命令，得到新的栅格图层 Dist_of_Point； 
打开“空间分析”工具栏，执行命令：<空间分析>－<栅格计算器>，基于栅格图层 Dist_of_point 进行栅格计算，分别提取（[Dist_of_point]≤0.1）和区域、及（[dist_of_point] >= 0.1 & [dist_of_point] <= 0.15）的区域，分别得到栅格图层：“计算”和“计算 2”；
 
          
 
 
 
 
3.对栅格图层：“计算 2”进行重分类运算（执行命令<空间分析>－<重分类>命令），使得原来的 True（1）值为 0，False（0）值为 1，得到新的栅格图层：Reclass of 计算 2。
 
 
4.将图层“计算”与“Reclass of 计算 2”相加（执行命令：<空间分析>－<栅格计算器>），得到中间结果：计算 3：执行命令<空间分析>－<重分类>对栅格图层“计算 3” 按分等定级的要求进行分类得到结果栅格图层“Reclass of 计算 3”。  (<=0.1 的区域污染级别定为 1，)=0.1 且<=0.15 的区域级别定为 2, >0.15 的区域级别定为 3)
 
 
 
 
 
  2.3 城市化的影响范围 
 
假定：urban 图层表示的是城市化进程中的一些工业小城镇，还包括一个自然生态保护区。这些小城镇的城市化会对周边地区产生一些扩张影响，但自然生态保护区周围 0.05 的范围内不能有污染性的工业，因此其城市化的范围就受到限制。
 
 
(1) 在 ArcMap 中，新建地图文档，添加图层：urban.shp、UrbanArea，对 urban 图层中的自然保护区图斑（属性 Type=”保护区”）执行[<Spatial Analyst Tools>-<Distance>-< 欧几里德距离>]命令，得到Dist_Nature）；对 urban 图层中的除了自然保护区的所有图斑执行[<Spatial Analyst Tools>-<Distance>-<欧几里德距离>]命令，得到Dist_Res
 
 
 
（2）对图层Dist_Nature执行栅格计算（使用空间分析工具中的栅格计算器），提取<=0.05 的区域，并进行重分类，使得原来的 True（1）值为 0，False（0）值为 1，得到“Reclass of 计算”。
 
 
（3）对图层 Dist_Res 进行栅格计算提取<=0.06 的区域，得到“计算 2”；
 
 
（4）将图层“Reclass of 计算”与“计算 2”相乘的栅格计算（使用空间分析工具栏中的栅格计算器），得到城市化范围“计算 3”图层。
 

趁着这次总结的机会，就想着把这几个星期编的小功能都写成博客记录下来，省的之后遗忘，我写这些东西也算是给自己加深印象吧，鱼的记忆哈哈。
 
一、AE缓冲区接口和类
 
还是老生常谈，先来了解实现缓冲区功能都需要那些接口和类：
 GP工具：也就是Geoprocessor类，这其中主要是封装了输入路径、输出路径还有就是可以调用各种工具的函数方法。它给我的感觉类似于原来我是一亩地一亩地的打农药，突然有一天来了一架飞机直接在天上散农药，一下子活就干完了。这GP工具就是那飞机，它大大的提高了效率。
 Buffer类：这里注意了这是AnalysisTools下的Buffer类，我记得有两个buffer类，下图是Buffer类中的参数设置表。
 图有点不清晰，我就只能截一点相对重要的信息。
 
 
二、功能的实现
 
2.1功能实现代码
 
   private IActiveView pActiveView;
   private AxMapControl mapControl;

   //通过构造函数获取MapControl控件
   public BufferAnalysis(AxMapControl mainmapControl)
       : this()
   {
       mapControl = mainmapControl;
       pActiveView = mainmapControl.ActiveView;
   }

   //输出路径
   private void button1_OutputPath_Click(object sender, EventArgs e)
   {
       SaveFileDialog saveDlg = new SaveFileDialog();
       saveDlg.Filter = "Shapefile(*.shp)|*.shp";
       if (saveDlg.ShowDialog() != DialogResult.OK) return;

       textBox2_OutputPath.Text = saveDlg.FileName;
   }

   //生成窗体函数
   private void BufferAnalysis_Load(object sender, EventArgs e)
   {
       if (mapControl == null || pActiveView.FocusMap.LayerCount == 0)
       {
           return;
       }
       IEnumLayer layers = pActiveView.FocusMap.get_Layers();
       layers.Reset();
       ILayer layer = layers.Next();
       while (layer != null)
       {
           comboBox1_ChooseLayer.Items.Add(layer.Name);
           layer = layers.Next();
       }
   }

   //确认分析
   private void button2_OK_Click(object sender, EventArgs e)
   {
       double bufferDistance = Convert.ToDouble(textBox1_BufferDistance.Text.Trim());
       if (bufferDistance == 0.0)
       {
           MessageBox.Show("缓冲区距离有误！");
           return;
       }

       if (comboBox1_ChooseLayer.Text == string.Empty)
       {
           MessageBox.Show("输入图层不能为空！");
           return;
       }

       if (textBox2_OutputPath.Text == string.Empty)
       {
           MessageBox.Show("输出路径不能为空！");
           return;
       }
       int index = comboBox1_ChooseLayer.SelectedIndex;
       string name = getLayerPath(pActiveView.FocusMap.get_Layer(index));
       string outPath = textBox2_OutputPath.Text.Trim();

       Geoprocessor pGp = new Geoprocessor();
       pGp.OverwriteOutput = true; //允许运算结果覆盖现有文件，可无
       ESRI.ArcGIS.AnalysisTools.Buffer pBuffer = new ESRI.ArcGIS.AnalysisTools.Buffer();
       //获取缓冲区分析图层
       ILayer pLayer = pActiveView.FocusMap.get_Layer(index);
       IFeatureLayer featLayer = pLayer as IFeatureLayer;
       //IFeatureCursor cursor = featLayer.Search(null, false);
       //IFeature feaClass = cursor.NextFeature();

       pBuffer.in_features = featLayer;

       pBuffer.out_feature_class = outPath; //输出路径
       pBuffer.buffer_distance_or_field = bufferDistance; //缓冲区参数
       pBuffer.dissolve_option = "NONE"; //融合缓冲区重叠交叉部分，如果不融合填"ALL"
       pGp.Execute(pBuffer, null); //执行

       string pFolder = System.IO.Path.GetDirectoryName(outPath); //得到字符串中文件夹位置
       string pFileName = System.IO.Path.GetFileName(outPath); //得到字符串中文件名字
       mapControl.AddShapeFile(pFolder, pFileName); //往地图控件里添加文件
       mapControl.ActiveView.Refresh(); //激活窗口刷新
       this.Close();
   }

   //获取图层源路径
   private string getLayerPath(ILayer pLayer)
   {
       IDatasetName pDatasetName = (pLayer as IDataLayer2).DataSourceName as IDatasetName;
       IWorkspaceName pWorkspaceName = pDatasetName.WorkspaceName;
       return pWorkspaceName.PathName +"\\" +pLayer.Name+".shp";
   }
 
为了省事，我直接将代码整个复制过来了哈哈。
 
2.2功能实现效果
 
 
三、过程总结
 
GP工具的出现大大的提升了编程的速度，但是这也使得我们变得“肤浅”，不过现在对我这种小白来讲现在的主要矛盾其实还是学习怎么使用GP工具，其他的日后再说吧哈哈。

 
Cesium缓冲区分析
 
前言
效果图
关键代码

 
前言
 
在Cesium中做缓冲区分析其实与在二维地图中做缓存区分析是一样的道理，缓冲后需要将二维坐标转到三维坐标，然后在场景中显示结果。
 
效果图
 
 
关键代码
 
//初始化点缓冲
    initPointBuffer() {
        let point = [106.422638966289, 29.5698367125623];
        this.addPoint(point);

        let pointF = turf.point(point);
        let buffered = turf.buffer(pointF, 60, { units: 'meters' });
        let coordinates = buffered.geometry.coordinates;
        let points = coordinates[0];
        let degreesArray = this.pointsToDegreesArray(points);
        this.addBufferPolyogn(Cesium.Cartesian3.fromDegreesArray(degreesArray));
    },

    //添加点
    addPoint(point) {
        this.viewer.entities.add({
            position: Cesium.Cartesian3.fromDegrees(point[0], point[1], 0),
            point: {
                pixelSize: 10,
                color: Cesium.Color.YELLOW,
                outlineWidth: 3,
                outlineColor: Cesium.Color.YELLOW.withAlpha(0.4),
            }
        });
    },

    //初始化线缓冲
    initPolylineBuffer() {
        let points = [
            [106.425203158107, 29.5694914480581],
            [106.428808047023, 29.569230166027],
            [106.431661917416, 29.5692674920729],
            [106.434708906857, 29.5693048181049]
        ];
        let degreesArray = this.pointsToDegreesArray(points);
        this.addPolyline(Cesium.Cartesian3.fromDegreesArray(degreesArray));

        let polylineF = turf.lineString(points);
        let buffered = turf.buffer(polylineF, 30, { units: 'meters' });
        let coordinates = buffered.geometry.coordinates;
        points = coordinates[0];
        degreesArray = this.pointsToDegreesArray(points);
        this.addBufferPolyogn(Cesium.Cartesian3.fromDegreesArray(degreesArray));
    },

    //添加线
    addPolyline(positions) {
        this.viewer.entities.add({
            polyline: {
                positions: positions,
                width: 2,
                material: Cesium.Color.YELLOW,
            }
        })
    },

    //初始化面缓冲
    initPolygonBuffer() {
        let points = [
            [106.438549830166, 29.5701073244566],
            [106.440695597377, 29.5701073244566],
            [106.440738512722, 29.5688755679036],
            [106.438700033871, 29.5687262630581],
            [106.438034846035, 29.5690248725284],
            [106.438549830166, 29.5701073244566]
        ];

        let degreesArray = this.pointsToDegreesArray(points);
        this.addPolygon(Cesium.Cartesian3.fromDegreesArray(degreesArray));

        let polygonF = turf.polygon([points]);
        let buffered = turf.buffer(polygonF, 60, { units: 'meters' });
        let coordinates = buffered.geometry.coordinates;
        points = coordinates[0];
        degreesArray = this.pointsToDegreesArray(points);
        this.addBufferPolyogn(Cesium.Cartesian3.fromDegreesArray(degreesArray));
    },

    //添加面
    addPolygon(positions) {
        this.viewer.entities.add({
            polygon: {
                hierarchy: new Cesium.PolygonHierarchy(positions),
                material: Cesium.Color.YELLOW.withAlpha(0.6),
                classificationType: Cesium.ClassificationType.BOTH
            },
            polyline: {
                positions: positions,
                width: 2,
                material: Cesium.Color.YELLOW.withAlpha(0.4),
            }
        });
    },

    //添加缓冲面
    addBufferPolyogn(positions) {
        this.viewer.entities.add({
            polygon: {
                hierarchy: new Cesium.PolygonHierarchy(positions),
                material: Cesium.Color.RED.withAlpha(0.6),
                classificationType: Cesium.ClassificationType.BOTH
            },
        });
    },

    //格式转换
    pointsToDegreesArray(points) {
        let degreesArray = [];
        points.map(item => {
            degreesArray.push(item[0]);
            degreesArray.push(item[1]);
        });
        return degreesArray;
    },

 
在前端我们可以借助turf插件进行缓冲区的生成，这样不用依赖后的GIS服务器就能做分析。
 
详情参见 Cesium实战专栏

 
实验平台：ArcGIS 9.3
 
实验目的：熟练掌握A rcGIS缓冲区分析和叠置分析操作，综合利用各项空间分析工具解决实际问题。
 
实验要求：对每个条件进行缓冲区分析，运用空间叠置分析对多个图层叠加，并分等级，确定合适的区域。
 
实验数据：ArcEx8
 
实验步骤
 
打开ArcMap,加载数据ArcEx8，如下所示：
 
 
◆ 选择view\data frame properties,在打开的对话框里将map,display选为meters，如下所示：