ImageFusion.zip

实现图像融合处理操作，以及融合影像质量的客观评价。包含有主成分分析PCA、高通滤波HPF、乘积性变换和其改进型、Brovey变换及其改进型、相关系数的加权平均、IHS变换（IHS、三角IHS、圆柱体IHS及其变种）等算法，以及信息熵、平均梯度、空间频率、相关关系等客观评价指标。

#include “Function.h“
#include 
#include 
#include “windows.h“
#include 
#include 
#include 
#include 
#include 
#include 
#include 


using namespace std;
///////////////////////////////////////////////////////////////////////////////
/*                        图像的读取与保存                                   */
///////////////////////////////////////////////////////////////////////////////

unsigned char* OpenTif（const char *FilePath）
{
	GDALAllRegister（）;//注册、读取图像
	CPLSetConfigOption（“GDAL_FILENAME_IS_UTF8““NO“）;//使之支持中文路径
	GDALDataset *poDataset = NULL;
	poDataset = （GDALDataset*）GDALOpen（FilePathGA_ReadOnly）;
	if（poDataset == NULL）
	{
		cout<<“无法打开影像！“<		GDALDestroyDriverManager（）;
	}
	//获取图像数据的参数
	int width = poDataset->GetRasterXSize（）;
	int height = poDataset->GetRasterYSize（）;
	int nRastercount = poDataset->GetRasterCount（）;
	//开辟内存
	unsigned char *pImageData = new unsigned char[width * height];
	int bandList = {1};
	poDataset->RasterIO（GF_Read00widthheightpImageDatawidthheightGDT_Byte1&bandList000）;
	//cout<<“单波段影像读入完成！“<	//关闭GDAL库相关驱动和释放内存
	GDALClose（poDataset）;
	return pImageData;
}

unsigned char** GdalOpenTif（const char *FilePath）
{
	GDALAllRegister（）;//注册、读取图像
	CPLSetConfigOption（“GDAL_FILENAME_IS_UTF8““NO“）;//使之支持中文路径
	GDALDataset *poDataset = NULL;
	poDataset = （GDALDataset*）GDALOpen（FilePathGA_ReadOnly）;
	if（poDataset == NULL）
	{
		cout<<“无法打开影像！“<		GDALDestroyDriverManager（）;
	}
	//获取图像数据的参数
	int width = poDataset->GetRasterXSize（）;
	int height = poDataset->GetRasterYSize（）;
	int nRastercount = poDataset->GetRasterCount（）;
	//开辟内存
	unsigned char **pImageData = new unsigned char *[nRastercount];
	if （nRastercount == 1）//单波段影像
	{
		int bandList = {1};
		*pImageData= new unsigned char[width*height];
		poDataset->RasterIO（GF_Read00widthheight*pImageDatawidthheightGDT_Byte1&bandList000）;
		cout<<“单波段影像读入完成！“<		//GDALClose（poDataset）;
		//return pImageData;
	} 
	else if（nRastercount >= 3 ）//多波段影像
	{
		for （int j=0;j< nRastercount;j++）
		{
			pImageData[j] =  new unsigned char[width*height];
		}
		for （int i = 1;i <= nRastercount;i++）
		{
			int bandList = {i};
			poDataset->RasterIO（GF_Read00widthheightpImageData[i-1]widthheightGDT_Byte1&bandList000）;
			//GDALRasterBand *pBand;
			//pBand = poDataset->GetRasterBand（i）;
			//CPLErr error;
			//error = pBand->RasterIO（GF_Read00widthheightpImageData[i-1]widthheightGDT_Byte00）;
			//   if （error == CE_Failure）
			//   {
			//	cout<<“读取图像数据时失败！“<			//	GDALDestroyDriverManager（）;
			//   }
		}
		cout<<“多光谱影像读入完成！“<	}
	//关闭GDAL库相关驱动和释放内存
	GDALClose（poDataset）;
	return pImageData;
}

unsigned char** GDALreadImage（int &widthint &heightint &nRastercountconst char *filepath）
{
	GDALAllRegister（）;//注册、读取图像
	CPLSetConfigOption（“GDA

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-05-18 12:08  ImageFusion_8bit\ImgFusion\
     目录           0  2019-01-21 15:23  ImageFusion_8bit\ImgFusion\Debug\
     目录           0  2019-05-17 18:46  ImageFusion_8bit\ImgFusion\GDAL\
     目录           0  2019-05-17 18:46  ImageFusion_8bit\ImgFusion\GDAL\gdal200\
     目录           0  2019-05-17 18:46  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\
     文件        8760  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\compdcs.csv
     文件        7921  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\coordinate_axis.csv
     文件       11977  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\cubewerx_extra.wkt
     文件      166471  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\datum_shift.csv
     文件      363706  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\ecw_cs.wkt
     文件       13866  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\ellipsoid.csv
     文件          50  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\epsg.wkt
     文件      209164  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\esri_extra.wkt
     文件      332546  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\esri_StatePlane_extra.wkt
     文件       75845  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\esri_Wisconsin_extra.wkt
     文件       50304  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gcs.csv
     文件         513  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gcs.override.csv
     文件      208777  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gdal_datum.csv
     文件        2021  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gdalicon.png
     文件       13022  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\GDALLogoBW.svg
     文件       12305  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\GDALLogoColor.svg
     文件       12305  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\GDALLogoGS.svg
     文件       16914  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gdalvrt.xsd
     文件       11929  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\geoccs.csv
     文件        3642  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gml_registry.xml
     文件       15804  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gt_datum.csv
     文件        1719  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\gt_ellips.csv
     文件        9120  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\header.dxf
     文件        1740  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\inspire_cp_BasicPropertyUnit.gfs
     文件        1650  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\inspire_cp_CadastralBoundary.gfs
     文件        2450  2015-06-15 03:06  ImageFusion_8bit\ImgFusion\GDAL\gdal200\data\inspire_cp_CadastralParcel.gfs
............此处省略341个文件信息