python处理FNL数据的grib文件和nc文件(纬度存储的问题)

2023-10-30 20:59

本文主要是介绍python处理FNL数据的grib文件和nc文件(纬度存储的问题),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

在使用python处理FNL数据时,2007年及之前的数据存储为grib文件,2007年以后为grib2文件,在数据处理时,将grib2转化为nc文件处理Windows上python读取grib2文件(不用Linux),对grib文件则安装了Anaconda环境处理Windows下xarray+cfgrib读取grib文件。但在根据经纬度读取变量的时候,发现不同文件的存储不同,而变量是根据经纬度存储的下标去获取值,所以需要注意一下不同的文件的存储内容。

1.python获取FNL的nc文件对应的变量

nc_obj =Dataset('F:/Zhu/download/2000-2015/NC/NC2008/fnl_20080113_06_00.nc')
nc_obj.variables[name][0][lat][lon]  根据变量名,经纬度提取变量值

从FNL的nc文件中获取变量的方式如上,需要根据变量名,第一个是time,这里是一维存储,下标只有0即可。然后是纬度,经度获取变量的值。

先看一下,FNL的nc文件中纬度的存储:

print(nc_obj.variables['latitude'][:])

在这里插入图片描述

从结果可以看到,纬度是整数存储,从-90到90,共181维。

再看一下FNL的nc文件中经度的存储:

print(nc_obj.variables['longitude'][:])

在这里插入图片描述

可以从结果看到经度是从0到359度,整数存储,共360维。

所有数据都存储在181x360的整数网格点上。

看一下存储array下标与纬度的关系:

lat = 12
lon = 128
print(nc_obj.variables['latitude'][lat])
print(nc_obj.variables['longitude'][lon])

赋予纬度经度特定的值,观察值在array中的取值

-78.0
128.0

可以发现下标与真实的纬度值之间差了90

print(nc_obj.variables['latitude'][lat + 90])
print(nc_obj.variables['longitude'][lon])
12.0
128.0

经度的下标则与真实的值对应。

所有变量的提取都是根据对应的经纬度下标
以读取UGRD_200mb示例:

print(nc_obj.variables['UGRD_200mb'][:])

在这里插入图片描述
结果则是全部网格点上的数据,提取latitude:-90,longitude:0的U200则如下:

print(nc_obj.variables['UGRD_200mb'][0][0][0])
-3.8

因此在计算真实的纬度后,从FNL的nc文件中获取对应变量的的值时,纬度需要加90即真实纬度对应的下标。

2.python获取FNL的grib文件对应的变量

grib文件变量的提取需要根据grib文件以及对应需要提取的level获取数据,

def readgrib_level(filein, lev):data = xr.open_dataset(filein, engine='cfgrib',backend_kwargs={'filter_by_keys': {'typeOfLevel': 'isobaricInhPa', 'level': lev}})return data

获取grib文件的纬度存储

filein = 'F:/Zhu/download/2000-2015/FNL_Data/fnl_20000101_00_00'
nc_obj = readgrib_level(filein, 200)
print(nc_obj.variables['latitude'])

在这里插入图片描述可以发现FNL的grib文件与nc文件纬度存储不同,从90到-90整数存储,维度为181。

再看一下经度存储:

print(nc_obj.variables['longitude'])

在这里插入图片描述
可以看到存储与nc文件相同,从0到359,360维。

以读取200hpa的u为例:

print(nc_obj.variables['u'][:])

在这里插入图片描述
根据经纬度下标获取变量

print((nc_obj.variables['u'][0][0]).values)
4.4

与nc文件提取变量相同,也是根据经纬度下标提取变量。

因此对于grib文件的数据,在计算真实维度后,需要使用90-真实纬度,得到提取变量对应的下标位置:

lat = 12
print((nc_obj.variables['latitude'][lat]).values)
78.0
print((nc_obj.variables['latitude'][90-lat]).values)
12.0

因此对于FNL的nc文件和grib文件存储不同,需要注意。

这篇关于python处理FNL数据的grib文件和nc文件(纬度存储的问题)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/310558

相关文章

C#使用SQLite进行大数据量高效处理的代码示例

《C#使用SQLite进行大数据量高效处理的代码示例》在软件开发中,高效处理大数据量是一个常见且具有挑战性的任务,SQLite因其零配置、嵌入式、跨平台的特性,成为许多开发者的首选数据库,本文将深入探... 目录前言准备工作数据实体核心技术批量插入:从乌龟到猎豹的蜕变分页查询:加载百万数据异步处理:拒绝界面

Python使用自带的base64库进行base64编码和解码

《Python使用自带的base64库进行base64编码和解码》在Python中,处理数据的编码和解码是数据传输和存储中非常普遍的需求,其中,Base64是一种常用的编码方案,本文我将详细介绍如何使... 目录引言使用python的base64库进行编码和解码编码函数解码函数Base64编码的应用场景注意

Python基于wxPython和FFmpeg开发一个视频标签工具

《Python基于wxPython和FFmpeg开发一个视频标签工具》在当今数字媒体时代,视频内容的管理和标记变得越来越重要,无论是研究人员需要对实验视频进行时间点标记,还是个人用户希望对家庭视频进行... 目录引言1. 应用概述2. 技术栈分析2.1 核心库和模块2.2 wxpython作为GUI选择的优

C# WinForms存储过程操作数据库的实例讲解

《C#WinForms存储过程操作数据库的实例讲解》:本文主要介绍C#WinForms存储过程操作数据库的实例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、存储过程基础二、C# 调用流程1. 数据库连接配置2. 执行存储过程(增删改)3. 查询数据三、事务处

Springboot处理跨域的实现方式(附Demo)

《Springboot处理跨域的实现方式(附Demo)》:本文主要介绍Springboot处理跨域的实现方式(附Demo),具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不... 目录Springboot处理跨域的方式1. 基本知识2. @CrossOrigin3. 全局跨域设置4.

Java利用JSONPath操作JSON数据的技术指南

《Java利用JSONPath操作JSON数据的技术指南》JSONPath是一种强大的工具,用于查询和操作JSON数据,类似于SQL的语法,它为处理复杂的JSON数据结构提供了简单且高效... 目录1、简述2、什么是 jsONPath?3、Java 示例3.1 基本查询3.2 过滤查询3.3 递归搜索3.4

Python如何使用__slots__实现节省内存和性能优化

《Python如何使用__slots__实现节省内存和性能优化》你有想过,一个小小的__slots__能让你的Python类内存消耗直接减半吗,没错,今天咱们要聊的就是这个让人眼前一亮的技巧,感兴趣的... 目录背景:内存吃得满满的类__slots__:你的内存管理小助手举个大概的例子:看看效果如何?1.

Python+PyQt5实现多屏幕协同播放功能

《Python+PyQt5实现多屏幕协同播放功能》在现代会议展示、数字广告、展览展示等场景中,多屏幕协同播放已成为刚需,下面我们就来看看如何利用Python和PyQt5开发一套功能强大的跨屏播控系统吧... 目录一、项目概述:突破传统播放限制二、核心技术解析2.1 多屏管理机制2.2 播放引擎设计2.3 专

Python中随机休眠技术原理与应用详解

《Python中随机休眠技术原理与应用详解》在编程中,让程序暂停执行特定时间是常见需求,当需要引入不确定性时,随机休眠就成为关键技巧,下面我们就来看看Python中随机休眠技术的具体实现与应用吧... 目录引言一、实现原理与基础方法1.1 核心函数解析1.2 基础实现模板1.3 整数版实现二、典型应用场景2

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很