# -*- coding: utf-8 -*- """ *************************************************************************** doMerge.py --------------------- Date : June 2010 Copyright : (C) 2010 by Giuseppe Sucameli Email : brush dot tyler at gmail dot com *************************************************************************** * * * 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. * * * *************************************************************************** """ __author__ = 'Giuseppe Sucameli' __date__ = 'June 2010' __copyright__ = '(C) 2010, Giuseppe Sucameli' # This will get replaced with a git SHA1 when you do a git archive __revision__ = '$Format:%H$' from PyQt4.QtCore import * from PyQt4.QtGui import * from qgis.core import * from qgis.gui import * from ui_widgetMerge import Ui_GdalToolsWidget as Ui_Widget from widgetPluginBase import GdalToolsBasePluginWidget as BasePluginWidget import GdalTools_utils as Utils class GdalToolsDialog(QWidget, Ui_Widget, BasePluginWidget): def __init__(self, iface): QWidget.__init__(self) self.iface = iface self.setupUi(self) BasePluginWidget.__init__(self, self.iface, "gdal_merge.py") self.inSelector.setType( self.inSelector.FILE ) self.outSelector.setType( self.outSelector.FILE ) self.recurseCheck.hide() # use this for approx. previous UI #self.creationOptionsTable.setType(QgsRasterFormatSaveOptionsWidget.Table) self.outputFormat = Utils.fillRasterOutputFormat() self.extent = None self.setParamsStatus( [ (self.inSelector, SIGNAL("filenameChanged()")), (self.outSelector, SIGNAL("filenameChanged()")), (self.noDataSpin, SIGNAL("valueChanged(int)"), self.noDataCheck), (self.inputDirCheck, SIGNAL("stateChanged(int)")), (self.recurseCheck, SIGNAL("stateChanged(int)"), self.inputDirCheck), ( self.separateCheck, SIGNAL( "stateChanged( int )" ) ), ( self.pctCheck, SIGNAL( "stateChanged( int )" ) ), ( self.intersectCheck, SIGNAL( "stateChanged( int )" ) ) ] ) self.connect(self.inSelector, SIGNAL("selectClicked()"), self.fillInputFilesEdit) self.connect(self.outSelector, SIGNAL("selectClicked()"), self.fillOutputFileEdit) self.connect(self.intersectCheck, SIGNAL("toggled(bool)"), self.refreshExtent) self.connect(self.inputDirCheck, SIGNAL("stateChanged( int )"), self.switchToolMode) self.connect(self.inSelector, SIGNAL("filenameChanged()"), self.refreshExtent) def switchToolMode(self): self.recurseCheck.setVisible( self.inputDirCheck.isChecked() ) self.inSelector.clear() if self.inputDirCheck.isChecked(): self.inFileLabel = self.label.text() self.label.setText( QCoreApplication.translate( "GdalTools", "&Input directory" ) ) QObject.disconnect(self.inSelector, SIGNAL("selectClicked()"), self.fillInputFilesEdit) QObject.connect(self.inSelector, SIGNAL("selectClicked()"), self.fillInputDir) else: self.label.setText( self.inFileLabel ) QObject.connect(self.inSelector, SIGNAL("selectClicked()"), self.fillInputFilesEdit) QObject.disconnect(self.inSelector, SIGNAL("selectClicked()"), self.fillInputDir) def fillInputFilesEdit(self): lastUsedFilter = Utils.FileFilter.lastUsedRasterFilter() files = Utils.FileDialog.getOpenFileNames(self, self.tr( "Select the files to Merge" ), Utils.FileFilter.allRastersFilter(), lastUsedFilter ) if files.isEmpty(): return Utils.FileFilter.setLastUsedRasterFilter(lastUsedFilter) self.inSelector.setFilename(files) def refreshExtent(self): files = self.getInputFileNames() self.intersectCheck.setEnabled( files.count() > 1 ) if not self.intersectCheck.isChecked(): self.someValueChanged() return if files.count() < 2: self.intersectCheck.setChecked( False ) return self.extent = self.getIntersectedExtent( files ) if self.extent == None: QMessageBox.warning( self, self.tr( "Error retrieving the extent" ), self.tr( 'GDAL was unable to retrieve the extent from any file. \nThe "Use intersected extent" option will be unchecked.' ) ) self.intersectCheck.setChecked( False ) return elif self.extent.isEmpty(): QMessageBox.warning( self, self.tr( "Empty extent" ), self.tr( 'The computed extent is empty. \nDisable the "Use intersected extent" option to have a nonempty output.' ) ) self.someValueChanged() def fillOutputFileEdit(self): lastUsedFilter = Utils.FileFilter.lastUsedRasterFilter() outputFile = Utils.FileDialog.getSaveFileName(self, self.tr( "Select where to save the Merge output" ), Utils.FileFilter.allRastersFilter(), lastUsedFilter ) if outputFile.isEmpty(): return Utils.FileFilter.setLastUsedRasterFilter(lastUsedFilter) self.outputFormat = Utils.fillRasterOutputFormat( lastUsedFilter, outputFile ) self.outSelector.setFilename( outputFile ) def fillInputDir( self ): inputDir = Utils.FileDialog.getExistingDirectory( self, self.tr( "Select the input directory with files to Merge" )) if inputDir.isEmpty(): return self.inSelector.setFilename( inputDir ) def getArguments(self): arguments = QStringList() if self.intersectCheck.isChecked(): if self.extent != None: arguments << "-ul_lr" arguments << str( self.extent.xMinimum() ) arguments << str( self.extent.yMaximum() ) arguments << str( self.extent.xMaximum() ) arguments << str( self.extent.yMinimum() ) if self.noDataCheck.isChecked(): arguments << "-n" arguments << str(self.noDataSpin.value()) if self.separateCheck.isChecked(): arguments << "-separate" if self.pctCheck.isChecked(): arguments << "-pct" if self.creationGroupBox.isChecked(): for opt in self.creationOptionsTable.options(): arguments << "-co" arguments << opt outputFn = self.getOutputFileName() if not outputFn.isEmpty(): arguments << "-of" arguments << self.outputFormat arguments << "-o" arguments << outputFn arguments << self.getInputFileNames() return arguments def getOutputFileName(self): return self.outSelector.filename() def getInputFileName(self): if self.inputDirCheck.isChecked(): return self.inSelector.filename() return self.inSelector.filename().split(",", QString.SkipEmptyParts) def getInputFileNames(self): if self.inputDirCheck.isChecked(): return Utils.getRasterFiles( self.inSelector.filename(), self.recurseCheck.isChecked() ) return self.inSelector.filename().split(",", QString.SkipEmptyParts) def addLayerIntoCanvas(self, fileInfo): self.iface.addRasterLayer(fileInfo.filePath()) def getIntersectedExtent(self, files): res = None for fileName in files: if res == None: res = Utils.getRasterExtent( self, fileName ) continue rect2 = Utils.getRasterExtent( self, fileName ) if rect2 == None: continue res = res.intersect( rect2 ) if res.isEmpty(): break return res