gdal 多面转换和图层融合（dissolve）

潘达小新

已于 2022-06-23 08:12:13 修改

阅读量1.4k

点赞数

分类专栏： GDAL 文章标签： gdal 多部件转换

于 2022-01-24 10:00:34 首次发布

本文链接：https://blog.csdn.net/fangyu723/article/details/122661616

版权

GDAL 专栏收录该内容

11 篇文章 0 订阅

订阅专栏

        private void multipolygon2polygon(Geometry geo, ref List<Geometry> polyList)
        {
            Geometry temGeo = null;
            int iGeoCount = geo.GetGeometryCount();
            for (int i = 0; i < iGeoCount; i++)
            {
                temGeo = geo.GetGeometryRef(i);
                wkbGeometryType geoType = temGeo.GetGeometryType();
                if (geoType == wkbGeometryType.wkbLineString)
                {
                    polyList.Add(geo);
                    return;
                }
                if (geoType == wkbGeometryType.wkbMultiPolygon)
                {
                    multipolygon2polygon(temGeo, ref polyList);

                }
                else
                {
                    int iCount = temGeo.GetGeometryCount();
                    Geometry mzPoly = geo.Clone();
                    mzPoly.RemoveGeometry(0);
                    int iTempCount = mzPoly.GetGeometryCount();
                    while (iTempCount > 0)
                    {
                        mzPoly.RemoveGeometry(0);
                        iTempCount = mzPoly.GetGeometryCount();
                    }
                    mzPoly.AddGeometry(temGeo); ;
                    polyList.Add(mzPoly);

                }

            }
        }

        private void polygon2multipolygon(Geometry SingleGeo, ref Geometry MultiGeo)
        {
            MultiGeo = null;
            MultiGeo = new Geometry(wkbGeometryType.wkbMultiPolygon);
            MultiGeo.AddGeometry(SingleGeo);

        }

        /// <summary>
        /// 图层融合
        /// </summary>
        /// <param name="input">输入图层</param>
        /// <param name="dissolve">结果图层</param>
        private void dissolve(Layer input, Layer dissolve)
        {
            try
            {

                FeatureDefn defn = dissolve.GetLayerDefn();
                Geometry multi = new Geometry(wkbGeometryType.wkbMultiPolygon);
                Feature inputFea = null;
                while ((inputFea = input.GetNextFeature()) != null)
                {
                    Geometry geo = inputFea.GetGeometryRef();
                    geo.CloseRings();
                    string wkt = "";
                    geo.ExportToWkt(out wkt);

                    multi.AddGeometryDirectly(Geometry.CreateFromWkt(wkt));
                }
                Geometry union = multi.UnionCascaded();
                int iGeoCount = union.GetGeometryCount();
                for (int i = 0; i < iGeoCount; i++)
                {
                    string wkt = "";
                    //union.GetGeometryRef(i).ExportToWkb();
                    union.GetGeometryRef(i).ExportToWkt(out wkt);
                    Geometry poly = Geometry.CreateFromWkt(wkt);

                    Feature fea = new Feature(defn);
                    fea.SetGeometry(poly);
                    dissolve.CreateFeature(fea);
                }

            }
            catch (Exception ex)
            {
                throw new Exception("图层融合异常" + ex.Message);
            }
        }

参考：

https://www.cnpython.com/qa/148110