Exemplo n.º 1
0
    def __init__(self, iface):
        """Constructor.

        :param iface: An interface instance that will be passed to this class
            which provides the hook by which you can manipulate the QGIS
            application at run time.
        :type iface: QgsInterface
        """
        # Save reference to the QGIS interface
        self.iface = iface
        # initialize plugin directory
        self.plugin_dir = os.path.dirname(__file__)
        # initialize locale
        locale = QSettings().value('locale/userLocale')[0:2]
        locale_path = os.path.join(
            self.plugin_dir,
            'i18n',
            'BlockModel_{}.qm'.format(locale))

        if os.path.exists(locale_path):
            self.translator = QTranslator()
            self.translator.load(locale_path)

            if qVersion() > '4.3.3':
                QCoreApplication.installTranslator(self.translator)

        # Create the dialog (after translation) and keep reference
        self.dlg = BlockModelDialog()

        # Declare instance attributes
        self.actions = []
        self.menu = self.tr(u'&Block Model')
        # TODO: We are going to let the user set this up in a future iteration
        self.toolbar = self.iface.addToolBar(u'BlockModel')
        self.toolbar.setObjectName(u'BlockModel')
        
        ## Initialize select block model button
        self.dlg.lineEdit.clear()
        self.dlg.pushButton.clicked.connect(self.select_blockmodel_file)
        self.dlg.outputFilePath.clear()
        self.dlg.outputFileButton.clicked.connect(self.select_output_path)
        self.dlg.buildExposedSurfaceButton.clicked.connect(self.buildExposedSurface)
Exemplo n.º 2
0
class BlockModel:
    """QGIS Plugin Implementation."""

    def __init__(self, iface):
        """Constructor.

        :param iface: An interface instance that will be passed to this class
            which provides the hook by which you can manipulate the QGIS
            application at run time.
        :type iface: QgsInterface
        """
        # Save reference to the QGIS interface
        self.iface = iface
        # initialize plugin directory
        self.plugin_dir = os.path.dirname(__file__)
        # initialize locale
        locale = QSettings().value('locale/userLocale')[0:2]
        locale_path = os.path.join(
            self.plugin_dir,
            'i18n',
            'BlockModel_{}.qm'.format(locale))

        if os.path.exists(locale_path):
            self.translator = QTranslator()
            self.translator.load(locale_path)

            if qVersion() > '4.3.3':
                QCoreApplication.installTranslator(self.translator)

        # Create the dialog (after translation) and keep reference
        self.dlg = BlockModelDialog()

        # Declare instance attributes
        self.actions = []
        self.menu = self.tr(u'&Block Model')
        # TODO: We are going to let the user set this up in a future iteration
        self.toolbar = self.iface.addToolBar(u'BlockModel')
        self.toolbar.setObjectName(u'BlockModel')
        
        ## Initialize select block model button
        self.dlg.lineEdit.clear()
        self.dlg.pushButton.clicked.connect(self.select_blockmodel_file)
        self.dlg.outputFilePath.clear()
        self.dlg.outputFileButton.clicked.connect(self.select_output_path)
        self.dlg.buildExposedSurfaceButton.clicked.connect(self.buildExposedSurface)
        

    # noinspection PyMethodMayBeStatic
    def tr(self, message):
        """Get the translation for a string using Qt translation API.

        We implement this ourselves since we do not inherit QObject.

        :param message: String for translation.
        :type message: str, QString

        :returns: Translated version of message.
        :rtype: QString
        """
        # noinspection PyTypeChecker,PyArgumentList,PyCallByClass
        return QCoreApplication.translate('BlockModel', message)


    def add_action(
        self,
        icon_path,
        text,
        callback,
        enabled_flag=True,
        add_to_menu=True,
        add_to_toolbar=True,
        status_tip=None,
        whats_this=None,
        parent=None):
        """Add a toolbar icon to the toolbar.

        :param icon_path: Path to the icon for this action. Can be a resource
            path (e.g. ':/plugins/foo/bar.png') or a normal file system path.
        :type icon_path: str

        :param text: Text that should be shown in menu items for this action.
        :type text: str

        :param callback: Function to be called when the action is triggered.
        :type callback: function

        :param enabled_flag: A flag indicating if the action should be enabled
            by default. Defaults to True.
        :type enabled_flag: bool

        :param add_to_menu: Flag indicating whether the action should also
            be added to the menu. Defaults to True.
        :type add_to_menu: bool

        :param add_to_toolbar: Flag indicating whether the action should also
            be added to the toolbar. Defaults to True.
        :type add_to_toolbar: bool

        :param status_tip: Optional text to show in a popup when mouse pointer
            hovers over the action.
        :type status_tip: str

        :param parent: Parent widget for the new action. Defaults None.
        :type parent: QWidget

        :param whats_this: Optional text to show in the status bar when the
            mouse pointer hovers over the action.

        :returns: The action that was created. Note that the action is also
            added to self.actions list.
        :rtype: QAction
        """

        icon = QIcon(icon_path)
        action = QAction(icon, text, parent)
        action.triggered.connect(callback)
        action.setEnabled(enabled_flag)

        if status_tip is not None:
            action.setStatusTip(status_tip)

        if whats_this is not None:
            action.setWhatsThis(whats_this)

        if add_to_toolbar:
            self.toolbar.addAction(action)

        if add_to_menu:
            self.iface.addPluginToMenu(
                self.menu,
                action)

        self.actions.append(action)

        return action

    def initGui(self):
        """Create the menu entries and toolbar icons inside the QGIS GUI."""

        icon_path = ':/plugins/BlockModel/icon.png'
        self.add_action(
            icon_path,
            text=self.tr(u'Block Model'),
            callback=self.run,
            parent=self.iface.mainWindow())


    def unload(self):
        """Removes the plugin menu item and icon from QGIS GUI."""
        for action in self.actions:
            self.iface.removePluginMenu(
                self.tr(u'&Block Model'),
                action)
            self.iface.removeToolBarIcon(action)
        # remove the toolbar
        del self.toolbar
        
    def select_output_path(self):
        fileEnding = '.asc'
        filename = QFileDialog.getSaveFileName(self.dlg, "Select output file ","", '*'+fileEnding)
        if filename == '':
        	return
        if filename[-4:] != '.asc':
            filename = filename+'.asc'
        self.dlg.outputFilePath.setText(filename)
    
    def select_blockmodel_file(self):
        filename = QFileDialog.getOpenFileName(self.dlg, "Select block model file ","", '*.csv')
        if filename == '':
        	return
        self.dlg.lineEdit.setText(filename)
        f=open(filename)
        lines = list()
        numCols = 0
        for i in range(11):
            line = f.next().strip().split(',')
            if len(line) > numCols:
                numCols = len(line)
            lines.append(line)
        self.dlg.tableWidget.setColumnCount(numCols)
        self.dlg.tableWidget.setRowCount(len(lines)-1)
        self.dlg.tableWidget.setHorizontalHeaderLabels(lines[0])
        for i in range(1,len(lines)):
            for j in range(len(lines[i])):
                item = QTableWidgetItem(lines[i][j])
                self.dlg.tableWidget.setItem(i-1,j,item)
        self.dlg.tableWidget.resizeColumnsToContents()
        self.dlg.tableWidget.resizeRowsToContents()
        
        self.dlg.eastingColumnSelection.clear()
        self.dlg.northingColumnSelection.clear()
        self.dlg.elevationColumnSelection.clear()
        self.dlg.valueColumnSelection.clear()

        self.dlg.eastingColumnSelection.addItems(lines[0])
        self.dlg.northingColumnSelection.addItems(lines[0])
        self.dlg.elevationColumnSelection.addItems(lines[0])
        self.dlg.valueColumnSelection.addItems(lines[0])
        
        self.dlg.eastingColumnSelection.setCurrentIndex(cni.guessEasting(lines[0]))
        self.dlg.northingColumnSelection.setCurrentIndex(cni.guessNorthing(lines[0]))
        self.dlg.elevationColumnSelection.setCurrentIndex(cni.guessElevation(lines[0]))
        self.dlg.valueColumnSelection.setCurrentIndex(3)
        
    def buildExposedSurface(self):
        self.dlg.statusBlockModelExtents.setValue(0)
        self.dlg.statusBlockModelLoading.setValue(0)
        self.dlg.statusIntersections.setValue(0)
        self.dlg.statusOutput.setValue(0)
        xCol = self.dlg.eastingColumnSelection.currentIndex()
        yCol = self.dlg.northingColumnSelection.currentIndex()
        zCol = self.dlg.elevationColumnSelection.currentIndex()
        valueCol = self.dlg.valueColumnSelection.currentIndex()
        blockModelFile = self.dlg.lineEdit.text()
        outputFile = self.dlg.outputFilePath.text()
        layers = self.iface.legendInterface().layers()
        for layer in layers:
            if str(layer.name()) == str(self.dlg.comboBox.currentText()):
				selectedLayer = layer
        sNODATA_value = -1
        print('Loading ASC surface.')
        sncols,snrows,sxllcorner,syllcorner,scellsize = cni.dtmProperties(selectedLayer)
        print('Finding block model limits.')
        bminX,bminY,bminZ,bnCols,bnRows,bnLvls,bcolSize,browSize,blvlSize = cni.blockModelLimits(blockModelFile,self.dlg.statusBlockModelExtents)
        print('Initializing 3d block model.')
        blockModel = cni.make3dList(bnCols,bnRows,bnLvls)
        print('Populating block model.')
        blockModelFileLength = cni.file_len(blockModelFile)
        with open(blockModelFile, 'r') as f:
            counter = 0
            status = 0
            for line in f:
                counter += 1
                if counter == 1:
                    continue
                incStatus = cni.percentStatus(counter,blockModelFileLength)
                if incStatus > status:
                    status = incStatus
                    self.dlg.statusBlockModelLoading.setValue(status)
                lineArray = line.split(",")
                xValue = float(lineArray[xCol])
                yValue = float(lineArray[yCol])
                zValue = float(lineArray[zCol])
                icol,irow,ilvl = cni.xyz2index(xValue,yValue,zValue,bminX,bminY,bminZ,bcolSize,browSize,blvlSize)
                try:
                    blockModel[icol][irow][ilvl] = float(lineArray[valueCol])
                except:
                    print(icol,irow,ilvl,float(lineArray[valueCol]))
        print('Building ASC')
        blankASC = cni.make2dList(sncols,snrows)
        ascCounter = 0
        ascArrayLength = len(blankASC)*len(blankASC[0])
        for i in range(len(blankASC)):
            for j in range(len(blankASC[0])):
                ascCounter += 1
                status = 0
                incStatus = cni.percentStatus(ascCounter,ascArrayLength)
                if incStatus > status:
                    status = incStatus
                    self.dlg.statusIntersections.setValue(status)
                x,y = cni.index2xy(i,j,sxllcorner,syllcorner,scellsize)
                value = selectedLayer.dataProvider().identify(QgsPoint(x,y),QgsRaster.IdentifyFormatValue).results()[1]
                try:
                    bi_row,bi_col,bi_lvl = cni.xyz2index(x,y,value,bminX,bminY,bminZ,bcolSize,browSize,blvlSize)
                    blockValue = blockModel[bi_row][bi_col][bi_lvl]
                    blankASC[i][j] = blockValue
                except:
                    blankASC[i][j] = sNODATA_value
        script = cni.generateASC(blankASC,sncols,snrows,sxllcorner,syllcorner,scellsize,sNODATA_value)
        f = open(outputFile, 'w')
        f.write(script)
        f.close()
        self.dlg.statusOutput.setValue(100)

    def run(self):
        self.dlg.comboBox.clear()
        self.dlg.lineEdit.clear()
        self.dlg.outputFilePath.clear()
        self.dlg.tableWidget.setRowCount(0)
        self.dlg.tableWidget.setColumnCount(0)
        self.dlg.eastingColumnSelection.clear()
        self.dlg.northingColumnSelection.clear()
        self.dlg.elevationColumnSelection.clear()
        self.dlg.valueColumnSelection.clear()
        self.dlg.statusBlockModelExtents.setValue(0)
        self.dlg.statusBlockModelLoading.setValue(0)
        self.dlg.statusIntersections.setValue(0)
        self.dlg.statusOutput.setValue(0)
        layers = self.iface.legendInterface().layers()
        layer_list = []
        if len(layers) == 0:
			return
        for layer in layers:
            if layer.type() == layer.RasterLayer:
                layer_list.append(layer.name())
        self.dlg.comboBox.addItems(layer_list)
        for layer in layers:
            if layer.name() == self.dlg.comboBox.currentText():
                selectedLayer = layer
        """Run method that performs all the real work"""
        # show the dialog
        self.dlg.show()