python将nc文件转成tif

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 多平台发布