msg_header_image_navig.~cpp

HRIT读取,用于在LINUX下显示高束数据图像

//-----------------------------------------------------------------------------
//
//  File        : MSG_header_image_navig.cpp
//  Description : MSG HRIT-LRIT format interface
//  Project     : Meteosatlib
//  Author      : Graziano Giuliani
//  References  : MSG/SPE/057 LRIT-HRIT Mission Specific Implementation,
//                V. 4.1 9 Mar 2001
//
//  This program is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2 of the License, or
//  (at your option) any later version.
//
//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License
//  along with this program; if not, write to the Free Software
//  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
//-----------------------------------------------------------------------------

#include <cstring>

#include <MSG_header_image_navig.h>

MSG_header_image_navig::MSG_header_image_navig( )
{
  satellite_h = 42164.0;
}

MSG_header_image_navig::MSG_header_image_navig( unsigned const char_1 *buff )
{
  read_from(buff);
}

MSG_header_image_navig::~MSG_header_image_navig( ) { }

void MSG_header_image_navig::read_from( unsigned const char_1 *buff )
{
  char_1 projection_name[32];
  memcpy(projection_name, buff+3, 32);
  column_scaling_factor = get_i4(buff+35);
  CFAC                  = column_scaling_factor;
  line_scaling_factor   = get_i4(buff+39);
  LFAC                  = line_scaling_factor;
  column_offset         = get_i4(buff+43);
  COFF                  = column_offset;
  line_offset           = get_i4(buff+47);
  LOFF                  = line_offset;
  if (projection_name[0] == 'G')
  {
    projection_name[4] = 0;
    pname = projection_name;
    char_1 *pnt = strchr(projection_name+5, ')');
    if (pnt != NULL) *pnt = 0;
    sscanf(projection_name+5, "%f", &subsatellite_longitude);
  }
  else if (projection_name[0] == 'P' && projection_name[1] == 'O')
  {
    char_1 *tmp;
    projection_name[5] = 0;
    pname = projection_name;
    tmp = strchr(projection_name+6, ',');
    *(tmp) = 0;
    sscanf(projection_name+6, "%f", &projection_plane);
    *(strchr(tmp, ')')) = 0;
    sscanf(tmp, "%f", &projection_longitude);
  }
  else if (projection_name[0] == 'P' && projection_name[1] == 'S')
    pname = "PSD";
  else if (projection_name[0] == 'M')
    pname = "MERCATOR";
  else
    std::cerr << "Unknown Projection in Navigation Header." << std::endl;
  return;
}

t_enum_MSG_projection MSG_header_image_navig::get_projection_code( )
{
  if (pname == "GEOS")
    return MSG_PROJECTION_GEOS;
  else if (pname == "POLAR")
    return MSG_PROJECTION_POLAR;
  else if (pname == "MERCATOR")
    return MSG_PROJECTION_MERCATOR;
  else
    return MSG_UNDEFINED;
}

std::ostream& operator<< ( std::ostream& os, MSG_header_image_navig &h )
{
  os << "------------------------------------------------------" << std::endl
     << "-            MSG IMAGE NAVIGATION HEADER             -" << std::endl
     << "------------------------------------------------------" << std::endl;
  os << "Projection Name     : " << h.pname << std::endl;
  switch (h.get_projection_code())
  {
    case MSG_PROJECTION_GEOS:
      os << "Subsatellite Lat.   : " << h.subsatellite_longitude << std::endl;
      break;
    case MSG_PROJECTION_POLAR:
      os << "Projection plane    : " << h.projection_plane << std::endl
         << "Projection long.    : " << h.projection_longitude << std::endl;
      break;
    default:
      break;
  }
  os << "Column scale factor : " << (int_8) h.column_scaling_factor << std::endl
     << "Line scale factor   : " << (int_8) h.line_scaling_factor << std::endl
     << "Column offset       : " << (int_8) h.column_offset << std::endl
     << "Line Offset         : " << (int_8) h.line_offset << std::endl;
  return os;
}