GLIMSView Source Code Documentation

extensions/GtoSConverter.cpp

Go to the documentation of this file.

#include "GtoSConverter.h"

#include <qmessagebox.h>
#include <iostream>
#include <sstream>


GtoSConverter::GtoSConverter( GLIMSProject *proj ){
  mProj = proj;
  mImage = ImageInterface::getInterface(mProj);
  //mImage = new ImageInterface(mProj);
}

GtoSConverter::~GtoSConverter(){
  ImageInterface::destroyInterface();
}

void GtoSConverter::convert(std::string path){
    outputFilePath = path ;
    getSelectedShapeList();
    getSelectedShapes();
    std::string str = getConversion();
//    writeConversion(str, path);

   // ImageInterface::destroyInterface();
}



void GtoSConverter::convert(std::string path, std::string delim){
  outputFilePath = path ;
  getSelectedShapeList();
  getSelectedShapes();
  std::string str = getConversion(delim);
//  writeConversion(str, path);
}


void GtoSConverter::getSelectedShapeList(){


  //get the lindata
  GLIMSLineData* linedata = &(mProj->getDataset())->getLineData();

  //ensure it's not a point type
  if(linedata->getType() == "POINT"){
    QMessageBox::information( 0, "GLIMSview", "You must have a line or polygon selected." );

    return;
  }

  //get the vector of shapes
  std::vector<Shape *> shapes = linedata->getLLSet();

  //get the list of selected shapes
  SelectionSet& selSetRef = linedata->getSelSet();

  if(selSetRef.size()<1){
    QMessageBox::information( 0, "GLIMSview", "You must have at least one shape selected." );
    return;
  }

  //for each selected shape, add it to our list
  for (int i=0; i< (int)(selSetRef.size()); i++){
    selectedShapes.push_back(shapes[selSetRef[i].obj]);
  }

}

void GtoSConverter::getSelectedShapes(){

  if( (selectedShapes.size()>0)){

    for(unsigned int i = 0 ; i < selectedShapes.size() ; i++) {

      Shape &curShape = *selectedShapes[i];
      Line& curLine = dynamic_cast<Line&>(curShape);
      double hx = 0 ;
      double lx = 0 ; 
      double hy = 0 ; 
      double ly = 0 ;
      double hlat, llat, hlon, llon;

      std::vector<ImageCube::ThreeDPoint> points;

      for(unsigned int j = 0 ; j < curLine.size() ; j++) {

        double x = curLine[j].x;
        double y = curLine[j].y;

        ((mProj->getImage())->getFormat())->getXY(0, x, y);

        //a hack...
        //since GLIMS uses the upper left corner as a reference,
        //and I had assumed the middle of a pixel, we have to
        //shift everything over by half a pixel
        x-=0.5;
        y-=0.5;


        //check that we haven't gone over the edge
        //just set to the edge if we do
        if(x < 0) x=0.0;
        else if(x> mImage->getWidth()) x= (mImage->getWidth())-1;

        if(y <0) y=0.0;
        else if(y> mImage->getHeight()) y = (mImage->getHeight())-1;


        if(j==0){//if it's the first run, prime max and mins
          hx = lx = hlat = llat = x;
          hy = ly = hlon = llon = y;
        } else {//not the first run, compare max and mins
          if(x>hx) hx = x;
          else if(x<lx) lx = x;

          if(y>hy) hy = y;
          else if(y<ly) ly = y;

        }

        ImageCube::ThreeDPoint tpoint;
        tpoint.x = x;
        tpoint.y = y;
        tpoint.z = 0;
        points.push_back(tpoint);

      }

      //if we're close to the low end of a pixel, round down
      //highs are always rounded up
      if(ly!=0) ly -=0.25;
      if(lx!=0) lx -=0.25;

      ImageSelectionCube *tcube = new ImageSelectionCube((int)(lx), (int)(ly), (((int)(hx+0.5))-((int)lx))+1, (((int)(hy+0.5))-((int)ly))+1);

      llat = lx;
      llon = ly;
      hlat = hx;
      hlon = hy;

      //translate to lat/lon, use the first layer as a reference
      ((mProj->getImage())->getFormat())->getLL(0, llat, llon);
      ((mProj->getImage())->getFormat())->getLL(0, hlat, hlon);

      //add the layers, and setup the selection
      tcube->addBands(mImage->getCubeLatLong(llat, llon, 0, hlat, hlon, (((mProj->getImage())->getFormat())->getNumBands())-1));
      //tcube->addBands(mImage->getCube((int)(lx), (int)(ly), 0, (((int)(hx+0.5))-((int)lx))+1, (((int)(hy+0.5))-((int)ly))+1, (((mProj->getImage())->getFormat())->getNumBands())-1));

      tcube->setSelection(points);

      selectedCubes.push_back(tcube);

    }
  }
}

std::string GtoSConverter::getConversion(){
  std::string str = "";

  //just grab the string for each cube
  for(unsigned int i = 0 ; i < selectedCubes.size() ; i++) {
    str.append(selectedCubes[i]->getAsString());
    str.append("\n");
  }
  return str;
}

std::string GtoSConverter::getConversion(std::string delim){
  std::string str = "";
  std::stringstream ss;

  std::vector<string> names = mImage->getBandNames();

  ss << "X" << delim << " Y" << delim << " ";

  for(int j=0; j<mImage->getBands(); j++){
    ss << names[j]  << delim << " ";
  }

  ss << "\n";

  str = ss.str();

  std::ofstream f( outputFilePath.c_str() );
    if ( f ) {
      QString qstr = QString(str.c_str());

    f.write( (const char*)qstr, qstr.length() );
    f.flush();

    //just grab the string for each cube
    for(unsigned int i = 0 ; i < selectedCubes.size() ; i++) {
      //str.append(selectedCubes[i]->getAsString(delim));
      selectedCubes[i]->getAsString(delim, f) ;
      //str.append("\n");
    }
  }
  f.close();

  return str;
}

void GtoSConverter::writeConversion(std::string str, std::string path){
    std::ofstream f( path.c_str() );
    if ( f ) {
      QString qstr = QString(str.c_str());

    f.write( (const char*)qstr, qstr.length() );
    f.flush();
  }
  f.close();
}

---

Home | Search | Disclaimers & Privacy | Contact Us GLIMSView Maintainer: dsoltesz@usgs.gov