PIE SDK波段运算

 

1.算法功能简介

    波段运算(Band Math)工具能够方便的执行图像中的各个波段的加减乘除、三角函数、指数、对数等数学函数计算,也可以使用IDL编写的函数。

    由于每个用户都有独特的需求,利用此工具用户可以自己定义处理算法,应用到某个波段或者整个图像中。波段运算实质上是对每个像素点对应的的像素值进行数学运算,运算表达式中的每一个变量不仅可以对应于单一波段,也可以是一个多波段的栅格文件。例如,在表达式 b1﹢b2 中,如果 b1 是一个多波段的图像, b2 为单一波段,则结果为 b1 所对应图像的所有波段分别与 b2 进行求和。

    PIE SDK支持算法功能的执行,下面对波段运算算法功能进行介绍。

2.算法功能实现说明

2.1 实现步骤

第一步

算法参数设置

第二步

算法执行

第三步

结果显示

2.2 算法参数

C#算法DLL

PIE.CommonAlgo.dll

C#算法名称

PIE.CommonAlgo.BandOperAlgo

参数结构体

BandOper_Exchange_Info

参数说明

StrExp

String

波段运算公式,如"(b4-b3)/(b4+b3)";

SelectFileBands

IList<Int>

选中的图像对应的波段band编号(是波段编号而非波段索引),根据波段运算公式的波段大小先后顺序b3、b4确定顺序。new List<int> {3, 4}

SelectFileNames

IList<String>

选中的图像对应的文件file名称,根据波段运算公式的波段大小先后顺序b3、b4确定顺序。new List<String> {b3对应文件路径, b4对应文件路径}

OutputFilePath

String

输出文件路径

FileTypeCode

String

根据输出类型获得文件编码类型

.tif/.tiff——GTiff

.img—————HFA

其他—————ENVI

FuncName

String

C#算法名称

PixelDataType

int(针对不同的类型设置对应的编号)

输出数据类型:

Byte8——1(编号)(字节型)

UInt16——2   (整型16位)

Int16——3      (无符号整型16位)

UInt32——4  (长整型32位)

Int32——5     (无符号长整型32位)

Float32——6 (浮点型32位)

Float64——7  (双精度浮点型64位)

2.3. 示例代码

项目路径

百度云盘地址下/PIE示例程序/10.算法调用/多功能工具/FundamentalToolDemo.BandOperDemo

数据路径

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

视频路径

百度云盘地址下/PIE视频教程/10.算法调用/多功能工具/波段运算算法.avi

示例代码

 1 /// <summary> 
 2 ///  波段运算算法测试,本算法实现了将World.tif的波段3和波段2波作波段差得到world3.tif
 3 /// </summary>
 4 public void BandOper()
 5 {
 6     #region 1、参数设置
 7     PIE.CommonAlgo.BandOper_Exchange_Info info = new PIE.CommonAlgo.BandOper_Exchange_Info();
 8     info.StrExp = "b2-b1";
 9     info.SelectFileBands = new List<int> { 3,2 };//band3和band2 根据运算公式的波段大小先后顺序确定 b1的选择波段就是3;
10     info.SelectFileNames = new List<string> { @"D:dataChina1worldWorld.tif", @"D:dataChina1worldWorld.tif" };//分别为band3和band2数据路径
11     info.OutputFilePath = @"D:dataworld3.tif";
12     info.FileTypeCode = "GTiff";
13 
14     PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.BandOperAlgo");
15     if (algo == null) return;
16     #endregion
17 
18     //2、算法执行
19     PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents;
20     algo.Name = "波段运算";
21     algo.Params = info;       
22 
23     //3、结果显示
24     bool result = PIE.SystemAlgo.AlgoFactory.Instance().ExecuteAlgo(algo);
25     if (result)
26     {
27         MessageBox.Show("波段算法执行成功");
28         ILayer layer = LayerFactory.CreateDefaultLayer(info.OutputFilePath);
29         if (layer == null) return;
30         m_HookHelper.ActiveView.FocusMap.AddLayer(layer);
31         m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);
32     }           
33

2.4.示例截图

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