zoukankan      html  css  js  c++  java
  • 删除重叠的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;
                }
            }
        }
  • 相关阅读:
    用PHPMailer在本地win环境,可以接收到邮件和附件,但在linux环境只能接收邮件信息接不到附件,是我的路
    linux 下 用phpmailer类smtp发送邮件始终不成功,提示:ERROR: Failed to co
    linux 下 用phpmailer类smtp发送邮件始终不成功,提示:ERROR: Failed to co
    phpmailer的SMTP ERROR: Failed to connect to server: 10
    SDK是什么?什么是SDK
    查看php的配置文件Php.ini的位置
    紧急求助!配置SMTP插件出错,SMTP connect() failed
    PHP move_uploaded_file() 函数
    HTML 5 video 视频标签全属性详解
    RAID级别与规范
  • 原文地址:https://www.cnblogs.com/zhangjun1130/p/2110181.html
Copyright © 2011-2022 走看看