Mengimpor Pemrosesan QGIS dalam skrip python yang berdiri sendiri?

10

Saya ingin menulis beberapa skrip yang berdiri sendiri yang menggunakan kotak peralatan pemrosesan Qgis.

Saya telah membaca beberapa utas (di sini dan di sini , misalnya) tetapi belum dapat menemukan solusi yang berfungsi.

Menggunakan Qgis 2.16.1 di Ubuntu Xenial 16.04 LTS

Bagian impor skrip saya terlihat seperti ini:

# Python modules
import sys
import time
import os

# Qgis modules
from qgis.core import *
import qgis.utils
from PyQt4.QtCore import QFileInfo, QSettings

Adakah yang tahu apa yang hilang bagi saya untuk dapat mengimpor modul pemrosesan?

Dengan pemrosesan impor sederhana , saya mendapatkan ini:

Original exception was:
Traceback (most recent call last):
 File "/home/steph/Documents/Projets/20141227-CIM_Bishkek/Scripts/python/00-projets/20160811-AnalysesUAVs/20160811-UAVAnalyse.py", line 36, in <module>
import processing
 File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 607, in _import
mod = _builtin_import(name, globals, locals, fromlist, level)
ImportError: No module named processing

EDIT (setelah komentar Joseph)

Saya sudah mencoba seperti ini:

# Python modules
import sys
import time
import os

# Qgis modules
from qgis.core import *
import qgis.utils
from PyQt4.QtGui import *
app = QApplication([])
QgsApplication.setPrefixPath("/usr", True)
QgsApplication.initQgis()
from PyQt4.QtCore import QFileInfo, QSettings
#from PyQt4.QtGui import *

# Prepare processing framework 
sys.path.append('/home/steph/.qgis2/python/plugins')
from processing.core.Processing import Processing
Processing.initialize()
from processing.tools import *

tetapi perilakunya aneh: skrip saya berjalan sampai akhir tanpa kesalahan tetapi tampaknya "melompati" tugas yang seharusnya dilakukan :-) Dengan kata lain, ia berjalan sampai akhir tetapi tanpa melakukan apa-apa.

Saya akui bahwa penjelasan saya tidak begitu jelas ... Saya akan menyelidiki lebih lanjut tetapi jika ada yang punya solusi ajaib (bukan solusi) untuk mengimpor modul ini, silakan!

EDIT 2: menambahkan seluruh skrip saya. Maaf kalau agak panjang ....

# -*- coding: cp1252 -*-
########################################################
## Name: Performs various analyses on UAV imagery using Qgis
## Source Name: UAVanalyse.py
## Version: Python 2.7
## Author: Stephane Henriod
## Usage: Performs a set of analyses on UAV imagery
## Date 11.08.2016
## Modified: 
########################################################


# Import required modules

# Python modules
import sys
import time
import os

# Qgis modules
from qgis.core import *
import qgis.utils
from PyQt4.QtCore import QFileInfo, QSettings

# Custom modules
from config_work import *
import display_msg as disp
import clean_time as cl

def make_raster_layer(raster, log_file):
    """Creates a raster layer from the path to a raster, if the path exists and if the raster is valid

    Param_in:
        raster (string) -- The path to the raster to be transformed into a layer
        log_file (string) -- The path to the log file to write in

    Param_out:
        list: 
            [0] = full path to the raster 
            [1] = raster layer

    """

    if os.path.exists(raster):
        fileName = raster
        fileInfo = QFileInfo(fileName)
        baseName = fileInfo.baseName()
        rlayer = QgsRasterLayer(fileName, baseName)
        if rlayer.isValid():
            return [raster, rlayer]
    else:
        return False

def study_raster(rlayer, log_file):
    """Returns properties of a raster, if this one exists and is valid

    Param_in:
        rlayer (bin) -- A raster layer
        log_file (string) -- The path to the log file to write in

    """

    infos = {}

    if rlayer:
        infos['a - Width'] = rlayer.width()
        infos['b - Height'] = rlayer.height()
        infos['c - Extent'] = rlayer.extent().toString()
        infos['d - # bands'] = rlayer.bandCount()
        infos['e - X resolution'] = rlayer.rasterUnitsPerPixelX()
        infos['f - Y resolution'] = rlayer.rasterUnitsPerPixelY()
        return infos
    else:
        return False


def project_raster(raster, to_crs, log_file):
    """Projects a raster into another crs

    Param_in:
        raster (string) -- The path to the raster to be transformed into a layer
        to_crs (string) -- The coordinate reference system to which the layer must be projected
        log_file (string) -- The path to the log file to write in

    """

    img_out_name = os.path.splitext(os.path.basename(raster))[0] + '_proj' + os.path.splitext(os.path.basename(raster))[1]
    img_out = os.path.join(output_folder, img_out_name)
    #processing.runalg("gdalwarp -overwrite -s_srs EPSG:32642 -t_srs " + to_crs + " " + rlayer[0] + " " + img_out)

    msg = img_out    
    disp.display_msg(log_file, msg, 'a')

    return img_out_name

if __name__ == "__main__":
    t_start_script = time.localtime()
    t_start_script_clean = time.strftime("%Y%m%d-%H%M", t_start_script)

    # Checking the folders
    if not os.path.exists(input_folder_path):
        os.makedirs(input_folder_path)
    if not os.path.exists(temp_folder_path):
        os.makedirs(temp_folder_path)
    if not os.path.exists(output_folder_path):
        os.makedirs(output_folder_path)

    # Creating the output and temp folders
    output_folder = os.path.join(output_folder_path, t_start_script_clean + '-UAVanalyse')
    temp_folder = os.path.join(temp_folder_path, t_start_script_clean + '-UAVanalyse')

    if not os.path.exists(output_folder):
        os.makedirs(output_folder)
    if not os.path.exists(temp_folder):
        os.makedirs(temp_folder)

    # Creating the log file
    log_file_name = t_start_script_clean + '-UAVanalyse.log'
    log_file = os.path.join(output_folder, log_file_name)

    # Heading of the log file
    msg = "Performs a set of analyses on UAV imagery" + os.linesep
    msg += "Input folder: " + input_folder_path
    msg += "\n RGB image: " + img_rgb_name
    msg += "\n NIR image: " + img_nir_name
    msg += "\n RGBIR image: " + img_rgbir_name
    msg += "\n DSM file: " + img_dsm_name
    disp.display_msg(log_file, msg, 'w')

    #msg = "Script started on " + cl.clean_time(t_start_script)
    #disp.display_msg(log_file, msg, 'a')


    # Initialize Qgis (source: http://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/intro.html)
    msg = 'Initialize Qgis'
    disp.display_msg(log_file, msg, 'a')
    # supply path to qgis install location
    QgsApplication.setPrefixPath("/usr", True)

    # create a reference to the QgsApplication, setting the
    # second argument to False disables the GUI
    qgs = QgsApplication([], False)

    # load providers
    qgs.initQgis()


    # Write your code here to load some layers, use processing algorithms, etc.

    # Make raster layers
    rlayer_rgb = make_raster_layer(img_rgb, log_file)
    rlayer_nir = make_raster_layer(img_nir, log_file)
    rlayer_rgbir = make_raster_layer(img_rgbir, log_file)
    rlayer_dsm = make_raster_layer(img_dsm, log_file)

    all_valid_layers = []
    if rlayer_rgb: all_valid_layers.append(rlayer_rgb)
    if rlayer_nir: all_valid_layers.append(rlayer_nir)
    if rlayer_rgbir: all_valid_layers.append(rlayer_rgbir)
    if rlayer_dsm: all_valid_layers.append(rlayer_dsm)




    # (I) Infos about the layers
    msg = os.linesep + frm_separator + os.linesep + '(I) Infos about the layers' + os.linesep + frm_separator + os.linesep
    disp.display_msg(log_file, msg, 'a')

    i = 1
    for layer in all_valid_layers:
        infos = study_raster(layer[1], log_file)
        msg = '\n (' + str(i) + ') ' + layer[0] + os.linesep
        for item in sorted(infos):
            msg += '\n ' + str(item) + ': ' + str(infos[item]) + os.linesep

        i+=1
        disp.display_msg(log_file, msg, 'a')

    msg = 'List of valid layers:' + os.linesep
    for layer in all_valid_layers:
        msg += layer[0]+ os.linesep
    disp.display_msg(log_file, msg, 'a')


    # (II) Projects the layers into the national coordinate system or any desired system
    msg = os.linesep + frm_separator + os.linesep + '(II) Projecting of the layers' + os.linesep + frm_separator + os.linesep
    disp.display_msg(log_file, msg, 'a')

    i = 1
    for layer in all_valid_layers:
        project_raster(layer[0], to_crs, log_file)




    # When script is complete, call exitQgis() to remove the provider and
    # layer registries from memory
    qgs.exitQgis()
    msg = 'Qgis has been closed'
    disp.display_msg(log_file, msg, 'a')

    #raw_input("Press Enter to continue...")
Stéphane Henriod
sumber
Di tautan kedua Anda, apakah Anda menggunakan kode yang disediakan oleh @ GermánCarrillo? Kode nya adalah apa yang saya juga gunakan untuk menjalankan skrip mandiri (saya mengedit jalur saya sedikit ketika saya menggunakan Windows).
Joseph
Apa tugas yang Anda coba lakukan? :) Bisakah Anda memasukkan kode dalam pertanyaan Anda? Ini dapat membantu orang lain mengidentifikasi apa yang salah.
Joseph
Saya mencoba membuat skrip beberapa fungsi untuk memproses gambar UAV. Hal pertama yang perlu saya lakukan adalah memproyeksikan ulang, itulah sebabnya saya membutuhkan modul pemrosesan. Izinkan saya menambahkan skrip lengkap saya, meskipun saya sedikit malu: Saya bukan pengembang "nyata" dan saya yakin kode saya bukan yang terbersih atau paling lurus yang pernah Anda lihat :-)
Stéphane Henriod
Jangan malu! Tidak ada yang saya poskan adalah yang paling bersih atau lurus ke depan =)
Joseph
Jadi saya menambahkan seluruh skrip saya. Silakan melihat-lihat :-)
Stéphane Henriod

Jawaban:

6

Linux QGIS 2.18.1

Dengan kode ini, jalankan pemrosesan dari skrip mandiri:

#!/usr/bin/env python
import qgis
from qgis.core import *
import sys

app = QgsApplication([],True, None)
app.setPrefixPath("/usr", True)
app.initQgis()
sys.path.append('/usr/share/qgis/python/plugins')
from processing.core.Processing import Processing
Processing.initialize()
from processing.tools import *

print Processing.getAlgorithm("qgis:creategrid")
Juanma Font
sumber
Ini adalah satu-satunya kombinasi yang berfungsi untuk lingkungan tempat saya bekerja (Pycharm C.Ed pada mesin Ubuntu 14.04 dengan Python 2.7). Sebelumnya, saya mencoba kombinasi gis.stackexchange.com/questions/129513/... , dan gis.stackexchange.com/questions/176821/… dan sayangnya, tidak ada satu pun dari mereka yang mengimpor "processing.core.Processing". Saya tidak tahu mengapa mengimpor modul membingungkan. Sebagai catatan: ada paket Python murni juga disebut "pemrosesan" yang membayangi yang QGIS.
Irene
3

Jadi saya berhasil membuatnya bekerja, terima kasih @ Joseph untuk petunjuk Anda:

# Import required modules

# Python modules
import sys
import time
import datetime
import os
from getpass import getuser

# Qgis modules and environment
from qgis.core import *
import qgis.utils
from PyQt4.QtCore import QFileInfo, QSettings

from PyQt4.QtGui import QApplication
app = QApplication([])

QgsApplication.setPrefixPath("/usr", True)
QgsApplication.initQgis()

# Prepare processing framework
sys.path.append('/home/' + getuser() + '/.qgis2/python/plugins')
from processing.core.Processing import Processing

Processing.initialize()

Dan saya bisa mengujinya dengan

print Processing.getAlgorithm("qgis:creategrid")

Masalah saya, saya kira, berasal dari fakta bahwa saya mengimpor modul di awal skrip saya dan kemudian mencoba membuat objek Qgi dari dalam suatu fungsi. Saya kira itu mungkin untuk dilakukan juga, tetapi mungkin ada kelemahan dalam keterampilan Python saya.

Sekarang saya akan mencoba menggunakan modul pemrosesan :-)

Stéphane Henriod
sumber
Sesuatu telah berubah ? Saya mencoba kode ini (dari Joseph) untuk menggunakan pemrosesan menjadi skrip mandiri dan saya mendapatkan: ImportError: Tidak ada modul bernama processing.core.Processing Menggunakan linux QGIS 2.18.1
Juanma Font