首页 > 空间数据查询

空间数据查询

发布时间：2018-07-02 06:54:19 来源：文档文库

小中大

字号：

手机查看

空间数据查询按查询方式来分，主要分为两类：根据属性条件查询对象和根据几何条件查询对象。无论是根据属性还是根据几何来查询对象，都必须设置相应的查询条件。在ArcGIS Engine中，FeatureLayer对象和IFeatureClass对象有一个Search()方法专门用于空间数据查询，同时ArcGIS Engine中提供了一个IQueryFilter接口来设置查询的条件（属性条件和几何条件）。首先来看一下IFeatureLayer接口的Search()方法，其定义如下：
public IFeatureCursor Search (IQueryFilter queryFilter, bool recycling);
该方法有两个参数：第一个参数queryFilter为一个IQueryFilter对象，用于指定查询的条件了；第二个参数recycling为一个Bool数据类型，表示查询结果游标是否循环。该方法的返回值为一个IFeatureCursor对象，该游标对象表示查询的结果对象。下面就分别根据属性条件查询对象和根据几何条件查询对象来介绍空间数据查询的过程：
一、根据属性条件查询对象
首先新建一个项目，设置项目名称为“空间数据查询”，并在窗体上添加如下图所示控件（axMapControl1、axTOCControl1、axLicenseControl1、axToolbarControl1）：
同时为了便于设置条件来查询对象，添加如下图所示控件（layerComboBox、queryFiltertextBox和根据属性条件查询对象）：

其中，layerComboBox控件用于表示要查询的图层，queryFiltertextBox控件表示查询的条件。
1、定义刷新图层函数
在每次添加新数据或打开新地图时，都要更新layerComboBox控件中的显示项来表示当前地图中的图层。定义刷新图层函数RefreshLayer()，其代码如下：

private void RefreshLayer()
{
layerComboBox.Items.Clear();

for (int i = 0; i < axMapControl1.LayerCount; i++)
{
layerComboBox.Items.Add(axMapControl1.get_Layer(i).Name);
}
layerComboBox.Text = layerComboBox.Items[0].ToString();
}
2、根据属性条件查询对象
根据属性条件查询对象主要通过“根据属性条件查询对象”按钮的Click()事件来实现，下面介绍的代码都在“根据属性条件查询对象”按钮的Click()事件中编写。
1）获取查询的图层
首先要获取到查询的图层对象，如下定义：

ESRI.ArcGIS.Carto.IFeatureLayer pFeatureLayer=axMapControl1.get_Layer(layerComboBox.SelectedIndex) as ESRI.ArcGIS.Carto.IFeatureLayer;

if (pFeatureLayer == null)

{

MessageBox.Show("选择图层不是Feature图层！");

return;

}
如上代码中，定义一个IFeatureLayer对象，该图层对象即为layerComboBox控件中选择的图层，如果该图层不是FeatureLayer，则返回。
2）定义查询条件
首先定义一个查询过滤对象，如下代码所示（由于IQueryFilter接口位于Geodatabase库中，所以还需要添加引用Geodatabase库）：

ESRI.ArcGIS.Geodatabase.IQueryFilter queryFilter = new ESRI.ArcGIS.Geodatabase.QueryFilterClass();

从queryFiltertextBox控件中获取查询条件,如下:

queryFilter.WhereClause = queryFiltertextBox.Text;
3）查询对象
定义一个游标对象IFeatureCursor表示查询到的结果，使用IFeatureLayer的Search()方法来查询，如下代码：

ESRI.ArcGIS.Geodatabase.IFeatureCursor featureCursor = pFeatureLayer.Search(queryFilter, false);

4）高亮显示数据
高亮显示数据可以使用axMapControl1的FlashShape()方法来实现，如下代码：
         ESRI.ArcGIS.Geodatabase.IFeature pFeature ;
         while((pFeature= featureCursor.NextFeature())!=null)
         {
            axMapControl1.FlashShape(pFeature.Shape);
         }
添加一些相关错误处理方式，最后“根据属性条件查询对象”按钮的Click()事件代码如下：
      private void 根据属性条件查询对象_Click(object sender, EventArgs e)
      {
         ESRI.ArcGIS.Carto.IFeatureLayer pFeatureLayer=axMapControl1.get_Layer(layerComboBox.SelectedIndex) as ESRI.ArcGIS.Carto.IFeatureLayer;
         if (pFeatureLayer == null)
         {
            MessageBox.Show("选择图层不是Feature图层！");
            return;
         }
         ESRI.ArcGIS.Geodatabase.IQueryFilter queryFilter = new ESRI.ArcGIS.Geodatabase.QueryFilterClass();
         queryFilter.WhereClause = queryFiltertextBox.Text;

         try
         {
            ESRI.ArcGIS.Geodatabase.IFeatureCursor featureCursor = pFeatureLayer.Search(queryFilter, false);

            ESRI.ArcGIS.Geodatabase.IFeature pFeature;
            while ((pFeature = featureCursor.NextFeature()) != null)
            {
                  axMapControl1.FlashShape(pFeature.Shape);
            }
         }
         catch (Exception pException)
         {
            MessageBox.Show(pException.Message);
         }
      }
值得注意的是：queryFiltertextBox控件中输入的查询条件需符号要求，如查询属性字段Elevation为800的数据对象，如果属性字段Elevation的数据类型为数值型（浮点型、整型、双精度），则直接输入“Elevation = 800”即可。如果属性字段Elevation的数据类型为Text，则输入“Elevation = ‘800’”，需要在值左右侧加上单引号。
运行程序，其结果如下图所示：

二、根据几何条件查询对象
根据几何条件查询对象同样主要是通过FeatureLayer对象或IFeatureClass对象的Search()方法来实现，在Search()方法的第一个参数为IQueryFilter接口，可以将其改为ISpatialFilter接口，通过ISpatialFilter接口来指定空间查询的几何范围，首先来看一下ISpatialFilter接口的一些常用属性和方法：

它本身就继承IQueryFilter接口，所以具有IQueryFilter接口的一些属性和方法，在ISpatialFilter接口的这些属性和方法中只需要注重两个属性：Geometry和SpatialRel。Geometry属性用于指定查询的范围，如果是点击查询则为一个IPoint对象，如果为面范围查询则为一个IPolygon对象。SpatialRel属性用于指定查询数据与查询范围的空间关系，如相交、覆盖、接触等关系，都是通过设置SpatialRel属性来完成。
首先添加如下图所示控件（pointRadioButton、polygonRadioButton、根据几何条件查询对象）：

由于点击或面范围查询，都需要绘制图形（点或面），这些都要在AxMapControl控件中完成，如点击查询可以在AxMapControl控件的OnMouseDown事件中完成，而绘制面则要在AxMapControl控件的OnMouseDown事件中绘制面，在OnDoubleClick()事件中完成绘制同时执行查询。当然这些都不是绝对的，如果已有几何图形范围则不需要去绘制几何图形了。
1、点击查询
1）全部变量定义
首先需要在全局中定义变量DoQueryIndex来表示当前要执行的查询模式：点击查询还是面范围查询。如果值为1表示点击查询，如果值为2表示面范围查询，默认其值为0，不执行任何操作。该变量定义如下：
      private int DoQueryIndex = 0;
2）设置查询模式
如果选择点击模式则设置变量DoQueryIndex值为1，如果选择面范围模式则设置变量DoQueryIndex值为2，如果不查询则设置DoQueryIndex值为0。“根据几何条件查询对象”按钮的Click()事件代码如下：
      private void 根据几何条件查询对象_Click(object sender, EventArgs e)
      {
         if (根据几何条件查询对象.Text == "根据几何条件查询对象")
         {
            if (pointRadioButton.Checked)
            {
                  DoQueryIndex = 1;
            }
            else if (polygonRadioButton.Checked)
            {
                  DoQueryIndex = 2;
            }
            根据几何条件查询对象.Text = "停止几何条件查询对象";
         }
         else
         {
            根据几何条件查询对象.Text = "根据几何条件查询对象";
            DoQueryIndex = 0;
         }
      }
3）点击查询
点击范围可以通过AxMapControl的OnMouseDown()事件来执行，在AxMapControl中按下鼠标即可获得点击位置，通过查询该位置的数据对象即可，如下代码：
      private void axMapControl1_OnMouseDown(object sender, ESRI.ArcGIS.Controls.IMapControlEvents2_OnMouseDownEvent e)
      {
         if (DoQueryIndex == 1)//点击查询
         {
            ESRI.ArcGIS.Carto.IFeatureLayer pFeatureLayer = axMapControl1.get_Layer(layerComboBox.SelectedIndex) as ESRI.ArcGIS.Carto.IFeatureLayer;
            ESRI.ArcGIS.Geometry.IPoint point = new ESRI.ArcGIS.Geometry.PointClass();
            point.PutCoords(e.mapX, e.mapY);

            ESRI.ArcGIS.Geodatabase.ISpatialFilter spatialFilter = new ESRI.ArcGIS.Geodatabase.SpatialFilterClass();
            spatialFilter.Geometry = point;
            spatialFilter.SpatialRel = ESRI.ArcGIS.Geodatabase.esriSpatialRelEnum.esriSpatialRelIntersects;
            ESRI.ArcGIS.Geodatabase.IFeatureCursor featureCursor = pFeatureLayer.Search(spatialFilter, false);

            ESRI.ArcGIS.Geodatabase.IFeature pFeature;
            while ((pFeature = featureCursor.NextFeature()) != null)
            {
                  axMapControl1.FlashShape(pFeature.Shape);
            }
         }
      }
其中，如果为点击模式（即DoQueryIndex = 1）则执行点击查询的代码。首先仍然是点击查询的图层对象或FeatureClass，在此使用IFeatureLayer，然后定义一个IPoint对象，通过PutCoords方法设置其值为地图上点击的位置。然后定义一个ISpatialFilter对象用于设置空间查询的条件，设置其查询的几何范围为前面定义的IPoint对象，同时设置几何对象的几何关系，最后通过IFeatureLayer的Search()方法来执行查询，查询返回结果为一个IFeatureCursor游标对象，通过该游标的NextFeature()方法可以获取游标中的每一个IFeature对象，这些IFeature就是要查询的结果。在此只是闪烁显示了这些对象，如果要获取该对象的属性等信息，直接使用IFeature的Fields属性即可获得。
运行程序，其结果如下图所示：

2、面范围查询
面范围查询首先还是要绘制面的几何图形，然后通过IFeatureClass或IFeatureLayer的Search()方法来执行查询，其方式和点击查询完全相同，只是绘制几何图形的方式不同。如果采用面范围查询方式，在AxMapControl控件的OnMouseDown()事件中记录鼠标点击位置，同时在OnMouseDown()事件中绘制这样一个面区域。
1）全局变量定义
在全局中定义变量pointCollection，其为IPointCollection对象，它用于保存每次在AxMapControl控件的OnMouseDown()事件中点击的鼠标位置，如下：
      private ESRI.ArcGIS.Geometry.IPointCollection pointCollection;
2）定义绘制面函数
由于在点击鼠标的同时，还要显示已经绘制的面的范围以便用户更好的选择查询区域，所以需要定义一个绘制面的函数DrawPolygon()，其代码如下：
      private void DrawPolygon(ESRI.ArcGIS.Geometry.IPointCollection pPointCollection, ESRI.ArcGIS.Controls.AxMapControl axMapControl)
      {
         ESRI.ArcGIS.Geometry.IPolygon pPolygon;
         pPolygon = (ESRI.ArcGIS.Geometry.IPolygon)pPointCollection;
         axMapControl.DrawShape(pPolygon);
      }
3）记录面的边界点
由于每次在AxMapControl上点击鼠标表示该为面创建一个边界点，所以在OnMouseDown()事件中要添加该点到变量pointCollection中去，同时绘制已经形成的面。故在axMapControl1的OnMouseDown()事件中继续添加如下代码：
         else if(DoQueryIndex == 2)//面范围查询
         {
            ESRI.ArcGIS.Geometry.IPoint point = new ESRI.ArcGIS.Geometry.PointClass();
            point.PutCoords(e.mapX, e.mapY);

            pointCollection.AddPoints(1, ref point);

            if (pointCollection.PointCount > 2)
            {
                  DrawPolygon(pointCollection, axMapControl1);
            }
         }

4）查询范围内对象
设计当鼠标在axMapControl1控件上双击时即完成面范围的绘制，所以执行面范围的查询就在axMapControl1的OnDoubleClick()事件中，如下代码所示：
      private void axMapControl1_OnDoubleClick(object sender, ESRI.ArcGIS.Controls.IMapControlEvents2_OnDoubleClickEvent e)
      {
         if(DoQueryIndex == 2)//面范围查询
         {
            ESRI.ArcGIS.Geometry.IPoint point = new ESRI.ArcGIS.Geometry.PointClass();
            point.PutCoords(e.mapX, e.mapY);
            pointCollection.AddPoints(1, ref point);
            if (pointCollection.PointCount > 2)
            {
                  axMapControl1.Refresh();
                  ESRI.ArcGIS.Carto.IFeatureLayer pFeatureLayer = axMapControl1.get_Layer(layerComboBox.SelectedIndex) as ESRI.ArcGIS.Carto.IFeatureLayer;
                  ESRI.ArcGIS.Geometry.IPolygon pPolygon;
                  pPolygon = (ESRI.ArcGIS.Geometry.IPolygon)pointCollection;

                  ESRI.ArcGIS.Geodatabase.ISpatialFilter spatialFilter = new ESRI.ArcGIS.Geodatabase.SpatialFilterClass();
                  spatialFilter.Geometry = pPolygon;
                  spatialFilter.SpatialRel = ESRI.ArcGIS.Geodatabase.esriSpatialRelEnum.esriSpatialRelIntersects;
                  ESRI.ArcGIS.Geodatabase.IFeatureCursor featureCursor = pFeatureLayer.Search(spatialFilter, false);

                  ESRI.ArcGIS.Geodatabase.IFeature pFeature;
                  int featureCount = 0;
                  while ((pFeature = featureCursor.NextFeature()) != null)
                  {
                     axMapControl1.FlashShape(pFeature.Shape);
                     featureCount++;
                  }
                  MessageBox.Show("选择范围内共有" + featureCount.ToString() + "个对象!");
            }
         }
      }
调用几何条件查询在“根据几何条件查询对象”按钮的Click()事件中完成，所以在调用面范围查询时还需要对全局变量pointCollection进行实例化，所以“根据几何条件查询对象”按钮Click()事件的代码如下：
      private void 根据几何条件查询对象_Click(object sender, EventArgs e)
      {
         if (根据几何条件查询对象.Text == "根据几何条件查询对象")
         {
            if (pointRadioButton.Checked)//点击查询
            {
                  DoQueryIndex = 1;
            }
            else if (polygonRadioButton.Checked)//面范围查询
            {
                  DoQueryIndex = 2;
                  pointCollection = new ESRI.ArcGIS.Geometry.PolygonClass();
            }
            根据几何条件查询对象.Text = "停止几何条件查询对象";
         }
         else
         {
            根据几何条件查询对象.Text = "根据几何条件查询对象";
            DoQueryIndex = 0;
         }