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:DataWorld.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             UInt16[] buf = new UInt16[500 * 500 * bandCount];
20             bool ok = rDataset.Read(0, 0, 500, 500, buf, 500, 500, type, bandCount, bandMap);
21 
22             //获取投影
23             string rPathGLL = @"D:DataWorld.tif";
24             IRasterDataset rDatasetGLL = DatasetFactory.OpenRasterDataset(rPathGLL, OpenMode.ReadOnly);
25             ISpatialReference spGLL = rDatasetGLL.SpatialReference;
26             string stGLL = spGLL.ExportToWkt();
27 
28             //写入数据并赋值投影
29             string rPathSaveGLL = @"D:DataWorld_new.tif";
30             IRasterDataset rasterDatasetSave = DatasetFactory.CreateRasterDataset(rPathSaveGLL, 500, 500, bandCount, type, "GTIFF", null);
31             ok = rasterDatasetSave.Write(0, 0, 500, 500, buf, 500, 500, type, bandCount, bandMap);
32             rasterDatasetSave.SpatialReference = SpatialReferenceFactory.CreateSpatialReference(stGLL);
33 
34             ((IDisposable)rasterDatasetSave).Dispose();
35             rasterDatasetSave = null;
36             ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(rPathSaveGLL);
37             mapControlMain.FocusMap.AddLayer(layer);            mapControlMain.ActiveView.PartialRefresh(PIE.Carto.ViewDrawPhaseType.ViewAll);
38         } 
View Code

2.4. 示例截图

 

原文地址:https://www.cnblogs.com/PIESat/p/10169156.html