很多遥感卫星数据使用的量化级别都要比8bit高，比如常用的WorldView用的是12bit的量化，对于一些图像处理软件，是不能直接处理12bit量化的图像，所以需要将12bit的数据转为16bit数据或者8bit数据来进行处理。

下面写了一个简单的函数来进行处理，具体原理很简单，就是使用GDAL将12bit的数据读进来，然后再使用线性拉伸为8bit存出去，或者直接保存为16bit数据。注意12bit的数据在GDAL中读取的时候会显示为16bit数据，就好比2bit的数据在GDAL中是8bit一样，因为在C或者C++中很难找到一个类型来表示2bit或者12bit的东西，最小的char是8bit，short是16bit。代码如下，首先是头文件：

/***************************************************************************
*
* Time: 2012-02-23
* Project: 遥感平台
* Purpose: 将12bit数据转换为8bit或者16bit
* Author:  李民录
* Copyright (c) 2012, liminlu0314@gmail.com
* Describe:将12bit数据转换为8bit或者16bit
*
****************************************************************************/
#ifndef DATARESCALE_H
#define DATARESCALE_H/*! 8U */
typedef unsigned char   DT_8U;
/*! 16U */
typedef unsigned short  DT_16U;
/**
* @brief 释放数组
*/
#define RELEASE(x)  if(x!=NULL) {delete []x; x = NULL;}/**
* @brief 图像转换，将图像存为16bit，前提确保输入的数据是12bit的
* @param pszSrcFile        输入文件路径
* @param pszDstFile            输出文件路径
* @param bTo8          是否转为8位，false为专为16bit数据，true表示转为8bit数据
* @param pszFormat         输出文件格式，详细参考GDAL支持数据类型
* @return 返回值，表示计算过程中出现的各种错误信息
*/
int ImageDataRescale216(const char* pszSrcFile, const char* pszDstFile, bool bTo8 = true, const char* pszFormat = "GTiff");#endif /* DATARESCALE_H */

下面是函数实现代码：

/***************************************************************************
*
* Time: 2012-02-23
* Project: 遥感平台
* Purpose: 将12bit数据转换为8bit或者16bit
* Author:  李民录
* Copyright (c) 2012, liminlu0314@gmail.com
* Describe:将12bit数据转换为8bit或者16bit
*
****************************************************************************/
#include "DataRescale.h"
#include "gdal_priv.h"/**
* @brief 图像转换，将图像存为16bit，前提确保输入的数据是12bit的
* @param pszSrcFile        输入文件路径
* @param pszDstFile            输出文件路径
* @param bTo8          是否转为8位，false为专为16bit数据，true表示转为8bit数据
* @param pszFormat         输出文件格式，详细参考GDAL支持数据类型
* @return 返回值，表示计算过程中出现的各种错误信息
*/
int ImageDataRescale(const char* pszSrcFile, const char* pszDstFile, bool bTo8 = true, const char* pszFormat = "GTiff")
{//判断输入路径是否为空if(pszSrcFile == NULL || pszDstFile == NULL)return -1;GDALAllRegister();GDALDataset *poSrcDS = (GDALDataset *) GDALOpen( pszSrcFile, GA_ReadOnly );if( poSrcDS == NULL ){//图像打开失败return -2;}GDALDriver *poDriver = GetGDALDriverManager()->GetDriverByName(pszFormat);if( poDriver == NULL ){//不能创建制定类型的文件，请检查该文件类型GDAL是否支持创建GDALClose( (GDALDatasetH) poSrcDS );return -3;}//获取图像宽高和波段数int iXSize = poSrcDS->GetRasterXSize();int iYSize = poSrcDS->GetRasterYSize();int iBandCount = poSrcDS->GetRasterCount();//确定输出图像的位数GDALDataType eDT = GDT_UInt16;if(bTo8)eDT = GDT_Byte;elseeDT = GDT_UInt16;//创建16bit的数据GDALDataset *poDstDS = poDriver->Create(pszDstFile, iXSize, iYSize, iBandCount, eDT, NULL);double dGeoTrans[6] = {0};//设置仿射变换参数poSrcDS->GetGeoTransform(dGeoTrans);poDstDS->SetGeoTransform(dGeoTrans);//设置图像投影信息poDstDS->SetProjection(poSrcDS->GetProjectionRef());//用于保存读取的12bit数据DT_16U *pSrcData = new DT_16U[iXSize];if(bTo8) //转换为8bit{//定义结果数据存储空间DT_8U *pDstData = new DT_8U[iXSize];//循环波段for(int iBand=1; iBand<=iBandCount; iBand++){GDALRasterBand *pSrcBand = poSrcDS->GetRasterBand(iBand);GDALRasterBand *pDstBand = poDstDS->GetRasterBand(iBand);for(int i=0; i<iYSize; i++)    //循环图像高{//将数据读出来pSrcBand->RasterIO(GF_Read, 0, i, iXSize, 1, pSrcData, iXSize, 1, GDT_UInt16, 0, 0);//循环，将12bit数据专为8bit数据，使用线性拉伸方式for(int j=0; j<iXSize; j++){double dTemp = pSrcData[j] * 255.0 / 4095.0;if(dTemp > 255.0)pDstData[j] = 255;elsepDstData[j] = (DT_8U)dTemp;}pDstBand->RasterIO(GF_Write, 0, i, iXSize, 1, pDstData, iXSize, 1, GDT_Byte, 0, 0);}}RELEASE(pDstData);    //释放内存}else //转换为16bit数据{//循环波段for(int iBand=1; iBand<=iBandCount; iBand++){GDALRasterBand *pSrcBand = poSrcDS->GetRasterBand(iBand);GDALRasterBand *pDstBand = poDstDS->GetRasterBand(iBand);for(int i=0; i<iYSize; i++)  //循环图像高{//将数据读出来，然后写入结果数据pSrcBand->RasterIO(GF_Read, 0, i, iXSize, 1, pSrcData, iXSize, 1, GDT_UInt16, 0, 0);pDstBand->RasterIO(GF_Write, 0, i, iXSize, 1, pSrcData, iXSize, 1, GDT_UInt16, 0, 0);}}}RELEASE(pSrcData);//关闭原始图像和结果图像GDALClose( (GDALDatasetH) poDstDS );GDALClose( (GDALDatasetH) poSrcDS );return 0;
}

对于上面的实现做一个简单的说明，12bit的数据读进来，对于16bit的直接写到结果图像里面，没有拉伸到16bit的范围，这样就是完全保留了原始数据的所有图像信息，对于保存为8bit的数据，肯定会造成部分信息的丢失，使用最简单的线性方程进行拉伸，将0～4095的范围拉伸到0～255的范围。