python读取nc数据
 
读取nc数据主要是用到netCDF4库。
 
#最基础的读取
nc_obj = nc.Dataset(filename)
#这里filename自己定义即可
 
这里已经将nc文件读取进来了，之后我们可以查看其中的各个维度以及数据。
 
#读取经度纬度值
lat=(nc_obj.variables['lat'][:])
lon=(nc_obj.variables['lon'][:])
 
这里variables方法即是读取变量的功能，一般nc文件中经度和纬度会命名为‘lon’和‘lat’，不确定的话可以阅读一下相关产品的说明书进行查询。
 或者可以直接调用variables.keys方法：
 
print(nc_obj.variables.keys())
 
输出结果：
 
 这里的变脸分别是经度、纬度、时间和蒸散发量。
 那么我们想读取蒸散发数值的话，就是和前面类似的操作：
 
Data=(nc_obj.variables[str(flag)][:])
 
看一下读出来的结果的形状：
 
print(Data.shape)
 

 这个产品的记录是以月份为单位，一共是40年的。第二个是纬度，第三位是经度，都是以0.5°为单位。
 到这里，数据已经存入数组Data中了，我们可以对它进行后续的操作。
 
补充一个题外话，如果想要灵活一点，让我们的程序可以每次读取我们选定的文件，而不需要我们每次更改路径，则可以使用win32ui库。
 
import win32ui
# 0代表另存为对话框，1代表打开文件对话框
dlg = win32ui.CreateFileDialog(1)
 
# 默认目录
dlg.SetOFNInitialDir('D:\Recently\Water_recycle\Data') 
 
# 显示对话框
dlg.DoModal()
 
# 获取用户选择的文件全路径
filename = dlg.GetPathName()

nc_obj = nc.Dataset(filename)

 
这部分文件读取的内容一开始是在CSDN上找到的，但是后来找不到了，所幸我的程序里还留着，所以就搬运了，如果有朋友找到了之前那个大佬的帖子，麻烦指一下路，我补充一下应用，非常感谢！！

本篇博客主要介绍采用MATLAB读取nc数据并进行显示。
 
首先是显示经纬度：
 
示例代码：
 
lon = ncread('met_em.d02.2018-09-12_00_00_00.nc', 'XLONG_C');
lat = ncread('met_em.d02.2018-09-12_00_00_00.nc', 'XLAT_C');
plot(lon, lat, 'r')
 
运行结果：
 
 
其次是显示地表温度：
 
示例代码：
 
lon = ncread('met_em.d02.2018-09-12_00_00_00.nc', 'XLONG_C');
lat = ncread('met_em.d02.2018-09-12_00_00_00.nc', 'XLAT_C');
skintemp = ncread('met_em.d02.2018-09-12_00_00_00.nc', 'SKINTEMP');
disp(skintemp)
skintemp1=skintemp(:,:);
pcolor(lat(1: 111, 1: 96),lon(1: 111, 1: 96),skintemp1);
 
运行结果：
 
 
欢迎关注我的公众号：
 
                     编程技术与生活（ID：hw_cchang）
 
 
 

 
netcdf读取NC数据，C语言+VS2017环境部署
 
简单说明
前期准备
环境、软件
  
环境搭建
环境配置
测试代码
 

 
简单说明
 
本来主要使用C#开发后台的，对C也只是简单会用，因C#调用各种第三方组件都多少存在一些无法解决的问题，尝试使用C来自己封装，所以编写的啰嗦一些，主要针对我这种小白，写的不对的地方，请大神告知。
 
前期准备
 
简单介绍环境软件介绍、下载地址等等。
 
环境、软件
 
使用环境、软件：
 
Panoply：也是使用netcdf编写的，读取NC文件，用来查看数据，也可以使用其他的；
vs：版本不要太古老，我用的是vs2017；
netcdf：我使用的是 4.6.1；
 
下载链接：
 
Panoply：https://www.giss.nasa.gov/tools/panoply/download/PanoplyWin-4.10.4.zip (官方地址)
vs：https://visualstudio.microsoft.com/zh-hans/downloads/ (官方地址)
 http://msdn.itellyou.cn/ (资源比较全，不保证安全)
netcdf：https://www.unidata.ucar.edu/downloads/netcdf/ftp/netCDF4.6.1-NC3-32.exe (官方地址)
 https://www.unidata.ucar.edu/downloads/netcdf/index.jsp (其他选择)
 
环境搭建
 
环境配置
 
项目属性\VC++目录\包含目录 添加 C:\Program Files (x86)\netCDF 4.6.1\include （安装位置的include目录）
项目属性\VC++目录\库目录 添加 C:\Program Files (x86)\netCDF 4.6.1\lib （安装位置的lib目录）
项目属性\链接器\附加库目录 添加 C:\Program Files (x86)\netCDF 4.6.1\bin\netcdf.dll （安装位置的bin\netcdf.dll文件）
项目属性\链接器\输入\附加依赖项 添加 netcdf.lib
 
必须以Release方式编译
 必须以Release方式编译
 必须以Release方式编译
 
测试代码
 
test.c
 
#include <stido.h>
int ncid;
//获取字段数据
int getColData(char *colname, float *dataArray)
{
	int varid;
	nc_inq_varid(ncid, colname, &varid);
	nc_get_var_float(ncid, varid, dataArray);
	return 0;
}
//获取字段内容数量
int getColLen(char *colname)
{
	int dimid;
	size_t size;
	//获取字段维度
	nc_inq_dimid(ncid, colname, &dimid);
	nc_inq_dimlen(ncid, dimid, &size);
	return (int)size;
}
//获取nc数据一维数据
void getNcData(char *ncFilePath,char *Name)
{
	int varid,i,lon_length;
	float *lon_data = NULL;
	//打开nc文件
	nc_open(ncFilePath,0,&ncid);
	lon_length = getColLen(Name);
	//声明存储数据空间
	lon_data = (float *)malloc(lon_length * sizeof(float));
	if (NULL == lon_data)
	{
		return NULL;
	}
	getColData(Name, lon_data);
	for(i = 0; i < lon_length; i++)
	{
		printf("%.4f\t",lon_data[i]);
	}
	
}
int main()
{
	getNcData("d:\******.nc","longitude");
	system("pause");
	return 0;
}

file_ID = NCDF_OPEN(theseFiles[fidx],/NOWRITE)   ;open netCDF file for READ only
    Tag = NCDF_INQUIRE(file_ID)
    timeid = NCDF_VARID(file_ID,'time')             ;initial_time0initial_time0_hours  initial_time0
    NCDF_VARGET, file_ID, timeid, time              ;获取 时间维变量
    nt = N_ELEMENTS(time)
    ;print,time
    latid = NCDF_VARID(file_ID,'latitude')               ;g0_lat_2
    NCDF_VARGET, file_ID, latid, latitude           ;获取 时间维变量
    nlat = N_ELEMENTS(latitude)                     ;纬向维长度
    ;print,'dsklfjdsjfkldjsfkljl'+latitude
    lonid = NCDF_VARID(file_ID,'longitude')                ;g0_lon_3
    NCDF_VARGET, file_ID, lonid, longitude           ;获取 时间维变量
    nlon = N_ELEMENTS(longitude)                     ;径向维长度

    t2mid = NCDF_VARID(file_ID,'t2m')                ;g0_lon_3
    NCDF_VARGET, file_ID, t2mid, t2m           ;获取 温度维变量
    nt2m = N_ELEMENTS(t2m)
    
    ncdf_attget,file_ID,t2mid,'scale_factor',a
    ncdf_attget,file_ID,t2mid,'add_offset',b
    year=STRMID(thesefiles[fidx],0,4)
    month=STRMID(thesefiles[fidx],5,2)
    outtmp=fltarr(480,15,band[fidx]*4)
    outtmp[*,*,*]=t2m[*,*,*]*a+b
    NCDF_CLOSE, file_ID
    print,'  All parameters of file ' + thesefiles[fidx] + ' have been writen out!'
    outfilename = STRCOMPRESS(outfilepath +string(year,FORMAT='(I04)')+'_'+string(month,FORMAT='(I02)'),/remove_all)
    openw,var_lun,outfilename,/get_lun
    writeu,var_lun,outtmp;
    close,var_lun
    free_lun,var_lun
 
scale_factor','add_offset一定要注意获取，也要注意最后读取出来的单位。