使用gdal ogr创建shapefile文件（c++）

使用ogr库创建点状要素的shapefile文件以及将经纬度坐标转为投影坐标。实例如下：

#include "ogrsf_frmts.h"#include "gdal.h"#include "gdal_priv.h"#include "cpl_string.h" #include <string>#include <iostream>#include <strstream>using namespace std;void convertToShp(double longitude, double latitude, char *outshp){       //使属性表字段支持中文    CPLSetConfigOption("SHAPE_ENCODING","");    OGRRegisterAll();//注册所有的驱动    //创建ESRI shp文件    char *pszDriverName = "ESRI Shapefile";    //调用对Shape文件读写的Driver    OGRSFDriver *poDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName(pszDriverName);    if (poDriver == NULL)    {        cout<<pszDriverName<<"驱动不可用！"<<endl;        return;    }    //创建数据源    OGRDataSource *poDs = poDriver->CreateDataSource(outshp, NULL);    if (poDs == NULL)    {        cout<<"DataSource Creation Error"<<endl;        return;    }    //创建图层Layer    string outShapName = outshp;    string layerName = outShapName.substr(0, outShapName.length()-4);    //layerName.c_str()表示将string转为char *类型    //参数说明：新图层名称，坐标系，图层的几何类型，创建选项，与驱动有关    OGRLayer *poLayer = poDs->CreateLayer(layerName.c_str(), NULL, wkbPoint, NULL);    if (poLayer == NULL)    {        cout<<"Layer Creation Failed"<<endl;        OGRDataSource::DestroyDataSource(poDs);        return;    }    //下面创建属性表，我们在属性表中创建两列数据即可    //先创建一个“ID”整型属性    OGRFieldDefn oFieldId("ID", OFTInteger);    oFieldId.SetWidth(5);    poLayer->CreateField(&oFieldId);    //再创建一个"X"double属性    OGRFieldDefn oFieldX("X", OFTString);    oFieldX.SetWidth(30);    poLayer->CreateField(&oFieldX);    //再创建一个"Y"double属性    OGRFieldDefn oFieldY("Y", OFTString);    oFieldY.SetWidth(30);    poLayer->CreateField(&oFieldY);    //创建一个feature    OGRFeature *poFeature;      poFeature = OGRFeature::CreateFeature(poLayer->GetLayerDefn());//GetLayerDefn()获取当前图层的属性表结构    //给属性表中我们刚创建的列赋值    int i = 0;    poFeature->SetField("ID", i);    poFeature->SetField("X", longitude);    poFeature->SetField("Y", latitude);    i++;    //创建一个OGRPoint对象    OGRPoint point;    point.setX(longitude);    point.setY(latitude);    //point.setZ(0);    poFeature->SetGeometry(&point);    if(poLayer->CreateFeature(poFeature) != OGRERR_NONE )    {        printf( "Failed to create feature in shapefile.\n" );        exit( 1 );    }    OGRFeature::DestroyFeature(poFeature);    OGRDataSource::DestroyDataSource(poDs);     }//======= 经纬度转化为投影坐标=============double* transfer(double longitude, double latitude){    OGRSpatialReference oSourceSRS;    //EPSG code 代表特定的椭球体、单位、地理坐标系或投影坐标系等信息    //This method will initialize the spatial reference based on the passed in EPSG GCS or PCS code    oSourceSRS.importFromEPSG(4326);//EPSG:4326代表地理坐标系WGS1984    OGRSpatialReference oTargetSRS;    oTargetSRS.importFromEPSG(2029);// utm zone 17N    OGRCoordinateTransformation *poTransform;    poTransform = OGRCreateCoordinateTransformation(&oSourceSRS, &oTargetSRS);    if (poTransform == NULL)    {        cout<<"poTransform is null"<<endl;        exit(1);    }       if (!poTransform->Transform(1, &longitude, &latitude, NULL))    {        cout<<"transform failed"<<endl;        exit(1);    }    //poTransform->Transform(1, &longitude, &latitude, NULL);    double *inout = new double[2];    inout[0] = longitude;    inout[1] = latitude;    return inout;}int main(int argc, char *argv[]){           char *xCoordinate = argv[1];    char *yCoordinate = argv[2];    char *outShp = argv[3];    double x = atof(xCoordinate);    double y = atof(yCoordinate);    convertToShp(x, y, outShp);    cout<<"success! file is saved to "<<outShp<<endl;    double *transferedLongLat = transfer(116.246742, 40.022211);    cout<<"转换后的投影坐标为："<<transferedLongLat[0]<<","<<transferedLongLat[1]<<endl;    getchar();}

程序运行结束后即在我们指定的目录下生成点状要素的shp文件。
关于坐标转换，可以参考：http://blog.csdn.net/jingss_3/article/details/8240328