geo/tcLandsatData.cpp

   1 /***************************************************************************
   2  *   This file is part of Tecorrec.                                        *
   3  *   Copyright 2008 James Hogan <james@albanarts.com>                      *
   4  *                                                                         *
   5  *   Tecorrec is free software: you can redistribute it and/or modify      *
   6  *   it under the terms of the GNU General Public License as published by  *
   7  *   the Free Software Foundation, either version 2 of the License, or     *
   8  *   (at your option) any later version.                                   *
   9  *                                                                         *
  10  *   Tecorrec is distributed in the hope that it will be useful,           *
  11  *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
  12  *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
  13  *   GNU General Public License for more details.                          *
  14  *                                                                         *
  15  *   You should have received a copy of the GNU General Public License     *
  16  *   along with Tecorrec.  If not, write to the Free Software Foundation,  *
  17  *   Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.   *
  18  ***************************************************************************/
  19
  20 /**
  21  * @file tcLandsatData.cpp
  22  * @brief Landsat image data.
  23  */
  24
  25 #include "tcLandsatData.h"
  26 #include "tcLandsatMetaData.h"
  27
  28 #include <QFile>
  29 #include <QImage>
  30
  31 /*
  32  * Constructors + destructor
  33  */
  34
  35 /// Load landsat data given the path.
  36 tcLandsatData::tcLandsatData(const QString& path)
  37 : tcPhotographyData()
  38 , m_meta(0)
  39 {
  40   /// @todo Find name of txt file automatically
  41   const QString metaFilePath = path + "L71195028_02819990725_MTL.txt";
  42
  43   QFile metaFile(metaFilePath);
  44   if (metaFile.open(QIODevice::ReadOnly))
  45   {
  46     m_meta = new tcLandsatMetaData(metaFile);
  47
  48     tcLandsatMetaData& productMetadata = (*m_meta)("L1_METADATA_FILE")("PRODUCT_METADATA");
  49
  50     // Load coordinates.
  51     m_coordinates[0][1] = tcGeo(productMetadata["PRODUCT_UL_CORNER_LON"].toDouble() * M_PI/180.0,
  52                                 productMetadata["PRODUCT_UL_CORNER_LAT"].toDouble() * M_PI/180.0);
  53     m_coordinates[1][1] = tcGeo(productMetadata["PRODUCT_UR_CORNER_LON"].toDouble() * M_PI/180.0,
  54                                 productMetadata["PRODUCT_UR_CORNER_LAT"].toDouble() * M_PI/180.0);
  55     m_coordinates[0][0] = tcGeo(productMetadata["PRODUCT_LL_CORNER_LON"].toDouble() * M_PI/180.0,
  56                                 productMetadata["PRODUCT_LL_CORNER_LAT"].toDouble() * M_PI/180.0);
  57     m_coordinates[1][0] = tcGeo(productMetadata["PRODUCT_LR_CORNER_LON"].toDouble() * M_PI/180.0,
  58                                 productMetadata["PRODUCT_LR_CORNER_LAT"].toDouble() * M_PI/180.0);
  59
  60     QString bands = productMetadata["BAND_COMBINATION"].toString();
  61
  62     QString band[8];
  63     band[0] = productMetadata["BAND1_FILE_NAME"].toString();
  64   }
  65 }
  66
  67 /// Destructor.
  68 tcLandsatData::~tcLandsatData()
  69 {
  70 }