首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

xarray走向netCDF处理(一):数据结构及数据读取

提取物理量 文件中读取数据ds = xarray.open_dataset() 假如数据中含有一个名为var的物理量可以通过ds.var或ds[var]来获取 实例 此处使用的是ERA-Interim...xarray封装了matplotlib的部分绘图函数,一行代码就可以数据画出来,不过作为一个负责的公众号,还是用cartopy顺带加载了地图。...import numpy as np import xarray as xr import cartopy.crs as ccrs import cartopy.feature as cfeat from...gl.xlabels_top = False #关闭顶端标签 gl.ylabels_right = False #关闭右侧标签 gl.xformatter = LONGITUDE_FORMATTER #x设为经度格式...本文数据 https://apps.ecmwf.int/datasets/data/interim-full-moda/levtype=sfc/ 本文自MeteoAI微信公众号,点击文末阅读原文按钮即可跳转原文

3K112

xarray走向netCDF处理(一):数据结构及数据读取

提取物理量 文件中读取数据ds = xarray.open_dataset() 假如数据中含有一个名为var的物理量可以通过ds.var或ds[var]来获取 实例 此处使用的是ERA-Interim...xarray封装了matplotlib的部分绘图函数,一行代码就可以数据画出来,不过作为一个负责的公众号,还是用cartopy顺带加载了地图。...import numpy as np import xarray as xr import cartopy.crs as ccrs import cartopy.feature as cfeat from...cfeat.RIVERS.with_scale('50m'), zorder=1) ax.add_feature(cfeat.LAKES.with_scale('50m'), zorder=1) # 设置网格属性...gl.xlabels_top = False #关闭顶端标签 gl.ylabels_right = False #关闭右侧标签 gl.xformatter = LONGITUDE_FORMATTER #x设为经度格式

24K1712
您找到你想要的搜索结果了吗?
是的
没有找到

Python气象绘图教程—(十九)剖面图

其中x表示经度全球东西360经度分为了10800刻度,相当于一个经度被分为30份;y表示纬度,全球南北180纬度分为了5400份,也是一个纬度分为30份。...以经度为例,前面已经讲到一个经度分为30份,那么我们要画东经70-140的图,那就需要对经度数据切片,原理如下(纬度同理): 起始:(180+70)×30=7500(在前面属性可知,切片是需加上西经...显然是不需要的,我们只想知道沿某个经度(或纬度)的地形变化如何,用ax.plot命令结合fill_between命令即可。而这两个命令,只需要传入一个一维的横坐标,和一维的纵坐标即可。...关键就在怎么把z二维的变为一维的。...我还没有画过,但是猜测应当是这个数据为四维数据,经度、纬度做降维处理,图上可以看出,这张图代表(30.28°E,108.93°N)这一个的整层数据随时间的变化。

12.5K75

用python在地球投影中轻松添加图形标注

circle = patches.Circle((2.5, 2.5), radius=1, linewidth=2, edgecolor='black', facecolor='none') # 圆形添加到子图中...ax.add_patch(circle) # 设置坐标范围 ax.set_xlim([0, 5]) ax.set_ylim([0, 5]) # 设置坐标纵横比例 ax.set_aspect('...1、非极地投影 以lat-lon投影为例,令纬度0°、经度0°为圆心,15°为半径,绘制圆形: code import matplotlib.pyplot as plt from cartopy import...figsize=(10, 10)) ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree()) # 设置地图的显示范围为-180到180经度...极地投影 以极射赤平投影为例,绘制75°N的纬度圈(黄色标注),gridlines有很多可选参数,大家可以官网了解一下: code import matplotlib.pyplot as plt from cartopy

48920

火星表面...

3.双经度畸变矫正方法介绍 传统经纬矫正算法因其无需外部设备标定而被广泛使用,算法基本思想是通过建立鱼眼图像坐标点与半球面模型坐标点之间的映射关系,球面点坐标换为球面经纬度坐标,并采用等距投影或正交投影原理球面上点经度映射为矫正图像水平坐标...模型拟通过横向经度和纵向经度对球面进行分割,建立鱼眼图像坐标点与球面双经度坐标的对应关系,然后纵向和横向经度值映射为平面坐标的横坐标和纵坐标,如图2 所示。 ?...球面模型侧视图, 与图 相对应, 其坐标为 分别为 点在 面和 面上的投影, 、 分别为 与 正向的夹角和 与 正向的夹角...由图3(c)可知, 竖直方向 为 与 负半的夹角, 水平方向经度 为 与 负半的夹角,因此可得: 继续进行化简,则: 此外, 根据球面特征以及角度与坐标之间关系, 可建立角度...我曾经花了4个月,跨专业双非上岸华五软工硕士,也从不会编程到进入到百度与腾讯实习。

1.7K20

【附jupyter代码】经验正交分解EOF详解及案例

由于这些主分量中前几个就可能占有原空间(即多个变量)总方差的很大部分,即通过EOF分解可以很容易地原始要素场的变化信息浓缩在前几个主分量及其对应的空间函数上,所以,EOF分解常被用于气象要素场时空变化特征规律的提取和研究...而我们的EOF分解方法,就是时空数据集 分解为空间函数 和时间函数 两部分,即 也可表示为 其中下标 表示空间,下标 表示时间。...color='k', lw=0.5) #添加海岸线 ax1.add_feature(cfeature.LAND, facecolor='white') #添加陆地 # 置数据...facecolor='white') def eof_contourf(EOFs, PCs, pers): ''' 绘制EOF填充图 ''' plt.close # 方差转换为百分数的形式...dpi=300) # 指定投影为经纬度投影,中心经纬度为180° projection = ccrs.PlateCarree(central_longitude=180) # 坐标换为日期格式

81610

如何在一张图上同时绘制云图和降水

单位转换为m mapx = ds1['x'].to_numpy() * 1000. mapy = ds1['y'].to_numpy() * 1000....通过gamma校正等方法这种线性关系转换为非线性,使较暗的区域变亮,较亮的区域保持不变。这样可以增加整个图像的对比度,使颜色更加饱和丰富 为什么修改单位km为m 图投影坐标系一般使用的是米为单位。...直接拿千米单位的影像坐标去绘制地图,会造成非常严重的坐标错位。 因此需要提前影像的坐标单位换算为与地图投影匹配的米单位,然后再传入投影变换,进行坐标转换到地图上。...gl.xlabels_top = False # 关闭顶端标签 gl.ylabels_right = False # 关闭右侧标签 gl.xformatter = LONGITUDE_FORMATTER # x设为经度格式...gl.yformatter = LATITUDE_FORMATTER # y设为纬度格式 ax.coastlines(color="w") ax.set_extent([100, 130,10,

7010

python绘图 | 气象雷达入门级讲解&多种雷达图像可视化方法

一般来说,业务中雷达的常规扫描方式是VPPI,在扫描时,雷达天线自体扫模式中最低仰角启动,并以固定仰角零度方位角(在多普勒天气雷达工作过程中,规定正北方为0°方位角,正东方为90°方位角,天线与水平面平行为...利用PyCINRAD,我们可以很方便的雷达基数据读取为xarray.Dataset #加载所需包 import math import time from itkwidgets import view...在这里,PyCinrad库同时计算出了圆锥面上每个的具体经纬度值和高度值,有了这些值,可以帮助我们更方便的在二维和三维的笛卡尔坐标下进行可视化。...gl.top_labels = False #关闭顶端标签 gl.right_labels = False #关闭右侧标签 gl.xformatter = LONGITUDE_FORMATTER #x设为经度格式...三维可视化 雷达数据并非分布在一个曲面上,在经度-维度-高度的笛卡尔坐标系下,一层的PPI数据在三维空间中呈圆锥面分布,因此可以对其进行三维的可视化 matmatplotlib三维静态可视化 fig =

19.2K89

xarray走向netCDF处理(二):数据索引

以下文章来源于MeteoAI ,作者学前班大队长 xarray专题再次开讲,错过第一部分的可以先去补个课xarray走向netCDF处理(一):数据结构及数据读取。...今天要介绍的就是xarray的索引功能,通过索引你可以对数据进行切片,整体中提取你所关注的区域、高度或者时间。...索引及可视化实战 import arrow import numpy as np import xarray as xr import cartopy.crs as ccrs import cartopy.feature...False #关闭顶端的经纬度标签 gl.ylabels_right = False #关闭右侧的经纬度标签 gl.xformatter = LONGITUDE_FORMATTER #x设为经度的格式...Spectral_r', extend='both', cbar_kwargs=cbar_kwargs,transform=ccrs.PlateCarree()) fig.show() 本文

1.7K122

PCL中SHOT1344描述子

SHOT特征的构造过程分成两大步: 建立以关键为中心的局部坐标系 1、云中提取关键p,关键提取的方法这里不作介绍; 2、以关键p为中心,搜索半径为R范围内的所有邻近Pk; 3、为了提高对噪声的鲁棒性...5、X、Z的方向确定需要通过下面公式判断: ? 根据关键建立的局部坐标系提取形状特征 以关键的局部坐标系建立球形坐标,如下: ?...2、球形坐标半径、经度和纬度三个维度划分成32个区域,半径分为内外球2份、维度分成2份、经度分成8份(图1中仅显示经度分成4份)。...3、对于落入每个区域内的计算该点法线方向nq与坐标系Z的方向余弦值,余弦值划分成11的区间,统计落入每个区域的点数。...XYZLAB的公式如下 ? 其中f函数如下: ? 2、计算关键与近邻域内的LAB距离值。 ? 然后γ量化到纹理特征直方图中。 完整的SHOT示意图 ?

89820

Python绘制垂直剖面流线图教程

核心是定义一个名为myStreamPlot的函数,它将经纬度和风速数据转换为流线图,利用三重网格插值确保准确性。 首先设置好坐标范围与刻度,以等高线形式呈现风切变率。...然后调用myStreamPlot函数,对经度、纬度、东西风向分量u、南北风向分量v以及垂直风速w_clm进行预处理,并调整单位便于理解。通过streamplot绘制出清晰易懂的全球风场流线图。...import xarray as xr import numpy as np import matplotlib.pyplot as plt import datetime as dt import cartopy.mpl.ticker...lon = lon.flatten() # 展平经度网格 lat = lat.flatten() # 展平纬度网格 # 确保风分量的方向正确...(标记的刻度是刻意修改的,以配合翻转的数组) 当然还有另一种方法就是气压转为高度,如此就不必对数组进行翻转 有更好的方法欢迎私信交流

8110

【教程】详解相机模型与坐标转换

; WGS84坐标系:与大地直角坐标系一样,只是采用纬度(M)、经度(L)和大地高(H)表示空间中任一位置。...相机坐标系:单位m,原点是光心,x和y与像素坐标系u和v平行,z为相机的光轴。光心到像素平面的距离为焦距f。相机坐标系上的和成像平面坐标系上的存在透视投影关系。...【我们这里是 => 右手坐标系+旋转坐标系本身】 旋转顺序:外旋(z->y->x)、内旋(x->y->z) 根据每次旋转是绕旋转之后的旋转,还是固定旋转,欧拉角分为内旋(intrisic roatation...世界坐标系到像素坐标系之间的转换关系可知,已知世界坐标系下的三维坐标,只要已知内外参矩阵,就可以求得像素坐标。...坐标格式:经度区纬度区以东以北,其中以东表示经度区的中心子午线的投影距离,而以北表示距离赤道的投影距离。单位为米。

12700
领券