编程实现BIP、BSQ、BIL三种格式转换范例.doc

编程实现BIP、BSQ、BIL三种格式转换 实验目的 理解遥感图像BIP、BIL、BSQ三种格式数据的组织方式，以及它们互相转换原理和方法。 实验原理 ENVI栅格图像文件，是以字节数据为单位，再按照指定顺序组织、排列而成，具体有如下三种方式： 1.BIP(像元波段交叉式):以一对像元为基本单位进行记录 ; 2.BIL(波段按行交叉格式 图像按顺序第一个像元所有的波段，接着是第二个像元的所有波段，然后是第三个像元的所有波段，等等直到像元总数为止 3.BSQ(波段顺序格式 每行数据后面紧接着同一波谱波段的下一行数据 实验步骤 新建项目文件：本实验选用的语言是C#，开发工具是VisualStudio2010，通过“文件—新建—项目—C#—Windows窗体应用程序”,命名“格式转换”即可； 编写代码：通过实验思路的流程图，我们确定代码也应该分这三个功能来写，对应可以组织为三个大方法，再细分，调用即可。 (1).编写读取头文件代码：ENVI中栅格图像被分为两个文件存储，一个“.HDR”头文件和一个数据文件，头文件存储了图像描述信息，它是一个文本文件，我们可以用System.IO.SreamReader字符读取流按行读取，再用字符串匹配方法，以读取字段的方式获取图像的行数、列数、波段数、基本数据类型，即可，具体代码可见后文代码附录； (2).编写数据读取数据文件代码：图像数据组织方式有BIP、BIL、BSQ三种，那么对应的也应该有三个方法按照行、列、波段交叉顺序对图像字节数据进行读取，具体代码见后文； (3).编写字节数据重排列存储代码：在读取了图像数据文件后，可以通过FileStream类按照目标类型格式文件组织顺序，把字节数据写入到文件中保存即可，具体代码见后文。 编译与调试：通过用预先知道内部字节数据的实验文件，进行测试，辅助断点调试，可排除一般问题，接着编译生成目标程序； 运行程序查看结果：通过对多种格式、多种基本数据类型、多种像素大小的文件进行多种方式转换，再利用ENVI打开两个图像查看是否一样，再辅以“Link”功能看同一像素点值是否相同，这样可以深层确定转换的正确性。 结果与分析 结果 程序运行界面 转换后结果 分析总结 (1).读取头文件时，不能通过各字段的现有位，按位读取，应该按行读取，再利用字符串匹配方式，找到指定字段，这样避免字段位变化时读取出错误信息； (2).再读取数据文件和写入存储时，应该要以基本数据类型长度个字节为单位进行，不能单纯的以一个字节操作，不然会不具有普遍性，遇到以“int、float”为基本类型的文件时，就会出错； (3).要想节省代码，可先将各种格式数据读取存储为一种格式字节数组中，再从这个格式往其他格式进行转换，这样可以节省不少工作量，但是引入了中间数组，耗费了内存。 源代码附录 1.读取头文件信息 ? ///?summary ??? ?///?读取头文件信息 ???? ///?/summary ???? ///?param?name=strFileName/param头文件路径和名称 ???? ///?param?name=iColumnsCount/param像素列数 ???? ///?param?name=iLinesCount/param像素行数 ???? ///?param?name=iBandsCount/param波段数 ???? ///?param?name=iType/param基本数据类型代码 ???? ///?param?name=strInterLeave/param文件组织格式 ???? ///?returns/returns返回是否读取成功 ???? public?static?bool?ReadHDR(String?strFileName,out爄nt爄ColumnsCount,out爄nt爄LinesCount,out爄nt爄BandsCount,out爄nt爄Type,out燬tring爏trInterLeave) 牋牋{ 牋牋牋燽ool燽lnSuccess?爁alse; 牋牋牋爄ColumnsCount??1; 牋牋牋爄LinesCount??1; 牋牋牋爄BandsCount??1; 牋牋牋爄Type??1; 牋牋牋爏trInterLeave??; 牋牋牋?//初始化各个变量 ???????StreamReader?hdrfile=null; ???????try ???????{ ??????????hdrfile?=?new?StreamReader(