您的位置：首页 > 其它

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

2011-01-06 22:02 239 查看

今天在使用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());
}

希望对大家有用！

内容来自用户分享和网络整理，不保证内容的准确性，如有侵权内容，可联系管理员处理

标签：

相关文章推荐

新的分享

章节导航