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
示例代码
/// <summary> /// 波段运算算法测试,本算法实现了将World.tif的波段3和波段2波作波段差得到world3.tif /// </summary> public void BandOper() { #region 1、参数设置 PIE.CommonAlgo.BandOper_Exchange_Info info = new PIE.CommonAlgo.BandOper_Exchange_Info(); info.StrExp = "b2-b1"; info.SelectFileBands = new List<int> { , };//band3和band2 根据运算公式的波段大小先后顺序确定 b1的选择波段就是3； info.SelectFileNames = new List<string> { @"D:\data\China1\world\World.tif", @"D:\data\China1\world\World.tif" };//分别为band3和band2数据路径 info.OutputFilePath = @"D:\data\world3.tif"; info.FileTypeCode = "GTiff"; PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.BandOperAlgo"); if (algo == null) return; #endregion //2、算法执行 PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents; algo.Name = "波段运算"; algo.Params = info; //3、结果显示 bool result = PIE.SystemAlgo.AlgoFactory.Instance().ExecuteAlgo(algo); if (result) { MessageBox.Show("波段算法执行成功"); ILayer layer = LayerFactory.CreateDefaultLayer(info.OutputFilePath); if (layer == null) return; m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); } }

2.4.示例截图