reproject.cpp File Reference

#include <iostream>
#include <iterator>
#include "base/data_types/ossimFilename.h"
#include "base/data_types/ossimString.h"
#include "imaging/factory/ossimImageHandlerRegistry.h"
#include "imaging/formats/ossimImageHandler.h"
#include "imaging/formats/ossimImageFileWriter.h"
#include "imaging/factory/ossimImageWriterFactoryRegistry.h"
#include "base/data_types/ossimGpt.h"
#include "projections/ossimProjection.h"
#include "projections/map_projections/ossimUtmProjection.h"
#include "projections/factory/ossimProjectionFactoryRegistry.h"
#include "imaging/tile_sources/ossimImageRenderer.h"
#include "init/ossimInit.h"

Go to the source code of this file.

Functions
void	usage ()
void	printOutputTypes ()
ossimProjection *	newUtmView (const ossimGpt &centerGround, const ossimDpt &metersPerPixel)
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 53 of file reproject.cpp.

References ossimConnectableObject::connectMyInputTo(), ossimProjectionFactoryRegistry::createProjection(), ossimImageWriterFactoryRegistry::createWriter(), ossimImageFileWriter::execute(), ossimImageHandler::getBoundingRect(), ossimImageHandler::getImageGeometry(), ossimProjection::getMetersPerPixel(), ossimInit::initialize(), ossimProjectionFactoryRegistry::instance(), ossimImageWriterFactoryRegistry::instance(), ossimInit::instance(), ossimImageHandlerRegistry::instance(), ossimProjection::lineSampleToWorld(), ossimIrect::midPoint(), newUtmView(), ossimImageHandlerRegistry::open(), ossimImageFileWriter::setFilename(), ossimImageRenderer::setView(), and usage().

{
   ossimInit::instance()->initialize(argc, argv);

   if(argc!=4)
   {
      usage();
   }
   else
   {
      // try to open up the passed in image
      //
      ossimImageHandler *handler   = ossimImageHandlerRegistry::instance()->open(ossimFilename(argv[2]));

      // try to create a writer for the output image type.
      //
      ossimImageFileWriter* writer = ossimImageWriterFactoryRegistry::instance()->createWriter(ossimString(argv[1]));


      if(!handler)
      {
         cout << "Unable to open input image: "<< argv[2] << endl;
         return 1;
      }
      if(!writer)
      {
         cout << "Unable to create writer of type: " << argv[1] << endl;
         delete handler;
         return 1;
      }

      ossimKeywordlist geom;
      handler->getImageGeometry(geom);

      // grab a projection if it exists
      //
      ossimProjection* inputProjection = ossimProjectionFactoryRegistry::instance()->createProjection(geom);

      if (!inputProjection)
      {
         cout << "the input image has no input projection and can't be reprojected" << endl;
         delete handler;
         delete writer;
         return 1;
      }

      // --------------------- SETUP The Resampleing process -------------
      //
      // the renderer is the resampler that fits into the chain.
      // we will set it to a UTM view. We first get the inputs geometry and
      // then get the scale from its projector.
      //
      // now lets set up the renderer
      ossimImageRenderer* renderer = new ossimImageRenderer;

      //
      // get the center ground by running the center pixel
      // through the handler's projection object
      ossimIrect bounds = handler->getBoundingRect();
      ossimGpt centerGround;
      inputProjection->lineSampleToWorld(bounds.midPoint(), centerGround);

      // now set the view to a new UTM view.
      // I pass true in to tell the renderer that it owns the
      // projection and will be responsible for deleting
      //
      renderer->setView(newUtmView(centerGround,
                                   inputProjection->getMetersPerPixel()),
                        true);

      // connect the renderer to the handler
      renderer->connectMyInputTo(handler);
      
      // specify the output file name
      writer->setFilename(ossimFilename(argv[3]));

      // now connect the writer to the renderer
      writer->connectMyInputTo(0, renderer);

      // execute the reprojection process.
      writer->execute();

      delete writer;
      delete handler;
   }

   return 0;
}

newUtmView()

ossimProjection * newUtmView	(	const ossimGpt &	centerGround,
		const ossimDpt &	metersPerPixel
	)

Definition at line 160 of file reproject.cpp.

References ossimMapProjection::setMetersPerPixel(), ossimUtmProjection::setZone(), ossimDpt::x, and ossimDpt::y.

Referenced by main().

{
   ossimUtmProjection* utm = new ossimUtmProjection;

   // we will make it a square pixel in meters
   double averageGsd = (metersPerPixel.x + metersPerPixel.y)*.5;
   utm->setZone(centerGround);
   utm->setMetersPerPixel(ossimDpt(metersPerPixel));

   
   return utm;
}

printOutputTypes()

void printOutputTypes

(

)

Definition at line 150 of file reproject.cpp.

References ossimImageWriterFactoryRegistry::getImageTypeList(), and ossimImageWriterFactoryRegistry::instance().

Referenced by usage().

{
    std::vector<ossimString> outputType;
   
    ossimImageWriterFactoryRegistry::instance()->getImageTypeList(outputType);
    std::copy(outputType.begin(),
              outputType.end(),
              std::ostream_iterator<ossimString>(std::cout, "\n"));
}

usage()

void usage

(

)

Definition at line 143 of file reproject.cpp.

References printOutputTypes().

Referenced by main().

{
   cout << "repoject <output_type> <input filename> <output filename>" << endl
        << "where output types are: " << endl;
   printOutputTypes();
}