1

我正在尝试从特定 shapefile 中读取所有特征数据。在这种情况下，我使用 DotSpatial 打开文件，并且我正在迭代功能。这个特殊的 shapefile 只有 9mb 大小，而 dbf 文件是 14mb。大约有 75k 个特征需要循环。

请注意，这一切都是通过控制台应用程序以编程方式完成的，因此没有渲染或任何涉及的内容。

加载形状文件时，我重新投影，然后进行迭代。加载重新投影非常快。然而，一旦代码到达我的 foreach 块，加载数据需要将近 2 分钟的时间，并且在 VisualStudio 中调试时使用大约 2GB 的内存。对于相当小的数据文件来说，这似乎非常非常过分。

我已经从命令行在 Visual Studio 之外运行了相同的代码，但是时间仍然大约是 2 分钟，并且该进程需要大约 1.3GB 的内存。

有没有办法加快速度？

下面是我的代码：

// Load the shape file and project to GDA94
Shapefile indexMapFile = Shapefile.OpenFile(shapeFilePath);
indexMapFile.Reproject(KnownCoordinateSystems.Geographic.Australia.GeocentricDatumofAustralia1994);

// Get's slow here and takes forever to get to the first item
foreach(IFeature feature in indexMapFile.Features)
{
    // Once inside the loop, it's blazingly quick.
}

有趣的是，当我使用 VS 即时窗口时，它超级超级快，完全没有延迟......

4

2 回答 2

4

这对于非常大的文件（120 万个特征）也有同样的问题，填充 .Features 集合永远不会结束。

但是，如果您要求该功能，则您没有内存或延迟开销。

        int lRows = fs.NumRows();
        for (int i = 0; i < lRows; i++)
        {

            // Get the feature
            IFeature pFeat = fs.GetFeature(i); 

            StringBuilder sb = new StringBuilder();
            {
                sb.Append(Guid.NewGuid().ToString());
                sb.Append("|");
                sb.Append(pFeat.DataRow["MAPA"]);
                sb.Append("|");
                sb.Append(pFeat.BasicGeometry.ToString());
            }
            pLinesList.Add(sb.ToString());
            lCnt++;

            if (lCnt % 10 == 0)
            {
                pOld = Console.ForegroundColor;
                Console.ForegroundColor = ConsoleColor.DarkGreen;
                Console.Write("\r{0} de {1} ({2}%)", lCnt.ToString(), lRows.ToString(), (100.0 * ((float)lCnt / (float)lRows)).ToString());
                Console.ForegroundColor = pOld;
            }

        }

寻找 GetFeature 方法。

于 2017-06-06T21:30:33.920 回答
4

我已经设法弄清楚了...

出于某种原因，在功能上调用 foreach 非常缓慢。

但是，由于这些文件具有与特征的 1-1 映射 - 数据行（每个特征都有一个相关的数据行），因此我将其稍微修改为以下内容。现在非常快.. 开始迭代不到一秒钟。

// Load the shape file and project to GDA94
Shapefile indexMapFile = Shapefile.OpenFile(shapeFilePath);
indexMapFile.Reproject(KnownCoordinateSystems.Geographic.Australia.GeocentricDatumofAustralia1994);

// Get the map index from the Feature data
for(int i = 0; i < indexMapFile.DataTable.Rows.Count; i++)
{

    // Get the feature
    IFeature feature = indexMapFile.Features.ElementAt(i);

    // Now it's very quick to iterate through and work with the feature.
}

我想知道为什么会这样。我想我需要查看 IFeatureList 实现上的迭代器。

干杯，贾斯汀

于 2016-02-12T09:31:35.353 回答