BSQ格式数据转换为RSD缺省的BIP格式数据_tiff格式是bsq吗-优快云博客

李国春

RSD内部统一以BIP格式排列数据，并且文件格式（非TFS）数据倒放（North Down）。早期是为了和设备无关位图（DIB）一致节省一点处理时间。现在设备处理能力增强了这点时间已经无关紧要，但是这种数据排列方式还是延续了下来。通常这些内部结构不需要用户关心，RSD自动处理成习惯的数据和图像。

高分卫星遥感数据运营部门分发的1级或者2级TIFF格式数据内部多保存为BSQ格式。RSD在读取时自动转换成BIP格式。但是也有一些部门（软件）开始就将其保存为BIP，RSD读取时再进行转换反而成了BSQ，读出来正射后结果如图1（左），4个波段的数据变成了16个小块。

下面介绍一个函数将其纠正为正常显示方式，图1（右）。

图1

下面的函数是一个读TIFF数据到缓存的函数。

int TifReadData(STRING name,WORD* buffer,int interleave);

这个函数将name 指定的TIFF数据以BIP格式读入到数据缓存buffer中。第三个参数int interleave是文件中的数据排列格式。如果你知道TIFF文件是BSQ格式排列的，将这个参数指定为BSQ,则读出数据如图1（右）。

下面是转换一个TIFF格式数据的完整代码。如果你正射校正一个数据集后发现结果如图一（左），将这个正射结果保存为一个TIFF格式文件，再用下述代码处理。代码如下：

/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
//	99脚本语言示例代码  —— TifReadData(name,buffer,BSQ)，将用BIP（INTERLEAVE==PIX）格式数据按 BSQ格式读出  
//
//	李国春 2022年 11月 6日
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
main()
{
	STRING name = OpenFileDialog(TRUE,"*.tif;*.tiff");//打开一个tif文件

	Print("开始打开 %s 文件，读BSQ格式的数据",name);
	DWORD tc = GetTickCount();

	STRING projstr;
	double box[4];		//投影后的四角坐标，次序为：minN，minE，maxN，maxE
	double quad[8];		//四角的经纬度，书写次序。
	double pixsize;
	int height,width,bands,datatype;
	int b2 = TifInqGTiffInfo(name,projstr,box,quad,pixsize,height,width,bands,datatype);	//读取 Geo TIFF 文件信息
	
	//申请内存，读入TIFF数据///////////////////////////////////////////////////////////////////////////////////////////////
	WORD buffer[height][width][bands];
	int e = TifReadData(name,buffer,BSQ);//最后一个参数进行 BSQ 转换。（读正常数据最后一个参数写 BIP，或者空着不写）

	//根据 PAN 数据范围创建任务，向任务框架添加正射校正后的 PAN 和 MSS 数据////////////////////////////////////////////////
	STRING misname = FileNameMain(name);
	int misID = CreateMission2(misname,projstr,box,pixsize);//创建任务，但是不显示窗口！！！
	FmShowFrameWindow(misID,TRUE);//显示主框架窗口

	//转换后数据创建新层
	int lyr = GdMakeBufferGridLayer(buffer,datatype,0,0,height,width,bands,BIP,0,32767,misname);//将栅格数据加载到任务框架
	int lyrs = InqNumLayers();//任务中层数

	//窗口刷新显示
	FmUpdateFrameWindow(lyrs-1);//用最后一层数据刷新任务主窗口

	double tm = (GetTickCount() - tc) / 1000.0;
	Print("转换并创建RSD任务 %s 并合成图像用时 %f 秒。",name,tm);

	Print("主程序结束");
	return TRUE;
}

运行后得到图1右侧的结果。

或者选中一个BSQ的GRID层，使用图2的菜单命令也可以将BSQ格式的数据转换为RSD缺省的BIP格式。