PIESDK二次开发——栅格数据的读取与保存

PIESDK二次开发——栅格数据的读取与保存

1、栅格数据的读取

• 使用IMapControl接口，

// 函数原型：
bool AddLayerFromFile(String filePath,int toIndex); // 在当前地图显示控件中通过文件路径添加图层，参数为目标文件路径及添加到的索引值。

// 示例：
IMapControl mapControl = null;
mapControl = mapControl1;
mapControl.AddLayerFromFile("aaa.tif",0)

• 使用PIE.DataSource.DatasetFactory类

//函数原型
IDataset OpenDataset(String strPath, OpenMode mode); // PIE.DataSource.OpenMode
//示例
PIE.DataSource.IRasterDataset rDataset = PIE.DataSource.DatasetFactory.OpenDataset("aaa.tif", PIE.DataSource.OpenMode.ReadOnly) as PIE.DataSource.IRasterDataset;

//函数原型
IRasterDataset OpenRasterDataset (String strPath, OpenMode mode); // PIE.DataSource.OpenMode
//示例
PIE.DataSource.IRasterDataset rDataset = PIE.DataSource.DatasetFactory.OpenRasterDataset("aaa.tif", PIE.DataSource.OpenMode.ReadOnly)

• 使用 PIE.Carto.LayerFactory类

// 函数原型：
public static ILayer CreateDefaultLayer(string strFile)
//示例：
 PIE.Carto.ILayer layer = null;
 layer = PIE.Carto.LayerFactory.CreateDefaultLayer("aaa.tif");

2、PIE.Carto.ILayer与PIE.DataSource.IRasterDataset之间的转换

从前面读取栅格数据的几个方法中，我们发现有的方法读取出来返回的结果是PIE.Carto.ILayer类型的，有的是PIE.DataSource.IRasterDataset类型的，那么这两种类型直接如何转换呢？

• PIE.Carto.ILayer 到 PIE.DataSource.IRasterDataset的转化

//方法1
IMapControl mapControl = null;
mapControl = mapControl1;
mapControl.AddLayerFromFile("aaa.tif",0)
IRasterDataset rasterDataset = (mapControl.FocusMap.GetLayer(0) as IRasterLayer).Dataset;
//方法2
PIE.Carto.ILayer layer = null;
layer = PIE.Carto.LayerFactory.CreateDefaultLayer("aaa.tif");
IRasterDataset rasterDataset = (layer  as IRasterLayer).Dataset;

• PIE.DataSource.IRasterDataset 到 PIE.Carto.ILayer的转化

//使用 PIE.Carto.LayerFactory.CreateDefaultRasterLayer方法
//函数原型：
public static IRasterLayer CreateDefaultRasterLayer(
	IRasterDataset ptrDataset
)

// 示例代码：
PIE.Carto.ILayer layer = PIE.Carto.LayerFactory.CreateDefaultRasterLayer(info.OutputRasterDataset) as PIE.Carto.ILayer;
if (m_HookHelper.FocusMap == null) AddNewMapPage(layer.Name);
m_HookHelper.FocusMap.AddLayer(layer);
m_HookHelper.ActiveView.PartialRefresh(Carto.ViewDrawPhaseType.ViewAll);

3、将栅格数据集保存为tif文件

// 我的函数：
     bool cml_saveRasterToFile(PIE.DataSource.IRasterDataset OutputRasterDataset,string outPutFile)
        {
            // 获取栅格数据范围信息
            PIE.Geometry.IEnvelope envelop = OutputRasterDataset.GetExtent();
            int xSize = OutputRasterDataset.GetRasterXSize();//行
            int ySize = OutputRasterDataset.GetRasterYSize();//列
                                                              // 	获取栅格数据集波段数量
            int BandCount = OutputRasterDataset.GetBandCount();
            // 获取数据类型
            PIE.DataSource.IRasterBand rasterband = OutputRasterDataset.GetRasterBand(BandCount - 1);
            PIE.DataSource.PixelDataType type = rasterband.GetRasterDataType();
            UInt16[] array = new UInt16[xSize * ySize * BandCount];
            int[] bandMap = new int[BandCount];
            for (int i = 0; i < BandCount; i++)
            {
                bandMap[i] = i + 1;
            }
            // 读取数据
            bool OK = OutputRasterDataset.Read(0, 0, xSize, ySize, array, xSize, ySize, type, BandCount, bandMap);
            if (!OK) return false;
            // 构建新的栅格数据集
            PIE.DataSource.IRasterDataset newRasterDataset = PIE.DataSource.DatasetFactory.CreateRasterDataset(outPutFile, xSize, ySize, BandCount, type, "GTiff", null);
            // 设设置空间参考
            newRasterDataset.SpatialReference = OutputRasterDataset.SpatialReference;
            // 栅格数据的写入
            bool OK1 = newRasterDataset.Write(0, 0, xSize, ySize, array, xSize, ySize, type, BandCount, bandMap);

            //7设置仿射变换
            //double[] geoTrans = OutputRasterDataset.GetGeoTransform();
            //newRasterDataset.SetGeoTransform(geoTrans);
            return OK1;
        }