去除碎斑

        private bool EliminateSpot(IFeatureClass pFeatureClass,double Min)
        {
            try
            {
                List<int> dl = new List<int>();
                int cout = pFeatureClass.FeatureCount(null);
                IFeatureCursor pcur = pFeatureClass.Update(null, false);
                IFeature pf = pcur.NextFeature();
                while (pf != null)
                {
                    IArea ar = pf.ShapeCopy as IArea;
                    if (ar.Area < Min)
                    {
                        dl.Add(pf.OID);
                        IFeature unfeature = pf;
                        for (int i = 0; i < cout; i++)
                        {
                            if (i == pf.OID)
                                continue;
                            IFeature listfeature = pFeatureClass.GetFeature(i);
                            IRelationalOperator ro = pf.ShapeCopy as IRelationalOperator;
                            ITopologicalOperator lt = listfeature.ShapeCopy as ITopologicalOperator;
                            ITopologicalOperator pt = pf.ShapeCopy as ITopologicalOperator;
                            IRelationalOperator rx = lt.Boundary as IRelationalOperator;
                            if (ro.Overlaps(listfeature.ShapeCopy))
                            {
                                IArea area = listfeature.ShapeCopy as IArea;
                                if (area.Area > (unfeature.ShapeCopy as IArea).Area)
                                    unfeature = listfeature;
                                break;
                            }
                            else if (rx.Contains(pt.Boundary))
                            {
                                IArea area = listfeature.ShapeCopy as IArea;
                                if (area.Area > (unfeature.ShapeCopy as IArea).Area)
                                    unfeature = listfeature;
                                break;
                            }
                            else if (rx.Overlaps(pt.Boundary))
                            {
                                IArea area = listfeature.ShapeCopy as IArea;
                                if (area.Area > (unfeature.ShapeCopy as IArea).Area)
                                    unfeature = listfeature;
                            }
                        }
                        IGeometry bag = new GeometryBagClass();
                        IGeometryCollection gc = bag as IGeometryCollection;
                        object m1 = Type.Missing;
                        gc.AddGeometry(pf.ShapeCopy, ref m1, ref m1);
                        gc.AddGeometry(unfeature.ShapeCopy, ref m1, ref m1);
                        ITopologicalOperator to = new PolygonClass();
                        to.ConstructUnion(bag as IEnumGeometry);
                        unfeature.Shape = to as IGeometry;
                        unfeature.Store();
                    }
                    pf = pcur.NextFeature();
                }

                string sql = dl[0].ToString();
                for (int i = 1; i < dl.Count; i++)
                {
                    sql += "," + dl[i].ToString();
                }
                sql = "(" + sql + ")";

                if (sql != "()")
                {
                    IQueryFilter pQueryFilter = new QueryFilterClass();
                    pQueryFilter.WhereClause = "FID in" + sql;
                    ITable pTable = pFeatureClass as ITable;
                    pTable.DeleteSearchedRows(pQueryFilter);
                    System.Runtime.InteropServices.Marshal.ReleaseComObject(pQueryFilter);
                }
                return true;
            }
            catch
            {
                return false;
            }

        }