gdal将点数据生成矢量文件（shape）

GDALAllRegister();
    OGRRegisterAll();
    const char *pszDriverName = "ESRI Shapefile";
    GDALDriver *poDriver;
    GDALAllRegister();
    poDriver = GetGDALDriverManager()->GetDriverByName(pszDriverName);
    if (poDriver == NULL)
    {
        return 0;
    }
    GDALDataset *poDS;
   poDS = poDriver->Create("d:/SHP/0331/0331Shp.shp", 0, 0, 0, GDT_Float64, NULL); //创建shp文件

if (poDS == NULL)
    {
        printf("Creation of output file failed.\n");
       return 0;
    }
    OGRLayer *poLayer;
    poLayer = poDS->CreateLayer("point_out", NULL, wkbPointZM, NULL);  //第二个参数为投影
    if (poLayer == NULL)
    {
        printf("Layer creation failed.\n");
        return 0;
   }

idField.SetPrecision(8);    //设置字段精度
    fir