ROI_PAC学习笔记-print.docVIP

ROI_PAC学习笔记(一)----总论ROI_PAC的全称是Repeat Orbit Interferometry Package. 由美国的JPL和Caltech开发. 软件运行于Linux平台, SGI或SUN平台. 是一个开源软件. 可以在openchannelfoundation的网站上申请下载. 相对于动辄数万元的INSAR商业软件而言, 免费的ROI_PAC还是比较有吸引力的. 对于研究而言, 更重要的还在于它的开源. 软件的使用者可以根据自己的需要改编, 或是在其基础上添加新功能. 这就提高了我们在INSAR软件实验上的较率. 虽然ROI_PAC的规模不是很大, 但要想使用它也非得先认真学习一番不可. 本人也是刚开始学习这个软件. 现在博客这样方便就不妨写上来, 如果有同人也有兴趣研究之, 那么一同学习讨论则是最好不过的了. 这个学习笔记的参考资料是由Buckley在2000年的博士论文改编的软件文档和ROI_PAC3.0. 整个软件的总体流程是由Perl脚本控制的. 各专门的功能模块则是用Fortran或C编写的. ??? ?? 处理步骤与模块的对应如下. 处理步骤 模块 流程控制 process.pl SAR原始数据整理 make_raw.pl SAR成像 roi SAR图像配准 ampcor 干涉图生成 resamp_roi 确定基线 baseest 干涉图去平 cecpxcc, rilooks 确定相关性 Makecc; icu 干涉图滤波 icu 相位解缠与绝对相位确定 Icu; baseline 确定形变 diffnsim 从解缠相位确定地形 inverse3d ??? ?? process.pl是总的流程控制脚本. 它控制着从原始数据一直到生成形变图的全过程ROI_PAC学习笔记(二)----Perl基础ROI_PAC是用Perl脚本控制的。笔者以为，了解一个软件应该从顶层入手，一开始先了解大体，然后再逐渐进入细部，而PERL恰恰就是ROI_PAC的最顶层的控制语言。因此在开始之前应该先学一点PERL，以便阅读源码。所谓“工欲善其事，必先利其器”也。 一、变量 有四种变量：标量，数组，哈希变量，句柄。 标量就是整形或浮点数和字符串，这一点类似C的int float double外加一个字符串类型。变量不用声明就可以使用，但是不声明的变量就是全局变量。局部变量在变量前加上my修饰。还有，标量要用$开头。例：$ a = 1; $b = 1.2; my $c = 22;$s=sdfjls; 数组也称列表。创建：@array=(1,2,3);访问元素：$array[0];要注意的是，数组变量用@开头，元素则是用$开头。如果把数组 直接赋值给标量，那么标量实际得到的是数据长度。PERL的数组是变长数组，可以随便访问元素，比如$array[100]=1;这是可以的，数组被自动 扩大了，下标为4－99的元素为空值。 哈希变量就是哈希表。可以用字符串索引元素，开头用%。创建：%hash=(ab=1,cd=2);访问：$hash {ab};注意访问时用的是花括号（{}）。使用时可以随时以赋值的形式插入新元素，如$hash{ef}=3;删除元素则只能用delete ($hash{ef}); 句柄有点像C里的文件指针，但它肯定不是，只是用着有点像。比如以读方式（另有覆盖和追加两种写方式，分别在文件名前加和 ）打开一个方件：open(fhd,test.txt);fhd就是一个文件句柄。然后：$line=fhd;这样就把文件 里的第一行读到$line里了，很简洁吧，注意，一定不能少。STDOUT是屏幕输出的句柄，print STDOUT dfjlsf;就是把字符串输出到屏幕上，一般STDOUT就略掉了，如果是输出到文件那就一定要有文件的句柄了。close(fhd);则是关闭了 句柄。 二、运算符 常用的加减乘除、逻辑与C的用法一样，也有自增自减运算符。另外比C多了一个乘方符，**。2的立方可以写作2**3; 三、函数 函数定义 sub function1 { ?? $arg1 = shift; ?? $arg2 = $_[0]; ??? #statements 　$arg1; } 函数调用 $ret = function1(arg1,arg2); 函数定义没有参数列表。传入的参数存在数组@_里（这个东西比较诡异，不过别慌，Perl里面诡异的系统变量还多着呢。）可以用上面例子里的 shift取出其中的第一个参数并将其从@_中删除，其实这是一个简写，完整的写法是shift(@_);由于@_是系统的特殊变量（数组），在函数的默 认参数就是它。顺便提一下，$_也是一个默认函数参数（标量）