icbklt.c

来自「spiht for linux this is used to decod an」· C语言 代码 · 共 104 行

/*
 *
 * QccPack: Quantization, compression, and coding utilities
 * Copyright (C) 1997-2009  James E. Fowler
 * 
 * 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., 675 Mass Ave, Cambridge, MA 02139, USA.
 *
 */


#include "icbklt.h"

#define USG_STRING "[-p %: %d:num_pcs] %s:kltfile %s:icbfile1 %s:icbfile2"

QccHYPklt Transform;
QccIMGImageCube Image;
QccString InputFilename;
QccString OutputFilename;

int NumPCsSpecified = 0;
int NumPCs;


int main(int argc, char *argv[])
{
  QccInit(argc, argv);

  QccHYPkltInitialize(&Transform);
  QccIMGImageCubeInitialize(&Image);
  
  if (QccParseParameters(argc, argv,
                         USG_STRING,
                         &NumPCsSpecified,
                         &NumPCs,
                         Transform.filename,
                         InputFilename,
                         OutputFilename))
    QccErrorExit();

  if (QccHYPkltRead(&Transform))
    {
      QccErrorAddMessage("%s: Error calling QccHYPkltRead()",
                         argv[0]);
      QccErrorExit();
    }

  if (NumPCsSpecified)
    {
      if ((NumPCs <= 0) || (NumPCs > Transform.num_bands))
        {
          QccErrorAddMessage("%s: Invalid number of PCs (%d) specified for a KLT with %d bands",
                             argv[0],
                             NumPCs,
                             Transform.num_bands);
          QccErrorExit();
        }
    }
  else
    NumPCs = Transform.num_bands;

  QccStringCopy(Image.filename, InputFilename);
  if (QccIMGImageCubeRead(&Image))
    {
      QccErrorAddMessage("%s: Error calling QccIMGImageCubeRead()",
                         argv[0]);
      QccErrorExit();
    }

  if (QccHYPkltTransform(&Image,
                         &Transform,
                         NumPCs))
    {
      QccErrorAddMessage("%s: Error calling QccHYPkltTransform()",
                         argv[0]);
      QccErrorExit();
    }

  QccStringCopy(Image.filename, OutputFilename);
  if (QccIMGImageCubeWrite(&Image))
    {
      QccErrorAddMessage("%s: Error calling QccIMGImageCubeWrite()",
                         argv[0]);
      QccErrorExit();
    }

  QccHYPkltFree(&Transform);
  QccIMGImageCubeFree(&Image);

  QccExit;
}