1、一般格点数据的绘图

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import cartopy.mpl.ticker as cticker
import matplotlib.pyplot as plt
import meteva.base as meb
def paint_contour_quiver(u,v,lon,lat,cmap,title,output_dir,shici):  #这里的u和v是dataArray格式的spd = (u**2+v**2)**0.5######画图尝试2fig = plt.figure(figsize=(12,8))proj = ccrs.PlateCarree(central_longitude=113)leftlon, rightlon, lowerlat, upperlat = (106,121,15,24)img_extent = [leftlon, rightlon, lowerlat, upperlat]ax = fig.add_axes([0.1, 0.1, 0.8, 0.6],projection = proj)# ax.set_extent(img_extent, crs=ccrs.PlateCarree())ax.add_feature(cfeature.COASTLINE) ax.set_xticks(np.arange(leftlon,rightlon,1), crs=ccrs.PlateCarree())ax.set_yticks(np.arange(lowerlat,upperlat,1), crs=ccrs.PlateCarree())lon_formatter = cticker.LongitudeFormatter()lat_formatter = cticker.LatitudeFormatter()ax.xaxis.set_major_formatter(lon_formatter)ax.yaxis.set_major_formatter(lat_formatter)c = ax.contourf(lon,lat,spd,extend='both' ,cmap=cmap,transform=ccrs.PlateCarree())cbar = fig.colorbar(c, orientation='vertical', pad=0.02, aspect=20, shrink=0.6)cbar.set_label('Speed (m/s)')q = ax.quiver(lon,lat,u,v,transform=ccrs.PlateCarree(),scale=150,color='k',regrid_shape = 10,width=0.003,headwidth=4,headlength=7)ax.quiverkey(q, 0.95, 1.05, 10, '10m/s', labelpos='E',coordinates='axes')ax.set_title(title)plt.savefig(output_dir+shici+'-'+title+'.png',dpi=400)

其中u和v是二维的numpy.ndarray格式的。

2、卫星轨道资料

但是卫星轨道资料并不是二维格点的。
它在解码后是dataframe格式的。

3、插值到格点并绘图

所以，要画图，需要把卫星轨道资料插值到格点上，然后再画图。这里用meteva库的cressman插值函数进行站点到格点的插值。

def paint_obs(obs_data,type_of_obs,output_dir,shici,cmap):#####################站点数据####################################################o#usta_u = meb.sta_data(obs_data[['类型','经度','纬度','U风']],columns = ["id","lon","lat","data0"])meb.set_stadata_coords(sta_u,level = 0,time = datetime.datetime(2019,1,1,8,0),dtime = 0)#vsta_v = meb.sta_data(obs_data[['类型','经度','纬度','V风']],columns = ["id","lon","lat","data0"])meb.set_stadata_coords(sta_v,level = 0,time = datetime.datetime(2019,1,1,8,0),dtime = 0)##################################################################################################################网格106,121,15,24#################################grid0 = meb.grid([106,121,0.01],[15,24,0.01])#start = time.time()##ogrd_u = meb.interp_sg_cressman(sta_u,grid = grid0,r_list = [16,12,8,4],nearNum = 100)grd_v = meb.interp_sg_cressman(sta_v,grid = grid0,r_list = [16,12,8,4],nearNum = 100)used_time = time.time() - start############################################################################################画图#########################################################cmap,clev = meb.def_cmap_clevs(meb.cmaps.wind_speed) #设置色标#按照种类把obs分开，画一个总的，再分别画每个种类的obs#画总的obs、o-b、o-apaint_contour_quiver(grd_u.values[0,0,0,0,:,:],grd_v.values[0,0,0,0,:,:],grd_u.lon.values,grd_u.lat.values,cmap,type_of_obs+'Obs',output_dir,shici)

4、困惑

我有几个困惑：1是洋面风的色标好丑；2是轨道之外虽然没有数据，但还是有点，可不可以让轨道之外没有点？

卫星轨道资料洋面风的绘图相关推荐

MATLAB STK构建卫星轨道姿态控制仿真演示系统
MATLAB & STK构建卫星轨道姿态控制仿真演示系统杜耀珂郑科宇摘要:通过在MATLAB中建立卫星的轨道运动和姿态动力学仿真模型就可以分析在不同的控制作用下,卫星的轨道和姿态的变化情 ...
卫星轨道和两行数据TLE
最近由于Sino-2和北斗的关系,很多网友贴了表示卫星运行轨道的TLE数据.这里想对卫星轨道参数和TLE的格式做一个简单介绍.虽然实际上没有人直接读TLE数据,而都是借助软件来获得卫星轨道和位置信息, ...
卫星轨道推演计算相关知识点总结（含欧拉角、旋转矩阵、及各坐标系转化等）
来源:轨道机动算法的C++实现_shirro123的专栏-CSDN博客卫星轨道推演计算相关基础知识点总结一.卫星的运动特性二.卫星的空间坐标系 ...
卫星轨道的估计问题(Matlab)(三):标准重采样粒子滤波(SIR)对新问题的尝试
SIR滤波器关于粒子滤波的基本知识可以参加下面的博客: 粒子滤波,讲的很通俗易懂基本粒子滤波算法过程: SIR算法伪代码: Matlab代码实现求解如下所示的滤波问题: %PF的应用 clc,c ...
卫星轨道的估计问题(Matlab)(二):扩展卡尔曼滤波(EKF)对新问题的尝试
前言在前面的问题中我们已经考虑到了用微分方程来描述卫星运动轨迹的方法: r¨=rθ˙2−GMr−2θ¨=−2r−1r˙θ˙\ddot r = r\dot \theta^2-GMr^{-2}\\\dd ...
cesium学习之 Entity 画卫星轨道（一）
cesium学习之 Entity 画卫星轨道 (一) 最近公司想做个卫星轨迹展示大屏,原型图看上去确实比较逼真,然后我就惨了,开始学习cesium,满满的English文档,看着头大. B站上面有入 ...
matlab 画卫星轨道,Matlab仿真实例-卫星轨迹
卫星轨迹一．问题提出设卫星在空中运行的运动方程为: d2rkd 2 r() 22 dtrdt 2 d 2drd rdtdt dt2 其中是k重力系数(k=401408km3/s).卫星轨道采用极坐 ...
036卫星轨道及卫星在轨运动
一.卫星轨道六要素 1.升交点赤经Ω\OmegaΩ 春分点:太阳沿黄道从天赤道以南向北通过天赤道的点,称为春分点. 升交点:卫星由南向北运行时,与地球赤道面的交点. 升交点赤经为地心与春分点连线和 ...
时间复杂度o(n^0.5)_铁路轨道资料汇总（含轨道资料、论文、视频讲解）【更新时间2020.11.5】...
小编主要从以下几个方面整理轨道资料:1线路资料,2机车车辆,3轨道视频,4论文,5打磨资料,6探伤资料,7测量资料,8戴老师语言讲堂,9老魏笔记 1关于铁路线路资料,点击下面链接即可钢轨摩擦冻结接头 ...

卫星轨道资料洋面风的绘图

1、一般格点数据的绘图

2、卫星轨道资料

3、插值到格点并绘图

4、困惑

卫星轨道资料洋面风的绘图相关推荐

最新文章

热门文章