GIS地理数据处理(c++)栅格矢量化
2023-12-20 10:50:54
系列文章目录
栅格矢量化
前言
这段代码主要执行GDALFPolygonize
CPLErr GDALFPolygonize(GDALRasterBandH hSrcBand, GDALRasterBandH hMaskBand, OGRLayerH hOutLayer, int iPixValField, char **papszOptions, GDALProgressFunc pfnProgress, void *pProgressArg)
一、完整代码
示例:pandas 是基于NumPy 的一种工具,该工具是为了解决数据分析任务而创建的。
代码如下:
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
import warnings
warnings.filterwarnings('ignore')
import ssl
ssl._create_default_https_context = ssl._create_unverified_context
二、详细代码
代码如下:
#include<gdal_priv.h>
#include "ogrsf_frmts.h"
#include <errno.h>
#include<iostream>
#include<gdal_alg.h>
using namespace std;
int main()
{
GDALAllRegister();
//栅格路径
const char* raster_path = "E:\\gdal\\JL_41_.tif";
// 矢量路径
string shp_path = "E:\\gdal\\JL_41_64_5.shp";
// 读取栅格文件
// 将GDALDatasetH类强制转为 GDALDataset
GDALDataset* Rasterdata = (GDALDataset*)GDALOpen(raster_path, GA_ReadOnly);
// 判断数据是否读取成功
if (Rasterdata == NULL)
{
printf("Open failed.\n");
exit(1);
}
// 获取波段
GDALRasterBand* RatserBand = Rasterdata->GetRasterBand(1);
// 创建空间参考
OGRSpatialReference spatialReference;
spatialReference.importFromEPSG(4326);
// 同GDALOpen,创建矢量输出文件
string driverName = "ESRI Shapefile";
GDALDriver* driver = GetGDALDriverManager()->GetDriverByName(driverName.c_str());
if (driver == NULL)
{
cout << "Faild Getting Driver:" << driverName << endl;
}
GDALDataset* shpdata = driver->Create(shp_path.c_str(), 0, 0, 0, GDT_Unknown, NULL);
if (shpdata == NULL)
{
cout << "create dataset failed :" << shp_path << endl;
return -1;
}
//
/*GDALDataset* shpdata = static_cast<GDALDataset*>(
GDALOpenEx(shp_path, GDAL_OF_VECTOR, NULL, NULL, NULL));*/
//判定数据是否读取成功
if (shpdata == NULL)
{
printf("Open failed.\n");
exit(1);
}
// 创建图层,存储矢量化的栅格
OGRLayer* layer = shpdata->CreateLayer("polygon", &spatialReference, wkbPolygon, NULL);
if (layer == NULL)
{
cout << "Faild Create Layer" << endl;
return -1;
}
OGRFieldDefn fieldDef1("value", OFTInteger);
fieldDef1.SetWidth(100);
// 判断该字段创建成功
if (layer->CreateField(&fieldDef1) != OGRERR_NONE)
{
cout << "字段创建成功" << endl;
return -1;
}
// 执行矢量化函数,GDALTermProgress显示矢量化进程
if(GDALFPolygonize(RatserBand,NULL,layer,0, NULL, GDALTermProgress, NULL) != CE_None);
{
// 函数执行成功,关闭数据集释放内存
GDALClose((GDALDatasetH)Rasterdata);
GDALClose((GDALDatasetH)shpdata);
// delete spatialReference;
/*spatialReference = NULL;*/
cout << "创建成功" << endl;
return 0;
}
//释放内存
GDALClose((GDALDatasetH)Rasterdata);
GDALClose((GDALDatasetH)shpdata);
return 0;
}
GDALTermProgress是矢量化函数的进程参数,还是蛮智能的
文章来源:https://blog.csdn.net/d8dongdong/article/details/135095030
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。 如若内容造成侵权/违法违规/事实不符,请联系我的编程经验分享网邮箱:veading@qq.com进行投诉反馈,一经查实,立即删除!
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。 如若内容造成侵权/违法违规/事实不符,请联系我的编程经验分享网邮箱:veading@qq.com进行投诉反馈,一经查实,立即删除!