DataAssassin
全网同号,数据免费获取请关注公众号「DataAssassin」,如有问题请添加公众号联系作者。感谢您的关注。
之前我们学习了如何从国家青藏高原科学数据中心下载数据,但我们会发现下载下来的气象数据是nc格式。今天我们就来学习用python将nc数据转成tif数据。
nc数据
「NetCDF」全称为Network Common Data Format,即网络通用数据格式,它是由美国大学大气研究协会的Unidata项目科学家针对科学数据的特点开发的,是一种面向数组型并适于网络共享的数据描述和编码标准,特别适用于包含多维数组和元数据的数据集。NetCDF 文件可以存储各种类型的数据,如气象数据、海洋数据、地理信息系统数据等。它可以在不同平台上进行读取和写入,因此非常适合用于交换和共享科学数据。
数据组成
NetCDF 文件通常具有以下组成部分:
-
维度(Dimensions): 描述数据的不同维度,例如时间、经度、纬度等。 -
变量(Variables): 存储实际数据的多维数组,这些数组与维度相关联。 -
属性(Attributes): 与数据集或变量相关联的元数据信息,用于描述数据的含义、单位等。
python查看nc数据
import netCDF4 as nc
# 打开NetCDF文件
file_path = 'path_to_your_file.nc'
dataset = nc.Dataset(file_path, 'r')
# 获取维度和变量
latitudes = dataset.variables['lat'][:]
longitudes = dataset.variables['lon'][:]
pm25_data = dataset.variables['PM2.5'][:]
# 关闭文件
dataset.close()
在上述示例中,假设你的NetCDF文件包含名为'lat'、'lon'和'PM2.5'的变量,分别表示纬度、经度和PM2.5数据。你可以通过variables属性从数据集中提取这些变量的数据。
通过打印latitudes、longitudes和pm25_data,你可以查看这些变量的数据。latitudes和longitudes将是数组,而pm25_data也将是一个数组,其中存储了对应的PM2.5数据。
nc转tif
下面就看看用python怎么实现nc转tif
导入必要的库
import os
from osgeo import gdal
import netCDF4 as nc
import numpy as np
from glob import glob
from osgeo import osr
# 导入所需库
# os 用于基本操作系统操作
# gdal 用于处理地理空间数据
# netCDF4 用于读取 netCDF 文件
# numpy 用于数值操作
# glob 用于文件路径匹配
# osr 用于处理空间参考信息
设置参数
# 定义工作路径和输出路径
file_path = r'.\data\2021_NO2.nc'
OutPath = r'.\data\'
# 指定空气污染物类型
AP = 'PM2.5'
# 定义空间分辨率
SP = 0.01 # 度
# 如果输出路径不存在,则创建
if not os.path.exists(OutPath):
os.makedirs(OutPath)
转换格式
f = nc.Dataset(file_path) # 打开 netCDF 文件
data = np.array(f[AP][:]) # 提取空气污染物数据
data[data == 65535] = np.nan # 将值为 65535 的数据替换为 NaN
lon = np.array(f['lon'][:]) # 提取经度数据
lat = np.array(f['lat'][:]) # 提取纬度数据
LonMin, LatMax, LonMax, LatMin = lon.min(), lat.max(), lon.max(), lat.min() # 计算范围和极值
N_Lat = len(lat)
N_Lon = len(lon)
Lon_Res = SP
Lat_Res = SP
# 提取文件名并构建输出路径
fname = os.path.basename(file_path).split('.nc')[0]
outfile = OutPath + '/{}.tif'.format(fname)
driver = gdal.GetDriverByName('GTiff') # 获取 GeoTIFF 驱动
# 创建输出 GeoTIFF 文件
outRaster = driver.Create(outfile, N_Lon, N_Lat, 1, gdal.GDT_Float32)
# 设置地理变换参数
outRaster.SetGeoTransform([LonMin - Lon_Res / 2, Lon_Res, 0, LatMax + Lat_Res / 2, 0, -Lat_Res])
sr = osr.SpatialReference()
sr.SetWellKnownGeogCS('WGS84') # 使用 WGS84 坐标系
outRaster.SetProjection(sr.ExportToWkt()) # 设置投影信息
outRaster.GetRasterBand(1).WriteArray(data) # 将数据写入栅格波段
print(fname + '.tif', '处理完成')
# 释放内存
del outRaster
f.close()
总结
转换完成后就能得到nc文件的tif格式啦。感兴趣的可以下载数据自己试试。下面是批量处理nc转tif的python代码。
获取方法
如有需要,请关注微信公众号「DataAssassin」后,后台回复「016」领取。
资料收集不易,请及时保存,如有侵权请第一时间联系作者。切勿用作商业用途,谢谢大家的支持。感兴趣的请「关注、点赞、收藏」哦!
本文由 mdnice 多平台发布