ogr矢量图层的创建(根据栅格影像创建矢量边框).docx

总体来说这个过程就是构建数据源 -> 构建层 ->构建要素 -> 构建形状 ->
关闭数据源。
实例代码:
// Create  : 定义控制台应用程序的入口点。
#include ""
// 要包含的 GDAL头文件
#include <>
#include <>
#ment(lib,"")
bool Creatshape(const char* pszFileName ,int line,int row);
int _tmain(int argc, _TCHAR* argv[])
{
const char *pszFileName="";
/*int m,n;
char pathName[100];
cout<<"请输入影像(路径名加文件名) \n";
cin>>pathName;
cout<<"请输入要划分的行数: m ";
cin>>m;
cout<<"请输入要划分的列数: n ";
cin>>n;
pszFileName=pathName;*/
Creatshape(pszFileName,7,9);
return 0;
}
/**********************************************************
**************/
/* 创建过程:构建数据源 ->构建层 ->构建要素 ->构建形状 ->关闭数据
源
/* 参数 pszFileName 为输入的文件名
参数 linenum 为划分的行数
参数 rows 为划分的列数 */
/**********************************************************
**************/
bool Creatshape(const char* pszFileName ,int linenum,int rows)
{
// 获取影像信息
GDALDataset *poDataSet;
GDALAllRegister();
poDataSet=(GDALDataset*)GDALOpen(pszFileName,GA_ReadOnly);
// 打开数据集
if (poDataSet==NULL)
{
//
"Failed to open this dataset!"; 代开失败的话则给出提示并退
出
exit(1);
}
函
double Trans[6];//坐标转换参数数组 ,作为 GetGeoTransform() 数
的参数
int width,height; //影像的像行列数
//
width=poDataSet->GetRasterXSize();获取影像列数, 为后续划分
网格做准备
获
height=poDataSet->GetRasterYSize();//取影像行数, 为后续划分
网格做准备
poDataSet->GetGeoTransform(Trans);
//Xp = padfTransform[0] + P*padfTrans