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)#读取图像第一波段并转换为数组
dat