liminlu0314

GDAL库扩展Landsat系列MTL文件格式支持

Landsat系列卫星提供的数据，一般都是每个波段一个tif文件，然后外加一个MTL.txt的元数据文件，使用gdal可以直接打开每个波段的tif文件，但是有时候想在打开tif数据的同时能够自动读取MTL文件里面的元数据信息，这个时候就只能自己再解析这个文件了。比较麻烦，下面就是针对这种情况，直接在gdal库里面扩展一种支持landsat的mtl的格式，直接打开mtl文件的同时自动打开tif数据以及元数据。

关于扩展GDAL库可以参考我之前写的CNSDTF格式扩展，以及gdal的官网说明，网址为http://www.gdal.org/gdal_drivertut.html。

首先我们看一下Landsat的数据目录格式，以Landsat8数据为例，如下图所示：

Landsat8一共有11个波段，里面包含12个tif文件和一个txt文件，Landsat8包括两个传感器，分别是OLI和TIRS，其中OLI传感器有9个光谱波段，空间分辨率为15m（全色波段1个）和30m（多光谱波段8个），成像幅宽为185km。 TIRS有2个热红外波段，空间分辨率为30m，所以一共11个波段，此外还有一个BQA的tif不太清楚干啥的。

知道了数据的类别，那么我们可以参考HDF数据的样子，将Landsat8的数据归类为4个子数据集，第一个OLI全色数据，第二个OLI多光谱数据，第三个TIRS数据，第四个就是BQA数据。

下面是Landsat的实现源码，将其放到frmts\landsat目录，然后参考官网的说明编写make文件等，重新编译gdal库即可。

/******************************************************************************
* $Id: landsatdataset.cpp 27044 2014-03-16 23:41:27Z rouault $
*
* Project:  Landsat 7/8 _MTL.TXT Driver
* Purpose:  Implementation of the LandsatDataset class.
* Author:   Minlu Li, [email protected]
*
******************************************************************************
* Copyright (c) 2015, Minlu Li
*
* Permission is hereby granted, free of charge, to any person obtaining a
* copy of this software and associated documentation files (the "Software"),
* to deal in the Software without restriction, including without limitation
* the rights to use, copy, modify, merge, publish, distribute, sublicense,
* and/or sell copies of the Software, and to permit persons to whom the
* Software is furnished to do so, subject to the following conditions:
*
* The above copyright notice and this permission notice shall be included
* in all copies or substantial portions of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
* FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
* DEALINGS IN THE SOFTWARE.
****************************************************************************/

#include "gdal_pam.h"
#include "gdal_proxy.h"
#include "ogr_spatialref.h"
#include "cpl_string.h"
#include "vrtdataset.h"
#include "cpl_multiproc.h"
#include "cplkeywordparser.h"

enum Satellite  // Satellites:
{
	LANDSAT7,	// Landsat 7
	LANDSAT8	// Landsat 8
};

/************************************************************************/
/* ==================================================================== */
/*			                   	LandsatDataset		             		*/
/* ==================================================================== */
/************************************************************************/

class CPL_DLL LandsatDataset : public GDALPamDataset
{
	VRTDataset *poVRTDS;
	std::vector<GDALDataset *> apoTifDS;

protected:
	virtual int         CloseDependentDatasets();

public:
	LandsatDataset();
	~LandsatDataset();

	virtual char **GetFileList(void);

	static GDALDataset *Open( GDALOpenInfo * );
	static int Identify( GDALOpenInfo *poOpenInfo );
	static int ParserLandsat7( const char* pszMtlFile, char **papszMetaInfo, int nSubdataset, LandsatDataset **ppoDS );
	static int ParserLandsat8( const char* pszMtlFile, char **papszMetaInfo, int nSubdataset, LandsatDataset **ppoDS );
};

/************************************************************************/
/* ==================================================================== */
/*                          LandSatRasterBand                           */
/* ==================================================================== */
/************************************************************************/

class LandSatRasterBand : public GDALPamRasterBand
{
	friend class LandsatDataset;

	GDALRasterBand *poVRTBand;

public:
	LandSatRasterBand( LandsatDataset *, int, GDALRasterBand * );
	virtual       ~LandSatRasterBand() {};

	virtual CPLErr IReadBlock( int, int, void * );
	virtual CPLErr IRasterIO( GDALRWFlag, int, int, int, int,
		void *, int, int, GDALDataType,
		int, int );
};

/************************************************************************/
/*                         LandSatRasterBand()                          */
/************************************************************************/

LandSatRasterBand::LandSatRasterBand( LandsatDataset *poTILDS, int nBand, 
									 GDALRasterBand *poVRTBand )

{
	this->poDS = poTILDS;
	this->poVRTBand = poVRTBand;
	this->nBand = nBand;
	this->eDataType = poVRTBand->GetRasterDataType();

	poVRTBand->GetBlockSize( &nBlockXSize, &nBlockYSize );
}

/************************************************************************/
/*                             IReadBlock()                             */
/************************************************************************/

CPLErr LandSatRasterBand::IReadBlock( int iBlockX, int iBlockY, void *pBuffer )

{
	return poVRTBand->ReadBlock( iBlockX, iBlockY, pBuffer );
}

/************************************************************************/
/*                             IRasterIO()                              */
/************************************************************************/

CPLErr LandSatRasterBand::IRasterIO( GDALRWFlag eRWFlag,
									int nXOff, int nYOff, int nXSize, int nYSize,
									void * pData, int nBufXSize, int nBufYSize,
									GDALDataType eBufType,
									int nPixelSpace, int nLineSpace )

{
	if(GetOverviewCount() > 0)
	{
		return GDALPamRasterBand::IRasterIO( eRWFlag, nXOff, nYOff, nXSize, nYSize,
			pData, nBufXSize, nBufYSize, eBufType,
			nPixelSpace, nLineSpace );
	}
	else //if not exist TIL overviews, try to use band source overviews
	{
		return poVRTBand->RasterIO( eRWFlag, nXOff, nYOff, nXSize, nYSize,
			pData, nBufXSize, nBufYSize, eBufType,
			nPixelSpace, nLineSpace );
	}
}

/************************************************************************/
/*                           LandsatDataset()                           */
/************************************************************************/

LandsatDataset::LandsatDataset()

{
	poVRTDS = NULL;
}

/************************************************************************/
/*                          ~LandsatDataset()                           */
/************************************************************************/

LandsatDataset::~LandsatDataset()

{
	CloseDependentDatasets();
}

/************************************************************************/
/*                        CloseDependentDatasets()                      */
/************************************************************************/

int LandsatDataset::CloseDependentDatasets()
{
	int bHasDroppedRef = GDALPamDataset::CloseDependentDatasets();

	if( poVRTDS )
	{
		bHasDroppedRef = TRUE;
		delete poVRTDS;
		poVRTDS = NULL;
	}

	while( !apoTifDS.empty() )
	{
		GDALClose( (GDALDatasetH) apoTifDS.back() );
		apoTifDS.pop_back();
	}

	return bHasDroppedRef;
}

/************************************************************************/
/*                              Identify()                              */
/************************************************************************/

int LandsatDataset::Identify( GDALOpenInfo *poOpenInfo )

{
	if( EQUALN(poOpenInfo->pszFilename,"LANDSAT:",8) )	//subdataset
		return TRUE;

	const char* pszExt = CPLGetExtension(poOpenInfo->pszFilename);
	if( !EQUAL(pszExt,"TXT") && !EQUAL(pszExt,"MET"))
		return FALSE;

	if( strstr((const char *) poOpenInfo->pabyHeader,"GROUP = L1_METADATA_FILE") == NULL )
		return FALSE;
	else
		return TRUE;
}

int LandsatDataset::ParserLandsat7( const char* pszMtlFile, char **papszMetaInfo, int nSubdataset, LandsatDataset **ppoDS )
{
	LandsatDataset *poDS = (LandsatDataset*)(*ppoDS);

	double dfCellSize = 0.0;
	char **papszRasterName = NULL;
	if(nSubdataset == 0)
	{
		dfCellSize = CPLAtof(CSLFetchNameValueDef(papszMetaInfo, "PROJECTION_PARAMETERS.GRID_CELL_SIZE_PAN", "0"));
		poDS->nRasterXSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_SAMPLES_PAN","0"));
		poDS->nRasterYSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_LINES_PAN","0"));
		papszRasterName = CSLAddString(papszRasterName, "BAND8");
	}
	else if(nSubdataset ==1)
	{
		dfCellSize = CPLAtof(CSLFetchNameValueDef(papszMetaInfo, "PROJECTION_PARAMETERS.GRID_CELL_SIZE_REF", "0"));
		poDS->nRasterXSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_SAMPLES_REF","0"));
		poDS->nRasterYSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_LINES_REF","0"));
		papszRasterName = CSLAddString(papszRasterName, "BAND1");
		papszRasterName = CSLAddString(papszRasterName, "BAND2");
		papszRasterName = CSLAddString(papszRasterName, "BAND3");
		papszRasterName = CSLAddString(papszRasterName, "BAND4");
		papszRasterName = CSLAddString(papszRasterName, "BAND5");
		papszRasterName = CSLAddString(papszRasterName, "BAND7");
	}
	else
	{
		dfCellSize = CPLAtof(CSLFetchNameValueDef(papszMetaInfo, "PROJECTION_PARAMETERS.GRID_CELL_SIZE_THM", "0"));
		poDS->nRasterXSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_SAMPLES_THM","0"));
		poDS->nRasterYSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_LINES_THM","0"));
		papszRasterName = CSLAddString(papszRasterName, "BAND61");
		papszRasterName = CSLAddString(papszRasterName, "BAND62");
	}

	if (!GDALCheckDatasetDimensions(poDS->nRasterXSize, poDS->nRasterYSize))
	{
		CSLDestroy(papszRasterName);
		return FALSE;
	}

	OGRSpatialReference oSRS;
	const char* pszMapProjection = CSLFetchNameValue(papszMetaInfo, "PROJECTION_PARAMETERS.MAP_PROJECTION");
	const char* pszDatum = NULL;
	const char* pszEllipsoid = NULL;
	const char* pszZone = NULL;

	if(strstr(pszMapProjection, "UTM") != NULL)
	{
		pszDatum = CSLFetchNameValue(papszMetaInfo, "PROJECTION_PARAMETERS.REFERENCE_DATUM");
		pszEllipsoid = CSLFetchNameValue(papszMetaInfo, "PROJECTION_PARAMETERS.REFERENCE_ELLIPSOID");
		pszZone = CSLFetchNameValue(papszMetaInfo, "UTM_PARAMETERS.ZONE_NUMBER");

		// trim double quotes. 
		if( pszDatum[0] == '"' )
			pszDatum++;
		if( pszDatum[strlen(pszDatum)-1] == '"' )
			((char *) pszDatum)[strlen(pszDatum)-1] = '\0';

		oSRS.SetWellKnownGeogCS(pszDatum);
		oSRS.SetUTM(atoi(pszZone));
	}

	char *pszWkt = NULL;
	if(oSRS.exportToWkt(&pszWkt) == OGRERR_NONE)
		poDS->SetProjection(pszWkt);

	double dfUL_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_UL_CORNER_MAPX","0"));
	double dfUL_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_UL_CORNER_MAPY","0"));
	double dfUR_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_UR_CORNER_MAPX","0"));
	double dfUR_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_UR_CORNER_MAPY","0"));
	double dfLL_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_LL_CORNER_MAPX","0"));
	double dfLL_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_LL_CORNER_MAPY","0"));
	double dfLR_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_LR_CORNER_MAPX","0"));
	double dfLR_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_LR_CORNER_MAPY","0"));

	double adfGeoTransform[6] = {dfUL_X, dfCellSize, 0, dfUL_Y, 0, -1*dfCellSize};
	poDS->SetGeoTransform(adfGeoTransform);

	/* -------------------------------------------------------------------- */
	/*      We need to open one of the images in order to establish         */
	/*      details like the band count and types.                          */
	/* -------------------------------------------------------------------- */
	GDALDataset *poTemplateDS = NULL;

	CPLString ostrFileName;
	ostrFileName.Printf( "L1_METADATA_FILE.%s_FILE_NAME", papszRasterName[0] );
	const char *pszFilename = CSLFetchNameValue( papszMetaInfo, ostrFileName.c_str() );
	if( pszFilename == NULL )
	{
		CPLError( CE_Failure, CPLE_AppDefined, "Missing %s in .TXT file.", ostrFileName.c_str() );
		CSLDestroy(papszRasterName);
		return FALSE;
	}

	CPLString osDirname = CPLGetDirname(pszMtlFile);

	// trim double quotes. 
	if( pszFilename[0] == '"' )
		pszFilename++;
	if( pszFilename[strlen(pszFilename)-1] == '"' )
		((char *) pszFilename)[strlen(pszFilename)-1] = '\0';

	const char* pszSubFile = CPLFormFilename(osDirname, pszFilename, NULL);
	poTemplateDS = (GDALDataset *) GDALOpen( pszSubFile, GA_ReadOnly );
	if( poTemplateDS == NULL || poTemplateDS->GetRasterCount() == 0)
	{
		if (poTemplateDS != NULL)
			GDALClose( poTemplateDS );

		CSLDestroy(papszRasterName);
		return FALSE;
	}

	poDS->SetProjection(poTemplateDS->GetProjectionRef());
	poTemplateDS->GetGeoTransform(adfGeoTransform);
	poDS->SetGeoTransform(adfGeoTransform);

	GDALRasterBand *poTemplateBand = poTemplateDS->GetRasterBand(1);
	GDALDataType eDT = poTemplateBand->GetRasterDataType();
	int          nBandCount = CSLCount(papszRasterName);

	poTemplateBand = NULL;
	GDALClose( poTemplateDS );

	/* -------------------------------------------------------------------- */
	/*      Create and initialize the corresponding VRT dataset used to     */
	/*      manage the tiled data access.                                   */
	/* -------------------------------------------------------------------- */
	poDS->poVRTDS = new VRTDataset(poDS->nRasterXSize, poDS->nRasterYSize);
	int iBand = 0;
	for( iBand = 0; iBand < nBandCount; iBand++ )
		poDS->poVRTDS->AddBand( eDT, NULL );

	/* Don't try to write a VRT file */
	poDS->poVRTDS->SetWritable(FALSE);

	/* -------------------------------------------------------------------- */
	/*      Create band information objects.                                */
	/* -------------------------------------------------------------------- */
	for( iBand = 0; iBand < nBandCount; iBand++ )
	{
		ostrFileName.Printf( "L1_METADATA_FILE.%s_FILE_NAME", papszRasterName[iBand] );
		pszFilename = CSLFetchNameValue( papszMetaInfo, ostrFileName.c_str() );
		if( pszFilename == NULL )
		{
			CPLError( CE_Failure, CPLE_AppDefined, "Missing %s in .TXT file.", ostrFileName.c_str() );
			CSLDestroy(papszRasterName);
			return FALSE;
		}

		// trim double quotes. 
		if( pszFilename[0] == '"' )
			pszFilename++;
		if( pszFilename[strlen(pszFilename)-1] == '"' )
			((char *) pszFilename)[strlen(pszFilename)-1] = '\0';

		pszSubFile = CPLFormFilename(osDirname, pszFilename, NULL);
		GDALDataset *poBandDS = (GDALDataset *) GDALOpen( pszSubFile, GA_ReadOnly );
		if( poBandDS == NULL || poBandDS->GetRasterCount() == 0)
		{
			if (poBandDS != NULL)
				GDALClose( poBandDS );

			CSLDestroy(papszRasterName);
			return FALSE;
		}

		poDS->apoTifDS.push_back( poBandDS );

		GDALRasterBand *pBand = poBandDS->GetRasterBand(1);
		LandSatRasterBand *pLandsatBand = new LandSatRasterBand( poDS, iBand+1, pBand);
		char** papszBandMeta = NULL;
		papszBandMeta = CSLAddString(papszBandMeta, pszFilename);
		pLandsatBand->SetMetadata(papszBandMeta);
		poDS->SetBand( iBand+1, pLandsatBand);
	}

	CSLDestroy(papszRasterName);
	return TRUE;
}

int LandsatDataset::ParserLandsat8( const char* pszMtlFile, char **papszMetaInfo, int nSubdataset, LandsatDataset **ppoDS )
{
	LandsatDataset *poDS = (LandsatDataset*)(*ppoDS);

	double dfCellSize = 0.0;
	char **papszRasterName = NULL;
	if(nSubdataset == 0)
	{
		dfCellSize = CPLAtof(CSLFetchNameValueDef(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_PANCHROMATIC", "0"));
		poDS->nRasterXSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.PANCHROMATIC_SAMPLES","0"));
		poDS->nRasterYSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.PANCHROMATIC_LINES","0"));
		papszRasterName = CSLAddString(papszRasterName, "BAND_8");
	}
	else if(nSubdataset ==1)
	{
		dfCellSize = CPLAtof(CSLFetchNameValueDef(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_REFLECTIVE", "0"));
		poDS->nRasterXSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.REFLECTIVE_SAMPLES","0"));
		poDS->nRasterYSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.REFLECTIVE_LINES","0"));
		papszRasterName = CSLAddString(papszRasterName, "BAND_1");
		papszRasterName = CSLAddString(papszRasterName, "BAND_2");
		papszRasterName = CSLAddString(papszRasterName, "BAND_3");
		papszRasterName = CSLAddString(papszRasterName, "BAND_4");
		papszRasterName = CSLAddString(papszRasterName, "BAND_5");
		papszRasterName = CSLAddString(papszRasterName, "BAND_6");
		papszRasterName = CSLAddString(papszRasterName, "BAND_7");
		papszRasterName = CSLAddString(papszRasterName, "BAND_9");
	}
	else if(nSubdataset ==2)
	{
		dfCellSize = CPLAtof(CSLFetchNameValueDef(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_THERMAL", "0"));
		poDS->nRasterXSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_SAMPLES","0"));
		poDS->nRasterYSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_LINES","0"));
		papszRasterName = CSLAddString(papszRasterName, "BAND_10");
		papszRasterName = CSLAddString(papszRasterName, "BAND_11");
	}
	else
	{
		dfCellSize = CPLAtof(CSLFetchNameValueDef(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_THERMAL", "0"));
		poDS->nRasterXSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_SAMPLES","0"));
		poDS->nRasterYSize = atoi(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_LINES","0"));
		papszRasterName = CSLAddString(papszRasterName, "BAND_QUALITY");
	}

	if (!GDALCheckDatasetDimensions(poDS->nRasterXSize, poDS->nRasterYSize))
	{
		CSLDestroy(papszRasterName);
		return FALSE;
	}

	OGRSpatialReference oSRS;
	const char* pszMapProjection = CSLFetchNameValue(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.MAP_PROJECTION");
	const char* pszDatum = NULL;
	const char* pszEllipsoid = NULL;
	const char* pszZone = NULL;

	if(strstr(pszMapProjection, "UTM") != NULL)
	{
		pszDatum = CSLFetchNameValue(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.DATUM");
		pszEllipsoid = CSLFetchNameValue(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.ELLIPSOID");
		pszZone = CSLFetchNameValue(papszMetaInfo, "L1_METADATA_FILE.PROJECTION_PARAMETERS.UTM_ZONE");

		// trim double quotes. 
		if( pszDatum[0] == '"' )
			pszDatum++;
		if( pszDatum[strlen(pszDatum)-1] == '"' )
			((char *) pszDatum)[strlen(pszDatum)-1] = '\0';

		oSRS.SetWellKnownGeogCS(pszDatum);
		oSRS.SetUTM(atoi(pszZone));
	}

	char *pszWkt = NULL;
	if(oSRS.exportToWkt(&pszWkt) == OGRERR_NONE)
		poDS->SetProjection(pszWkt);

	double dfUL_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_PROJECTION_X_PRODUCT","0"));
	double dfUL_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_PROJECTION_Y_PRODUCT","0"));
	double dfUR_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_PROJECTION_X_PRODUCT","0"));
	double dfUR_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_PROJECTION_Y_PRODUCT","0"));
	double dfLL_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_PROJECTION_X_PRODUCT","0"));
	double dfLL_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_PROJECTION_Y_PRODUCT","0"));
	double dfLR_X = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_PROJECTION_X_PRODUCT","0"));
	double dfLR_Y = CPLAtof(CSLFetchNameValueDef(papszMetaInfo,"L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_PROJECTION_Y_PRODUCT","0"));

	double adfGeoTransform[6] = {dfUL_X, dfCellSize, 0, dfUL_Y, 0, -1*dfCellSize};
	poDS->SetGeoTransform(adfGeoTransform);

	/* -------------------------------------------------------------------- */
	/*      We need to open one of the images in order to establish         */
	/*      details like the band count and types.                          */
	/* -------------------------------------------------------------------- */
	GDALDataset *poTemplateDS = NULL;
	
	CPLString ostrFileName;
	ostrFileName.Printf( "L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_%s", papszRasterName[0] );
	const char *pszFilename = CSLFetchNameValue( papszMetaInfo, ostrFileName.c_str() );
	if( pszFilename == NULL )
	{
		CPLError( CE_Failure, CPLE_AppDefined, "Missing %s in .TXT file.", ostrFileName.c_str() );
		CSLDestroy(papszRasterName);
		return FALSE;
	}

	CPLString osDirname = CPLGetDirname(pszMtlFile);

	// trim double quotes. 
	if( pszFilename[0] == '"' )
		pszFilename++;
	if( pszFilename[strlen(pszFilename)-1] == '"' )
		((char *) pszFilename)[strlen(pszFilename)-1] = '\0';

	const char* pszSubFile = CPLFormFilename(osDirname, pszFilename, NULL);
	poTemplateDS = (GDALDataset *) GDALOpen( pszSubFile, GA_ReadOnly );
	if( poTemplateDS == NULL || poTemplateDS->GetRasterCount() == 0)
	{
		if (poTemplateDS != NULL)
			GDALClose( poTemplateDS );

		CSLDestroy(papszRasterName);
		return FALSE;
	}

	poDS->SetProjection(poTemplateDS->GetProjectionRef());
	poTemplateDS->GetGeoTransform(adfGeoTransform);
	poDS->SetGeoTransform(adfGeoTransform);

	GDALRasterBand *poTemplateBand = poTemplateDS->GetRasterBand(1);
	GDALDataType eDT = poTemplateBand->GetRasterDataType();
	int          nBandCount = CSLCount(papszRasterName);

	poTemplateBand = NULL;
	GDALClose( poTemplateDS );

	/* -------------------------------------------------------------------- */
	/*      Create and initialize the corresponding VRT dataset used to     */
	/*      manage the tiled data access.                                   */
	/* -------------------------------------------------------------------- */
	poDS->poVRTDS = new VRTDataset(poDS->nRasterXSize, poDS->nRasterYSize);
	int iBand = 0;
	for( iBand = 0; iBand < nBandCount; iBand++ )
		poDS->poVRTDS->AddBand( eDT, NULL );

	/* Don't try to write a VRT file */
	poDS->poVRTDS->SetWritable(FALSE);

	/* -------------------------------------------------------------------- */
	/*      Create band information objects.                                */
	/* -------------------------------------------------------------------- */
	for( iBand = 0; iBand < nBandCount; iBand++ )
	{
		ostrFileName.Printf( "L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_%s", papszRasterName[iBand] );
		pszFilename = CSLFetchNameValue( papszMetaInfo, ostrFileName.c_str() );
		if( pszFilename == NULL )
		{
			CPLError( CE_Failure, CPLE_AppDefined, "Missing %s in .TXT file.", ostrFileName.c_str() );
			CSLDestroy(papszRasterName);
			return FALSE;
		}

		// trim double quotes. 
		if( pszFilename[0] == '"' )
			pszFilename++;
		if( pszFilename[strlen(pszFilename)-1] == '"' )
			((char *) pszFilename)[strlen(pszFilename)-1] = '\0';

		pszSubFile = CPLFormFilename(osDirname, pszFilename, NULL);
		GDALDataset *poBandDS = (GDALDataset *) GDALOpen( pszSubFile, GA_ReadOnly );
		if( poBandDS == NULL || poBandDS->GetRasterCount() == 0)
		{
			if (poBandDS != NULL)
				GDALClose( poBandDS );

			CSLDestroy(papszRasterName);
			return FALSE;
		}

		poDS->apoTifDS.push_back( poBandDS );

		GDALRasterBand *pBand = poBandDS->GetRasterBand(1);
		LandSatRasterBand *pLandsatBand = new LandSatRasterBand( poDS, iBand+1, pBand);
		char** papszBandMeta = NULL;
		papszBandMeta = CSLAddString(papszBandMeta, pszFilename);
		pLandsatBand->SetMetadata(papszBandMeta);
		poDS->SetBand( iBand+1, pLandsatBand);
	}

	CSLDestroy(papszRasterName);
	return TRUE;
}

/************************************************************************/
/*                                Open()                                */
/************************************************************************/

GDALDataset *LandsatDataset::Open( GDALOpenInfo * poOpenInfo )

{
	if( !Identify( poOpenInfo ) )
		return NULL;

	/* -------------------------------------------------------------------- */
	/*      Confirm the requested access is supported.                      */
	/* -------------------------------------------------------------------- */
	if( poOpenInfo->eAccess == GA_Update )
	{
		CPLError( CE_Failure, CPLE_NotSupported, 
			"The LANDSAT driver does not support update access to existing"
			" datasets.\n" );
		return NULL;
	}

	CPLString osFilename;
	int iSubdatasetIndex = -1;//0:PAN 1:REF 2:THM

	if( EQUALN(poOpenInfo->pszFilename, "LANDSAT:",8) )
	{
		const char *pszRest = poOpenInfo->pszFilename+8;

		iSubdatasetIndex = atoi(pszRest);
		while( *pszRest != '\0' && *pszRest != ':' )
			pszRest++;

		if( *pszRest == ':' )
			pszRest++;

		osFilename = pszRest;
	}
	else
		osFilename = poOpenInfo->pszFilename;

	if(iSubdatasetIndex <-1 || iSubdatasetIndex > 3)
	{
		CPLError( CE_Failure, CPLE_IllegalArg, "The LANDSAT driver does not support %d subdatasets.\n", iSubdatasetIndex );
		return NULL;
	}

	/* -------------------------------------------------------------------- */
	/*      Try to load and parse the .MTL .TXT file.                       */
	/* -------------------------------------------------------------------- */
	VSILFILE *fp = VSIFOpenL( osFilename, "r" );
	if( fp == NULL )
	{
		return NULL;
	}

	CPLKeywordParser oParser;
	if( !oParser.Ingest( fp ) )
	{
		VSIFCloseL( fp );
		return NULL;
	}

	VSIFCloseL( fp );

	char **papszMTL = oParser.GetAllKeywords();
	Satellite sat = LANDSAT7;

	const char* pszSatID = CSLFetchNameValue(papszMTL, "L1_METADATA_FILE.PRODUCT_METADATA.SPACECRAFT_ID");
	if(pszSatID == NULL)
	{
		return NULL;
	}

	if(strstr(pszSatID,"Landsat7") != NULL)
		sat = LANDSAT7;
	else if(strstr(pszSatID,"LANDSAT_8") != NULL)
		sat = LANDSAT8;
	else
	{
		return NULL;
	}

	if(sat == LANDSAT7 && iSubdatasetIndex ==3)	//landsat7 not quality
	{
		return NULL;
	}

	/* -------------------------------------------------------------------- */
	/*      Create a corresponding GDALDataset.                             */
	/* -------------------------------------------------------------------- */
	LandsatDataset 	*poDS;
	poDS = new LandsatDataset();
	poDS->SetMetadata(papszMTL);

	char *papszName[4] = {"PAN", "REF", "THE", "QUA"};
	char *papszDesc[4] = {"Panchromatic", "Reflective", "Thermal", "Quality"};
	if( iSubdatasetIndex == -1 )	//not a subdataset
	{
		int nSubdataset = (sat==LANDSAT7) ? 3 : 4;
		for (int i=0; i<nSubdataset; i++)
		{
			CPLString osKey, osValue;
			osKey.Printf( "SUBDATASET_%d_NAME", i+1 );
			osValue.Printf( "LANDSAT:%d:%s", i, osFilename.c_str() );
			poDS->SetMetadataItem( osKey, osValue, "SUBDATASETS" );

			osKey.Printf( "SUBDATASET_%d_DESC", i+1 );
			osValue.Printf( "LANDSAT:%s", papszDesc[i] );
			poDS->SetMetadataItem( osKey, osValue, "SUBDATASETS" );
		}

		return( poDS );
	}

	int nSuccess = FALSE;
	if(sat == LANDSAT7)
		nSuccess = ParserLandsat7(osFilename.c_str(), papszMTL, iSubdatasetIndex, &poDS);
	else
		nSuccess = ParserLandsat8(osFilename.c_str(), papszMTL, iSubdatasetIndex, &poDS);

	if(!nSuccess)
	{
		delete poDS;
		return NULL;
	}

	const char* pszDirname = CPLGetDirname(osFilename);
	const char* pszBasename = CPLGetBasename(osFilename);
	CPLString osTemp;
	osTemp.Printf("%s_%s", pszBasename, papszName[iSubdatasetIndex]);
	const char* pszSubFile = CPLFormFilename(pszDirname, osTemp, "");

	/* -------------------------------------------------------------------- */
	/*      Initialize any PAM information.                                 */
	/* -------------------------------------------------------------------- */
	poDS->SetDescription( pszSubFile );
	poDS->TryLoadXML();

	/* -------------------------------------------------------------------- */
	/*      Check for overviews.                                            */
	/* -------------------------------------------------------------------- */
	poDS->oOvManager.Initialize( poDS, pszSubFile );

	return( poDS );
}

/************************************************************************/
/*                            GetFileList()                             */
/************************************************************************/

char **LandsatDataset::GetFileList()

{
	unsigned int  i;
	char **papszFileList = GDALPamDataset::GetFileList();

	for( i = 0; i < apoTifDS.size(); i++ )
		papszFileList = CSLAddString( papszFileList, apoTifDS[i]->GetDescription() );

	return papszFileList;
}

/************************************************************************/
/*                         GDALRegister_Landsat()                       */
/************************************************************************/

void GDALRegister_Landsat()

{
	GDALDriver	*poDriver;

	if( GDALGetDriverByName( "LANDSAT" ) == NULL )
	{
		poDriver = new GDALDriver();

		poDriver->SetDescription( "LANDSAT" );
		poDriver->SetMetadataItem( GDAL_DMD_LONGNAME, "LANDSAT 7/8 GeoTiff with Metadata" );
		poDriver->SetMetadataItem( GDAL_DMD_HELPTOPIC, "frmt_landsat.html" );

		poDriver->SetMetadataItem( GDAL_DCAP_VIRTUALIO, "YES" );

		poDriver->pfnOpen = LandsatDataset::Open;
		poDriver->pfnIdentify = LandsatDataset::Identify;

		GetGDALDriverManager()->RegisterDriver( poDriver );
	}
}

编译之后，使用gdalinfo命令查看是否编译成功，处理输出如下，可以看到最后一行就是新加的landsat的驱动。

F:\gdal\bin>gdalinfo.exe --formats
Supported Formats:
  ECW (rov): ERDAS Compressed Wavelets (SDK 4.2)
  JP2ECW (rov): ERDAS JPEG2000 (SDK 4.2)
  FITS (rw+): Flexible Image Transport System
  GeoRaster (rw+s): Oracle Spatial GeoRaster
  HDF4 (ros): Hierarchical Data Format Release 4
  HDF4Image (rw+): HDF4 Dataset
  HDF5 (ros): Hierarchical Data Format Release 5
  HDF5Image (ro): HDF5 Dataset
  MG4Lidar (ro): MrSID Generation 4 / Lidar (.sid)
  MrSID (rov): Multi-resolution Seamless Image Database (MrSID)
  JP2MrSID (rov): MrSID JPEG2000
  netCDF (rw+s): Network Common Data Format
  GMT (rw): GMT NetCDF Grid Format
  PostGISRaster (rws): PostGIS Raster driver
  VRT (rw+v): Virtual Raster
  GTiff (rw+vs): GeoTIFF
  NITF (rw+vs): National Imagery Transmission Format
  RPFTOC (rovs): Raster Product Format TOC format
  ECRGTOC (rovs): ECRG TOC format
  HFA (rw+v): Erdas Imagine Images (.img)
  SAR_CEOS (rov): CEOS SAR Image
  CEOS (rov): CEOS Image
  。。。。太多了，删除了一些。。。
  IRIS (rov): IRIS data (.PPI, .CAPPi etc)
  WMTS (rwv): OGC Web Mab Tile Service
  CNSDTF (rwv): China Geospatial Data Transfer Grid Format (.grd)
  LANDSAT (rov): LANDSAT 7/8 GeoTiff with Metadata

再使用gdalinfo --format LANDSAT看看驱动详情

F:\gdal\bin>gdalinfo.exe --format LANDSAT
Format Details:
  Short Name: LANDSAT
  Long Name: LANDSAT 7/8 GeoTiff with Metadata
  Help Topic: frmt_landsat.html
  Supports: Virtual IO - eg. /vsimem/

如果输出上面的信息，说明gdal扩展成功，下面使用gdalinfo打开landsat8的数据试试，上面的驱动在实现的时候同时考虑了Landsat7的数据，所以这两个数据应该都可以支持。使用gdalinfo输出mtl.txt文件之后，输出内容如下：

F:\gdal\bin>gdalinfo.exe F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
Driver: LANDSAT/LANDSAT 7/8 GeoTiff with Metadata
Files: F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
Size is 512, 512
Coordinate System is `'
Metadata:
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.CLOUD_COVER=11.95
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.EARTH_SUN_DISTANCE=0.9925442
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_MODEL=8.958
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_MODEL_X=6.205
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_MODEL_Y=6.460
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_VERIFY=6.144
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GROUND_CONTROL_POINTS_MODEL=141
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GROUND_CONTROL_POINTS_VERIFY=36
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.IMAGE_QUALITY_OLI=9
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.IMAGE_QUALITY_TIRS=9
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.ROLL_ANGLE=-0.001
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.SUN_AZIMUTH=109.94537007
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.SUN_ELEVATION=59.45969108
  L1_METADATA_FILE.METADATA_FILE_INFO.FILE_DATE=2014-03-08T06:54:18Z
  L1_METADATA_FILE.METADATA_FILE_INFO.LANDSAT_SCENE_ID="LC81270562014067LGN00"
  L1_METADATA_FILE.METADATA_FILE_INFO.ORIGIN="Image courtesy of the U.S. Geological Survey"
  L1_METADATA_FILE.METADATA_FILE_INFO.PROCESSING_SOFTWARE_VERSION="LPGS_2.3.0"
  L1_METADATA_FILE.METADATA_FILE_INFO.REQUEST_ID="0501403074153_00047"
  L1_METADATA_FILE.METADATA_FILE_INFO.STATION_ID="LGN"
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_10=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_11=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_1=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_2=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_3=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_4=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_5=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_6=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_7=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_8=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_9=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_10=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_11=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_1=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_2=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_3=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_4=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_5=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_6=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_7=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_8=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_9=1
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_10=22.00180
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_11=22.00180
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_1=771.52448
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_2=790.05048
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_3=728.02478
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_4=613.91150
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_5=375.68323
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_6=93.42900
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_7=31.49057
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_8=694.77887
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_9=146.82552
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_10=0.10033
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_11=0.10033
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_1=-63.71274
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_2=-65.24262
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_3=-60.12052
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_4=-50.69701
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_5=-31.02404
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_6=-7.71540
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_7=-2.60050
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_8=-57.37506
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_9=-12.12490
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_1=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_2=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_3=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_4=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_5=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_6=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_7=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_8=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_9=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_1=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_2=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_3=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_4=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_5=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_6=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_7=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_8=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_9=-0.099980
  L1_METADATA_FILE.PRODUCT_METADATA.BPF_NAME_OLI="LO8BPF20140308030904_20140308035315.01"
  L1_METADATA_FILE.PRODUCT_METADATA.BPF_NAME_TIRS="LT8BPF20140308030510_20140308035408.01"
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_LAT_PRODUCT=4.73612
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_LON_PRODUCT=100.95098
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_PROJECTION_X_PRODUCT=716400.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_PROJECTION_Y_PRODUCT=523800.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_LAT_PRODUCT=4.72733
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_LON_PRODUCT=102.99592
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_PROJECTION_X_PRODUCT=943500.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_PROJECTION_Y_PRODUCT=523800.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_LAT_PRODUCT=6.83546
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_LON_PRODUCT=100.95819
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_PROJECTION_X_PRODUCT=716400.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_PROJECTION_Y_PRODUCT=756000.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_LAT_PRODUCT=6.82274
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_LON_PRODUCT=103.01065
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_PROJECTION_X_PRODUCT=943500.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_PROJECTION_Y_PRODUCT=756000.000
  L1_METADATA_FILE.PRODUCT_METADATA.CPF_NAME="L8CPF20140101_20140331.03"
  L1_METADATA_FILE.PRODUCT_METADATA.DATA_TYPE="L1T"
  L1_METADATA_FILE.PRODUCT_METADATA.DATE_ACQUIRED=2014-03-08
  L1_METADATA_FILE.PRODUCT_METADATA.ELEVATION_SOURCE="GLS2000"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_10="LC81270562014067LGN00_B10.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_11="LC81270562014067LGN00_B11.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_1="LC81270562014067LGN00_B1.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_2="LC81270562014067LGN00_B2.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_3="LC81270562014067LGN00_B3.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_4="LC81270562014067LGN00_B4.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_5="LC81270562014067LGN00_B5.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_6="LC81270562014067LGN00_B6.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_7="LC81270562014067LGN00_B7.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_8="LC81270562014067LGN00_B8.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_9="LC81270562014067LGN00_B9.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_QUALITY="LC81270562014067LGN00_BQA.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.METADATA_FILE_NAME="LC81270562014067LGN00_MTL.txt"
  L1_METADATA_FILE.PRODUCT_METADATA.NADIR_OFFNADIR="NADIR"
  L1_METADATA_FILE.PRODUCT_METADATA.OUTPUT_FORMAT="GEOTIFF"
  L1_METADATA_FILE.PRODUCT_METADATA.PANCHROMATIC_LINES=15481
  L1_METADATA_FILE.PRODUCT_METADATA.PANCHROMATIC_SAMPLES=15141
  L1_METADATA_FILE.PRODUCT_METADATA.REFLECTIVE_LINES=7741
  L1_METADATA_FILE.PRODUCT_METADATA.REFLECTIVE_SAMPLES=7571
  L1_METADATA_FILE.PRODUCT_METADATA.RLUT_FILE_NAME="L8RLUT20130211_20431231v09.h5"
  L1_METADATA_FILE.PRODUCT_METADATA.SCENE_CENTER_TIME=03:28:21.9117383Z
  L1_METADATA_FILE.PRODUCT_METADATA.SENSOR_ID="OLI_TIRS"
  L1_METADATA_FILE.PRODUCT_METADATA.SPACECRAFT_ID="LANDSAT_8"
  L1_METADATA_FILE.PRODUCT_METADATA.TARGET_WRS_PATH=127
  L1_METADATA_FILE.PRODUCT_METADATA.TARGET_WRS_ROW=56
  L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_LINES=7741
  L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_SAMPLES=7571
  L1_METADATA_FILE.PRODUCT_METADATA.WRS_PATH=127
  L1_METADATA_FILE.PRODUCT_METADATA.WRS_ROW=56
  L1_METADATA_FILE.PROJECTION_PARAMETERS.DATUM="WGS84"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.ELLIPSOID="WGS84"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_PANCHROMATIC=15.00
  L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_REFLECTIVE=30.00
  L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_THERMAL=30.00
  L1_METADATA_FILE.PROJECTION_PARAMETERS.MAP_PROJECTION="UTM"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.ORIENTATION="NORTH_UP"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.RESAMPLING_OPTION="CUBIC_CONVOLUTION"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.UTM_ZONE=47
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_10=0.10000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_11=0.10000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_1=-63.72549
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_2=-65.25568
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_3=-60.13255
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_4=-50.70715
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_5=-31.03025
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_6=-7.71694
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_7=-2.60102
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_8=-57.38654
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_9=-12.12732
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_10=3.3420E-04
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_11=3.3420E-04
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_1=1.2745E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_2=1.3051E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_3=1.2027E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_4=1.0141E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_5=6.2060E-03
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_6=1.5434E-03
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_7=5.2020E-04
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_8=1.1477E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_9=2.4255E-03
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_1=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_2=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_3=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_4=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_5=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_6=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_7=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_8=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_9=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_1=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_2=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_3=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_4=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_5=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_6=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_7=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_8=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_9=2.0000E-05
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K1_CONSTANT_BAND_10=774.89
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K1_CONSTANT_BAND_11=480.89
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K2_CONSTANT_BAND_10=1321.08
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K2_CONSTANT_BAND_11=1201.14
Subdatasets:
  SUBDATASET_1_NAME=LANDSAT:0:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
  SUBDATASET_1_DESC=LANDSAT:Panchromatic
  SUBDATASET_2_NAME=LANDSAT:1:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
  SUBDATASET_2_DESC=LANDSAT:Reflective
  SUBDATASET_3_NAME=LANDSAT:2:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
  SUBDATASET_3_DESC=LANDSAT:Thermal
  SUBDATASET_4_NAME=LANDSAT:3:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
  SUBDATASET_4_DESC=LANDSAT:Quality
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

从上面的输出内容可以看出，扩展的驱动可以将MTL中的所有的信息作为元数据进行输出，同时构造了四个subdataset。下面分别输出四个subdataset的信息。首先是subdataset1，也就是全色数据，子数据集的名称为LANDSAT:0:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt，下面以这个作为gdalinfo的参数，输出内容如下：

F:\gdal\bin>gdalinfo.exe LANDSAT:0:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
Driver: LANDSAT/LANDSAT 7/8 GeoTiff with Metadata
Files: F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B8.TIF
Size is 15141, 15481
Coordinate System is:
PROJCS["WGS 84 / UTM zone 47N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",99],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32647"]]
Origin = (716392.500000000000000,756007.500000000000000)
Pixel Size = (15.000000000000000,-15.000000000000000)
Metadata:
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.CLOUD_COVER=11.95
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.EARTH_SUN_DISTANCE=0.9925442
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_MODEL=8.958
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_MODEL_X=6.205
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_MODEL_Y=6.460
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GEOMETRIC_RMSE_VERIFY=6.144
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GROUND_CONTROL_POINTS_MODEL=141
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.GROUND_CONTROL_POINTS_VERIFY=36
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.IMAGE_QUALITY_OLI=9
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.IMAGE_QUALITY_TIRS=9
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.ROLL_ANGLE=-0.001
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.SUN_AZIMUTH=109.94537007
  L1_METADATA_FILE.IMAGE_ATTRIBUTES.SUN_ELEVATION=59.45969108
  L1_METADATA_FILE.METADATA_FILE_INFO.FILE_DATE=2014-03-08T06:54:18Z
  L1_METADATA_FILE.METADATA_FILE_INFO.LANDSAT_SCENE_ID="LC81270562014067LGN00"
  L1_METADATA_FILE.METADATA_FILE_INFO.ORIGIN="Image courtesy of the U.S. Geological Survey"
  L1_METADATA_FILE.METADATA_FILE_INFO.PROCESSING_SOFTWARE_VERSION="LPGS_2.3.0"
  L1_METADATA_FILE.METADATA_FILE_INFO.REQUEST_ID="0501403074153_00047"
  L1_METADATA_FILE.METADATA_FILE_INFO.STATION_ID="LGN"
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_10=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_11=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_1=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_2=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_3=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_4=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_5=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_6=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_7=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_8=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MAX_BAND_9=65535
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_10=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_11=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_1=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_2=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_3=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_4=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_5=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_6=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_7=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_8=1
  L1_METADATA_FILE.MIN_MAX_PIXEL_VALUE.QUANTIZE_CAL_MIN_BAND_9=1
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_10=22.00180
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_11=22.00180
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_1=771.52448
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_2=790.05048
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_3=728.02478
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_4=613.91150
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_5=375.68323
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_6=93.42900
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_7=31.49057
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_8=694.77887
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MAXIMUM_BAND_9=146.82552
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_10=0.10033
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_11=0.10033
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_1=-63.71274
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_2=-65.24262
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_3=-60.12052
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_4=-50.69701
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_5=-31.02404
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_6=-7.71540
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_7=-2.60050
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_8=-57.37506
  L1_METADATA_FILE.MIN_MAX_RADIANCE.RADIANCE_MINIMUM_BAND_9=-12.12490
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_1=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_2=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_3=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_4=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_5=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_6=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_7=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_8=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MAXIMUM_BAND_9=1.210700
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_1=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_2=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_3=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_4=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_5=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_6=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_7=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_8=-0.099980
  L1_METADATA_FILE.MIN_MAX_REFLECTANCE.REFLECTANCE_MINIMUM_BAND_9=-0.099980
  L1_METADATA_FILE.PRODUCT_METADATA.BPF_NAME_OLI="LO8BPF20140308030904_20140308035315.01"
  L1_METADATA_FILE.PRODUCT_METADATA.BPF_NAME_TIRS="LT8BPF20140308030510_20140308035408.01"
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_LAT_PRODUCT=4.73612
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_LON_PRODUCT=100.95098
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_PROJECTION_X_PRODUCT=716400.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LL_PROJECTION_Y_PRODUCT=523800.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_LAT_PRODUCT=4.72733
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_LON_PRODUCT=102.99592
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_PROJECTION_X_PRODUCT=943500.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_LR_PROJECTION_Y_PRODUCT=523800.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_LAT_PRODUCT=6.83546
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_LON_PRODUCT=100.95819
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_PROJECTION_X_PRODUCT=716400.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UL_PROJECTION_Y_PRODUCT=756000.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_LAT_PRODUCT=6.82274
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_LON_PRODUCT=103.01065
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_PROJECTION_X_PRODUCT=943500.000
  L1_METADATA_FILE.PRODUCT_METADATA.CORNER_UR_PROJECTION_Y_PRODUCT=756000.000
  L1_METADATA_FILE.PRODUCT_METADATA.CPF_NAME="L8CPF20140101_20140331.03"
  L1_METADATA_FILE.PRODUCT_METADATA.DATA_TYPE="L1T"
  L1_METADATA_FILE.PRODUCT_METADATA.DATE_ACQUIRED=2014-03-08
  L1_METADATA_FILE.PRODUCT_METADATA.ELEVATION_SOURCE="GLS2000"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_10="LC81270562014067LGN00_B10.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_11="LC81270562014067LGN00_B11.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_1="LC81270562014067LGN00_B1.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_2="LC81270562014067LGN00_B2.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_3="LC81270562014067LGN00_B3.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_4="LC81270562014067LGN00_B4.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_5="LC81270562014067LGN00_B5.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_6="LC81270562014067LGN00_B6.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_7="LC81270562014067LGN00_B7.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_8="LC81270562014067LGN00_B8.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_9="LC81270562014067LGN00_B9.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.FILE_NAME_BAND_QUALITY="LC81270562014067LGN00_BQA.TIF"
  L1_METADATA_FILE.PRODUCT_METADATA.METADATA_FILE_NAME="LC81270562014067LGN00_MTL.txt"
  L1_METADATA_FILE.PRODUCT_METADATA.NADIR_OFFNADIR="NADIR"
  L1_METADATA_FILE.PRODUCT_METADATA.OUTPUT_FORMAT="GEOTIFF"
  L1_METADATA_FILE.PRODUCT_METADATA.PANCHROMATIC_LINES=15481
  L1_METADATA_FILE.PRODUCT_METADATA.PANCHROMATIC_SAMPLES=15141
  L1_METADATA_FILE.PRODUCT_METADATA.REFLECTIVE_LINES=7741
  L1_METADATA_FILE.PRODUCT_METADATA.REFLECTIVE_SAMPLES=7571
  L1_METADATA_FILE.PRODUCT_METADATA.RLUT_FILE_NAME="L8RLUT20130211_20431231v09.h5"
  L1_METADATA_FILE.PRODUCT_METADATA.SCENE_CENTER_TIME=03:28:21.9117383Z
  L1_METADATA_FILE.PRODUCT_METADATA.SENSOR_ID="OLI_TIRS"
  L1_METADATA_FILE.PRODUCT_METADATA.SPACECRAFT_ID="LANDSAT_8"
  L1_METADATA_FILE.PRODUCT_METADATA.TARGET_WRS_PATH=127
  L1_METADATA_FILE.PRODUCT_METADATA.TARGET_WRS_ROW=56
  L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_LINES=7741
  L1_METADATA_FILE.PRODUCT_METADATA.THERMAL_SAMPLES=7571
  L1_METADATA_FILE.PRODUCT_METADATA.WRS_PATH=127
  L1_METADATA_FILE.PRODUCT_METADATA.WRS_ROW=56
  L1_METADATA_FILE.PROJECTION_PARAMETERS.DATUM="WGS84"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.ELLIPSOID="WGS84"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_PANCHROMATIC=15.00
  L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_REFLECTIVE=30.00
  L1_METADATA_FILE.PROJECTION_PARAMETERS.GRID_CELL_SIZE_THERMAL=30.00
  L1_METADATA_FILE.PROJECTION_PARAMETERS.MAP_PROJECTION="UTM"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.ORIENTATION="NORTH_UP"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.RESAMPLING_OPTION="CUBIC_CONVOLUTION"
  L1_METADATA_FILE.PROJECTION_PARAMETERS.UTM_ZONE=47
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_10=0.10000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_11=0.10000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_1=-63.72549
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_2=-65.25568
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_3=-60.13255
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_4=-50.70715
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_5=-31.03025
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_6=-7.71694
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_7=-2.60102
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_8=-57.38654
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_ADD_BAND_9=-12.12732
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_10=3.3420E-04
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_11=3.3420E-04
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_1=1.2745E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_2=1.3051E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_3=1.2027E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_4=1.0141E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_5=6.2060E-03
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_6=1.5434E-03
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_7=5.2020E-04
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_8=1.1477E-02
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.RADIANCE_MULT_BAND_9=2.4255E-03
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_1=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_2=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_3=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_4=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_5=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_6=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_7=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_8=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_ADD_BAND_9=-0.100000
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_1=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_2=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_3=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_4=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_5=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_6=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_7=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_8=2.0000E-05
  L1_METADATA_FILE.RADIOMETRIC_RESCALING.REFLECTANCE_MULT_BAND_9=2.0000E-05
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K1_CONSTANT_BAND_10=774.89
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K1_CONSTANT_BAND_11=480.89
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K2_CONSTANT_BAND_10=1321.08
  L1_METADATA_FILE.TIRS_THERMAL_CONSTANTS.K2_CONSTANT_BAND_11=1201.14
Corner Coordinates:
Upper Left  (  716392.500,  756007.500) (100d57'29.23"E,  6d50' 7.91"N)
Lower Left  (  716392.500,  523792.500) (100d57' 3.27"E,  4d44' 9.79"N)
Upper Right (  943507.500,  756007.500) (103d 0'38.59"E,  6d49'22.11"N)
Lower Right (  943507.500,  523792.500) (102d59'45.54"E,  4d43'38.14"N)
Center      (  829950.000,  639900.000) (101d58'44.08"E,  5d46'52.89"N)
Band 1 Block=15141x1 Type=UInt16, ColorInterp=Undefined
  Metadata:
    LC81270562014067LGN00_B8.TIF

从上面的输出可以看出，子数据集也将mtl中的元数据信息进行了输出，同时输出的全色数据的投影信息，四角坐标，分辨率等信息。

下面是多光谱数据，子数据集名称为LANDSAT:1:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt，由于输出的元数据信息一样，所以使用参数-nomd禁止输出元数据。输出内容如下：

F:\gdal\bin>gdalinfo.exe -nomd LANDSAT:1:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
Driver: LANDSAT/LANDSAT 7/8 GeoTiff with Metadata
Files: F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B1.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B2.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B3.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B4.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B5.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B6.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B7.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B9.TIF
Size is 7571, 7741
Coordinate System is:
PROJCS["WGS 84 / UTM zone 47N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",99],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32647"]]
Origin = (716385.000000000000000,756015.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Corner Coordinates:
Upper Left  (  716385.000,  756015.000) (100d57'28.99"E,  6d50' 8.15"N)
Lower Left  (  716385.000,  523785.000) (100d57' 3.03"E,  4d44' 9.55"N)
Upper Right (  943515.000,  756015.000) (103d 0'38.83"E,  6d49'22.36"N)
Lower Right (  943515.000,  523785.000) (102d59'45.78"E,  4d43'37.90"N)
Center      (  829950.000,  639900.000) (101d58'44.08"E,  5d46'52.89"N)
Band 1 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 2 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 3 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 4 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 5 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 6 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 7 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 8 Block=7571x1 Type=UInt16, ColorInterp=Undefined

从上面的信息中可以看出，输出的多光谱一共有8个波段，以及空间参考信息，分辨率，四至范围等。

下面输出TIRS的数据，命令行如上，子数据集的名称为LANDSAT:2:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt，输出的信息如下：

F:\gdal\bin>gdalinfo.exe -nomd LANDSAT:2:F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_MTL.txt
Driver: LANDSAT/LANDSAT 7/8 GeoTiff with Metadata
Files: F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B10.TIF
       F:\Data\LandSat\LC81270562014067LGN00\LC81270562014067LGN00_B11.TIF
Size is 7571, 7741
Coordinate System is:
PROJCS["WGS 84 / UTM zone 47N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",99],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32647"]]
Origin = (716385.000000000000000,756015.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Corner Coordinates:
Upper Left  (  716385.000,  756015.000) (100d57'28.99"E,  6d50' 8.15"N)
Lower Left  (  716385.000,  523785.000) (100d57' 3.03"E,  4d44' 9.55"N)
Upper Right (  943515.000,  756015.000) (103d 0'38.83"E,  6d49'22.36"N)
Lower Right (  943515.000,  523785.000) (102d59'45.78"E,  4d43'37.90"N)
Center      (  829950.000,  639900.000) (101d58'44.08"E,  5d46'52.89"N)
Band 1 Block=7571x1 Type=UInt16, ColorInterp=Undefined
Band 2 Block=7571x1 Type=UInt16, ColorInterp=Undefined

这里有个问题，可以看到输出的TIRS数据的分辨率为30m，但是网上的资料写的TIRS的分辨率都是100m，我直接使用gdalinfo输出了10.tif和11.tif的信息，发现里面的也是30m的分辨率，所以我觉得TIRS的分辨率应该也是30m才对。

有了上面的驱动，就可以直接使用gdal打开mtl文件从而自动将数据组合在一起了，在使用的时候也比较方便。希望对大家能有所帮助。

你可能感兴趣的:(GDAL,dirver,landsat8)

python tif转png Python与遥感 python 开发语言
importosfromosgeoimportgdalimportnumpyasnpfromPILimportImage#提取432三波段fromspectralimport*#输入文件夹路径defget_img(dataset_img):width=dataset_img.RasterXSize#获取行列数height=dataset_img.RasterYSizebands=dataset_i
Python实现TIFF 文件转换为 PNG 和 JPG 格式 sand&wich python 开发语言
在日常的图像处理工作中，可能会遇到需要将TIFF格式的图像转换为其他格式的情况，例如PNG和JPG。下面，本文将介绍如何使用Python和GDAL库实现这一功能。准备工作在开始之前，请确保已经安装了必要的库：GDAL（GeospatialDataAbstractionLibrary）可以使用以下命令安装GDAL：pipinstallgdal代码实现以下是一个将TIFF文件转换为PNG文件的示例代码
Python，gdal，osgeo，报错，ERROR 1: PROJ: proj_create_from_database: XX\Library\share\proj\proj.db lacks 飞天小女警出击 python 数据库开发语言
问题描述：Python运行使用gdal得程序，程序正常运行，但是出现了警告，我同时有Python3.9，和Python3.6的环境，分别安装了gdal。一开始Python3.6正常运行，无任何警告；Python3.9运行，会出现下方的警告；ERROR1:PROJ:proj_create_from_database:XX\Library\share\proj\proj.dblacksDATABASE
QT/C++中的GDAL多线程应用（读取）：发生的问题以及解决方案寻找你的25号底片 c++c语言 qt 图像处理
1.引言在使用GDAL库对TIF文件进行切割和创建瓦片金字塔时，为了提高创建效率，不得不考虑使用多线程处理。然而，在实际实现过程中，我遇到了许多问题。通过不断的尝试和优化，最终找到了有效的解决方案。本文将详细记录这一过程中的问题和解决方法。2.初始多线程尝试与问题2.1常规多线程实现最初，我尝试使用常规的多线程方法来处理TIF文件切割，但很快发现这一方法导致了大量错误：ERROR1:GF2-035
合并Shapefile文件 wavemenu python
在地理信息系统中，Shapefile（.shp）是一种非常常见的矢量数据格式，用于存储地理要素的空间信息和属性数据。在处理大量Shapefile时，有时需要将多个文件合并成一个，以便于统一管理和分析。本文将介绍如何使用Python的GDAL/OGR库实现多个Shapefile文件的合并操作，并讨论其中的一些技术细节和注意事项。背景介绍GDAL/OGR是一个开源的地理数据转换库，提供了丰富的地理数据
python gdal 地心坐标系转wgs84坐标系 hehehehe
fromosgeoimportosr,ogrimportnumpyasnpclassOsr_Coord_Convert(object):source=osr.SpatialReference()source.ImportFromProj4("+proj=geocent+datum=WGS84+units=m+no_defs")target=osr.SpatialReference()target.
GDAL下载地址砍柴人_ GDAL GDAL c++
gdal源码从下面地址下载：http://www.gdal.org/不想自己编译，编译好的gdal库，下载地址：http://www.gisinternals.com/index.html
配置anaconda虚拟环境并下载GDAL库白日做大梦 python
1.打开AnacondaPrompt2.创建虚拟环境condacreate-nmynamepython=3.63.激活新创建的虚拟环境condaactivatemyname4.下载gdal库（此步骤如果不成功用第5步的方法）condainstall-cconda-forgegdal用上面代码下载失败的话，可以加镜像网站下载：condainstall-cconda-forgegdal-ihttps:
简易配置C++ GDAL环境 VS2022 SatVision炼金士 C++地理数据处理 c++开发语言
提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档文章目录一、vcpkg安装gdal库二、使用步骤1.前置条件:2.VS2022环境配置2.1创建项目所需库文件2.2为项目添加/链接库文件一、vcpkg安装gdal库Vcpkg是一个命令行包管理工具，它可以极大地简化三方库的获取、安装、使用和部署流程。Vcpkg是微软团队在GitHub上的一个开源项目，它提供一系列简单的命令，自动下载
ModuleNotFoundError: No module named ‘osgeo‘或者ModuleNotFoundError: No module named ‘‘gdal” cjj_ python python linux conda jupyter
Linux系统下总是报错，Nomodulenamed‘osgeo‘或者ModuleNotFoundError:Nomodulenamed‘‘gdal”。尝试很多方法安装‘osgeo'或者’gdal'不成功最后安装成功，不报错方法如下：找linux下python的轮子，通过下载wheel文件（编译发布的.whl包）进行安装。1.根据python版本下载对应版本的.whl文件链接https://sou
【BI&AI】Lecture 12 - Emotion in Brain & AI 头发没了还会再长人工智能脑机 BCI 脑科学心理
EmotioninBrain&AI专业术语Limbicsystem边缘系统Amygdala杏仁核temporallobe颞叶hippocampus海马体CentralNucleus中央核medialamygdala内侧杏仁核ventralperiaqueductalgray腹侧中脑导水管周围灰质课程大纲Whatisemotion当大脑检测到某些具有挑战性的情况时，会引发一系列更多或更少无意识的生理
一个坐标系查询网站python获取所有坐标系数字化信息化智能化解决方案 python
技术路线选择我是使用的vue3开发的网页界面，element-plus构建网页组件，openlayer展示地图，express提供后端API，vercel进行在线部署。python获取所有坐标系想要展示所有坐标系，那需要先获取坐标系，怎么获取了？可以读取gdal，也可以读取其他库，这里我使用了地理与气象分析包gma导出坐标系的相关信息，代码如下：import jsonimport osfrom m
Spark SQL调优实战 sighting_info spark sql 大数据
1、新添参数说明//Driver和Executor内存和CPU资源相关配置--是否开启executor动态分配，开启时spark.executor.instances不生效spark.dynamicAllocation.enabled=false--配置Driver内存spark.dirver.memory=5g--driver最大结果大小，设置为0代表不限制，driver在拉取结果时，如果结果超
Springboot+GDAL3.x在windows上配置开发环境后端小肥肠 GDAL spring boot 后端 java
1.GDAL简介GDAL（GeospatialDataAbstractionLibrary）是一个开源的地理空间数据抽象库，用于处理各种栅格和矢量地理空间数据格式。它提供了一个统一的抽象数据模型，用于访问各种地理空间文件格式。GDAL通常用于GIS（地理信息系统）、遥感和地理空间分析等领域，可以用来转换不同的文件格式、处理栅格数据以及处理地理空间元数据。如果你在处理地理数据，GDAL是一个非常强大
基于python的CV变异系数计算软件（适用于长时间栅格数据） RS GIS遥感地信学习 python 开发语言 CV变异系数
1.下面是一个例子，计算下面22幅影像的CV变异系数，分析不同区域的波动性。每幅影像只有一个波段，且为tif格式。2.计算结果如下，CV值越高表明该地区的生态环境质量时间序列的稳定性越低，反之则稳定性越高。3.话不多说，直接上代码，有编程基础的可以自己调试，小白可以使用本人编好的exe程序，链接在下面。fromosgeoimportgdalimportnumpyasnpimportosimport
GEE入门篇|栅格数据集概述（二）：Image Collections（单幅图像的集合）不止GIS Google Earth Engine入门及进阶教程开发语言 javascript
在前面的内容中了解图像集合时，我们使用了Landsat8原始图像数据集，这些原始图像已经为您完成了一些重要的修正。然而，原始图像只是为Landsat8生成的多个图像集合之一，遥感界开发了额外的图像校正，有助于提高分析的准确性和一致性，每个不同图像处理路径的结果都存储在EarthEngine中不同的ImageCollection中。其中最突出的是ImageCollection，旨在最大限度地减少地球
麒麟v10 Ubuntu miniconda3 gdal 发布flask web项目 qq_33739712 ubuntu flask 数据库
麒麟v10Ubuntuminiconda3gdal发布flaskweb项目一、需要安装的软件VMwareWorkstationPro，虚拟机推荐17.5.x及以上版本，不容易蓝屏死机Xshell7、Xftp7，用于连接linux服务器终端，传输文件等Miniconda3,用于配置python脚本运行的虚环境麒麟v10sp1镜像文件二、新建虚拟机1.虚拟机配置2.系统安装不会格式化本地磁盘，只操作虚
Python环境安装GDAL @王坤 python python 开发语言
一、去官网下载对应版本的GDALhttps://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal很多包有windows-386.tar.gz和windows-amd64.tar.gz的字样-i386i386首先可以简化一个概念，i386=Intel80386。其实dui386通常被用来作为对Intel(英zhi特尔)32位微dao处理器的统称。-AMD64AMD64
python中写出遥感影像、保存数组到geotiff文件（geotifwrite）汤姆波利波 python python 图像处理矩阵 arcgis
python读入遥感影像后进行数组和矩阵的运算，在计算完成后还需要导出为geotiff格式。总共需要五个参数，其中：“坐标系”、"仿射参数"与输入影像的相同；“数据矩阵”、“数据类型”、“输出路径”自行设置。坐标系projection、仿射参数geo_tansform在读取影像的时候获取：python利用gdal读遥感卫星影像deftifwrite(savepath,data,geo_transf
python栅格遥感影像转矢量汤姆波利波 python python 数据库矩阵图像处理数据库开发 json
栅格影像转矢量其中，读取栅格影像，见另一篇文章：python利用gdal读遥感卫星影像fromosgeoimportgdal,ogrdefcreateShpfile_from_tiffile(shpfile,tiffile):fieldName="MASK"geotiff=tifread(tifpath)##读取栅格影像，具体代码见另一篇文章data=geotiff.dataarraydata_d
【实用的GDAL】使用 gdalwarp 影像裁剪和重采样单线程zz python gdal
目录主要功能全部参数常用参数介绍使用例子影像裁剪影像重采样2022.12.22补充主要功能Imagereprojectionandwarpingutility影像重采样和裁剪全部参数gdalwarp[--help-general][--formats][-s_srssrs_def][-t_srssrs_def][-ctstring][-to"NAME=VALUE"]*[-vshift|-novsh
【Python&RS】GDAL批量裁剪遥感影像/栅格数据 RS迷途小书童 Sensing python 开发语言
GDAL(GeospatialDataAbstractionLibrary)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。Python的GDAL库作为栅格数据的处理转换库，其支持几百种栅格数据格式，如常见的TIFF、ENVI、HFA、HDF4等。因为遥感影像大部分都是栅格数据，所以GDAL库非常适合处
python gdal使用shp裁剪遥感影像阿星_ gdal ogr GIS 开发语言
fromosgeoimportgdalimporttimet0=time.time()#dst=r"D:\temp\aus11.tif"#shp=r"D:\temp\shp\shp.shp"#src=r"D:\data\aus11.tif"dst=r"D:\temp\sz3.tif"shp=r"D:\data\gmy\shp\rect.shp"src=r"D:\data\gmy\building\
python中利用gdal.warp完成矢量裁剪栅格（遥感影像掩膜、shp裁剪tif、图像裁剪）汤姆波利波 python python 开发语言图像处理数据分析矩阵数据结构
gdal库能够很方便的完成栅格数据的裁剪，代码如下：importnumpyasnpfromosgeoimportgdaldefextract_by_shp(in_shp_path,inputpath,outputpath):input_raster=gdal.Open(inputpath)#利用gdal.Warp进行裁剪result=gdal.Warp(outputpath,input_raste
Python与ArcGIS系列（二十）GDAL之合并shp和geojson要素图层 Auto工程师 Python与Arcgis系列 python arcgis gdal arcpy gis自动化 gis二次开发
目录0简述1代码实现2结果展示0简述Shp格式是GIS中非常重要的数据格式，主要在Arcgis中使用，但在进行很多基于网页的空间数据可视化时，通常只接受GeoJSON格式的数据，众所周知JSON是利用键值对+嵌套来表示数据的一种格式，以其轻量、易解析的优点，被广泛使用与各种领域，而GeoJSON就是指在一套规定的语法规则下用JSON格式存储矢量数据。本篇将介绍如何利用Python的GDAL实现多个
M1 mac 编译java gdal 3.0.4和3.6.3到指定目录 Pisererer java macos arm
目录3.0.4说明设置环境变量以便能找到gdal依赖的库编译gdal编译java代码测试问题记录3.6.3环境安装步骤java其他参考文档此文档记录在mac上安装gdal到指定目录，3.0.4使用make构建，3.6.3使用cmake构建。3.0.4说明机器环境：m1pro芯片macos12。java环境：zulujdk8，通过dmg安装。软件环境：因为之前通过homebrew安装过postgis
Python 下的GDAL、rasterio的安装与学习游无穷 Python GDAL 地理空间数据分析
GDAL与rasterio的安装使用conda进行安装condainstallGDALcondainstallrsaterio使用pip进行安装pipinstallGDALpipinstallrasterio使用python扩展包进行安装扩展包下载链接：https://www.lfd.uci.edu/~gohlke/pythonlibs/#pip找到对应python版本的GDAL包和rasteri
用pip安装GDAL时出错解决办法 Circle-C Python gdal
用pip安装gdal会报错：pipinstall-ihttps://pypi.tuna.tsinghua.edu.cn/simplegdal错误代码如下：原因：GDAL不是纯净的python库，无法像pipinstallrequests这样安装库文件。Buildingwheelforgdal(setup.py)...errorERROR:Commanderroredoutwithexitstatu
Python pygame贪吃蛇小游戏 (200行完整代码+注释+可运行）星卯教育tony Python编程 pygame python
一、运行效果图二、完整代码#!/usr/bin/envpython#-*-coding:utf-8-*-#author：Wangdalitime:2021年1月20日16:08:44#python实现：贪吃蛇'''游戏玩法：回车开始游戏；空格暂停游戏/继续游戏；方向键/wsad控制小蛇走向''''''思路：用列表存储蛇的身体；用浅色表示身体，深色背景将身体凸显出来；蛇的移动：仔细观察，是：身体除头
python实现tif转为数组、数组转为栅格深度神经质患者 python python numpy
文章目录前言具体代码1.引入库2.读入数据总结前言在遥感应用领域，有时候需要对tif进行读取转为数组，有时候也需要将数组写为tif。基于此，我在想能否写一个python函数实现，方便以后调用。请看下文...具体代码1.引入库importosimportnumpyasnpfromosgeoimportgdal2.tif读取成数组defread_img(self,filename):dataset=g
ViewController添加button按钮解析。（翻译）张亚雄 c
<div class="it610-blog-content-contain" style="font-size: 14px"></div>// ViewController.m // Reservation software // // Created by 张亚雄 on 15/6/2.
mongoDB 简单的增删改查开窍的石头 mongodb
在上一篇文章中我们已经讲了mongodb怎么安装和数据库/表的创建。在这里我们讲mongoDB的数据库操作在mongo中对于不存在的表当你用db.表名他会自动统计下边用到的user是表明，db代表的是数据库添加(insert):
log4j配置 0624chenhong log4j
1) 新建java项目 2) 导入jar包，项目右击，properties—java build path—libraries—Add External jar，加入log4j.jar包。 3) 新建一个类com.hand.Log4jTest package com.hand; import org.apache.log4j.Logger; public class
多点触摸(图片缩放为例) 不懂事的小屁孩多点触摸
多点触摸的事件跟单点是大同小异的，上个图片缩放的代码，供大家参考一下 import android.app.Activity; import android.os.Bundle; import android.view.MotionEvent; import android.view.View; import android.view.View.OnTouchListener
有关浏览器窗口宽度高度几个值的解析换个号韩国红果果 JavaScript html
1 元素的 offsetWidth 包括border padding content 整体的宽度。 clientWidth 只包括内容区 padding 不包括border。 clientLeft = offsetWidth -clientWidth 即这个元素border的值 offsetLeft 若无已定位的包裹元素
数据库产品巡礼：IBM DB2概览蓝儿唯美 db2
IBM DB2是一个支持了NoSQL功能的关系数据库管理系统，其包含了对XML，图像存储和Java脚本对象表示（JSON）的支持。DB2可被各种类型的企业使用，它提供了一个数据平台，同时支持事务和分析操作，通过提供持续的数据流来保持事务工作流和分析操作的高效性。 DB2支持的操作系统 DB2可应用于以下三个主要的平台: 工作站，DB2可在Linus、Unix、Windo
java笔记5 a-john java
控制执行流程： 1，true和false 利用条件表达式的真或假来决定执行路径。例：（a==b）。它利用条件操作符“==”来判断a值是否等于b值，返回true或false。java不允许我们将一个数字作为布尔值使用，虽然这在C和C++里是允许的。如果想在布尔测试中使用一个非布尔值，那么首先必须用一个条件表达式将其转化成布尔值，例如if(a!=0)。 2，if-els
Web开发常用手册汇总 aijuans PHP
一门技术，如果没有好的参考手册指导,很难普及大众。这其实就是为什么很多技术，非常好，却得不到普遍运用的原因。正如我们学习一门技术，过程大概是这个样子： ①我们日常工作中，遇到了问题，困难。寻找解决方案，即寻找新的技术； ②为什么要学习这门技术？这门技术是不是很好的解决了我们遇到的难题，困惑。这个问题，非常重要，我们不是为了学习技术而学习技术，而是为了更好的处理我们遇到的问题，才需要学习新的
今天帮助人解决的一个sql问题 asialee sql
今天有个人问了一个问题，如下： type AD value A
意图对象传递数据百合不是茶 android 意图Intent Bundle对象数据的传递
学习意图将数据传递给目标活动; 初学者需要好好研究的 1,将下面的代码添加到main.xml中 <?xml version="1.0" encoding="utf-8"?> <LinearLayout xmlns:android="http:/
oracle查询锁表解锁语句 bijian1013 oracle object session kill
一.查询锁定的表如下语句，都可以查询锁定的表语句一： select a.sid, a.serial#, p.spid, c.object_name, b.session_id, b.oracle_username, b.os_user_name from v$process p, v$s
mac osx 10.10 下安装 mysql 5.6 二进制文件［tar.gz］征客丶 mysql osx
场景：在 mac osx 10.10 下安装 mysql 5.6 的二进制文件。环境：mac osx 10.10、mysql 5.6 的二进制文件步骤：[所有目录请从根“/”目录开始取，以免层级弄错导致找不到目录] 1、下载 mysql 5.6 的二进制文件，下载目录下面称之为 mysql5.6SourceDir；下载地址：http://dev.mysql.com/downl
分布式系统与框架 bit1129 分布式
RPC框架 Dubbo 什么是Dubbo Dubbo是一个分布式服务框架，致力于提供高性能和透明化的RPC远程服务调用方案，以及SOA服务治理方案。其核心部分包含: 远程通讯: 提供对多种基于长连接的NIO框架抽象封装，包括多种线程模型，序列化，以及“请求-响应”模式的信息交换方式。集群容错: 提供基于接
那些令人蛋痛的专业术语白糖_ spring Web SSO IOC
spring 【控制反转(IOC)/依赖注入(DI)】：由容器控制程序之间的关系，而非传统实现中，由程序代码直接操控。这也就是所谓“控制反转”的概念所在：控制权由应用代码中转到了外部容器，控制权的转移，是所谓反转。简单的说：对象的创建又容器(比如spring容器)来执行，程序里不直接new对象。 Web 【单点登录(SSO)】：SSO的定义是在多个应用系统中，用户
《给大忙人看的java8》摘抄 braveCS java8
函数式接口：只包含一个抽象方法的接口 lambda表达式：是一段可以传递的代码你最好将一个lambda表达式想象成一个函数，而不是一个对象，并记住它可以被转换为一个函数式接口。事实上，函数式接口的转换是你在Java中使用lambda表达式能做的唯一一件事。方法引用：又是要传递给其他代码的操作已经有实现的方法了，这时可以使
编程之美-计算字符串的相似度 bylijinnan java 算法编程之美
public class StringDistance { /** * 编程之美计算字符串的相似度 * 我们定义一套操作方法来把两个不相同的字符串变得相同，具体的操作方法为： * 1.修改一个字符（如把“a”替换为“b”）; * 2.增加一个字符（如把“abdd”变为“aebdd”）; * 3.删除一个字符（如把“travelling”变为“trav
上传、下载压缩图片 chengxuyuancsdn 下载
/** * * @param uploadImage --本地路径(tomacat路径) * @param serverDir --服务器路径 * @param imageType --文件或图片类型 * 此方法可以上传文件或图片.txt,.jpg,.gif等 */ public void upload(String uploadImage,Str
bellman-ford(贝尔曼-福特)算法 comsci 算法 F#
Bellman-Ford算法(根据发明者 Richard Bellman 和 Lester Ford 命名)是求解单源最短路径问题的一种算法。单源点的最短路径问题是指：给定一个加权有向图G和源点s，对于图G中的任意一点v，求从s到v的最短路径。有时候这种算法也被称为 Moore-Bellman-Ford 算法，因为 Edward F. Moore zu 也为这个算法的发展做出了贡献。与迪科
oracle ASM中ASM_POWER_LIMIT参数 daizj ASM oracle ASM_POWER_LIMIT 磁盘平衡
ASM_POWER_LIMIT 该初始化参数用于指定ASM例程平衡磁盘所用的最大权值，其数值范围为0~11，默认值为1。该初始化参数是动态参数，可以使用ALTER SESSION或ALTER SYSTEM命令进行修改。示例如下： SQL>ALTER SESSION SET Asm_power_limit=2;
高级排序:快速排序 dieslrae 快速排序
public void quickSort(int[] array){ this.quickSort(array, 0, array.length - 1); } public void quickSort(int[] array,int left,int right){ if(right - left <= 0
C语言学习六指针_何谓变量的地址一个指针变量到底占几个字节 dcj3sjt126com C语言
# include <stdio.h> int main(void) { /* 1、一个变量的地址只用第一个字节表示 2、虽然他只使用了第一个字节表示，但是他本身指针变量类型就可以确定出他指向的指针变量占几个字节了 3、他都只存了第一个字节地址，为什么只需要存一个字节的地址，却占了4个字节，虽然只有一个字节，但是这些字节比较多，所以编号就比较大，
phpize使用方法 dcj3sjt126com PHP
phpize是用来扩展php扩展模块的，通过phpize可以建立php的外挂模块,下面介绍一个它的使用方法,需要的朋友可以参考下安装（fastcgi模式）的时候，常常有这样一句命令：代码如下: /usr/local/webserver/php/bin/phpize 一、phpize是干嘛的？ phpize是什么？ phpize是用来扩展php扩展模块的，通过phpi
Java虚拟机学习 - 对象引用强度 shuizhaosi888 JAVA虚拟机
本文原文链接：http://blog.csdn.net/java2000_wl/article/details/8090276 转载请注明出处！无论是通过计数算法判断对象的引用数量，还是通过根搜索算法判断对象引用链是否可达，判定对象是否存活都与“引用”相关。引用主要分为：强引用(Strong Reference)、软引用(Soft Reference)、弱引用(Wea
.NET Framework 3.5 Service Pack 1（完整软件包）下载地址 happyqing .net 下载 framework
Microsoft .NET Framework 3.5 Service Pack 1（完整软件包） http://www.microsoft.com/zh-cn/download/details.aspx?id=25150 Microsoft .NET Framework 3.5 Service Pack 1 是一个累积更新，包含很多基于 .NET Framewo
JAVA定时器的使用 jingjing0907 java timer 线程定时器
1、在应用开发中，经常需要一些周期性的操作，比如每5分钟执行某一操作等。对于这样的操作最方便、高效的实现方式就是使用java.util.Timer工具类。 privatejava.util.Timer timer; timer = newTimer(true); timer.schedule( newjava.util.TimerTask() { public void run()
Webbench 流浪鱼 webbench
首页下载地址 http://home.tiscali.cz/~cz210552/webbench.html Webbench是知名的网站压力测试工具，它是由Lionbridge公司（http://www.lionbridge.com）开发。 Webbench能测试处在相同硬件上，不同服务的性能以及不同硬件上同一个服务的运行状况。webbench的标准测试可以向我们展示服务器的两项内容：每秒钟相
第11章动画效果（中） onestopweb 动画
index.html <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> <html xmlns="http://www.w3.org/
windows下制作bat启动脚本. sanyecao2314 java cmd 脚本 bat
java -classpath C:\dwjj\commons-dbcp.jar;C:\dwjj\commons-pool.jar;C:\dwjj\log4j-1.2.16.jar;C:\dwjj\poi-3.9-20121203.jar;C:\dwjj\sqljdbc4.jar;C:\dwjj\voucherimp.jar com.citsamex.core.startup.MainStart
Java进行RSA加解密的例子 tomcat_oracle java
加密是保证数据安全的手段之一。加密是将纯文本数据转换为难以理解的密文；解密是将密文转换回纯文本。　　数据的加解密属于密码学的范畴。通常，加密和解密都需要使用一些秘密信息，这些秘密信息叫做密钥，将纯文本转为密文或者转回的时候都要用到这些密钥。　　对称加密指的是发送者和接收者共用同一个密钥的加解密方法。　　非对称加密(又称公钥加密)指的是需要一个私有密钥一个公开密钥，两个不同的密钥的
Android_ViewStub 阿尔萨斯 ViewStub
public final class ViewStub extends View java.lang.Object android.view.View android.view.ViewStub 类摘要： ViewStub 是一个隐藏的，不占用内存空间的视图对象，它可以在运行时延迟加载布局资源文件。当 ViewSt