[转]使用GDAL创建Erdas格式的金字塔

   在使用Erdas或者ArcGIS打开栅格图像的时候，会创建一个后缀名为rrd的金字塔文件，用于快速显示图像。那么在使用GDAL编写自己的图像算法后，像快速的在Erdas或者ArcGIS中显示，就需要自己创建rrd格式的金字塔文件，这样在打开该图像文件时，打开速度会非常快，在我的电脑上一个2G的img不到一秒钟可以全部加载进来。    查看GDAL中，有个gdaladdo的工具，就是一个专门用于创建金字塔文件的，但是默认的创建的是一个后缀名为ovr的金字塔文件，该种格式的金字塔不能被Erdas或者ArcGIS使用，但是可以在QGIS等使用。为了能够在Erdas中使用，在创建的时候需要指定一个选项，那就是 USE_RRD=YES，使用该选项后，创建的金字塔格式是一个aux文件，虽然后缀名不是rrd，但是该文件是可以被Erdas或者ArcGIS中使用的。
    关于gdaladdo的使用帮助，可以参考网址：http://www.gdal.org/gdaladdo.html
    Erdas的金字塔是按照2的次方来采样，金字塔顶层的大小应该是小于等于64*64，创建金字塔的，于是按照gdaladdo中的说明，其命令行参数应该是： gdaladdo --config USE_RRD YES airphoto.img 2 4 8 16 ...

最后的...表示采样级别，一直到最顶层的像元个数小于等于64*64结束。有了上面的知识，下面就给出我写的一个函数，用于创建金字塔。

/**
* @brief 创建金字塔
* @param pszFileName        输入的栅格文件
* @param pProgress          进度条指针
* @return 成功返回0
*/
int CreatePyramids(const char* pszFileName, LT_Progress *pProgress)
{
if (pProgress != NULL)
{
pProgress->SetProgressCaption("创建金字塔");
pProgress->SetProgressTip("正在创建金字塔...");
}

GDALAllRegister();
CPLSetConfigOption("USE_RRD","YES");    //创建Erdas格式的字塔文件

/* -------------------------------------------------------------------- */
/*      Open data file.                                                 */
/* -------------------------------------------------------------------- */

GDALDatasetH     hDataset;
hDataset = GDALOpen( pszFileName, GA_ReadOnly );

GDALDriverH hDriver = GDALGetDatasetDriver(hDataset);
const char* pszDriver = GDALGetDriverShortName(hDriver);
if (EQUAL(pszDriver, "HFA") || EQUAL(pszDriver, "PCIDSK"))
{
GDALClose(hDataset);    //如果文件是Erdas的img或者PCI的pix格式，创建内金字塔，其他的创建外金字塔
hDataset = GDALOpen( pszFileName, GA_Update );
}

if( hDataset == NULL )
{
if (pProgress != NULL)
pProgress->SetProgressTip("打开图像失败，请检查图像是否存在或文件是否是图像文件！");

return RE_NOFILE;
}

/* -------------------------------------------------------------------- */
/*      Get File basic infomation                                       */
/* -------------------------------------------------------------------- */
int iWidth = GDALGetRasterXSize(hDataset);
int iHeigh = GDALGetRasterYSize(hDataset);

int iPixelNum = iWidth * iHeigh;    //图像中的总像元个数
int iTopNum = 4096;                 //顶层金字塔大小，64*64
int iCurNum = iPixelNum / 4;

int anLevels[1024] = { 0 };
int nLevelCount = 0;                //金字塔级数

do    //计算金字塔级数，从第二级到顶层
{
anLevels[nLevelCount] = static_cast<int>(pow(2.0, nLevelCount+2));
nLevelCount ++;
iCurNum /= 4;
} while (iCurNum > iTopNum);

const char      *pszResampling = "nearest"; //采样方式
GDALProgressFunc pfnProgress = GDALProgress;//进度条

/* -------------------------------------------------------------------- */
/*      Generate overviews.                                             */
/* -------------------------------------------------------------------- */
if (nLevelCount > 0 &&
GDALBuildOverviews( hDataset,pszResampling, nLevelCount, anLevels,
0, NULL, pfnProgress, pProgress ) != CE_None )
{
if (pProgress != NULL)
pProgress->SetProgressTip("创建金字塔失败！");

return RE_FAILD;
}

/* -------------------------------------------------------------------- */
/*      Cleanup                                                         */
/* -------------------------------------------------------------------- */
GDALClose(hDataset);
GDALDestroyDriverManager();

if (pProgress != NULL)
pProgress->SetProgressTip("创建金字塔完成！");

return RE_SUCCESS;
}

    需要说明的是，这段代码创建img格式和pix格式的金字塔会创建内金字塔，Erdas的img格式和PCI的pix格式可以把金字塔存放在文件内部。
PS：在给img文件创建内金字塔后，使用除ArcGIS9.2以外的软件打开后，都正常，但是使用ArcGIS9.2打开后会出现图层偏移的问题，不知道是否ArcGIS9.2的bug。ArcGIS10正常！ArcGIS9.3没有测试。
测试代码：void main()
{
LT_ConsoleProgress *pProgress = new LT_ConsoleProgress();

int f = CreatePyramids("C://Work//Data//ttttt.img", pProgress);

if (f == RE_SUCCESS)
printf("计算成功/n");
else
printf("计算失败/n");

delete pProgress;
}测试：在使用ArcGIS10打开没有金字塔的文件时，提示：



运行测试代码：



再次打开刚才的文件，没有上面的提示对话框了，而且很快加载进来，说明已经有金字塔了，如果使用Erdas打开的话，可以看到详细的金字塔信息。不过可以试用QGIS查看金字塔信息。右侧的列表显示的是金字塔的级别，Erdas的金字塔是从第二级开始建立的，所以看到第一级的图标上有个红色的小叉叉。见下图：



转自 http://blog.csdn.net/liminlu0314/article/details/6127755