问题要求：

编写一个函数完成以下任务：截取经度在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参数全解

实验总结与体会

本次实验是我第一次使用python来处理数据集。通过本次实验，

加深了对nc文件的认识。

对python处理数据的方法有了一定了解。

对python用于数据处理计算的库例如netCDF4和numpy有了一定了解。

对于python的文件操作有了一定认识。import pandas

import netCDF4 as ne

import numpy as np

import csv

import glob

import sys

def 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=0

writer.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')