多光谱影像与全色影像的融合代码

由于多种软件和编程语言都可以用于多光谱影像与全色影像融合，以下提供几种可能的代码示例：

1. 使用ENVI软件的IDL语言

ENVI软件提供了IDL语言的接口，可以通过编写IDL代码实现多光谱影像和全色影像的融合。以下是一个简单的IDL代码示例：

; 读取多光谱影像和全色影像
file1 = 'multispectral.img'
file2 = 'panchromatic.img'
multispectral = ENVI_OPEN_DATA(file1, /READ)
panchromatic = ENVI_OPEN_DATA(file2, /READ)

; 将全色影像插值到多光谱影像的空间分辨率
panchromatic_resampled = ENVIRESAMPLE(panchromatic, multispectral)

; 将全色影像融合到多光谱影像中
fused = ENVIMULTIFUSE(multispectral, panchromatic_resampled)

; 保存融合后的影像
outfile = 'fused.img'
ENVISAVE, fused, outfile, /INTERLEAVE

2. 使用Python的GDAL库

GDAL是一个开源的地理信息系统库，可以用Python语言调用其API实现多光谱影像和全色影像的融合。以下是一个简单的Python代码示例：

import gdal

# 读取多光谱影像和全色影像
multispectral_file = 'multispectral.tif'
panchromatic_file = 'panchromatic.tif'
multispectral = gdal.Open(multispectral_file, gdal.GA_ReadOnly)
panchromatic = gdal.Open(panchromatic_file, gdal.GA_ReadOnly)

# 将全色影像插值到多光谱影像的空间分辨率
panchromatic_resampled = gdal.Warp('', panchromatic, format='VRT', outputBounds=multispectral.GetGeoTransform(),
                                   xRes=multispectral.GetGeoTransform()[1], yRes=multispectral.GetGeoTransform()[5])

# 将全色影像融合到多光谱影像中
fused = gdal.BandMerge(multispectral.GetRasterBand(1), multispectral.GetRasterBand(2), multispectral.GetRasterBand(3),
                       panchromatic_resampled.GetRasterBand(1))

# 保存融合后的影像
driver = gdal.GetDriverByName('GTiff')
outfile = 'fused.tif'
outdataset = driver.Create(outfile, multispectral.RasterXSize, multispectral.RasterYSize, 4, gdal.GDT_Byte)
outdataset.SetProjection(multispectral.GetProjection())
outdataset.SetGeoTransform(multispectral.GetGeoTransform())
outdataset.GetRasterBand(1).WriteArray(fused[:,:,0])
outdataset.GetRasterBand(2).WriteArray(fused[:,:,1])
outdataset.GetRasterBand(3).WriteArray(fused[:,:,2])
outdataset.GetRasterBand(4).WriteArray(fused[:,:,3])
outdataset.FlushCache()

3. 使用MATLAB的图像处理工具箱

MATLAB的图像处理工具箱提供了多种图像融合算法，可以用MATLAB语言实现多光谱影像和全色影像的融合。以下是一个简单的MATLAB代码示例：

% 读取多光谱影像和全色影像
multispectral_file = 'multispectral.tif';
panchromatic_file = 'panchromatic.tif';
multispectral = imread(multispectral_file);
panchromatic = imread(panchromatic_file);

% 将全色影像插值到多光谱影像的空间分辨率
panchromatic_resampled = imresize(panchromatic, size(multispectral), 'bilinear');

% 将全色影像融合到多光谱影像中
fused = imlincomb(0.3, multispectral(:,:,1), 0.59, multispectral(:,:,2), 0.11, multispectral(:,:,3), 0.4, panchromatic_resampled);

% 保存融合后的影像
outfile = 'fused.tif';
imwrite(fused, outfile);

以上代码示例仅为参考，实际应用中需要根据具体数据和需求进行修改和优化。