目录

简述

GDALBuildVRT工具和函数默认都是不支持异构波段的（即要建立虚拟影像目录的文件中有波段与其它有异的，少于指定波段等特征的），公司一个项目中需要支持这个特性，不足的以第1波段替代（除非波段数为0，那这在影像质检的时候就已经检查完了），这里记录一下修改方式。

修改源码

使用的GDAL版本是2.2.3，修改的源码文件为gdal-2.2.3\apps\gdalbuildvrt_lib.cpp。

1、修改DatasetProperty结构体

修改这个结构体，增加用于记录每个影像波段数的成员变量，修改结果如下

大致在源文件的75-90行

typedef struct{    int    isFileOK;    int    nRasterXSize;    int    nRasterYSize;    double adfGeoTransform[6];    int    nBlockXSize;    int    nBlockYSize;    GDALDataType firstBandType;    int*         panHasNoData;    double*      padfNoDataValues;    int    bHasDatasetMask;    int    nMaskBlockXSize;    int    nMaskBlockYSize;    int    nRasterBandCount;   /*新增，记录影像波段数*/} DatasetProperty;

2、修改VRTBuilder::AnalyseRaster函数

大致在源文件的406-789行

修改这个函数主要是去除对异构波段的检测，同时加上对来源影像波段数的记录。

将以下代码

大致在源文件的543行

int _nBands = GDALGetRasterCount(hDS);

修改为如下，记录来源影像的波段数

int _nBands = psDatasetProperties->nRasterBandCount =GDALGetRasterCount(hDS);

将以下代码注释掉，跳过波段数的检测

大致在源文件的685-690行

if (!bSeparate)        {            // if (nMaxBandNo > _nBands)            // {            //     CPLError(CE_Warning, CPLE_NotSupported,            //                 "gdalbuildvrt does not support heterogeneous band numbers. Skipping %s",            //             dsFileName);            //     return FALSE;            // }

将以下代码注释，去掉对波段颜色特征信息的匹配

大致在源文件的695行

for(j=0;j
    
      _nBands){ bandNo = 1;}                GDALRasterBandH hRasterBand = GDALGetRasterBand( hDS, bandNo );                /*去掉对波段颜色特征信息的比对*/                if (/* pasBandProperties[j].colorInterpretation != GDALGetRasterColorInterpretation(hRasterBand) || */                    pasBandProperties[j].dataType != GDALGetRasterDataType(hRasterBand))                {
    

3、修改VRTBuilder::CreateVRTNonSeparate函数

因为我们不做每个影像单独一个波段的版本，所以这里只修改NonSeparate版本。

这里主要修改VRT波段配置源(poVRTBand->ConfigureSource)时候的来源影像波段的选择，修改的地方不多，这里贴出稍微完整点的带注释的代码。

大致在源文件的950-978行

for(int j=0;j
    
      psDatasetProperties->nRasterBandCount ? 0 : j;            int nSelBand = panBandList[sourceBandIndex ] - 1;  // 修改的位置也就这两行            VRTSourcedRasterBand* poVRTBand = (VRTSourcedRasterBand*)hVRTBand;           /* 检测源影像的Nodata信息 */            VRTSimpleSource* poSimpleSource;            if (bAllowSrcNoData && psDatasetProperties->panHasNoData[nSelBand])            {                poSimpleSource = new VRTComplexSource();                poSimpleSource->SetNoDataValue( psDatasetProperties->padfNoDataValues[nSelBand] );            }            else                poSimpleSource = new VRTSimpleSource();            if( pszResampling )                poSimpleSource->SetResampling(pszResampling);            /*对VRT波段配置正确的源信息*/            poVRTBand->ConfigureSource( poSimpleSource,                                        (GDALRasterBand*)GDALGetRasterBand((GDALDatasetH)hProxyDS, nSelBand + 1),                                        FALSE,                                        dfSrcXOff, dfSrcYOff,                                        dfSrcXSize, dfSrcYSize,                                        dfDstXOff, dfDstYOff,                                        dfDstXSize, dfDstYSize );            /*添加源*/            poVRTBand->AddSource( poSimpleSource );        }