Membangun peta panas yang efektif?

65

Menggunakan ArcGIS, QGIS, Rumput, dan / atau GVSIG:

  • Apa saja alat dan proses yang terlibat dalam membangun peta panas yang efektif?
  • Apa saja plugin yang terlibat?
  • Apa persyaratan data utama?
  • Apa saja kelemahan dengan peta panas yang ada?
  • Apa saja masalah yang tidak bisa dibahas oleh peta panas secara efektif?
  • Bagaimana tidak melakukan peta panas?
  • Apakah ada alternatif yang lebih baik (dalam konteks yang sama) daripada peta panas untuk representasi data?
dassouki
sumber
Meskipun tidak menggunakan salah satu alat yang Anda tentukan, Anda mungkin ingin melihat skrip Python ini juga jjguy.com/heatmap
radek
1
Dassouki, bisakah Anda menjelaskan apa yang Anda maksud dengan "peta panas"? Wikipedia tampaknya berpikir itu adalah terjemahan choropleth dari sebuah array nilai. Ini secara halus, tetapi penting, berbeda dari semua balasan di utas ini, yang menganggapnya hanya berarti peta dari setiap atribut berbasis grid (atau gambar-), terutama yang telah diinterpolasi ke grid dari data titik. Jawaban untuk setiap pertanyaan berpoin Anda akan berbeda untuk peta panas yang sebenarnya.
whuber

Jawaban:

74

Pertanyaan ini telah dikonversi ke Wiki Komunitas dan wiki dikunci karena ini adalah contoh pertanyaan yang mencari daftar jawaban dan tampaknya cukup populer untuk melindunginya dari penutupan. Itu harus diperlakukan sebagai kasus khusus dan tidak boleh dipandang sebagai jenis pertanyaan yang didorong pada ini, atau situs Stack Exchange, tetapi jika Anda ingin berkontribusi lebih banyak konten untuk itu maka jangan ragu untuk melakukannya dengan mengedit jawaban ini .


Setidaknya ada dua jenis peta panas:

  1. Heatmaps mewakili konsentrasi titik, dan
  2. Heatmaps mewakili distribusi nilai atribut

Setiap metode memiliki kelebihan dan masalah, saya khawatir masuk ke detail jauh melampaui T&J ini.

Saya akan mencoba membuat daftar beberapa metode dan fungsi untuk QGIS dan GRASS.

Konsentrasi poin

Jika Anda melacak pergerakan satwa liar, kendaraan, dll. Akan berguna untuk menilai daerah dengan konsentrasi pesan lokasi yang tinggi.

Alat: mis. Plugin QGIS Heatmap (tersedia dalam versi> 1.7.x) atau GRASS v.neighbors atau v.kernel

Distribusi nilai atribut

Di sini, kita pada dasarnya berbicara lebih atau kurang tentang metode interpolasi. Metode meliputi:

  1. IDW

    Bergantung pada implementasinya, ini bisa bersifat global (menggunakan semua titik yang tersedia di set) atau lokal (dibatasi oleh jumlah titik atau jarak maksimum antara titik dan posisi yang diinterpolasi).

    Alat: Plugin interpolasi QGIS (global), GRASS v.surf.idw atau r.surf.idw (lokal)

  2. Splines

    Sekali lagi, sejumlah besar kemungkinan implementasi. B-Splines populer.

    Alat: GRASS v.surf.bspline

  3. Kriging

    Metode statistik dengan berbagai sub-tipe.

    Alat: GRASS v.krige (terima kasih kepada om_henners untuk tipnya) atau menggunakan R.

underdark
sumber
1
Ada antarmuka melalui GRASS untuk kriging, v.krige ( grass.osgeo.org/wiki/V.krige_GSoC_2009 ), tetapi masih membutuhkan R, dan berbagai paket R dan binding yang disebutkan pada halaman Wiki GRASS.
om_henners
Sejak QGIS 2.8 ada renderer Heatmap untuk layer titik. Tidak perlu membuat data baru.
Alexandre Neto
34

Secara statistik, inilah cara Anda melakukan peta panas:

1) Mengintegrasikan fitur titik. Ide integrasi adalah untuk mengambil poin yang harus dianggap kebetulan dan menggabungkannya sebagai satu lokasi. Saya suka menggunakan analisis tetangga terdekat dan menggunakan nilai yang sesuai dari sana. (Misalnya, ketika melakukan peta panas kejahatan, saya menggunakan tetangga terdekat rata - rata 1 untuk dataset paket yang mendasari terhadap mana kejahatan tersebut di-geocode).

2) Kumpulkan acara . Ini menciptakan bobot spasial untuk semua titik terintegrasi Anda. Misalnya, jika Anda memiliki 5 acara di satu lokasi, itu akan menjadi satu poin dengan bobot 5. Ini penting untuk dua langkah berikutnya. Jika Anda perlu menjumlahkan atribut dalam acara yang dikumpulkan, yaitu acara yang berbeda memiliki bobot lebih tinggi, maka Anda dapat menggunakan gabungan spasial satu-ke-satu . Gunakan output 'kumpulkan acara' sebagai target dan acara terintegrasi asli Anda sebagai fitur bergabung. Tetapkan aturan gabungan peta bidang yang secara statistik menggabungkan atribut pada acara terintegrasi (biasanya dengan SUM, meskipun Anda dapat menggunakan statistik lainnya).

3) Tentukan autokorelasi spasial puncak menggunakan global Moran saya . Seperti yang dikatakannya, jalankan global moran's I pada interval yang berbeda untuk menentukan pita puncak autokorelasi spasial dalam skala yang sesuai dengan analisis yang Anda lakukan. Anda mungkin ingin menjalankan tetangga terdekat lagi pada acara yang Anda kumpulkan untuk menentukan rentang awal untuk pengujian I moran Anda. (mis. gunakan nilai maks untuk tetangga terdekat pertama)

4) Jalankan Getis-Ord Gi * . Gunakan pita jarak tetap berdasarkan analisis I moran Anda, atau gunakan pita jarak tetap sebagai zona ketidakpedulian. Bobot spasial Anda dari acara kumpulkan adalah bidang jumlah numerik Anda. Ini akan memberi Anda skor-z untuk setiap titik acara di set Anda.

5) Jalankan IDW terhadap hasil Anda dari Getis-Ord Gi *.

Hasil ini sangat berbeda dari yang Anda dapatkan dengan kepadatan kernel. Ini akan menunjukkan kepada Anda di mana nilai tinggi dan nilai rendah dikelompokkan bersama daripada hanya di mana nilai tinggi, tanpa memperhatikan pengelompokan, seperti dalam kepadatan kernel.

tuan-castillo
sumber
20

Sementara saya suka peta panas, saya menyadari mereka sering salah digunakan.

Biasanya yang saya lihat adalah proses di mana warna setiap piksel didasarkan pada hasil dari fungsi berbobot jarak terbalik yang diterapkan pada kumpulan titik. Setiap kali peta memiliki banyak penanda titik yang tumpang tindih, saya pikir perlu mempertimbangkan peta panas.

Ini api berbasis web .

GeoChalkboard memiliki tutorial yang bagus untuk itu .

Anda dapat menggunakan IDW di ArcGIS.

Kirk Kuykendall
sumber
7
Perlu diketahui bahwa IDW sangat sensitif terhadap lokasi pengumpulan data. Jika data berkerumun, misalnya, Anda berpotensi mendapatkan anomali matematika yang buruk.
Reed Copsey
@Reed Copsey Alternatif apa yang akan Anda sarankan?
fmark
2
@ tanda: Ada banyak rutin interpolasi yang dapat Anda gunakan sebagai ganti IDW, termasuk pendekatan bertetangga / triangulasi alami, Kriging, splining / tegangan minimum, dll.
Reed Copsey
@Reed Aku tidak pernah benar-benar khawatir tentang kebenaran matematis dari heatmaps (mungkin aku harus). Tapi saya pikir mereka berkomunikasi dengan baik dalam banyak situasi. Berikut adalah contoh dari peta yang saya pikir dapat bermanfaat sebagai peta panas
Kirk Kuykendall
2
Saya pikir mereka alat yang hebat. Ketepatan matematis / statistik mungkin hanya penting jika Anda menggunakan hasil untuk pengambilan keputusan, tetapi jika itu untuk menyampaikan pengertian umum tentang distribusi, IDW mungkin baik-baik saja. (Ini lebih merupakan masalah cluster yang menyebabkan "kemiringan" besar dalam hasil peta panas, terutama antara cluster, karena anomali matematika.)
Reed Copsey
12

Untuk peta panas sederhana dan menghasilkan garis countour, saya telah menggunakan QGis dengan integrasi Rumput:

  1. Memuat titik data
  2. Memuat bentuk pembatas - misalnya batas wilayah
  3. Buat mapset Grass
  4. Buka kotak alat Rumput dan klik pada daftar modul untuk mencari setiap alat
  5. Muat modul v.in.ogr.qgis dan muat data titik dan bentuk batas, setiap kali ingat untuk mengeklik tampilan keluaran untuk masing-masing - beri masing-masing nama yang berguna seperti pointdata dan maskshape
  6. Ubah maskshape menjadi raster untuk menggunakannya sebagai mask dengan v.to.rast dan tambahkan ke mapset - sebut saja sesuatu seperti maskraster - ini bisa memakan waktu untuk poligon yang kompleks.
  7. Muat modul r.mask untuk memaksa tindakan selanjutnya dibatasi pada wilayah buffer.
  8. Jalankan v.surf.rst untuk menghasilkan kisi interpolasi dari pointdata - pilih kolom yang sesuai sebagai bidang atribut untuk melakukan interpolasi, dan beri nama seperti rasterurface . Ini adalah bit yang membutuhkan waktu dan menghasilkan raster yang dapat digunakan sebagai peta panas atau dapat diarsir 3D.
  9. Tutup kotak alat Rumput
  10. Gunakan plug-in Kontur GDAL Raster untuk memilih raster GRASS sebagai input; biarkan nilai level default pada 10, dan pilih direktori keluaran tempat contour shapefile akan disimpan. Periksa "Nama Atribut" dan ketikkan nama.

NB: Agar ini berfungsi, dataset harus dalam proyeksi yang sama!

Adrian
sumber
4

Saya pikir pertanyaan ini sebagian besar telah dijawab kecuali beberapa poin tentang masalah tersebut.

Peta panas bisa bagus, tetapi kelemahan klasik dan masalah terletak pada interpretasi. Ambil perbedaan antara peta panas peristiwa kejahatan dibandingkan dengan peta (panas atau lainnya) tingkat kejahatan / proporsi. Sementara peta panas peristiwa mungkin berguna dalam hal mengidentifikasi kepadatan acara keseluruhan, itu buta sebagai perkiraan risiko, tetapi sering ditafsirkan atau disalahgunakan dengan cara ini. Pertimbangkan jumlah kejadian yang sama di suatu wilayah dengan ukuran dan bentuk yang sama, tetapi dengan populasi yang berbeda, sementara kejahatan mungkin terkonsentrasi di suatu daerah, itu bisa saja karena ada lebih banyak orang di ruang itu. Selain itu, angka untuk data kejadian, seperti untuk kejahatan, bisa sulit untuk dimodelkan, karena untuk menghasilkan peta panas raster, mereka dapat memerlukan acara seperti model populasi, tetapi orang-orang cenderung tidak tinggal diam.

Masalah kedua adalah bahwa peta panas terbatas untuk mempertimbangkan skala ruang tunggal, dan memilih skala ruang ini, yaitu ukuran kernel atau tingkat peluruhan, dapat rumit dan tergantung pada tujuan penelitian, tetapi harus dibenarkan. . Jika intinya adalah untuk mengidentifikasi pusat dari kelompok terkuat, dan skala di mana ia terjadi (mungkin untuk mengidentifikasi sumber wabah penyakit, dan faktor penyebarannya) pilihan yang lebih baik mungkin dengan mempertimbangkan berbagai skala. Dengan bobot yang sesuai dengan skala / area untuk menghasilkan raster 3 dimensi, di mana maksimum lokal dalam raster skala ruang 3D menunjukkan lokasi pusat cluster dan ukurannya masing-masing, dan persistensi antara skala.

Craig Stobbs
sumber
1
Anda membuat beberapa poin yang sangat valid. Kedua masalah ini sebenarnya adalah masalah klasik geografi. Masalah pertama terkait dengan interpretasi yang mendasari non-keseragaman ruang, yaitu karena distribusi orang tidak seragam (dengan beberapa daerah berpenduduk jarang), peluang untuk kejahatan juga tidak seragam. Pola yang satu dipaksakan oleh pola yang lain. Masalah skala kedua adalah bagian dari Masalah Unit yang Dapat Dimodifikasi (MAUP) yang akan memengaruhi setiap ukuran yang bergantung pada area untuk diukur, misalnya kepadatan. Ini adalah masalah klasik dalam sebagian besar pekerjaan geografis.
WhiteboxDev
Pertimbangkan juga poligon Thiessen sebagai tampilan pertama pada pola spasial. Mereka cukup sederhana, tetapi juga mampu menunjukkan banyak detail, adalah interpolator yang tepat, dan tampaknya dapat mengatasi beberapa masalah skala yang mengganggu metode interpolasi yang lebih canggih lainnya - gislandscapeecology.blogspot.com/2016/04/…
Tom Dilts