python画南极的风场

最新推荐文章于 2024-10-03 22:30:59 发布

原创

最新推荐文章于 2024-10-03 22:30:59 发布 · 1.3k 阅读

13 ·

CC 4.0 BY-SA版权

文章标签：

#python

# ----------------------------------1.这个地区长期的风向图
# 读数据
import netCDF4 as nc
import pandas as pd
import numpy as np

file = 'adaptor.mars.internal-1645705255.1744275-4126-10-f3cf0eb9-b755-45ec-9de1-92b8694846c7.nc'
dataset = nc.Dataset(file)

uwind = dataset.variables['u'][:]
vwind = dataset.variables['v'][:]

# 读取经纬度数据lat=46, lon=72
lon = dataset.variables['longitude'][:]
lat = dataset.variables['latitude'][:]
# 对风向进行平均
uwind_mean=np.mean(uwind,axis=0)
vwind_mean=np.mean(vwind,axis=0)

# 绘图可视化
import matplotlib.pyplot as plt
import matplotlib.path as mpath
import matplotlib.ticker as mticker
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
import cartopy.mpl.ticker as cticker
import cartopy.io.shapereader as shpreader
import seaborn as sns
sns.set(style='whitegrid')
font = {
   
   'family' : 'Arial',
        'weight' : 'normal',
        'size'   : 30,
        }
fig1 = plt.figure(figsize=(18,14))
######################################################################整体投影图
leftlon, rightlon, lowerlat, upperlat = (-180,180,-60,-90)
img_extent =