NetCDF (Network Common Data Form) adalah format file yang digunakan untuk menyimpan data ilmiah multidimensi (misal: latitude, longitude, waktu), seperti data hujan dan data remote sensing lainnya (seperti citra satelit atau data radar). Format ini memungkinkan pengaturan data secara efisien serta memudahkan untuk dianalisis untuk kebutuhan riset dan dan visualisasi.
NetCDF menyediakan keuntungan dalam mengatur data secara efisien, mempertahankan metadata, dan mendukung operasi ilmiah seperti pemotongan area, subset data, dan penggabungan data yang berbeda.
Dalam artikel ini saya akan berikan panduan untuk membuka data netCDF menggunakan pemrograman python serta memvisualisasikannya. Anda tidak membutuhkan program lainnya untuk membuka dan memvisualisasi file netCDF. Python (dan R) sudah sangat powerful untuk melakukan tugas ini, dengan cara yang cukup mudah!
Untuk mendapatkan tutorial install dan setting python dapat dilihat disini.
1. Install dan import library
Pastikan Anda telah menginstal paket-paket Python yang diperlukan, seperti numpy, netCDF4, dan matplotlib. Anda dapat menginstalnya menggunakan Anaconda (jika anda menggunakan environment Anaconda) dengan menjalankan perintah berikut:
conda install -c anaconda xarray
Impor modul-modul yang diperlukan:
import numpy as np
import xarray as xr
import matplotlib.pyplot as plt2. Buka netCDF
Buka file NetCDF menggunakan Dataset dari netCDF4.
Misal setelah anda mendownload file GPM IMERG DF pada tutorial sebelumnya, kemudian anda akan mendapatkan file dengan format berikut: “3B-DAY.MS.MRG.3IMERG.20210930-S000000-E235959.V06.nc4.nc4”
File diatas menunjukkan tanggal dan spesifikasi data lainnya. Buka file tersebut dengan mengetikkan code berikut:
nc_file ="C:/Users/vempi/Downloads/3B-DAY.MS.MRG.3IMERG.20210930-S000000-E235959.V06.nc4.nc4"d =xarray.open_dataset(nc_file)
Identifikasi variabel yang ingin Anda visualisasikan dalam file NetCDF. Anda dapat menggunakan ‘dataset’ untuk melihat daftar variabel yang tersedia.
print(d)
Untuk data netCDF GPM dengan variabel hujan, nama variabel adalah “precipitationCal“. Lihat gambar berikut.
Baca data variabel yang ingin Anda visualisasikan dari file NetCDF:
variable = d['precipitationCal']
3. Visualisasi
Visualisasikan data menggunakan metode plot dari Xarray:
variable.plot()
plt.title('Judul Plot')
plt.xlabel('Label sumbu x')
plt.ylabel('Label sumbu y')
plt.show()
Itulah langkah-langkah dasar untuk membuka dan memvisualisasi file NetCDF menggunakan Python. Anda dapat menyesuaikan kode tersebut sesuai dengan kebutuhan Anda, misalnya dengan menggunakan jenis plot atau tampilan visual yang berbeda.
Selamat mencoba!
Data diatas (GPM IMERG) masih membutuhkan setting tambahan jika dibuka menggunakan xarray. Salah satunya adalah penyesuaian latitude longitude, dimana terlihat plot perlu diputar (dengan menukar kordinat lat lon). Analisis netCDF Lebih detail akan diberikan selanjutnya.