QGIS » QGIS Application

from qgis.PyQt.QtCore import QCoreApplication

from qgis.core import (QgsProcessing,QgsFeatureSink,QgsProcessingException,QgsProcessingAlgorithm,#QgsRasterLayer,QgsMapLayerRegistry,QgsProcessingAlgorithm,

QgsProcessingParameterFeatureSource,QgsProcessingParameterFeatureSink,QgsFeatureRequest,QgsFeature,QgsProcessingParameterField,QgsProcessingParameterString,QgsRasterLayer,QgsProject,QgsProcessingParameterRasterDestination)

from qgis.utils import iface,QgsMessageBar

import processing

class ExampleProcessingAlgorithm(QgsProcessingAlgorithm):

        #Declaration de nos shapes comme des variables text

        Vector_layer_of_communes = 'Vector_layer_of_communes'

        OUTPUT_raster = 'OUTPUT_raster'

        Commune_name='Commune_name'

        INSEE_code='INSEE_code'

        EPSG_code='EPSG_code'

        def tr(self, string):

                """

                Returns a translatable string with the self.tr() function.

                """

                return QCoreApplication.translate('Processing', string)

        def createInstance(self):

                return ExampleProcessingAlgorithm()

        def name(self):

                return 'Cadastre FR - WMS - Add a cadastral map'

        def displayName(self):

                return self.tr('Cadastre FR - WMS - Add a cadastral map')

        #Creation de group pour organiser laffichage des scripts

        def group(self):

                return self.tr(self.groupId())

        def groupId(self):

                return 'Circet'

        #Creation de description du script

        def shortHelpString(self):

                #Partie pour afficher la description du script

                return self.tr("Cadastre FR - WMS - Add a cadastral mapt")

        #Partie diniatilisation du script

        def initAlgorithm(self, config=None):

                #Initialisation de nos variables shapes dentree comme des shapes

                self.addParameter(

                        QgsProcessingParameterFeatureSource(

                                self.Vector_layer_of_communes,

                                self.tr('Input Vector_layer_of_communes'),

                                [QgsProcessing.TypeVectorAnyGeometry]

                        )

                )

        #Initialisation de nos variables shapes de sorti comme des shapes

                self.addParameter(QgsProcessingParameterRasterDestination(

                        self.OUTPUT_raster,

                        self.tr("OUTPUT_raster"),

                        None, False))

                #Initialisation de nos variables dattribut des shapes dentree

                self.addParameter(

        QgsProcessingParameterField(

            self.Commune_name,#"SpeciesAttribute",

            description=self.Commune_name,#"Species attribute",

            parentLayerParameterName=self.Vector_layer_of_communes,#"SpeciesRangePolygons",INSEE_code

            type=QgsProcessingParameterField.Any#String

        )

    )

                self.addParameter(

                        QgsProcessingParameterField(

                                self.INSEE_code,#"SpeciesAttribute",

                description=self.INSEE_code,#"Species attribute",

                parentLayerParameterName=self.Vector_layer_of_communes,#"SpeciesRangePolygons",INSEE_code

                                type=QgsProcessingParameterField.Any

                        )

                )

                self.addParameter(

                        QgsProcessingParameterString(

                                self.EPSG_code,#name="SpeciesAttribute",

                description='Code du projection, exemple: 2154'#self.EPSG_code

                                #self.EPSG_code,#description="Species attribute",

                                #self.Vector_layer_of_communes,#parentLayerParameterName="SpeciesRangePolygons",

                                #type=QgsProcessingParameterField.Any

                        )

                )

        #Partie dimplentation du code python qui fait le traitement

        def processAlgorithm(self, parameters, context, feedback):

                #Parametrage des shapes dentree

                shape_commune = self.parameterAsSource(parameters, self.Vector_layer_of_communes, context)

                #Parametrage du resultat en fichier de sorti

                output_path_raster = self.parameterAsOutputLayer(parameters, self.OUTPUT_raster, context)

                attribut_INSEE_code = self.parameterAsString(parameters, self.INSEE_code, context)

                attribut_Commune_name = self.parameterAsString(parameters, self.Commune_name, context)

                attribut_EPSG_code = self.parameterAsString(parameters, self.EPSG_code, context)

                def recuperation_cadastre_commune(layer,INSEE_code,Commune_name,EPSG_code):

                        if EPSG_code == '2154' or EPSG_code == '3942' or EPSG_code == '3943' or EPSG_code == '3944' or EPSG_code == '3945' or EPSG_code == '3946' or EPSG_code == '3947' or EPSG_code == '3948' or EPSG_code == '3949' or EPSG_code == '3950' or EPSG_code == '32630' or EPSG_code == ' 32631' or EPSG_code == '32632' or EPSG_code == '3857' or EPSG_code == '4326' or EPSG_code == '4258' or EPSG_code == '32620' or EPSG_code == '2970' or EPSG_code == '2972' or EPSG_code == '2973' or EPSG_code == '2975' or EPSG_code == '32622' or EPSG_code == '32740' or EPSG_code == '32738' or EPSG_code == '4471' or EPSG_code == '32621' :

                                #progress.setText(u'EPSG code : ' + EPSG_code)

                                tab = []

                                for f in layer.getFeatures():

                                        col_select =str(f[INSEE_code]),f[Commune_name]

                                        tab.append(col_select)

                                        #Permet la suppression des doublons

                                        Lt= list(set(tab))

                                        Lt.sort()

                                for c_insee, n_couche in Lt  :

                    #AMORCES_CAD,LIEUDIT,CP.CadastralParcel,SUBFISCAL,CLOTURE,DETAIL_TOPO,HYDRO,VOIE_COMMUNICATION,BU.Building,BORNE_REPERE

                                        urlWithParams = "url=http://inspire.cadastre.gouv.fr/scpc/"+c_insee+".wms?contextualWMSLegend=0&crs=EPSG:"+EPSG_code+"&dpiMode=7&featureCount=10&format=image/png&layers=AMORCES_CAD&layers=LIEUDIT&layers=CP.CadastralParcel&layers=SUBFISCAL&layers=CLOTURE&layers=DETAIL_TOPO&layers=HYDRO&layers=VOIE_COMMUNICATION&layers=BU.Building&layers=BORNE_REPERE&styles=&styles=&styles=&styles=&styles=&styles=&styles=&styles=&styles=&styles=&maxHeight=1024&maxWidth=1280"

                                        rlayer = QgsRasterLayer(urlWithParams, 'Cadastre_'+n_couche+'_'+c_insee, 'wms')

                                        #My problems is to insert this result to the raster created "output_path_raster"

                                        output_path_raster=rlayer

                                        #progress.setText(u'Commune name : ' + n_couche+' - '+c_insee)

                                        #progress.setText(u'Validity of WMS : %s' % rlayer.isValid())

                                        #QgsProject.instance().addMapLayer(rlayer)

                                        '''if  rlayer.isValid() == True :

                                                iface.messageBar().pushMessage("Information :", "Ajout du cadastre dans Qgis_map : "+n_couche, Qgis.Info, duration=5)

                                                iface.mapCanvas().refresh()

                                        else :

                                                iface.messageBar().pushMessage("Warning :", "WMS du Cadastre est invalide : "+n_couche, Qgis.Warning, duration=10)

                        else :

                                iface.messageBar().pushMessage("Warning :", "Projection Inconnue ", Qgis.Warning, duration=10)

                                progress.setText(u'EPSG is unknown ')'''

                #Fonction de calcul du besoin

                recuperation_cadastre_commune(shape_commune,attribut_INSEE_code,attribut_Commune_name,attribut_EPSG_code)

                results = {}

                results[self.OUTPUT_raster] = output_path_raster

                #return {self.OUTPUT_PT: dest_id_boite_sans_PT, self.OUTPUT_CABLE: dest_id_boite_sans_cables}

                return results#{}