0

我编写了这段代码是为了在 opengl 窗口中显示线形文件。
实际上，它是代码项目中的 OpenGL 中渲染 Shapefile 的代码，但我对其进行了一些更改，以便使用 OGR 库而不是 shapelib 读取形状。

#include "ogrsf_frmts.h" 
//#include "shapelib\shapefil.h"

void OpenShapeFile(char* filename)
{
   int i = 0;
   int j = 0;
   OGRErr error;
   OGRDataSource *poDataSource;
   poDataSource = OGRSFDriverRegistrar::Open(filename,false);
   OGRLayer *poLayer;
   poLayer = poDataSource ->GetLayer(0);
   OGREnvelope *poEnvelope = new OGREnvelope();
   error = poLayer ->GetExtent(poEnvelope,true);
   sBoundingBox.fMaxX = poEnvelope ->MaxX;
   sBoundingBox.fMaxY = poEnvelope ->MaxY;
   sBoundingBox.fMinX = poEnvelope ->MinX;
   sBoundingBox.fMinY = poEnvelope ->MinY;
   delete poEnvelope;
   OGRwkbGeometryType GeometryType = poLayer ->GetGeomType();
   int NumberOfFeatures = poLayer ->GetFeatureCount(true);
   poLayer ->ResetReading();

   //Line Shapefile
   if ( wkbFlatten ( GeometryType ) == wkbLineString )
   {
       OGRFeature *poFeature;
       MyLineString2D lineString;

       //temporary pointer in order to store coordinates of individual line vertexes
       OGRPoint *poPointTemp = new OGRPoint();
       for ( i = 0; i < NumberOfFeatures; i++ )
       {
           poFeature = poLayer ->GetNextFeature();
           OGRGeometry *poGeometry;
           poGeometry = poFeature ->GetGeometryRef();
           if ( poGeometry != NULL )
           {
               OGRLineString *poLineString = (OGRLineString *)poGeometry;
               int NumberOfVertexes = poLineString ->getNumPoints();
               lineString.vPointList.resize(NumberOfVertexes);
               for ( j = 0; j < NumberOfVertexes ; j++ )
               {
                   poLineString ->getPoint(j,poPointTemp);
                   MyPoint2D ptTemp;
                   ptTemp.dX = poPointTemp ->getX();
                   ptTemp.dY = poPointTemp ->getY();
                   lineString.vPointList.at(j) = ptTemp;
               }
               vLines.push_back(lineString);
           }
           OGRFeature::DestroyFeature( poFeature );
       }
   delete poPointTemp;
   } 
}  

void initializeGL()
{
    //glClearColor (0.0, 0.0, 0.0, 0.0);
    glClearColor (1.0, 1.0, 1.0, 1.0);  

int main(int argc, char** argv)
{
//OpenShapeFile("Shapefiles\\poi.shp");//Point Shapefile
OpenShapeFile("Shapefiles\\strassen.shp");//Line Shapefile
//OpenShapeFile("Shapefiles\\gruenflaechen.shp");//Polygon Shapefile

当我评论该行时： 

delete poPointTemp;

没有错误，并且正确地在 OpenGL 窗口上绘制了形状。但正如你所知，这是一种内存泄漏，所以当我不需要时，poPointTemp我应该删除它。
但在这条线上： 

delete poPointTemp;

我收到运行时错误：

在此处输入图像描述

我已经调试了我的代码，运行此行之前的地址和内容poPointTemp如下： 

poPointTemp 0x00503a90 {x=3435936.3300000001 y=5790327.5999999996 z=0.00000000000000000 }

我的意思poPointTemp是不是 NULL 指针，而且我们无权访问的地址与poPointTemp's 地址不同！！！

你认为我的代码有什么问题？

我试图尽可能地减少代码。因此，如果您想自己调试代码，您可以在代码项目中下载它，然后根据我上面编写的代码对其进行更改。

  • 请使用我的OpenShapeFile方法代替代码对应的方法，并对代码进行其他更改。
    另外我已经在 4shared中上传了我的代码，你可以把它拿在这里，然后链接GDAL 库并运行它。

当然，在约翰的回答的帮助下，我的问题得到了解决。但我仍然会感谢有人告诉我上述代码有什么问题？

4

2 回答 2

1

这不是一个真正的答案，但它可能会有所帮助，并且发表评论有点太长了。还要记住，我不知道这个库，所以我的建议可能不正确。

毫无疑问，这个错误掩盖了代码中其他地方的真正错误。正如您所说，此时您应该删除内存。但是，在我看来，您确实在不必要地分配内存。例如这个

OGREnvelope *poEnvelope = new OGREnvelope();
error = poLayer ->GetExtent(poEnvelope,true);
sBoundingBox.fMaxX = poEnvelope ->MaxX;
sBoundingBox.fMaxY = poEnvelope ->MaxY;
sBoundingBox.fMinX = poEnvelope ->MinX;
sBoundingBox.fMinY = poEnvelope ->MinY;
delete poEnvelope;

可以在没有任何内存分配的情况下重写

OGREnvelope poEnvelope;
error = poLayer ->GetExtent(&poEnvelope,true);
sBoundingBox.fMaxX = poEnvelope.MaxX;
sBoundingBox.fMaxY = poEnvelope.MaxY;
sBoundingBox.fMinX = poEnvelope.MinX;
sBoundingBox.fMinY = poEnvelope.MinY;

可以对poPointTemp. 通常，如果您发现自己为一个对象分配内存，然后在同一范围内将其删除，则表明您根本不需要分配内存。

这些都不能修复仍然隐藏的错误，但关键是通过简化代码并分配更少的内存，您将更接近真正的错误，这无疑是某种内存损坏。

于 2013-09-11T07:36:42.223 回答
0

我有同样的问题。你的代码是正确的。

在我的情况下，问题是PATH中的其他gdal.dll。确保您的可执行文件加载正确的库。

（gdal.dll 附带了很多关闭应用程序，并将它们的文件夹添加到 PATH。像 QGis 或 GeoConcept ......）

于 2014-04-18T14:18:52.373 回答