python nc basemap

最新推荐文章于 2024-09-11 11:00:10 发布

啊泽不爱香菜

最新推荐文章于 2024-09-11 11:00:10 发布

阅读量2.1k

点赞数 1

分类专栏： basemap&nc 文章标签： nc-basemap

本文链接：https://blog.youkuaiyun.com/weixin_41108410/article/details/78729823

版权

本文介绍如何利用Python和netCDF4库读取NWS降水数据，并使用Basemap库进行可视化绘制。通过实例展示了从读取数据到绘制降水分布图的完整过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

官网链接

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Dec  5 23:10:46 2017

@author: chenze
"""
from mpl_toolkits.basemap import Basemap, cm
# requires netcdf4-python (netcdf4-python.googlecode.com)

from netCDF4 import Dataset as NetCDFFile

import numpy as np
import matplotlib.pyplot as plt
#import h5py

# plot rainfall from NWS using special precipitation
# colormap used by the NWS, and included in basemap.

nc = NetCDFFile('/home/chenze/practice/nws_precip_1day_20171203_netcdf/nws_precip_conus_20061222.nc', mode='r')
#nc= NetCDFFile('/home/chenze/practice/nws_precip_1day_20171203_netcdf/nws_precip_1day_20171203_ak.nc')
#nc2 = NetCDFFile('/home/chenze/practice/nws_precip_1day_20171203_netcdf/nws_precip_1day_20171203_pr.nc')
# data from http://water.weather.gov/precip/
print(nc)#,'1\n',nc1,'2\n',nc2)
print(nc.variables.keys())
prcpvar = nc.variables['amountofprecip']
data = 0.01*prcpvar[: