pihmgis.cpp

Go to the documentation of this file.

/***************************************************************************
  pihmgis.cpp 
  An Integrated GIS Interface to PennState Integrated Hydrologic Model (PIHM)

  -------------------
         begin                : [PluginDate]
         copyright            : [(C) Your Name and Date]
         email                : [Your Email]

 ***************************************************************************
 *                                                                         *
 *   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.                                   *
 *                                                                         *
 ***************************************************************************/
/*  $Id: plugin.cpp 5476 2006-05-19 09:53:08Z g_j_m $ */

//
// QGIS Specific includes
//

#include <qgisapp.h>
#include <qgisgui.h>
#include <qgsmaplayer.h>
#include <qgsrasterlayer.h>
#include "pihmgis.h"
//the gui subclass
#include "pihmgisgui.h"

//
// Qt4 Related Includes
//

#include <QToolBar>
#include <QMenuBar>
#include <QMessageBox>
#include <QMenu>
#include <QLineEdit>
#include <QAction>
#include <QApplication>
#include <QCursor>

#include <Q3PopupMenu>
#include <QToolButton>

//non qt includes
#include <iostream>


#include "RasterProcessing/FillPits/fillpits.h"
#include "RasterProcessing/FlowGrid/flowgrid.h"
#include "RasterProcessing/StreamGrid/streamgrid.h"
#include "RasterProcessing/LinkGrid/linkgrid.h"
#include "RasterProcessing/StreamPolyLine/streampolyline.h"
#include "RasterProcessing/CatchmentGrid/catchmentgrid.h"
#include "RasterProcessing/CatchmentPolygon/catchmentpolygon.h"

#include "VectorProcessing/shpFilesInputDialog/shpfilesinputdialog.h"
#include "VectorProcessing/polygonToPolyLineDialog/polygontopolylinedialog.h"
#include "VectorProcessing/simplifyLineDialog/simplifylinedialog.h"
#include "VectorProcessing/splitLineDialog/splitlinedialog.h"
#include "VectorProcessing/mergeFeaturesDialog/mergefeaturesdialog.h"

#include "DomainDecomposition/GenerateShapeTopology/generateshapetopology.h"
#include "DomainDecomposition/RunTriangle/runtriangle.h"
#include "DomainDecomposition/CreateTINs/createtins.h"

#include "ExtractData/mshFile/mshfile.h"
#include "ExtractData/attFile/attfile.h"
#include "ExtractData/rivFile/rivfile.h"
#include "ExtractData/paraFile/parafile.h"

#include "runPIHM/runpihm.h"
#include "sundials_types.h"   /* realtype, integertype, booleantype defination */
#include "cvode.h"           /* CVODE header file                             */
#include "cvode_spgmr.h"         /* CVSPGMR linear header file                    */
#include "sundials_smalldense.h"      /* use generic DENSE linear solver for "small"   */
#include "nvector_serial.h"  /* contains the definition of type N_Vector      */
#include "sundials_math.h"    /* contains UnitRoundoff, RSqrt, SQR functions   */
#include "cvode_dense.h"         /* CVDENSE header file                           */
#include "sundials_dense.h"

#include "Analysis/TimeSeries/timeseries.h"
#include "Analysis/SpatialPlot/spatialplot.h"

#include "Help/about/about.h"


#ifdef WIN32
#define QGISEXTERN extern "C" __declspec( dllexport )
#else
#define QGISEXTERN extern "C"
#endif

static const char * const sIdent = "$Id: plugin.cpp 5476 2006-05-19 09:53:08Z g_j_m $";
static const char * const sName = "PIHMgis";
static const char * const sDescription = "An Integrated GIS Interface to PennState Integrated Hydrologic Model (PIHM)";
static const char * const sPluginVersion = "Version 0.1";
static const QgisPlugin::PLUGINTYPE sPluginType = QgisPlugin::UI;

//
// THE FOLLOWING METHODS ARE MANDATORY FOR ALL PLUGINS
//

PIHMgis::PIHMgis(QgisApp * theQGisApp, 
                                       QgisIface * theQgisInterface):
                 QgisPlugin(sName,sDescription,sPluginVersion,sPluginType),
                 mQGisApp(theQGisApp), 
                 mQGisIface(theQgisInterface)
{
}

PIHMgis::~PIHMgis()
{

}

/*
 * Initialize the GUI interface for the plugin - this is only called once when the plugin is 
 * added to the plugin registry in the QGIS application.
 */
void PIHMgis::initGui()
{

  // Create the action for tool
  //mQActionPointer = new QAction(QIcon(":/pihmgis/pihmgis.png"),"PIHMgis", this);
  // Set the what's this text
  //mQActionPointer->setWhatsThis(tr("Replace this with a short description of the what the plugin does"));
  // Connect the action to the run
  //connect(mQActionPointer, SIGNAL(activated()), this, SLOT(run()));
  // Add the toolbar
  //mToolBarPointer = new QToolBar((QMainWindow *) mQGisApp, "PIHMgis");
  //mToolBarPointer->setLabel("PIHMgis");
  // Add the to the toolbar
  //mQGisIface->addToolBarIcon(mQActionPointer);
  //mQGisIface->addPluginMenu("&PIHMgis", mQActionPointer);

/*
QToolBar *mToolBarPointer2 = new QToolBar((QMainWindow *) mQGisApp, "&PIHMgis");
mToolBarPointer2->setLabel("PIHMgis");
  mToolBarPointer2->setMovable(1);

QMenu* rastMenu = new QMenu("&Raster Processing", mQGisIface->app());
  QAction *act1 = rastMenu->addAction(QIcon(":/pihmgis/pihmgis.png"), "Fill Pits", this, SLOT(runDoNothing()));
  QAction *act2 = rastMenu->addAction(QIcon(":/pihmgis/pihmgis.png"), "Flow Grid", this, SLOT(runDoNothing()));

mToolBarPointer2->addAction(act1);
mToolBarPointer2->addAction(act2);
*/
//mToolBar
/*
  QAction *action ;

  QToolBar* mToolBarPointer1 = new QToolBar((QMainWindow *) mQGisApp, "&PIHMgis");
  mToolBarPointer1->setLabel("PIHMgis");
  mToolBarPointer1->setMovable(1);

  QToolButton* rastTool = new QToolButton((QMainWindow *) mQGisApp);//,"  Raster Processing  ");
  rastTool->setText("&Raster Processing");
  rastTool->setPopupMode(QToolButton::MenuButtonPopup);
*/
  QMenu* rastMenu = new QMenu("&Raster-Processing", mQGisIface->app());
  rastMenu->addAction("Run All", this, SLOT(runFillPits()));
  rastMenu->addSeparator();
//      rastMenu->addAction(QIcon(":/pihmgis/pihmgis.png"), "Fill Pits", this, SLOT(runFillPits()));
  rastMenu->addAction( "Fill Pits", this, SLOT(runFillPits()));
  rastMenu->addAction( "Flow Grid", this, SLOT(runFlowGrid()));
  rastMenu->addAction( "Stream Grid", this, SLOT(runStreamGrid()));
  rastMenu->addAction( "Link Grid", this, SLOT(runLinkGrid()));
  rastMenu->addAction( "Stream Polyline", this, SLOT(runStreamPolyLine()));
  rastMenu->addAction( "Catchment Grid", this, SLOT(runCatchmentGrid()));
  rastMenu->addAction( "Catchment Polygon", this, SLOT(runCatchmentPolygon()));
  
  QAction *rastAction = new QAction(" Raster Processing ", this);
  rastAction->setMenu(rastMenu);
  mQGisIface->addToolBarIcon(rastAction);

  /*
  action=mToolBarPointer1->addWidget(rastTool);  


  action = mToolBarPointer1->insertSeparator(action);
  
  QToolButton* vectTool = new QToolButton((QMainWindow *) mQGisApp);//,"  Raster Processing  ");
  vectTool->setText("&Vector Processing");
  vectTool->setPopupMode(QToolButton::MenuButtonPopup);
*/
  QMenu* vectMenu = new QMenu("&Vector-Processing", mQGisIface->app());
  vectMenu->addAction( "Run All", this, SLOT(runShpFileInput()));
  vectMenu->addSeparator();
  vectMenu->addAction( "Input Files", this, SLOT(runShpFileInput()));
  vectMenu->addAction( "Polygon To Line", this, SLOT(runPolygonToPolyline()));
  vectMenu->addAction( "Simplify Line", this, SLOT(runSimplifyLine()));
  vectMenu->addAction( "Split Line", this, SLOT(runSplitLine()));
  vectMenu->addAction( "Vector Merge", this, SLOT(runVectorMerge()));

  QAction *vectAction = new QAction(" Vector Processing ", this);
  vectAction->setMenu(vectMenu);
  mQGisIface->addToolBarIcon(vectAction);


  QMenu* domainMenu = new QMenu("&Domain-Decomposition", mQGisIface->app());
  domainMenu->addAction( "Run All", this, SLOT(runGenerateShapeTopology()));
  domainMenu->addSeparator();
  domainMenu->addAction( "Read ShapeTopology", this, SLOT(runGenerateShapeTopology()));
  domainMenu->addAction( "Run Triangle", this, SLOT(runRunTriangle()));
  domainMenu->addAction( "TIN Generation", this, SLOT(runCreateTINs()));

  QAction *domainAction = new QAction("Domain Decomposition", this);
  domainAction->setMenu(domainMenu);
  mQGisIface->addToolBarIcon(domainAction);

  QMenu* xtractMenu = new QMenu("&Extract-Data", mQGisIface->app());
  xtractMenu->addAction( "Run All", this, SLOT(runMshFile()));
  xtractMenu->addSeparator();
  xtractMenu->addAction( "Generate MeshFile", this, SLOT(runMshFile()));
  xtractMenu->addAction( "Generate AttFile", this, SLOT(runAttFile()));
  xtractMenu->addAction( "Generate RivFile", this, SLOT(runRivFile()));
  xtractMenu->addAction( "Generate ParaFile", this, SLOT(runParaFile()));

  QAction *xtractAction = new QAction("&DataModel Loader", this);
  xtractAction->setMenu(xtractMenu);
  mQGisIface->addToolBarIcon(xtractAction);

  QMenu* pihmMenu = new QMenu("&PIHM", mQGisIface->app());
  pihmMenu->addAction( "Run PIHM", this, SLOT(runPIHM()));
  //pihmMenu->addAction( "Run PIHM", this, SLOT(runDoNothing()));
  QAction *pihmAction = new QAction("      PIHM      ", this);
  pihmAction->setMenu(pihmMenu);
  mQGisIface->addToolBarIcon(pihmAction);

  QMenu* analysisMenu = new QMenu("&Analysis", mQGisIface->app());
  //analysisMenu->addAction("TimeSeries Plots", this, SLOT(runTimeSeriesPlots()));
  analysisMenu->addAction("Spatial Plots", this, SLOT(runSpatialPlots()));

  QAction *analysisAction = new QAction("      Analysis      ", this);
  analysisAction->setMenu(analysisMenu);
  mQGisIface->addToolBarIcon(analysisAction);

  QMenu* helpMenu = new QMenu("&HELP", mQGisIface->app());
  helpMenu->addAction( "HelpContents", this, SLOT(runDoNothing()));
  helpMenu->addAction( "PIHMgis Homepage", this, SLOT(runDoNothing()));
  helpMenu->addAction( "About", this, SLOT(runAbout()));
  QAction *helpAction = new QAction("      Help      ", this);
  helpAction->setMenu(helpMenu);
  mQGisIface->addToolBarIcon(helpAction);

/*  
  vectTool->setMenu(vectMenu);

  action = mToolBarPointer1->insertWidget(action, vectTool);
*/
  //mToolBarPointer1->addSeparator();
  //mQGisIface->addToolBarIcon(mQActionPointer);
  //mQGisIface->addPluginMenu("&PIHMgis", mQActionPointer);

}
//method defined in interface

void PIHMgis::help()
{
  //implement me!
}

// Slot called when the menu item is activated
// If you created more menu items / toolbar buttons in initiGui, you should 
// create a separate handler for each action - this single run() method will
// not be enough
void PIHMgis::run()
{
  PIHMgisGui *myPluginGui=new PIHMgisGui(mQGisApp, QgisGui::ModalDialogFlags);
  //listen for when the layer has been made so we can draw it
  connect(myPluginGui, SIGNAL(drawRasterLayer(QString)), this, SLOT(drawRasterLayer(QString)));
  connect(myPluginGui, SIGNAL(drawVectorLayer(QString,QString,QString)), this, SLOT(drawVectorLayer(QString,QString,QString)));
  myPluginGui->show();
}

// Unload the plugin by cleaning up the GUI
void PIHMgis::unload()
{
  // remove the GUI
  mQGisIface->removePluginMenu("&PIHMgis",mQActionPointer);
  mQGisIface->removeToolBarIcon(mQActionPointer);
  delete mQActionPointer;
}

void PIHMgis::runDoNothing(){

}
void PIHMgis::runFillPits(){
        fillpitsDlg* Dlg = new fillpitsDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}

void PIHMgis::runFlowGrid(){
        FlowGridDlg* Dlg = new FlowGridDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}

void PIHMgis::runStreamGrid(){
        StreamGridDlg* Dlg = new StreamGridDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}

void PIHMgis::runLinkGrid(){
        LinkGridDlg* Dlg = new LinkGridDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}

void PIHMgis::runStreamPolyLine(){
        StreamPolyLineDlg* Dlg = new StreamPolyLineDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}

void PIHMgis::runCatchmentGrid(){
        CatchmentGridDlg* Dlg = new CatchmentGridDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}

void PIHMgis::runCatchmentPolygon(){
        CatchmentPolygonDlg* Dlg = new CatchmentPolygonDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}


void PIHMgis::runShpFileInput(){
        shpFilesInputDialogDlg* Dlg = new shpFilesInputDialogDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runPolygonToPolyline(){
        polygonToPolyLineDialogDlg* Dlg = new polygonToPolyLineDialogDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runSimplifyLine(){
        simplifyLineDialogDlg* Dlg = new simplifyLineDialogDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runSplitLine(){
        splitLineDialogDlg* Dlg = new splitLineDialogDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runVectorMerge(){
        mergeFeaturesDialogDlg* Dlg = new mergeFeaturesDialogDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runGenerateShapeTopology(){
        generateShapeTopologyDlg* Dlg = new generateShapeTopologyDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runRunTriangle(){
        runTriangleDlg* Dlg = new runTriangleDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runCreateTINs(){
        createTINsDlg* Dlg = new createTINsDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}
void PIHMgis::runMshFile(){
        mshFileDlg* Dlg = new mshFileDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runAttFile(){
        attFileDlg* Dlg = new attFileDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runRivFile(){
        rivFileDlg* Dlg = new rivFileDlg;
        Dlg->show();
}

void PIHMgis::runParaFile(){
        paraFileDlg* Dlg = new paraFileDlg; //(NULL,NULL,TRUE,0); //,"/opt3/helpDialog/poly.html","Help Caption");
        //Dlg->setApplicationPointer(mQGisIface->app());
        Dlg->show();
}

void PIHMgis::runPIHM(){
        runPIHMDlg* Dlg = new runPIHMDlg;
        Dlg->show();
}
void PIHMgis::runTimeSeriesPlots(){
//      timeSeriesDlg *Dlg = new timeSeriesDlg;
//      Dlg->show();
}

void PIHMgis::runSpatialPlots(){
        spatialPlotDlg* Dlg = new spatialPlotDlg;
        Dlg->setApplicationPointer(mQGisApp);
        Dlg->show();
}
void PIHMgis::runAbout(){
        AboutDlg* Dlg = new AboutDlg; 
        Dlg->show();
}
//
//                  END OF MANDATORY PLUGIN METHODS
//
//
// The following methods are provided to demonstrate how you can 
// load a vector or raster layer into the main gui. Please delete
// if you are not intending to use these. Note also that there are
// ways in which layers can be loaded.
//

//layer
void PIHMgis::drawRasterLayer(QString theQString)
{
  mQGisIface->addRasterLayer(theQString);
}

// dialog when it as finished creating a layer. It needs to be given 
// vectorLayerPath, baseName, providerKey ("ogr" or "postgres");
void PIHMgis::drawVectorLayer(QString thePathNameQString, QString theBaseNameQString, QString theProviderQString)
{
  mQGisIface->addVectorLayer( thePathNameQString, theBaseNameQString, theProviderQString);
}


//
//
//  THE FOLLOWING CODE IS AUTOGENERATED BY THE PLUGIN BUILDER SCRIPT
//    YOU WOULD NORMALLY NOT NEED TO MODIFY THIS, AND YOUR PLUGIN
//      MAY NOT WORK PROPERLY IF YOU MODIFY THIS INCORRECTLY
//
//


// Class factory to return a new instance of the plugin class
QGISEXTERN QgisPlugin * classFactory(QgisApp * theQGisAppPointer, QgisIface * theQgisInterfacePointer)
{
  return new PIHMgis(theQGisAppPointer, theQgisInterfacePointer);
}
// Return the name of the plugin - note that we do not user class members as
// the class may not yet be insantiated when this method is called.
QGISEXTERN QString name()
{
  return sName;
}

// Return the description
QGISEXTERN QString description()
{
  return sDescription;
}

// Return the type (either UI or MapLayer plugin)
QGISEXTERN int type()
{
  return sPluginType;
}

// Return the version number for the plugin
QGISEXTERN QString version()
{
  return sPluginVersion;
}

// Delete ourself
QGISEXTERN void unload(QgisPlugin * thePluginPointer)
{
  delete thePluginPointer;
}

---