地球科学常见数据的处理技术

地球科学指一切研究地球的科学，主要包括地质学、地理学，以及其它衍生学科。各学科通常会以物理、地理、地质、气象、数学、化学、生物的角度研究地球。它和人类的生活息息相关，人们手上所戴的黄金饰品和钻石，都是来自地球的矿产资源；盖房子所用的砂、石、水泥，其原料也是来自地球；所吃的鱼虾，大都取自海洋；气温的变化影响生活甚巨；天体的运行，也时时刻刻影响着我们。因此，地球科学是一门很基础、很重要的学科。

地球科学的范围很广，涵盖地质学、海洋学、气象学和天文学等领域。地质学在探讨地球的历史与各部分组成，包括其演化和各种矿学、岩石以及矿产的分布；海洋学在研究海水的运动、海水的物理与化学性质及海底地形；气象学在分析大气的组成、构造和运动；而有关地球起源、太阳系的形成和天体的运动变化，乃至宇宙的演化，均属天文学的研究范围。以陨石撞击地球为例：高温高压撞击地球的结果，势必引起地形与地质的变化；飞扬在大气中的粉尘微粒会遮蔽阳光，大气和海水温度因而降底。因此，看似简单的天文事件，却引起地质、气象和海洋的变化，可见各领域关系密切、环环相扣。

在地球科学中，不同数据根据具体学科的特色存储为多种数据格式。在科研工作中需要将多种数据进行综合使用分析，因此需要寻找学习通用的数据格式解决方法，把研究的精力聚焦到具体科学问题上。

针对上述问题，选取大气科学、水文学和生态学常见的数据进行讲解。

讲解的主要格式：

1. 全球大气再分析数据netCDF
2. 雪深ASCII
3. 积雪覆盖ASCII/TIFF
4. 海温数据netCDF
5. 植被指数数据NDVI
6. 土地利用数据HDF

需要的处理工具：

Anaconda 5.0+（python 3.6）、xarray==0.13、netcdf4==1.5.3、rasterio、pandas、pyhdf、fiona、shapely、gdal

学习的预备知识：

一定的Linux和Python的学习基础。

一、全球大气再分析数据-ERA5数据

ECMWF中心推出的ERA5全球大气再分析数据提供了大量大气、陆地和海洋气候变量的逐小时数据。这些数据在30km网格上覆盖了全球，在时间跨度上从1979至今。该数据能够提供全球范围的格点气象数据。

将针对该数据介绍以下内容：

1. 多年数据的读取

一般在下载过程中，会将按照最小的时间单位获取文件。因此，需要将数据进行后处理进行合并，之后开展后续的操作。在本案例中，将对气温变量进行读取。

1. 趋势分析

在多年气温数据的基础上，使用mann kendall趋势检验方法，判断趋势；使用线性回归方法计算变化倾向率。

1. 多时间尺度统计

利用小时数据，计算日、月、年尺度数据。

二、北半球雪深数据-加拿大雪深数据

加拿大气象中心（CMC）制作了全球每日雪深业务分析数据。数据时间分辨率为1天，空间分辨率为0.25°。数据包括雪深和估计雪水当量（SWE）2个变量。该数据可以为全球和区域尺度冰冻圈积雪相关水文和生态研究提供支持。

将针对该数据介绍以下内容：

1. 将ASCII数据转化为netCDF数据

将日数据进行整理。使用multiprocess工具并行转化netCDF格式数据。

1. 分析北半球雪深的气候态空间分布

使用Xarray工具计算多年、年和月的均值。

1. 分析时间序列的趋势变化及与气温的关系

使用线性回归分析月尺度上趋势；计算雪深变化序列与气温之间的相关系数。

三、积雪覆盖数据-IMS积雪覆盖数据

美国国家冰雪中心（NSIDC）从1997年2月至今的北半球雪盖和海冰的地图。这些数据以ASCII文本和GeoTIFF格式提供，并提供有三种不同的分辨率：1公里、4公里和24公里。该数据可以为全球和区域尺度冰冻圈积雪相关水文和生态研究提供支持。

将针对该数据介绍以下内容：

1. 将ASCII数据转化为netCDF数据

将日数据进行整理。使用multiprocess工具并行转化netCDF格式数据。

1. 分析新疆北疆、青藏高原和东北地区气候态积雪分布

使用Xarray工具计算多年、年和月的均值。使用basemap进行可视化。

1. 分析新疆北疆、青藏高原和东北地区积雪面积变化规律

使用Xarray工具进行汇总，绘制不同子区的面积图时间序列。

四、海温数据-NOAA OISST

NOAA 1/4°每日最佳内插海面温度（OISST）是一个长期的气候数据记录，它将来自不同平台（卫星、船舶、浮标和Argo浮标）的观测数据纳入一个定期的全球网格。该数据集经过插值处理，以填补网格上的空白，并创建一个空间上完整的海面温度图。卫星和船舶的观测数据以浮标为参照，以补偿平台的差异和传感器的偏差。OISST数据空间分辨率0.25°，时间跨度从1981至今。该数据可以为气候变化、海~气关系等研究提供支持。

将针对该数据介绍以下内容：

1. 读取多年月海温数据

使用Xarray工具读取海温数据

1. 分析海温变化趋势
2. 全球海温与长江流域夏季降水关系

    分析CN05中长江流域夏季降水与全球海温之间的相关关系

五、植被指数数据-GIMMS NDVI

植被指数是对植被冠层绿叶中叶绿素吸收的光合有效辐射的辐射测量，因此是衡量一个地区生理功能的植被覆盖水平的指数。最新版本的GIMMS NDVI数据集横跨1981年7月至2015年12月。

课程中将针对该数据介绍以下内容：

A、读取多年生长季NDVI数据

B、计算青藏高原生长季NDVI与降水关系

六、土地利用数据-MODIS土地利用数据

地覆盖动态产品包括标志着季节性周期的植被生长、成熟和衰老的时间层。MODIS Terra+Aqua组合土地覆盖产品包含了五个不同的土地覆盖分类方案，通过监督决策树分类方法得出。主要的土地覆盖方案确定了IGBP定义的17个等级，包括11个自然植被等级，3个人为改变的等级，以及3个非植被等级。

将针对该数据介绍以下内容：

A、读取多年土地利用数据

B、裁剪长江流域土地利用

C、分析长时间序列长江流域土地利用变化