Cara memvisualisasikan Set Data NetCDF
Misalkan saya memiliki NetCDFfile yang agak besar yang saya buka sebagai dataset menggunakan xarray.
import xarray
nc_in = r"....\D_P_WGS84_comp.nc"
nc_ds = xarray.open_dataset(nc_in)
Ini terdiri dari beberapa variabel raster, masing-masing dengan satu band, koordinat x, y dan crs (4326):
nc_ds
Misalkan saya ingin memvisualisasikan ini dengan sangat sederhana di Python Notebook saya. Saya tidak benar-benar ingin menambahkan peta dasar atau legenda atau semacamnya. Saya mencoba menjalankan show(nc_ds) but throws an errors AttributeError: objek 'Dataset' tidak memiliki atribut 'shape'`. Saya juga memperhatikan bahwa dimensinya membeku. Apa cara termudah untuk memvisualisasikan kumpulan data ini atau satu variabelnya dengan Python? Apakah ada mod atau ekstensi paket yang dapat melakukan ini?
Jawaban
Karena sebuah Datasetcan berisi lebih dari satu variabel, tidaklah mudah bagi perangkat lunak untuk mengetahui apa yang Anda ingin tampilkan. Namun, jika Anda memilih salah satu variabelnya DataArray, itu dapat dengan mudah ditampilkan (jika dalam format yang sesuai).
Variabel a Datasettersedia melalui data_varsatributnya. Ini mengekspos variabel sebagai DataArrays dalam sebuah dict dengan nama variabel sebagai kunci.
Selangkah demi selangkah:
variables = nc_ds.data_vars
f02_nbr_e0002 = variables['F02_NBR_E0002']
show(f02_nbr_e0002)
Satu perjalanan yang efisien:
show(nc_ds.data_vars['F02_NBR_E0002'])