PIE SDK栅格数据的创建

 

1. 功能简介

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

2. 功能实现说明

2.1. 实现思路及原理说明

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

2.2. 核心接口与方法

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

2.3. 示例代码

项目路径	百度云盘地址下/PIE示例程序/04数据操作/10栅格数据的创建
数据路径	百度云盘地址下/PIE示例数据/栅格数据/04.World/World.tif
视频路径	百度云盘地址下/PIE视频教程/04数据操作/10栅格数据的创建.avi
示例代码
1 /// <summary> 2 /// 创建栅格数据 3 /// </summary> 4 /// <param name="sender"></param> 5 /// <param name="e"></param> 6 private void toolStripButton_CreateRasterData_Click(object sender, EventArgs e) 7 { 8 //1获取当前的栅格图层 9 //可以自定义打开一个tiff栅格数据DatasetFactory.OpenRasterDataset(strPath,OpenMode.ReadOnly); 10 IRasterDataset rasterDataset = (mapControlMain.FocusMap.GetLayer(0) as IRasterLayer).Dataset; 11 int BandCount = rasterDataset.GetBandCount(); 12 //2栅格数据范围等信息的读取 13 IEnvelope envelop = rasterDataset.GetExtent(); 14 int xSize = rasterDataset.GetRasterXSize();//行 15 int ySize = rasterDataset.GetRasterYSize();//列 16 UInt16[] array = new UInt16[xSize * ySize * BandCount]; 17 int[] bandMap = new int[BandCount]; 18 for (int i = 0; i < BandCount; i++) 19 { 20 bandMap[i] = i + 1; 21 } 22 SaveFileDialog saveFileDialog = new SaveFileDialog(); 23 saveFileDialog.Title = "请选择要保存的路径："; 24 saveFileDialog.Filter = "Raster File|*.tiff"; 25 if (saveFileDialog.ShowDialog() != DialogResult.OK) return; 26 //3栅格数据读取 27 bool OK = rasterDataset.Read(0, 0, xSize, ySize, array, xSize, ySize, PixelDataType.UInt16, BandCount, bandMap); 28 //4构建栅格数据集 29 IRasterDataset newRasterDataset = DatasetFactory.CreateRasterDataset(saveFileDialog.FileName, xSize, ySize, BandCount, PixelDataType.UInt16, "GTiff", null); 30 //5设设置空间参考 31 newRasterDataset.SpatialReference = rasterDataset.SpatialReference; 32 //6栅格数据的写入 33 bool OK1 = newRasterDataset.Write(0, 0, xSize, ySize, array, xSize, ySize, PixelDataType.UInt16, BandCount, bandMap); 34 //7设置仿射变换 35 double[] geoTrans = rasterDataset.GetGeoTransform(); 36 newRasterDataset.SetGeoTransform(geoTrans); 37 //8创建图层添加到map中 38 ILayer layer = LayerFactory.CreateDefaultRasterLayer(newRasterDataset) as ILayer; 39 //第一种添加数据的方法 40 mapControlMain.FocusMap.AddLayer(layer, 0); 41 //第二种添加数据方法 42 // mapControl1.AddLayerFromFile(saveFileDialog.FileName, 0); 43 mapControlMain.PartialRefresh(ViewDrawPhaseType.ViewAll); 44 } View Code

2.4. 示例截图

转载于:https://www.cnblogs.com/PIESat/p/10168847.html