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 |
示例代码 | |
1 方法(一)打开自定义栅格数据 2 if (mapControlMain.FocusMap.LayerCount == 0)return; 3 IRasterLayer pRasterLayer = mapControlMain.ActiveView.CurrentLayer as IRasterLayer; 4 if (pRasterLayer == null)return; 5 IRasterLayer memLayer = OpenMemRasterLayer(pRasterLayer, 0, 0); 6 if(memLayer!=null) 7 { 8 mapControlMain.FocusMap.AddLayer(memLayer as ILayer); 9 mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); 10 } 11 12 /// <summary> 13 /// 自定义栅格数据读取(通过内存RasterData构建RasterLayer) 14 /// </summary> 15 /// <param name="rasterLayer">栅格数据路径</param> 16 /// <param name="startX">读取起始位置X</param> 17 /// <param name="startY">读取起始位置Y</param> 18 /// <param name="?"></param> 19 /// <returns></returns> 20 private IRasterLayer OpenMemRasterLayer(IRasterLayer rasterLayer, int startX, int startY) 21 { 22 IRasterLayer memRasterLayer = null; 23 object array = null;//存储读取出的数据 24 //栅格图层参数读取 25 IRasterDataset rasterDataset = rasterLayer.Dataset; 26 int bandCount = rasterDataset.GetBandCount(); 27 PixelDataType rasterDataType = rasterDataset.GetRasterBand(0).GetRasterDataType(); 28 29 //读取栅格数据块大小 30 int width = 2048; 31 int height = 1024; 32 switch (rasterDataType) 33 { 34 case PixelDataType.UInt16: 35 array = new UInt16[width * height * bandCount]; 36 break; 37 case PixelDataType.Byte: 38 array = new Byte[width * height * bandCount]; 39 break; 40 default: 41 break; 42 } 43 int[] bandMap = new int[bandCount]; 44 for (int i = 0; i < bandCount; i++) 45 { 46 bandMap[i] = i + 1; 47 } 48 //将pRasterLayer中的数据读取至数组中 49 bool oK = rasterDataset.Read(startX, startY, width, height, array, width, height, rasterDataType, bandCount, bandMap); //栅格读取 50 double[] geoTrans = rasterDataset.GetGeoTransform(); 51 52 //创建内存数据集 53 IRasterDataset newRDataset = DatasetFactory.CreateRasterDataset("", width, height, bandCount, rasterDataType, "MEM", null); 54 //数组中的数据写入内存数据集中 55 oK = newRDataset.Write(0, 0, width, height, array, width, height, rasterDataType, bandCount, bandMap); //栅格写入 56 //设置参数 57 newRDataset.SetGeoTransform(geoTrans);//设置仿射变换参数 58 newRDataset.SpatialReference = rasterDataset.SpatialReference;//设置空间参考 59 60 //创建RasterLayer 61 memRasterLayer = PIE.Carto.LayerFactory.CreateDefaultRasterLayer(newRDataset); 62 return memRasterLayer; 63 } |
5. 示例截图