前言

這兩天幫一個朋友處理了些 nc 數(shù)據(jù)，本以為很簡單的事情，沒想到里面涉及到了很多的細節(jié)和坑，無論是“知難行易”還是“知易行難”都不能充分的說明問題，還是“知行合一”來的更靠譜些，既要知道理論又要知道如何實現(xiàn)，于是經(jīng)過不太充分的研究后總結(jié)成此文，以記錄如何使用 python 處理 nc 數(shù)據(jù)。

一、nc 數(shù)據(jù)介紹

nc 全稱 netCDF（The Network Common Data Form），可以用來存儲一系列的數(shù)組，就是這么簡單（參考https://www.unidata.ucar.edu/software/netcdf/docs/netcdf_introduction.html）。

既然 nc 可以用來一系列的數(shù)組，所以經(jīng)常被用來存儲科學觀測數(shù)據(jù)，最好還是長時間序列的。

試想一下一個科學家每隔一分鐘采集一次實驗數(shù)據(jù)并存儲了下來，如果不用這種格式存儲，時間長了可能就需要創(chuàng)建一系列的 csv 或者 txt 等，而采用 nc 一個文件就可以搞定，是不是很方便。

更方便的是如果這個科學實驗與氣象、水文、溫度等地理信息稍微沾點邊的，完全也可以用 nc 進行存儲， GeoTiff 頂多能多存幾個波段（此處波段可以認為是氣象、水文等不同信號），而 nc 可以存儲不同波段的長時間觀測結(jié)果，是不是非常方便。

可以使用 gdal 查看數(shù)據(jù)信息，執(zhí)行：

gdalinfo name.nc

即可得到如下信息：

Driver: netCDF/Network Common Data FormatFiles: test.ncSize is 512, 512Coordinate System is `'Subdatasets: SUBDATASET_1_NAME=NETCDF:"test.nc":T2 SUBDATASET_1_DESC=[696x130x120] T2 (32-bit floating-point) SUBDATASET_2_NAME=NETCDF:"test.nc":PSFC SUBDATASET_2_DESC=[696x130x120] PSFC (32-bit floating-point) SUBDATASET_3_NAME=NETCDF:"test.nc":Q2 SUBDATASET_3_DESC=[696x130x120] Q2 (32-bit floating-point) SUBDATASET_4_NAME=NETCDF:"test.nc":U10 SUBDATASET_4_DESC=[696x130x120] U10 (32-bit floating-point) SUBDATASET_5_NAME=NETCDF:"test.nc":V10 SUBDATASET_5_DESC=[696x130x120] V10 (32-bit floating-point) SUBDATASET_6_NAME=NETCDF:"test.nc":RAINC SUBDATASET_6_DESC=[696x130x120] RAINC (32-bit floating-point) SUBDATASET_7_NAME=NETCDF:"test.nc":SWDOWN SUBDATASET_7_DESC=[696x130x120] SWDOWN (32-bit floating-point) SUBDATASET_8_NAME=NETCDF:"test.nc":GLW SUBDATASET_8_DESC=[696x130x120] GLW (32-bit floating-point) SUBDATASET_9_NAME=NETCDF:"test.nc":LAT SUBDATASET_9_DESC=[130x120] LAT (32-bit floating-point) SUBDATASET_10_NAME=NETCDF:"test.nc":LONG SUBDATASET_10_DESC=[130x120] LONG (32-bit floating-point)Corner Coordinates:Upper Left ( 0.0, 0.0)Lower Left ( 0.0, 512.0)Upper Right ( 512.0, 0.0)Lower Right ( 512.0, 512.0)Center ( 256.0, 256.0)

每一個 SUBDATASET 表示記錄的是一種格式的數(shù)據(jù)（氣象、水文等等），如果要想查看此 SUBDATASET 的具體信息，可以執(zhí)行：

gdalinfo NETCDF:name.nc:SUBDATASET_NAME