108 строки
4.0 KiB
C++
108 строки
4.0 KiB
C++
/******************************************************************************
|
|
* $Id$
|
|
*
|
|
* Project: MSG Driver
|
|
* Purpose: GDALDataset driver for MSG translator for read support.
|
|
* Author: Bas Retsios, retsios@itc.nl
|
|
*
|
|
******************************************************************************
|
|
* Copyright (c) 2004, ITC
|
|
*
|
|
* 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_priv.h"
|
|
#include "cpl_csv.h"
|
|
#include "ogr_spatialref.h"
|
|
#include "msgcommand.h"
|
|
|
|
#include <string>
|
|
#include <fstream>
|
|
|
|
/************************************************************************/
|
|
/* MSGRasterBand */
|
|
/************************************************************************/
|
|
|
|
class MSGDataset;
|
|
class ReflectanceCalculator;
|
|
|
|
class MSGRasterBand final : public GDALRasterBand
|
|
{
|
|
friend class MSGDataset;
|
|
|
|
public:
|
|
MSGRasterBand(MSGDataset *, int);
|
|
virtual ~MSGRasterBand();
|
|
virtual CPLErr IReadBlock(int, int, void *) override;
|
|
|
|
private:
|
|
double rRadiometricCorrection(unsigned int iDN, int iChannel, int iRow,
|
|
int iCol, MSGDataset *poGDS) const;
|
|
bool fScanNorth;
|
|
int iLowerShift; // nr of pixels that lower HRV image is shifted compared
|
|
// to upper
|
|
int iSplitLine; // line from top where the HRV image splits
|
|
int iLowerWestColumnPlanned;
|
|
int iSatellite; // satellite number 1,2,3,4 for MSG1, MSG2, MSG3 and MSG4
|
|
ReflectanceCalculator *m_rc;
|
|
static const double rRTOA[12];
|
|
};
|
|
|
|
/************************************************************************/
|
|
/* MSGDataset */
|
|
/************************************************************************/
|
|
class MSGDataset final : public GDALDataset
|
|
{
|
|
friend class MSGRasterBand;
|
|
|
|
public:
|
|
MSGDataset();
|
|
virtual ~MSGDataset();
|
|
|
|
static GDALDataset *Open(GDALOpenInfo *);
|
|
|
|
const OGRSpatialReference *GetSpatialRef() const override
|
|
{
|
|
return &m_oSRS;
|
|
}
|
|
|
|
virtual CPLErr GetGeoTransform(double *padfTransform) override;
|
|
|
|
private:
|
|
MSGCommand command;
|
|
double adfGeoTransform[6]; // Calculate and store once as GetGeoTransform
|
|
// may be called multiple times
|
|
OGRSpatialReference m_oSRS{};
|
|
OGRSpatialReference oLL;
|
|
OGRCoordinateTransformation *poTransform = nullptr;
|
|
double rCalibrationOffset[12];
|
|
double rCalibrationSlope[12];
|
|
int iCurrentSatellite =
|
|
0; // satellite number 1,2,3,4 for MSG1, MSG2, MSG3 and MSG4
|
|
static int iCurrentSatelliteHint; // hint for satellite number 1,2,3,4 for
|
|
// MSG1, MSG2, MSG3 and MSG4
|
|
static const double rCentralWvl[12];
|
|
static const double rVc[12];
|
|
static const double rA[12];
|
|
static const double rB[12];
|
|
static const int iCentralPixelVIS_IR;
|
|
static const int iCentralPixelHRV;
|
|
static const char *metadataDomain;
|
|
};
|