使用Python截取nc文件数据保存到CSV文件
问题要求：
编写一个函数完成以下任务：截取经度在23°N-40°N，纬度在118°E-131°E范围内各属性不同深度的数据，使用Python中合适的数据结构将截取的数据保存到同名CSV文件中。（nc文件数据格式参见笔者其他文章）
实验内容（附代码）
- 实验数据介绍（通过实验介绍你对NC数据的认识）
nc文件即NetCDF全称为network Common Data Format，中文译法为“网络通用数据格式”。 一个NetCDF文件的结构包括以下对象：变量(Variables)，维(dimension)，属性(Attribute) 。
本次查看的nc文件其中包括 (['salinity', 'time', 'depth', 'lat', 'lon', 'water_temp']) 一共5个变量。其中’salinity’和’water_temp’两个变量有四个维度，表示为
‘salinity(time, depth, lat, lon)’和’water_temp(time, depth, lat, lon)’。每个变量有自己的属性值，例如纬度’lat’,表示为
float64 lat(lat)
long_name: Latitude
standard_name: latitude
units: degrees_north
axis: Y
NAVO_code: 1
_CoordinateAxisType: Lat
- 实验步骤（记录实验过程，配图加以说明）
实验1：简单的安装netCDF4包，然后用这个包下的Dateset函数来读取文件信息。
实验2：
1.设一个函数，然后用Dateset函数读取nc文件
2.获取变量相应数组集合-纬度经度温度深度
3.文件路径
4.创建csv文件，设置异常处理
这里用嵌套循环将符合条件的值writerow到创建的csv文件中。
5.’__name__==__main__’模块
这个模块的作用是既能保证当前的.py文件能直接运行，也能保证其可以作为模块被其他.py文件导入。
- 实验结果
成功生成了对应的csv文件。其中共有34345行数据。
- 实验参考资料
参考博客：
如何在PyCharm中安装第三方库(包)？
http://www.360doc.com/content/18/0822/17/11881101_780384182.shtml
NC文件读写方式
https://blog.csdn.net/ennaymonkey/article/details/62886843
Python读取nc文件
https://blog.csdn.net/showpingzhang/article/details/83384780
Python读写csv文件的几种方法 及 pandas.read_csv参数全解
https://blog.csdn.net/secondlieutenant/article/details/79157956
- 实验总结与体会
本次实验是我第一次使用python来处理数据集。通过本次实验，
- 加深了对nc文件的认识。
- 对python处理数据的方法有了一定了解。
- 对python用于数据处理计算的库例如netCDF4和numpy有了一定了解。
- 对于python的文件操作有了一定认识。
import pandas
import netCDF4 as ne
import numpy as np
import csv
import glob
import sysdef to_csv(source_file):#读取nc数据dataset=ne.Dataset(source_file)print(dataset.variables.keys())#获取相应数组集合--纬度经度温度深度lat_set = dataset.variables['lat'][:]lon_set = dataset.variables['lon'][:]temp_set=dataset.variables['water_temp'] [:]dep_set=dataset.variables['depth'][:]index = [] # 把满足条件的经纬度放入列表index之中for j in range(len(lat_set)): # j为纬度for k in range(len(lon_set)): # k为经度if lat_set[j] > 23 and lat_set[j] < 40:if lon_set[k] > 118 and lon_set[k] < 131:index.append((j, k)) # 插入满足条件的数据print('输出index列表：')print(index)print('-------------------------------------------------------------------')#文件名不含扩展名source_file=source_file.split('.')file_name=source_file[0]#创建csv目标文件try:#打开目标csv文件with open(file_name+'.csv','a',newline='') as targetFile:# 创建写入流writer = csv.writer(targetFile)# 写入表头writer.writerow(('lat', 'lon', 'temperature', 'depth',))# 写入数据for j in range(len(lat_set)): # j为纬度for k in range(len(lon_set)): # k为经度if lat_set[j] > 23 and lat_set[j] < 40:if lon_set[k] > 118 and lon_set[k] < 131:i=0writer.writerow((lat_set[j], lon_set[k], temp_set[0][i][j][k], dep_set[i]))targetFile.close()#关闭文件print('Get'+file_name+'.csv Success!')except Exception as e:#抛异常print('Error :'+str(e))if '__name__ ==__main__':print("start transfrom!")to_csv('20150101.nc')print('Transform successfully')