PIE-SDK For C++打开自定义栅格数据

1、数据介绍

   信息提取和解译的过程中，经常会生成一部分中间临时栅格数据，这些数据在执行完对应操作后就失去了存在的价值，针对这种情况，PIE增加了内存栅格数据集，来协助用户完成对自定义栅格数据的读取和显示。下面以一副影像数据的抠图算法来演示自定义栅格数据的应用。

2、实现思路及原理说明

第一步	打开原始栅格数据；
第二步	读取部分数据至内存中，设置读取的起始点、读取范围等参数；
第三步	建立栅格内存数据集；
第四步	将内存中的数据写入内存数据集；
第五步	通过内存数据集创建栅格图层；
第六部	添加栅格图层到地图并刷新；

3、核心接口与方法

接口/类	方法	说明
SysDataSource::RasterDatasetPtr	GetBandCount()	获取波段个数
	GetRasterBand()	获取指定索引的波段
	Read()	将数据读取至内存中
	Write()	将内存数据写入栅格数据集中
	SetGeoTransform	设置六参数
	SpatialReference	空间参考
	GetGeoTransform	获取六参数
SysDataSource::DatasetFactory	CreateRasterDataset()	创建栅格数据集

4、示例代码

项目名称	百度云盘地址下/PIE示例程序/03.数据加载/11.打开自定义栅格数据
数据位置	百度云盘地址下/PIE示例数据/栅格数据/ 04.World/World.tif
视频位置	百度云盘地址下/PIE视频教程/03.数据加载/11.打开自定义栅格数据.avi
示例代码
方法（一）打开自定义栅格数据 void PIEMainWindow::On_ActionGetMemRaster_Triggered(bool checked) {     SysCarto::RasterLayerPtr rasterLayerPtr = m_pCurrentControl->GetActiveView()->GetCurrentLayer();     if (rasterLayerPtr == nullptr) return;     //1、获取源Dataset     SysDataSource::RasterDatasetPtr sourceRasterPtr = rasterLayerPtr->GetDataset();     int mWidth = 500;     int mHeigh = 600;     int count = sourceRasterPtr->GetBandCount();     double* geoTrans = new double[6];     sourceRasterPtr->GetGeoTransform(geoTrans);     QVector<int> bandMap;     int* bandList = new int[count];          for (int i=0;i<count;i++)     {         bandMap.insert(i, i + 1);         bandList[i] = i + 1;     }     //2、读取数据至内存     SysDataSource::PixelBufferPtr bufferPtr=sourceRasterPtr->Read(0, 0, mWidth, mHeigh, mWidth, mHeigh, bandMap);     //3、创建内存数据集     SysDataSource::RasterDatasetPtr memRasterPtr = SysDataSource::DatasetFactory::Instance()->CreateRasterDataset("", mWidth, mHeigh, count,         sourceRasterPtr->GetRasterBand(0)->GetRasterDataType(), "MEM");     if (memRasterPtr!=nullptr)     {         bool ok=memRasterPtr->Write(0, 0, mWidth, mHeigh, bufferPtr->GetData(), mWidth, mHeigh, sourceRasterPtr->GetRasterBand(0)->GetRasterDataType(),             count, bandList);         if (ok)         {             //4、设置空间参考和仿射变换参数             memRasterPtr->SetSpatialReference(sourceRasterPtr->GetSpatialReference());             memRasterPtr->SetGeoTransform(geoTrans);               SysCarto::LayerPtr memLayer = SysCarto::LayerFactory::CreateDefaultRasterLayer(memRasterPtr);             memLayer->SetName("内存图层");             m_pCurrentControl->GetMap()->AddLayer(memLayer);             m_pCurrentControl->GetActiveView()->Refresh();         }     } }

5、示例截图

前后对比图如下