删除重叠的rasterdataset

        public static bool DeleteDuplicatedRasterDataset(ref IRasterCatalog targetRasterCatalog)

        {

            try

            {

                IFeatureClass pFeatureClass = targetRasterCatalog as IFeatureClass;

                int featurecount = pFeatureClass.FeatureCount(null);

                int oidindex = pFeatureClass.FindField(pFeatureClass.OIDFieldName);

                int shpindex = pFeatureClass.FindField(pFeatureClass.ShapeFieldName);

                IList<int> oids = new List<int>();

                IFeatureCursor pFeatureCursor = pFeatureClass.Search(null, false);

                IFeature pFeature = pFeatureCursor.NextFeature();

                while (pFeature != null)

                {

                    oids.Add((int)pFeature.get_Value(oidindex));

                    pFeature = pFeatureCursor.NextFeature();

                }

                IFeature pfirstFeature;

                IFeature potherFeature;

                IRelationalOperator pRelationalOperator;

                ITopologicalOperator pTopologicalOperator;

                IGeometryBag pGeometryBag;

                IGeometryCollection pGeometryCollection;

                IPolygon4 pOutPolygon;

                object before = Type.Missing;

                object after = Type.Missing;

                foreach (int o in oids)

                {

                    if (IsFeatureExisted(pFeatureClass, o))

                    {

                        pfirstFeature = pFeatureClass.GetFeature(o);

                        pRelationalOperator = (IRelationalOperator)pfirstFeature.Shape.Envelope;

                        pGeometryBag = new GeometryBagClass();

                        pGeometryCollection = pGeometryBag as IGeometryCollection;

                        foreach (int i in oids)   //效率太低，可以优化成先查询获取相交的feature.shape

                        {

                            if (IsFeatureExisted(pFeatureClass, i))

                            {

                                potherFeature = pFeatureClass.GetFeature(i);

                                if (pRelationalOperator.Overlaps(potherFeature.Shape.Envelope))

                                {

                                    pGeometryCollection.AddGeometry(potherFeature.Shape.Envelope as IGeometry, ref before, ref after);

                                }

                            }

                        }

                        if (pGeometryCollection.GeometryCount > 0)

                        {

                            pTopologicalOperator = new PolygonClass();

                            pTopologicalOperator.ConstructUnion(pGeometryBag as IEnumGeometry);

                            pOutPolygon = pTopologicalOperator as IPolygon4;

                            if ((pTopologicalOperator as IRelationalOperator).Contains(pfirstFeature.Shape))

                            {

                                pfirstFeature.Delete();

                                UpdataFootprint(targetRasterCatalog);

                            }

                        }

                    }

                }

                return true;

            }

            catch (Exception ex)

            {

                MessageBox.Show(ex.ToString());

                return false;

            }

        }

//改进后

   public  bool DeleteDuplicatedRasterDataset(ref IRasterCatalog targetRasterCatalog,string shppath,string shpname)

        {

            try

            {

                IFeatureClass pFeatureClass = targetRasterCatalog as IFeatureClass;

                int featurecount = pFeatureClass.FeatureCount(null);

                int oidindex = pFeatureClass.FindField(pFeatureClass.OIDFieldName);

                int shpindex = pFeatureClass.FindField(pFeatureClass.ShapeFieldName);

                IList<int> oids = new List<int>();

                IFeatureCursor pFeatureCursor = pFeatureClass.Search(null, false);

                IFeature pFeature = pFeatureCursor.NextFeature();

                while (pFeature != null)

                {

                    oids.Add((int)pFeature.get_Value(oidindex));

                    pFeature = pFeatureCursor.NextFeature();

                }

                ITopologicalOperator pTopologicalOperator;

                IGeometryBag pGeometryBag;

                IGeometryCollection pGeometryCollection;

                object before = Type.Missing;

                object after = Type.Missing;

                IFeature pnFeature;

                ISpatialFilter pnspatialFilter;

                IFeatureCursor pnFeatureCursor;

                IFeatureClass pHelperFeatureClass;

                IFeature pHelperFeature;

                int count = oids.Count;

                int i = 0;

                this.prgBarAll.Maximum = count;

                foreach (int o in oids)

                {

                    this.Refresh();

                    lblPrgInfo.Text =  "共计" + count.ToString()+"条数据，" +"开始扫描第" + (i+1).ToString() + "个feature";

                    if (AEHelper.IsFeatureExisted(pFeatureClass, o))

                    {

                        pFeature = pFeatureClass.GetFeature(o);

                        pGeometryBag = new GeometryBagClass();

                        pGeometryCollection = (IGeometryCollection)pGeometryBag;

                        pnspatialFilter = new SpatialFilterClass();

                        pnspatialFilter.Geometry = pFeature.Shape;

                        pnspatialFilter.GeometryField = pFeatureClass.ShapeFieldName;

                        pnspatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;

                        pnFeatureCursor = pFeatureClass.Search(pnspatialFilter, false);

                        pnFeature = pnFeatureCursor.NextFeature();

                        while (pnFeature != null)

                        {

                            if (pnFeature.OID != o)

                            {

                                pGeometryCollection.AddGeometry(pnFeature.Shape, ref before, ref after);

                            }

                            pnFeature = pnFeatureCursor.NextFeature();

                        }

                        if (pGeometryCollection.GeometryCount > 0)

                        {

                            pTopologicalOperator = new PolygonClass();

                            pTopologicalOperator.ConstructUnion((IEnumGeometry)pGeometryBag);

                            //赋值给一个feature，然后再获取feature.shape

                            pHelperFeatureClass = AEHelper.GetFeatureClassFromShapefileOnDisk(shppath, shpname);

                            pHelperFeature = pHelperFeatureClass.CreateFeature();

                            pHelperFeature.Shape = (IGeometry)pTopologicalOperator;

                            pHelperFeature.Store();

                            pHelperFeature = pHelperFeatureClass.CreateFeature();

                            pHelperFeature.Shape = pFeature.Shape;

                            pHelperFeature.Store();

                            if (((IRelationalOperator)pTopologicalOperator).Contains(pFeature.Shape))

                            {

                                lblPrgInfo.Text =  "共计"+ count.ToString()+"条数据，"+"删除第" + (i + 1).ToString() + "个feature";

                                pFeature.Delete();

                            }

                            (pHelperFeatureClass as ITable).DeleteSearchedRows(null); //删除所有的feature

                        }

                    }

                    this.prgBarAll.Value = ++i;

                }

                MessageBox.Show("删除完成！");

                return true;

            }

            catch (Exception ex)

            {

                MessageBox.Show(ex.ToString());

                return false;

            }

        }

    }