0

我有一个带有 EPSG:32749 的 shapefile,它将被插入 Oracle 数据库并显示在 geoserver 中。在此之前,我想使用点空间库将我的 shapefile 重新投影到 ESPG:4326,这是我的代码

var EXTRACTED_NAME = Server.MapPath("~/upload/shp/example/");
string shapeFilePath = @"\example.shp";
shapeFilePath = EXTRACTED_NAME + shapeFilePath;
Shapefile indexMapFile = Shapefile.OpenFile(shapeFilePath);
indexMapFile.Reproject(KnownCoordinateSystems.Geographic.World.WGS1984);

但是当我在geoserver中预览时,我的shapefile显示是这样的,当原来的shapefile是这样的

在此处输入图像描述

在此处输入图像描述

我的问题是,为什么重新投影到 EPSG 4326 的 shapefile 与原始文件不同?

谢谢

4

1 回答 1

0

这有点晚了,但您应该能够从几何图形中访问内环。您可能需要将 IGeometry 转换为 IPolygon 以专门处理内环,而不仅仅是 getGeometryN。以下代码尚未经过测试,但至少应该让您指出正确的方向。

Shapefile file = Shapefile.OpenFile(@"D:\Data\Counties\Counties.shp");
foreach(Feature f in file.Features){
    if(f.Geometry is IPolygon){
        IPolygon p = (IPolygon)f.Geometry;
        Debug.WriteLine("Feature " + f.Fid + "\n");
        foreach(ILineString innerRing in p.InteriorRings){
            // Do stuff with your innerRing
            Debug.WriteLine("Ring length : " + innerRing.Length);
        }
    }
    if (f.Geometry is IMultiPolygon)
    {
        IMultiPolygon multi = (IMultiPolygon)f.Geometry;
        for (int i = 0; i < multi.NumGeometries; i++)
        {
            IGeometry g = multi.GetGeometryN(i);
            if (g is IPolygon)
            {
                IPolygon p = (IPolygon)g;
                foreach (ILineString innerRing in p.InteriorRings)
                {
                    // Do stuff with your innerRing
                }
            }

        }

    }

}
于 2019-03-06T19:01:07.673 回答