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#算法名称 |
2.3示例代码
项目路径 |
百度云盘地址下/PIE示例程序/ FundamentalToolDemo. BandOperDemo |
数据路径 |
百度云盘地址下/ PIE示例数据/栅格数/World/world.tiff |
视频路径 |
百度云盘地址下/PIE视频教程/波段运算算法avi |
示例代码 |
|
1 /// <summary> 2 ///波段运算算法测试,本算法实现了将World.tif的第一波段生成World3.tif文件 3 /// </summary> 4 private void Test_KrigingInterpolationAlgo() 5 { 6 #region 1、参数设置 7 PIE.CommonAlgo.BandOper_Exchange_Info info = new PIE.CommonAlgo.BandOper_Exchange_Info(); 8 info.StrExp = "b1"; 9 info.SelectFileBands = new List<int> { 1 }; 10 info.SelectFileNames = new List<string> { @"D:\Data\World.tif", @"D:\Data\World.tif" }; 11 info.OutputFilePath = @"D:\Data\World3.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; PIE.SystemAlgo.AlgoFactory.Instance().ExecuteAlgo(algo); 22 23 //3、结果显示 24 ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(@"D:\Data\World3.tif"); 25 m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); 26 } |