Как создать полигональный шейп-файл ESRI на С++?

Я использую opencv в С++ для извлечения набора 2D-точек из геопространственного изображения. Используя эти точки в качестве вершины, я хочу создать шейп-файл полигона ESRI.

Учебники по Gdal API (https://gdal.org/1.11/ogr/ogr_apitut.html) есть только примеры того, как создавать шейп-файлы точек и линий. Такие функции, как «OGR_G_AddPoint_2D», не работают с шейп-файлами полигонального типа.

У меня есть функция, которая добавляет одну точку в точечный слой и работает хорошо.

void addPoint(OGRLayerH hLayer, long x, long y, const char *szName){  
    OGRFeatureH hFeature;
    hFeature = OGR_F_Create( OGR_L_GetLayerDefn( hLayer ) );
    OGR_F_SetFieldString( hFeature, OGR_F_GetFieldIndex(hFeature, "Name"), szName );

    OGRGeometryH hPt;
    hPt = OGR_G_CreateGeometry(wkbPoint);
    OGR_G_SetPoint_2D(hPt, 0, x, y);
    OGR_F_SetGeometry( hFeature, hPt );
    OGR_G_DestroyGeometry(hPt);

    if( OGR_L_CreateFeature( hLayer, hFeature ) != OGRERR_NONE ){
       printf( "Failed to create feature in shapefile.\n" );
       exit(1);
    }
    OGR_F_Destroy( hFeature );
};

Попытался изменить его на форму многоугольника. Ошибка не печатается, но шейп-файл не сохраняется.

void addPolygon(OGRLayerH hLayer, double *x, double *y, const char *szName){  
    OGRFeatureH hFeature;
    hFeature = OGR_F_Create( OGR_L_GetLayerDefn( hLayer ) );
    OGR_F_SetFieldString(hFeature, OGR_F_GetFieldIndex(hFeature, "Name"), szName);

    OGRGeometryH hGeo;
    hGeo = OGR_G_CreateGeometry(wkbLineString);

    cout << "Performing for loop" << endl;
    for (int i=0; i < 4; i++){
        OGR_G_SetPoint_2D(hGeo, 0, x[i], y[i]);
    }
    OGR_G_CloseRings(hGeo);
    hGeo = OGR_G_ForceToPolygon(hGeo);

    if( OGR_L_CreateFeature( hLayer, hFeature ) != OGRERR_NONE ){
       printf( "Failed to create feature in shapefile.\n" );
       exit(1);
    }
    OGR_F_Destroy( hFeature );
    return;
};

c shapefile gdal ogr
person underfloor    schedule 21.05.2019    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Для функции OGR_G_SetPoint_2D требуется индекс i точки. Если все индексы равны 0, вместо многоугольника будет отображаться одна точка. Это работает для меня:

void addPolygon(OGRLayerH hLayer, double *x, double *y, int n, const char *szName)
{  OGRFeatureH hFeature=OGR_F_Create(OGR_L_GetLayerDefn(hLayer));
   OGR_F_SetFieldString(hFeature,OGR_F_GetFieldIndex(hFeature,"Name"),szName);
   OGRGeometryH hGeo=OGR_G_CreateGeometry(wkbLinearRing);
   printf("Performing LinearRing of %d vertices\n",n);
   for(int i=0;i<n;i++) { OGR_G_SetPoint_2D(hGeo,i,x[i],y[i]); }
   OGR_G_CloseRings(hGeo);
   hGeo=OGR_G_ForceToPolygon(hGeo);
   OGR_F_SetGeometry(hFeature,hGeo);
   OGR_G_DestroyGeometry(hGeo);
   if(OGR_L_CreateFeature(hLayer,hFeature) != OGRERR_NONE)
   {  printf("Failed to create feature in shapefile.\n"); exit(1); }
   OGR_F_Destroy(hFeature);
   return;
}
person fran4004    schedule 25.10.2019