Sampling Data Raster menggunakan Poin atau Poligon

Banyak dataset sains dan lingkungan berupa raster bergrid. Data elevasi (DEM) juga terdistribusi berupa file raster. Pada raster file ini, parameter yang tersaji terencoded sebagai nilai pixel raster. Seringkalim diperlukan untuk mengekstrak nilai pixel pada lokasi tertentu atau mengagregasi mereka pada sejumlah area. Fungsi ini tersedia di QGIS dalam 2 plugin - Point Sampling Tool dan Zonal Statistics plugin.

Tinjauan Tugas

Diberikan sebuah grid raster untuk temperatur maksimum di AS, kita perlu untuk mengekstrak temperatur pada semua area urban dan juga menghitung temperatur rata-rata untuk setiap county di AS.

Skill lain yang akan anda pelajari

  • mere-proyeksi sebuah layer vektor.

  • Pilih dan hapus beberapa layer dari Table of Content QGIS.

Mendapatkan data

NOAA’s Climate Prediction Center provides GIS data related to temperature and precipitation in the US. Download the latest grid filei for maximum temperatures. The file will be named us.tmax_nohads_ll_{YYYYMMDD}_float.tif

Kita akan menggunakan sebuah file CSV dari 2013 US Gazetteer merepresentasikan area urban di AS. Unduh Urban Areas Gazetteer File

Karena kita ingin mengagregasi temperatur pada tiap county, kita akan menggunakan 2013 TIGER/Line Shapefiles . Unduh Counties (and equivalents) shapefile.

For convenience, you may directly download a copy of the datasets from the links below:

us.tmax_nohads_ll_20140525_float.tif

2013_Gaz_ua_national.zip

tl_2013_us_county.zip

Sumber Data: [NOAACPC], [USGAZETTEER] [TIGER]

Prosedur

  1. Akses Layer ‣ Add Raster Layer dan jelajahi ke file unduhan us.tmax_nohads_ll_{YYYYMMDD}_float.tif dan klik Open.

../_images/1131.png
  1. Ketika layer terbuka, pilih took Identify dan klik di amna saja pada layer. Anda akan melihat nilai temperatur dalam celseius atau Band 1 di lokasi.

../_images/279.png
  1. Sekarang unzip file unduhan 2013_Gaz_ua_national.zip dan ekstrak file 2013_Gaz_ua_national.txt pada disk anda. Akses Layer ‣ Add Delimited Text Layer.

../_images/340.png
  1. Pada dialog Create a Layer from Delimited Text File , klik Browse dan 2013_Gaz_ua_national.txt . Pilih Tab pada Custom delimiters . Koordinat poin dalam lintang atau latitude dan bujur atau longitude, jadi pilih INTPTLONG sebagai X field dan INTPTLAT sebagai Y field . Beri tanda cek pada box Use spatial index dan klik OK.

../_images/430.png
  1. Sekarang kita siap untuk mengekstrak nilai temperatur dari layer raster,. Install plugin Point Sampling Tool . Lihat docs/using_plugins untuk deteil bagaimana menginstall plugins.

../_images/531.png
  1. Buka dialog plugin dari Plugins ‣ Analyses ‣ Point sampling tool.

../_images/629.png
  1. Pada dialog Point Sampling Tool , pilih 2013_Gaz_ua_national sebagai Layer containing sampling points . Kita harus secara eksplisit mengambil field dari layer input yang kita inginkan pada layer output. Pilih field GEOID dan NAME dari layer 2013_Gaz_ua_national . Kita dapat mensampel nilai dari beberapa band raster sekali bersamaan, tapi karena raster kita hanya memiliki 1 band, pilih us.tmax_nohads_ll_{YYYYMMDD}_float: Band 1 . Beri nama layer vektor output sebagai max_temparature_at_urban_locations.shp . Klik OK untuk memulai proses sampling. Klik Close ketika proses selesai.

../_images/728.png
  1. Anda akan melihat sebuah layer baru max_temparature_at_urban_locations terbuka di QGIS. Gunakan tool Identify untuk mengklik poin mana saja untuk melihat attribut. Anda akan melihat field us.tmax_no - yang berisi nilai pixel raster pada lokasi poin.

../_images/827.png
  1. Bagian pertama dari anlisis kita sudah selesai. Mari hapus layer yang tidak diperlukan. Tahan Shift dan pilih layer max_temparature_at_urban_locations dan 2013_Gaz_ua_national . Klik-kanan dan pilih Remove untuk menghapus mereka dari QGIS TOC atau Table of Content.

../_images/925.png
  1. Akses Layer ‣ Add Vector Layer . Jelajahi file unduhan tl_2013_us_county.zip dan klik Open . Pilih tl_2013_us_county.shp sebagai layer dan klik:guilabel:OK.

../_images/1026.png
  1. tl_2013_us_county akan ditambahkan ke QGIS. Layer ini dalam proyeksi EPSG:4269 NAD83 . Ini tidak sesuai dengan proyeksi dari layer raster. Kita akan mere-proyeksi layer ini ke proyeksi EPSG:4326 WGS84

../_images/1132.png
  1. Klik kanan layer tl_2013_us_county dan pilih Save As...

../_images/1227.png
  1. Pada dialog Save Vector layer as.., klik Browse dan beri nama output file sebagai counties.shp . Pilih Selected CRS dari menu dropdown CRS . Klik Browse dan pilih WGS 84 sebgai CRS. Beri tanda cek pada Add saved file to map dan klik OK.

../_images/1324.png
  1. Sebuah layer baru bernama counties akan ditambahkan ke QGIS.

../_images/1423.png
  1. Aktifkan Zonal Statistics Plugins . ini adalah plugin inti, jadi sudah terinstall. Lihat docs/using_plugins untuk mengetahui bagaimana mengaktifkan plugin inti.

../_images/1522.png
  1. Akses Raster ‣ Zonal statistics ‣ Zonal statistics.

../_images/1621.png
  1. Pilih us.tmax_nohads_ll_{YYYYMMDD}_float sebagai Raster layer dan counties sebagai Polygon layer containing the zones. Masukkan ZS_ sebagai Output column prefix . Klik OK.

../_images/1720.png
  1. Analisis mungkin membutuhkan beberapa waktu bergantung pada ukuran dataset.

../_images/1820.png
  1. Ketika proses selesai, pilih layer counties . Gunakan tool Identify dan klik pada poligon county mana saja. Anda akan melihat tiga attribut baru ditambahkan ke layer : ZS_count` , ZS_mean dan ZS_sum . Attribut ini berisi jumlah pixel raster, rata-rata nilai pixel raster dan jumlah nilai pixek raster secara berurutan. Karena kita tertarik pada tempertaur rata-rata, field ZS_mean akan menjadi yang kita cari.

../_images/1918.png
  1. Mari menstyle layer ini untuk membuat sebuah peta temperatur. Klik-kanan layer counties dan pilih Properties.

../_images/2015.png
  1. Pindah ke tab Style . Pilih style Graduated dan pilih ZS_mean sebagai Column . Pilih Color Ramp dan Mode sesuai pilihan anda. Klik Classify untuk membuat kelas-kelas. Klik OK. (lihat docs/pembuatan_vector_sederhana untuk detail lebih dalam styling )

../_images/2119.png
  1. Anda akan melihat poligon county yang sudah distyle menggunakan temperatur maksimum rata-rata dari grid raster.

../_images/2218.png
comments powered by Disqus

This work is licensed under a Creative Commons Attribution 4.0 International License