PIE SDK栅格数据集的读写

 

1. 功能简介

    栅格数据包含很多信息，在数据的运用中需要对数据的信息进行读取或写入，目前PIE SDK支持多种数据格式的数据读取和写入，下面对栅格数据格式的数据读写功能进行介绍。

2. 功能实现说明

2.1. 实现思路及原理说明

第一步

获取要读取的栅格数据

第二步

读取栅格数据

第三步

写入栅格数据并赋值投影

2.2. 核心接口与方法

接口/类

方法/属性

说明

 

 

IRasterDataset

GetBandCount()

获取波段数

GetRasterBand(int nIndex)

获取栅格波段对象

GetRasterDataType()

获取栅格数据类型

Read(int nx, int ny, int nWid, int nHei, int nBufXSize, int nBufYSize, IList<int> bandMap)

栅格数据集读取数据

2.3. 示例代码

项目路径

百度云盘地址下/PIE示例程序/04数据操作/12栅格数据集的读写

数据路径

百度云盘地址下/PIE示例数据/栅格数据/04.World/World.tif

视频路径

百度云盘地址下/PIE视频教程/04数据操作/12栅格数据集的读写.avi

示例代码

 
 1          /// <summary>
 2         /// 栅格数据集的读写
 3         /// </summary>
 4         /// <param name="sender"></param>
 5         /// <param name="e"></param>
 6         private void toolStripButton_ReadRasterData_Click(object sender, EventArgs e)
 7         {
 8             string rPathLBT = @"D:\Data\World.tif";
 9             IRasterDataset rDataset = DatasetFactory.OpenRasterDataset(rPathLBT, OpenMode.ReadOnly);
10             int bandCount = rDataset.GetBandCount();
11             IRasterBand rasterband = rDataset.GetRasterBand(bandCount - 1);
12             PixelDataType type = rasterband.GetRasterDataType();
13 
14             int[] bandMap = new int[bandCount];
15             for (int i = 0; i < bandCount; i++)
16             {
17                 bandMap[i] = i + 1;
18             }
19             int nWidth = 500;
20             int nHeight = 500;
21             UInt16[] buf = new UInt16[nWidtn * nHight * bandCount];
22             bool ok = rDataset.Read(0, 0, 500, 500, buf, 500, 500, type, bandCount, bandMap);
23 
24             //获取投影
25             ISpatialReference spGLL = rDatasetGLL.SpatialReference;
26 
27             //写入数据并赋值投影
29             string rPathSaveGLL = @"D:\Data\World_new.tif";
30             IRasterDataset rasterDatasetSave = DatasetFactory.CreateRasterDataset(rPathSaveGLL, nWidth, nHeight, bandCount, type, "GTIFF", null);
31             rasterDatasetSave.SpatialReference = spGLL;
32             ok = rasterDatasetSave.Write(0, 0, nWidth, nHeight, buf, nWidth, nHeight, type, bandCount, bandMap);
33 
33             double[] geoTrans = rDataset.GetGeoTransform();
34             rasterDatasetSave.SetGeoTransform(geoTrans);
35             ((IDisposable)rasterDatasetSave).Dispose();
36             rasterDatasetSave = null;
37             ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(rPathSaveGLL);
38             mapControlMain.FocusMap.AddLayer(layer);            mapControlMain.ActiveView.PartialRefresh(PIE.Carto.ViewDrawPhaseType.ViewAll);
39         }
View Code

2.4. 示例截图

 

 

posted @ 2018-12-24 16:07  PIESAT  阅读(2560)  评论(0)  编辑  收藏  举报