pix格式的摸索(二)

作者：朱金灿

来源：http://blog.csdn.net/clever101

PCI的系统格式pix是一个设计很巧妙的遥感图像格式，而且其设计巧妙之处不止一处两处，这些都有待我日后一一去摸索。今天就具体谈谈它的设计巧妙的一个地方——它是认真地考虑图像的大小端字节存储的问题并作出了针对性的设计。

关于计算机处理器的大小端字节序存储的基础知识不在本文的论述范围之列，有兴趣的读者可以阅读《深入理解计算机系统》这本书。首先要讨论的问题是图像作为一种二进制文件格式，是否有大小端存储之分？我认为是有的。可能有些读者对此有疑问：大小端不是对处理器而言的吗？图像怎么有大小端存储之分的呢？实际上是有的。我们想象这样一个场景：假如一个图像是在一个大端机器创建并处理了，毫无疑问它是按大端的字节序存储的，现在用户把它拷贝到一个小端机器上浏览，软件不知道它是按大端存储还是按小端存储的，就想当然地按小端字节序来浏览，自然就看到一堆乱码而不是一个预想中的有序图像。这个问题其实是一直存在的，为什么我们几乎很少碰到？一是我们日常接触的都是8位图像（8位图像不存在字节序存储问题，大于8位的图像在遥感领域很常见），二则很可能是因为PC领域一直是由intel的小端模式的x86架构所统治。

现在问题来了：如何解决这个问题？显然要根本解决问题是图像文件要提供记忆自己是哪种字节序来存储的功能。pix文件正是基于这种考虑，在每个波段头都存储了它是按哪种字节序存储的标记（绝妙的是它是每个波段都有这个标记，这就意味着它的每个波段是可以按照不同的字节序存储的）。现在我们看看pix文件是如何具体实现的。这个我们从PCI公司提供的读写pix文件的开源库中找到答案。首先源码中的PCIDSK::Create函数是负责创建pix文件的。在pix文件中每个波段头都有1024字节用于存储波段的描述信息，其中它的第201个字节就是用于保存是大端字节序存储还是小端字节序存储的，设为'S'为小端存储，设为'N'为大端存储。注意要直接指定ih.buffer[201]='S';而不能使用它的设置函数，如ih.Put('S',201,1);Put函数是一个坑爹的函数！

在读取和写入pix文件是这样处理的：

unsigned short test_value = 1;
// byte_order保存的是图像是按哪种字节序存储
//假如是小端机器，同时又是小端存储，那么就不需要交换高低位的数//据，否则需要交换
if( ((uint8 *) &test_value)[0] == 1 )
	needs_swap = (byte_order != 'S');
else  //假如是大端机器，办法和小端一样
     needs_swap = (byte_order == 'S');

交换的做法很简单，以交换32位浮点数为例，具体如下：

   BYTE  *pbBuffer = NULL;
   BYTE  *pbtPixel = NULL;
float fValue=(static_cast<float*>(pvBuffer))[i];
					   pbtPixel=(BYTE*)(&fValue);
					   *pbBuffer = pbtPixel[3];
					   pbBuffer++;
					   *pbBuffer = pbtPixel[2];
					   pbBuffer++;
					   *pbBuffer = pbtPixel[1];
					   pbBuffer++;
					   *pbBuffer = pbtPixel[0];
					   pbBuffer++;