Echo Coverage (Part 1)

Dari beberapa paper yang sudah saya pelajari (baca: sebagian saya pelajari), ada beberapa gambar menarik yang menggunakan data CAPPI radar. Salah satunya adalah Echo Coverage (EC), yang menunjukkan tutupan echo/sinyal balik dari presipitasi yang terdeteksi oleh radar. 

Salah satu paper menarik yg membahas soal EC adalah paper Mori-san (kebetulan saya jadi co-authornya .. ehm ) : Convective System Developed along the Coastline of Sumatera Island, Indonesia, Observed with an X-band Doppler Radar during the HARIMAU2006 Campaign

Analisis echo coverage ini biasanya ditunjukkan dengan beberapa plot, antara lain : 

1. Time-Height cross section. Pada grafik ini, sumbu X adalah waktu, sumbu Y adalah ketinggian, dan sumbu Z adalah echo coverage.
2. Time-Time cross section. Pada grafik ini, sumbu X adalah waktu (hari), sumbu Y adalah waktu (jam) dan sumbu Z adalah echo coverage.

Analisis EC sangat bermanfaat untuk mengetahui pola spasial (statiform atau cumuliform) dan temporal (diurnal atau semi-diurnal) dari curah hujan di suatu wilayah.

Saya sebenarnya sudah pernah melakukan analisis data yang mirip, dengan menggunakan data satelit TRMM produk TRMM3B43 dan GSMaP NRT/MWR, dengan menggunakan Time vs Time plot, tapi dimensi ketiga yang saya gunakan adalah intensitas rata-rata curah hujan (rainrate).

Kali ini rencananya saya akan melakukan analisis dengan data XDR, selama perionde IOP2011 (Desember 2011). Well, ini belum pernah saya lakukan sebelumnya, jadi waktu pagi ini habis dengan coretan di kertas untuk menghitung EC sebelum dikonversi menjadi kode program. 

Masalah : Data radar adalah data grid 641x641x20. Karena yang dibutuhkan adalah data dalam coverage 160 km, maka tidak semua data akan dimasukkan dalam perhitungan, hanya data yang lebih kecil dari diameter 160 km yg diambil.

Algoritma :

1. Buka data binary radar.
2. Assign data radar, kedalam array 3 dimensi (x,y,z).
3. Mulai hitung jarak tiap grid dari pusat radar. Grid yang jaraknya lebih kecil dari 160 km akan dihitung, yg lain akan diabaikan. Phytagoras theorm will do the trick ... 
4. Periksa grid yang diperoleh dari langkah 3, apakah memiliki echo atau tidak. Bila ada, hitung berapa banyak grid yang memiliki echo.
5. EC diperoleh dengan membandingkan jumlah grid yang memiliki echo dengan jumlah total grid dalam diameter 160 km.
6. Ulangi langkah 1-5 untuk tiap data radar (10 menitan).

to be continued ....