PIE-SDK For C++栅格数据的创建

1．功能简介

    目前在地理信息领域中数据包括矢量和栅格两种数据组织形式。每一种数据有不同的数据格式，目前PIE SDK支持多种数据格式的数据创建，下面对栅格数据格式的数据创建功能进行介绍。

2．功能实现说明

2.1 实现思路及原理说明

第一步	获取要读取的栅格数据
第二步	获取栅格数据的读取
第三步	构建新的栅格数据集
第四步	栅格数据写入
第五步	设置空间参考
第六步	设置仿射变换

2.2 核心接口与方法

接口/类	方法/属性	说明
SysDataSource::RasterDatasetPtr	GetBandCount()	获取波段数
	Read()	栅格读取
	Write()	栅格写入
	GetGeoTransform()	获取仿射变换参数
	SetGeoTransform()	设置仿射变换参数
	GetSpatialReference()	空间参考
SysDataSource::DatasetFactory	CreateRasterDataset()	空间栅格数据集

2.3 示例代码

项目路径	百度云盘地址下/PIE示例程序/04数据操作/10栅格数据的创建
数据路径	百度云盘地址下/PIE示例数据/栅格数据/04.World/World.tif
视频路径	百度云盘地址下/PIE视频教程/04数据操作/10栅格数据的创建.avi
示例代码
void PIEMainWindow::On_ActionCreateRaster_Triggered(bool checked) {     SysCarto::RasterLayerPtr rasterLayerPtr = m_pCurrentControl->GetActiveView()->GetCurrentLayer();     if (rasterLayerPtr == nullptr) return;       SysDataSource::RasterDatasetPtr datasetPtr = rasterLayerPtr->GetDataset();     int BandCount = datasetPtr->GetBandCount();     //2栅格数据范围等信息的读取     int xSize = datasetPtr->GetRasterXSize();//行     int ySize = datasetPtr->GetRasterYSize();//列     QVector<int> vectorBand;     int* bandMap = new int[BandCount];     for (int i = 0; i < BandCount; i++)     {         bandMap[i] = i + 1;         vectorBand.insert(i, i + 1);     }     QString filter = "Raster Files (*.tif *.tiff *.img *.bmp *.jpg *.ldf *.dat *.1bd *.1b)";     QString saveFilePath= QFileDialog::getSaveFileName(nullptr, "保存", "", filter);     if (saveFilePath.isEmpty()) return;     SysDataSource::PixelBufferPtr buffer=datasetPtr->Read(0, 0, xSize, ySize, xSize, ySize, vectorBand);       SysDataSource::RasterDatasetPtr newDataSet = SysDataSource::DatasetFactory::Instance()->CreateRasterDataset(saveFilePath, xSize, ySize, BandCount, datasetPtr->GetRasterBand(0)->GetRasterDataType(), "GTiff");     double* geo=new double[6];     datasetPtr->GetGeoTransform(geo);     newDataSet->SetSpatialReference(datasetPtr->GetSpatialReference());     newDataSet->SetGeoTransform(geo);     bool flag= newDataSet->Write(0, 0, xSize, ySize, buffer->GetData(), xSize, ySize, datasetPtr->GetRasterBand(0)->GetRasterDataType(), BandCount, bandMap);     if (flag)     {         SysCarto::LayerPtr newLayer = SysCarto::LayerFactory::Instance()->CreateDefaultRasterLayer(newDataSet);         if (newLayer!=nullptr)         {             m_pCurrentControl->GetMap()->AddLayer(newLayer);             m_pCurrentControl->GetActiveView()->Refresh();         }     } }

2.4 示例截图