GDAL:shp文件创建

创建一个只有点要素的shp文件

#include "gdal_priv.h"
#include "ogrsf_frmts.h"
#include <iostream>
using namespace std;
int main(int argc, char* argv[])
{
	OGRRegisterAll();
	const char *pszDriverName = "ESRI Shapefile";
	OGRSFDriver *poDriver;
	poDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName(pszDriverName);//注册OGRSFDriver驱动器

	if (poDriver == NULL)
	{
		cout << pszDriverName << " driver not available." << endl;
		return 0;
	}
	OGRDataSource *poDS;
	poDS = poDriver->CreateDataSource("D:/TEST/point_out2.shp", NULL);//创建shp数据集，如果已存在文件，会报错
	if (poDS == NULL)
	{
		cout << "Creation of point_out.shp file failed." << endl;
		return 0;
	}
	OGRLayer *poLayer;
	poLayer = poDS->CreateLayer("point_out", NULL, wkbPoint, NULL);//创建图层wkbPoint为点图层

	if (poLayer == NULL)
	{
		cout << "Layer creation failed." << endl;
		return 0;
	}
	OGRFieldDefn firstField("faci_code", OFTInteger);//创建字段
	OGRFieldDefn secondField("X", OFTReal);
	OGRFieldDefn thirdField("Y", OFTReal);
	firstField.SetWidth(32);//设置字段长度
	secondField.SetWidth(32);
	thirdField.SetWidth(32);
	poLayer->CreateField(&firstField);//给图层新增字段列表
	poLayer->CreateField(&secondField);
	poLayer->CreateField(&thirdField);

	double x, y;
	int code;
	for (int i = 0; i != 100; ++i)
	{
		code = i + 1;
		x = i;
		y = 100 + i;
		OGRFeature *poFeature;//创建要素指针
		poFeature = OGRFeature::CreateFeature(poLayer->GetLayerDefn());//构造函数中加入字段信息
		poFeature->SetField("faci_code", code);//将字段和对应数据进行绑定
		poFeature->SetField("X", x);
		poFeature->SetField("Y", y);
 		OGRPoint pt;//创建点几何类型
 		pt.setX(x);//给点几何类型赋值
 		pt.setY(y);
		poFeature->SetGeometry(&pt);//将当前点的几何类型和对应的要素绑定

		if (poLayer->CreateFeature(poFeature) != OGRERR_NONE)//将要素和图层进行绑定
		{
			cout << "Failed to create feature in shapefile." << endl;
			return 0;
		}
		OGRFeature::DestroyFeature(poFeature);//销毁要素指针
	}
	OGRDataSource::DestroyDataSource(poDS);//最后关闭数据源指针
	system("pause");
	return 0;
}

打开shp

#include "gdal_priv.h"
#include "ogrsf_frmts.h"
#include <iostream>
using namespace std;
int main(){
	GDALAllRegister();
	OGRRegisterAll();
	CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");
	CPLSetConfigOption("SHAPE_ENCODING", "");
	const char* shpFileName = "D:/分块影像/shp/shp1.shp";
	OGRDataSource *poSrcDS = OGRSFDriverRegistrar::Open(shpFileName);
	if (poSrcDS == NULL){
		cout << shpFileName << "文件打开失败！" << endl;
		return 0;
	}
	int iLayerNumber = poSrcDS->GetLayerCount();
	OGRLayer *poLayer = poSrcDS->GetLayer(0);
	if (poLayer==NULL)
	{
		cout << "获取图层失败！" << endl;
		return 0;
	}
	poLayer->ResetReading();//重置图层，清楚筛选格式
	OGRFeatureDefn *poDefn = poLayer->GetLayerDefn();//定义要素的指针
	int iFieldCount = poDefn->GetFieldCount();//获取要素的字段个数
	for (int i = 0; i < iFieldCount;i++)
	{
		OGRFieldDefn* poField = poDefn->GetFieldDefn(i);//定义字段指针
		cout << poField->GetNameRef() << "\n" << poField->GetFieldTypeName(poField->GetType()) << "\n" << poField->GetWidth() << "\n" << poField->GetPrecision() << endl;
		cout << "--------------------------------------" << endl;
	}
	cout << "要素个数 = "<<poLayer->GetFeatureCount() << endl;
	OGRFeature *poFeature = NULL;
	while ((poFeature = poLayer->GetNextFeature())!=NULL)
	{
		cout << "当前是第" << poFeature->GetFID() << "属性值" << endl;
		for (int i = 0; i < iFieldCount;i++)
		{
			OGRFieldDefn *POFIELD = poDefn->GetFieldDefn(i);
			const char *str = POFIELD->GetNameRef();
			OGRFieldType type_ = POFIELD->GetType();
			cout << poFeature->GetFieldAsString(i) << endl;
		}
		OGRGeometry *poGeometry = poFeature->GetGeometryRef();
		char*GeoJson = poGeometry->exportToJson();//导出到JSON文件
		char*GeoWkt;
		poGeometry->exportToWkt(&GeoWkt);//导出格式为WKT
		
		const char* nameGeo = poGeometry->getGeometryName();//获得几何体名称
		OGRwkbGeometryType GeoType = poGeometry->getGeometryType();//获得几何体类型
		//写出--------------------------------------------------
		FILE *fp;
		fp = fopen("D:/分块影像/shp/Jack.txt", "w");
		fwrite(GeoWkt, strlen(GeoWkt) + 1, 1, fp);
		fclose(fp);
	}
	system("pause");
	return 0;
}

#include "gdal_priv.h"
#include "ogrsf_frmts.h"
 
GDALAllRegister();
 
    const char * pFileName = "shp.tif";
    const char * pszDriverName = "ESRI Shapefile";
    GDALDriver * poDriver = GetGDALDriverManager()->GetDriverByName(pszDriverName);
    if(poDriver == NULL)
        return ;
    GDALDataset * poDs = poDriver->Create(pFileName, 0,0,0, GDT_Unknown, NULL);
    if(poDs == NULL)
        return ;
 
    //wkbPolygon表示绘制多边形；wkbLinearRing表示绘制线环
    OGRLayer * poLayer = poDs->CreateLayer("ring", NULL, wkbLineString, NULL);//不删
    if(poLayer == NULL)
    {
        return;
    }
 
 
    OGRFeature * poFeature = OGRFeature::CreateFeature(poLayer->GetLayerDefn());//放到DestroyFeature删
    //加入几何图形
 
    //例1，绘制线环，内部不标记颜色.线的绘制原则是，从第一个点开始，逐次绘制直到最后一个点。
    //对于闭合的线环，由于首点和尾点是一个点。所以4边形绘制需要5个点，首点被重复一次
    OGRLineString line;
    line.setNumPoints(5);
    line.setPoint(0, 0, 1000, 0);//第一个输入变量是序号，第二、三、四分别是线XYZ坐标
    line.setPoint(1, 0, 0, 0);
    line.setPoint(2, 1000, 0, 0);
    line.setPoint(3, 1000, 1000, 0);
    line.setPoint(4, 0, 1000, 0);
 
 
    poFeature->SetGeometry(&line);
 
 
 
    //例2，绘制多边形。多边形内部亦被颜色标记
    /*OGRPolygon poly;
    std::vector<QPoint> vec;
    vec.push_back(QPoint(0,1));
    vec.push_back(QPoint(1,1));
    vec.push_back(QPoint(1,0));
    vec.push_back(QPoint(0,0));
    OGRLinearRing ring;
    for(int i = 0; i < 4; i++)
    {
        ring.addPoint(vec[i].x(), vec[i].y());
    }
    ring.closeRings();
    poly.addRing(&ring);
    poFeature->SetGeometry(&poly);*/
 
    //例3 多个线
    /*OGRMultiLineString multiLine;
    OGRLineString line1, line2;
        line1.setNumPoints(5);
        line1.setPoint(0, 0, 1000, 0);//第一个输入变量是序号，第二、三、四分别是线XYZ坐标
        line1.setPoint(1, 0, 0, 0);
        line1.setPoint(2, 1000, 0, 0);
        line1.setPoint(3, 1000, 1000, 0);
        line1.setPoint(4, 0, 1000, 0);
        line2.setNumPoints(5);
        line2.setPoint(0, 500, 1500, 0);//第一个输入变量是序号，第二、三、四分别是线XYZ坐标
        line2.setPoint(1, 500, 500, 0);
        line2.setPoint(2, 1500, 500, 0);
        line2.setPoint(3, 1500, 1500, 0);
        line2.setPoint(4, 500, 1500, 0);
    multiLine.addGeometry((const OGRGeometry *)&line1);
    multiLine.addGeometry((const OGRGeometry *)&line2);
    poFeature->SetGeometry(&multiLine);*/
 
    if(poLayer->CreateFeature(poFeature) != OGRERR_NONE)
    {
        return;
    }
 
    OGRFeature::DestroyFeature(poFeature);
    GDALClose(poDs);
    OGRCleanupAll();