使用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(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的金字塔是从第二级开始建立的,所以看到第一级的图标上有个红色的小叉叉。见下图:

你可能感兴趣的:(GIS,RS,GDAL,C++编程技术,GDAL交流)