我已经编写了这段代码,以便在opengl窗口中显示线形图文件 实际上它是代码Rendering Shapefile in OpenGL here in code-project,但我正在对它进行一些更改,以便使用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
的地址不同!
您认为我的代码有什么问题?
我尽可能地尽量减少代码。因此,如果您想自己调试代码,可以在code-project下载代码,然后根据我上面编写的代码对其进行更改。
OpenShapeFile
方法代替代码
相应的方法并对代码进行其他更改
我还上传了我的代码here in 4shared,您可以将其带到此处,然后链接GDAL Library并运行它。 当然,在约翰回答的帮助下,我的问题得到解决。但我仍然会感谢有人告诉我上述代码的问题是什么?
答案 0 :(得分: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
进行类似的更改。一般情况下,如果您发现自己为一个对象分配内存,然后在同一范围内删除 ,则表明您根本不需要分配内存。
这些都不能解决仍然隐藏的错误,但重点是通过简化代码并分配更少的内存,您将更接近真正的错误,这无疑是某种内存损坏。
答案 1 :(得分:0)
我遇到了同样的问题。你的代码是正确的。
在我的情况下,问题是PATH中的其他gdal.dll。确保您的可执行文件加载了正确的库。
(gdal.dll附带了很多应用程序,并将其文件夹添加到PATH。就像QGis或GeoConcept ...)