ae_将面积小于1500的Feature同附近Feature进行合并

private void 合并1500图斑ToolStripMenuItem_Click(object sender, EventArgs e)
        {
            /*
             *将图层中面积小于1500的图斑与之相同BSM且相邻的最大的图斑合并。 
             */

            if(axMapControl2.LayerCount<=0)
            { MessageBox.Show("请先添加图层再合并！"); }
            else
            {
            pLayer = axMapControl2.get_Layer(0);//拿到层
            pFLayer = pLayer as IFeatureLayer;//要素层
            pFC = pFLayer.FeatureClass;//要素类
            
            IFeatureCursor pFCursor = pFC.Search(null, false);//拿到该要素下的指针

            IFeature pFeature = pFCursor.NextFeature();//递进光标的位置，拿到返回在该位置的要素对象
             


            //创建表，获得要素ID，标识码，面积，用于筛选相同表示码，面积较小的要素
            DataTable pTable = new DataTable();

            DataColumn colID = new DataColumn("要素ID");//新建一个要素列，并设置列名
            colID.DataType = System.Type.GetType("System.String");//设置该列类型
            pTable.Columns.Add(colID);//塞入表

            DataColumn colIdent = new DataColumn("标识码");
            colIdent.DataType = System.Type.GetType("System.String");
            pTable.Columns.Add(colIdent);

            DataColumn colArea = new DataColumn("面积");
            colArea.DataType = System.Type.GetType("System.Double");
            pTable.Columns.Add(colArea);

            int indexOfFID_HBTC_Z = pFC.FindField("FID_HBTC_Z");//返回要素类的指定字段的索引
            int indexOfBSM = pFC.FindField("BSM");//返回要素类的指定字段的索引
            int indexOfTBMJ = pFC.FindField("TBMJ");//返回要素类的指定字段的索引
           
            while (pFeature != null)
            {
                //string FID_HBTC_Z = pFeature.get_Value(indexOfFID_HBTC_Z).ToString();
                string BSM = pFeature.get_Value(indexOfBSM).ToString();
                double TBMJ = (double)pFeature.get_Value(indexOfTBMJ);


                //挑选出小于1500的图斑
                if (TBMJ < 1500)
                {
                    //查找标识码相同的图斑
                    MergeFeature(pFeature);

                    //IQueryFilter queryFilter1 = new QueryFilterClass();
                    //queryFilter1.WhereClause = "TBMJ=" + TBMJ.ToString();
                    //IFeatureCursor pFUpdateCursor = pFC.Update(queryFilter1, false);
                    //pFUpdateCursor.DeleteFeature();
                    
                }
                
                pFeature = pFCursor.NextFeature();
            }

            MessageBox.Show("合并成功");

            }
        }

        private void MergeFeature(IFeature tarFeature)
        {
            int indexOfBSM = pFC.FindField("BSM");
            string BSM = tarFeature.get_Value(indexOfBSM).ToString();//拿到当前feature的BSM值
            IQueryFilter queryFilter = new QueryFilterClass();
            queryFilter.WhereClause = "BSM=" + BSM;

            int tarindexOfTBMJ = pFC.FindField("TBMJ");
            double tarTBMJ = (double)tarFeature.get_Value(tarindexOfTBMJ);//拿到当前feature的TBMJ值

            IFeatureCursor updateCursor = pFC.Update(queryFilter, true);//按照查询，拿到更新要素的指针

            IFeature feature1 = updateCursor.NextFeature();//拿到BSM相同的下一个要素

            //feature1.get_Value(30);

            while (feature1 != null)
            {
                //是否相邻
                ITopologicalOperator topo = (feature1.Shape as IPolygon) as ITopologicalOperator;

                IRelationalOperator rela = (feature1.Shape as IPolygon)  as IRelationalOperator;


                //bool biii = rela.Touches(tarFeature.Shape);
                if (rela.Touches(tarFeature.Shape))//如果相邻
                {
                    int indexOfTBMJ = pFC.FindField("TBMJ");//从要素类中拿到该字段索引值
                    double TBMJ = (double)feature1.get_Value(indexOfTBMJ);//拿到该要素该字段的值

                    //合并图斑
                    if (TBMJ>=1500)
                    {
                        feature1.Shape = topo.Union(tarFeature.Shape) as IPolygon;
                        feature1.Store();
                        tarFeature.Delete();//删除被合并要素
                        break;
                    }
                }
                //ESRI.ArcGIS.SpatialStatisticsTools.CalculateAreas

                //updateCursor.Flush(); //是否。。。
                
                feature1 = updateCursor.NextFeature();

            }

           

            //TokayWorkspace.ComRelease(del_featcur);
            //updateCursor = null;  


            
            axMapControl2.ActiveView.Refresh();
        }