C++、GDAL创建shapefile文件

#include "ogrsf_frmts.h" #include <iostream> using namespace std; int main() { const char *pszDriverName = "ESRI Shapefile"; OGRSFDriver *poDriver; OGRRegisterAll(); poDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName( pszDriverName ); if( poDriver == NULL ) { printf( "%s driver not available.\n", pszDriverName ); exit( 1 ); } OGRDataSource *poDS; poDS = poDriver->CreateDataSource( "point_out.shp", NULL );//但文件夹中以前存在point_out.shp文件时会报错; if( poDS == NULL ) { printf( "Creation of output file failed.\n" ); exit( 1 ); } OGRLayer *poLayer; poDS->CreateLayer( "point_out_layer", NULL, wkbPoint, NULL ); poLayer=NULL; poLayer = poDS->GetLayer(0);//shp文件只有一个图层; if( poLayer == NULL ) { printf( "Layer creation failed.\n" ); exit( 1 ); } OGRFieldDefn oField( "Name", OFTString ); oField.SetWidth(32); if( poLayer->CreateField( &oField ) != OGRERR_NONE ) { printf( "Creating Name field failed.\n" ); exit( 1 ); } double x, y; char szName[33]; cout<<"输入：x,y,name"<<"(Example:30,30,lu)"<<"逗号为英文下的逗号，否则按回车便结束程序"<<endl; cout<<"要想结束输入，按ctrl+d,再按回车键"<<endl; while( !feof(stdin) // stdin文件流的结束符按ctrl+d,这样便结束while循环; && fscanf( stdin, "%lf,%lf,%32s", &x, &y, szName ) == 3 ) { OGRFeature *poFeature; poFeature=OGRFeature::CreateFeature(poLayer->GetLayerDefn());//必须用CreateFeature来生成对象，用new生成对象出错 //poFeature = new OGRFeature( poLayer->GetLayerDefn() );//必须用CreateFeature来生成对象，用new生成对象出错 poFeature->SetField( "Name", szName ); OGRPoint pt; pt.setX( x ); pt.setY( y ); poFeature->SetGeometry( &pt ); if( poLayer->CreateFeature( poFeature ) != OGRERR_NONE ) { printf( "Failed to create feature in shapefile.\n" ); exit( 1 ); } OGRFeature::DestroyFeature( poFeature ); } OGRDataSource::DestroyDataSource( poDS ); }