PIE SDK打开自定义栅格数据

 

1. 数据介绍

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

    下面以一副影像数据的抠图算法来演示自定义栅格数据的应用。

2. 实现思路及原理说明

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

3. 核心接口与方法

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

4. 示例代码

项目名称	百度云盘地址下/PIE示例程序/03.数据加载/11.打开自定义栅格数据
数据位置	百度云盘地址下/PIE示例数据/栅格数据/04.World/World.tif
视频位置	百度云盘地址下/PIE视频教程/03.数据加载/11.打开自定义栅格数据.avi
示例代码
方法（一）打开自定义栅格数据 if (mapControlMain.FocusMap.LayerCount == 0)return; IRasterLayer pRasterLayer = mapControlMain.ActiveView.CurrentLayer as IRasterLayer; if (pRasterLayer == null)return; IRasterLayer memLayer = OpenMemRasterLayer(pRasterLayer, 0, 0); if(memLayer!=null) { mapControlMain.FocusMap.AddLayer(memLayer as ILayer); mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); } /// <summary> /// 自定义栅格数据读取（通过内存RasterData构建RasterLayer） /// </summary> /// <param name="rasterLayer">栅格数据路径</param> /// <param name="startX">读取起始位置X</param> /// <param name="startY">读取起始位置Y</param> /// <param name="?"></param> /// <returns></returns> private IRasterLayer OpenMemRasterLayer(IRasterLayer rasterLayer, int startX, int startY) { IRasterLayer memRasterLayer = null; object array = null;//存储读取出的数据 //栅格图层参数读取 IRasterDataset rasterDataset = rasterLayer.Dataset; int bandCount = rasterDataset.GetBandCount(); PixelDataType rasterDataType = rasterDataset.GetRasterBand(0).GetRasterDataType(); //读取栅格数据块大小 int width = 2048; int height = 1024; switch (rasterDataType) { case PixelDataType.UInt16: array = new UInt16[width * height * bandCount]; break; case PixelDataType.Byte: array = new Byte[width * height * bandCount]; break; default: break; } int[] bandMap = new int[bandCount]; for (int i = 0; i < bandCount; i++) { bandMap[i] = i + 1; } //将pRasterLayer中的数据读取至数组中 bool oK = rasterDataset.Read(startX, startY, width, height, array, width, height, rasterDataType, bandCount, bandMap); //栅格读取 double[] geoTrans = rasterDataset.GetGeoTransform(); //创建内存数据集 IRasterDataset newRDataset = DatasetFactory.CreateRasterDataset("", width, height, bandCount, rasterDataType, "MEM", null); //数组中的数据写入内存数据集中 oK = newRDataset.Write(0, 0, width, height, array, width, height, rasterDataType, bandCount, bandMap); //栅格写入 //设置参数 newRDataset.SetGeoTransform(geoTrans);//设置仿射变换参数 newRDataset.SpatialReference = rasterDataset.SpatialReference;//设置空间参考 //创建RasterLayer memRasterLayer = PIE.Carto.LayerFactory.CreateDefaultRasterLayer(newRDataset); return memRasterLayer; }

5. 示例截图