Brovey图像融合方法+Python代码

Brovey图像融合方法+Python代码

1、方法原理

原理：Brovey变换融合又叫做比值变换融合，通过对多光谱数据进行归一化后，与高分辨率影像进行乘积来增加影像信息；

                                           

2、实现代码

# coding=utf-8
'''
Created on 2020-6-1
@author: jiangao
Project: Brovey图像融合方法
'''
import numpy as np
import cv2
import scipy.misc as smi
from osgeo import gdal
from PIL import Image
def gdal_open(path):
    """
    读取图像函数
    输入：图像路径
    返回：np.ndArray格式的三维数组
    """
    data = gdal.Open(path)
    col = data.RasterXSize#读取图像长度
    row = data.RasterYSize#读取图像宽度
    data_array_r = data.GetRasterBand(1).ReadAsArray(0,0,col,row).astype(np.float)#读取图像第一波段并转换为数组
    data_array_g = data.GetRa