关于GDAL打开hfa大文件的问题

     今天在使用GDAL打开大的img文件的时候,(这里所谓的大文件指的是img文件太大,会将数据文件存放到ige格式raw文件中)。在讲img文件和ige文件重命名后,使用GDAL打开文件后,只能读取到文件信息,不能读取图像的数据文件。仔细查看GDAL源代码发现,在img文件中记录了对应的ige文件的名称,重命名后img文件中的记录ige文件名还是原来的,找不到ige文件,所以就打不开了。但是在使用Erdas和ArcGIS打开该文件时,会正常打开,于是查看GDAL代码,修改部分代码,能够使GDAL正常打开。

    修改的代码位置如下,gdal源代码目录/frmts/hfa/hfaband.cpp中 367行处的代码修改为下面的:

/* -------------------------------------------------------------------- */
/*      Open raw data file.                                             */
/* -------------------------------------------------------------------- */
    const char *pszRawFilename = poDMS->GetStringField( "fileName.string" );
    const char *pszFullFilename;

    pszFullFilename = CPLFormFilename( psInfo->pszPath, pszRawFilename, NULL );

    if( psInfo->eAccess == HFA_ReadOnly )
        fpExternal = VSIFOpenL( pszFullFilename, "rb" );
    else
        fpExternal = VSIFOpenL( pszFullFilename, "r+b" );

    if( fpExternal == NULL )
    {
        CPLString strFileName = psInfo->pszFilename;
        strFileName = strFileName.substr(strFileName.find_last_of('.')+1) + "ige";
        pszFullFilename = CPLFormFilename( psInfo->pszPath, strFileName.c_str(), NULL );

        if( psInfo->eAccess == HFA_ReadOnly )
            fpExternal = VSIFOpenL( pszFullFilename, "rb" );
        else
            fpExternal = VSIFOpenL( pszFullFilename, "r+b" );

        if( fpExternal == NULL )
        {
            CPLError( CE_Failure, CPLE_OpenFailed, 
                      "Unable to open external data file:/n%s/n", 
                      pszFullFilename );
            return CE_Failure;
        }

        psInfo->pszIGEFilename = const_cast<char*>(strFileName.c_str());
    }
   
希望对大家有用!

posted on 2011-01-06 22:02  王大王  阅读(291)  评论(0编辑  收藏  举报

导航