QGIS DENGAN PYTHON

Pengertian

QGIS memiliki sebuah antarmuka pemrograman yang ampuh sehingga kita bisa memperlebar fungsi inti dari software dan juga menulis skrip untuk otomasi pelaksanaan tugas. QGIS mendukung Bahasa skript Phyton yang terkenal. Walaupun anda seorang pemula, mempelajari sedikit phyton dan antarmuka Pemrogramman QGIS akan membuat anda jauh lebih produktif dalam pekerjaan anda. Tutorial ini mengasumsikan anda tanpa pengalaman pemrograman dan ditujukan untuk memberikan perkenalan pada scripting Python di QGIS (PyQGIS).

Tinjauan Tugas 

Kita akan membjuka sebuah layer poin yang merepresentasikan semua airport besan dan menggunakan skript python untuk menghasilkan sebuah file teks dengan nama airport, kode airport, garis lintang dan garis bujur untuk setiap airport pada layer.

Mendapatkan data

kita akan menggunakan dataset dari Natural Earth Airports

Unduh Airports shapefile.

Sumber data [NATURALEARTH]

Prosedur

1. Dalam QGIS, akses Layers ‣ Add Vector Layer . Jelajah ke file ne_10m_airports.zip yang sudah terunduh dan klik ne_10m_airports.zip . Pilih layer ne_10m_airports.shp dan klik OK.

2. Anda akan melihat layer ne_10m_airports terbuka di QGIS

3. Pilih tool Identify dan klik pada poin mana saja untuk memeriksa attribut yang tersedia. Anda akan melihat bahwa nama airport dan 3 digit kode didalamnya adalah attribut name dan iata_code secara berurutan

4. QGIS menyediakan konsol built-in dimana anda bisa mengetik command phyton dan mendapatkan hasilnya. Konsol ini merupakan cara yang sangat baik untuk belajar scripting dan juga melakukan proses data cepat. Buka Python Console dengan mengakses Plugins ‣ Python Console.

5. Anda akan melihat sebuah panel baru di bawah kanvas QGIS. Anda akan melihat sebuah tanda seperti >>> di bawah dimana anda dapat mengetik command. Untuk interaksi dengan environment QGIS, anda harus menggunakan variabel iface . Untuk mengakses layer yang sedang aktif di QGIS, anda dapat mengetik sebagai berikut dan tekan iface . Command ini menarik referensi ke layer yang sedang dibuka dan menyimpannya sebagai variabel layer.

layer = iface.activeLayer()

6. Ada function berguna yang disebut dir() di phyton yang menunjukkan anda semua metode yang tersedia untuk objek mana saja. Ini berguna ketika anda tidak yakin function mana yang tersedia untuk obyek tertentu. Jalankan command berikut dan lihat operasi apa yang bisa kita lakukan pada variabel layer .

dir(layer)

7. Anda akan melihat sebuah daftar panjang dari function yang tersedia. Untuk sekarang, kita akan menggunakan sebuah fungsi bernama getFeatures() dimana akan memberikan anda referensi untuk semua fitur pada sebuah layer. Pada kasus kita, setiap fitur akan menjadi sebuah poin yang merepresentasikan airport. Anda dapat mengetik command untuk memgiterasi melalui setiap fitur pada layer yang sedang aktif. Pastikan untuk menambah 2 spasi sebelum mengetik garis kedua.

for f in layer.getFeatures():

  print f

8. Seperti yang anda lihat di output, setiap garis mengandung sebuah referensi ke sebuah fitur di layer. Referensi ke fitur disimpan dalam variabel f . Kita dapat menggunakan variabel f``untuk mengakses attribut pada setiap fitur. Ketik sebagai berikut untuk mencetak ``name dan iata_code untuk setiap fitur airport.

for f in layer.getFeatures():

  print f['name'], f['iata_code']

9. Jadi sekarang anda tahu bagaimana secara programatik mengakses attribut dari seitap fitur di layer. Sekarang, mari kita lihat bagaimana mengakses koordinat fitur. Koordinat dari sebuah fitur vektor dapat diakses dengan memanggil function geometry() . Function ini memberikan sebuah obyek geometri yang dapat kita simpan pada variable geom . Anda dapat menjalankan function asPoint() pada obyek geometri untuk mendapatkan koordinat x dan y poin tersebut. Jika fitur anda adalah sebuah garis atau sebuah poligon, anda dapat menggunakan fungsi asPolyline() atau asPolygon(). Ketik Kode berikut pada prompt dan tekan Enter untuk melihat koordinat x dan y dari setiap fitur.

for f in layer.getFeatures():

  geom = f.geometry()

  print geom.asPoint()

10. Bagaimana jika kita hanya ingin mendapatkan koordinat x dari fitur ? Anda dapat memanggil fungsi x() pada obyek poin dan mendapatkan koordinat x nya.

for f in layer.getFeatures():

  geom = f.geometry()

  print geom.asPoint().x()

11. Sekarang kita mempunyai semua kepingan yang akan kita satukan untuk menghasilkan output yang kita inginkan. Ketik kode berikut untuk mencetak name, iata_code, latitude dan longitude dari setiap fitur airport. Notasi %s dan %f adalah cara untuk memformat sebuah variabel string dan numerik.

for f in layer.getFeatures():

  geom = f.geometry()

  print '%s, %s, %f, %f' % (f['name'], f['iata_code'],

         geom.asPoint().y(), geom.asPoint().x())

12. Anda dapat melihat outpu tercetak pada konsol. Sebuah cara yang lebih berguna untuk menyimpan output adalah dengan sebuah file. Anda dapat mengetik kode berikut untuk membuat sebuah file dan menulis outputnya di sana, Ganti path atau alamat file dengan sebuah path dalam sistem anda sendiri. Catat bahwa kita menambahkan \n di akhir baris. Ini untuk menambah sebuah garis baru setelah kita menambahkan data untuk setiap fitur. Anda seharusnya juga memperhatikan garis unicode_line = line.encode('utf-8') . Karena layer kita mengandung sejumlah fitur dengan karakter unicode, kita tidak bisa begitu saja menulisnya ke dalam file teks. Kita encode teks menggunakan UTF-8 encoding dan kemudian menulis file teks tersebut.

output_file = open('c:/Users/Ujaval/Desktop/airports.txt', 'w')

for f in layer.getFeatures():

  geom = f.geometry()

  line = '%s, %s, %f, %f\n' % (f['name'], f['iata_code'],

          geom.asPoint().y(), geom.asPoint().x())

  unicode_line = line.encode('utf-8')

  output_file.write(unicode_line)

output_file.close()

13. Anda dapat mengakses lokasi file output yang sudah ditentukan dan membuka file teks. Anda akan melihat data dari shapefile airport yang kita ekstraksi menggunakan skripting phyton.