立即注册 登录
气象家园 返回首页

忽忽儿的个人空间 http://bbs.06climate.com/?2646 [收藏] [复制] [分享] [RSS]

留言板

facelist doodle涂鸦板

您需要登录后才可以留言 登录 | 立即注册


1048651077 2023-11-9 17:30
忽忽儿: 求助水汽后向追踪是怎么做的呢,可以加个好友咨询一下吗
https://zhuanlan.zhihu.com/p/34885340
qq469015280 2023-3-4 16:35
忽忽儿: 想问一下,后来OHP文件读取成功了吗
但是我不懂读出来ohp是什么单位? 数据很小,会不会是降水等级?
qq469015280 2023-3-4 16:33
忽忽儿: 想问一下,后来OHP文件读取成功了吗
#加载所需包 程序来自python气象数据处理与可视化202111209.pdf
####关键1:将文件名作为参数,传入metpy.io. Level3File(),
#生成io.nexrad.Level13File类型的对象。这个对象里包含站点信息、
#时间信息等基础信息。含有扫描探测数据的存储在f.sym block列表中,
#内有起始方位角'start_az'、结束方位角'end_az'、数据’data'、中心'center'、扫描宽度' gate_scale'、起始度数’first’。
import math
import time
# from itkwidgets import view
import matplotlib.pyplot as plt
import matplotlib.ticker as mticker
import numpy as np
from matplotlib import colors
from scipy.interpolate import griddata
import scipy.misc
import cartopy.crs as ccrs
import cartopy.feature as cfeat
import cinrad
import xarray as xr
from cartopy.mpl.gridliner import LATITUDE_FORMATTER, LONGITUDE_FORMATTER
from cinrad.io import CinradReader, StandardData
from cinrad.visualize import PPI, Section
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm
from mpl_toolkits.mplot3d import Axes3D
import pandas as pd
import plotly.express as px
from metpy.io import Level3File
from metpy.plots import add_metpy_logo,colortables
%matplotlib inline
##数据和绘图部分
fig,axes = plt.subplots(1,1,figsize=(15,8))
add_metpy_logo(fig,790,35,size='small')
f = Level3File("G:/雷达产品/2021/ohp2/Z9092/OHP/78/Z_RADR_I_Z9092_20210707153800_P_DOR_SA_OHP_20_230_NUL.092.bin")#### Level3File 对象里包含站点信息,时间信息等基础信息
#file=("D:/Python/ohp/20190420.030312.00.78.92")
#可以读取升级前的产品数据
ctable=('rainbow')

#print(f.sym_block)
datadict=f.sym_block[0][0]
data=np.ma.array(datadict['data'])#####需要将数据转换成数组###########
print(type(data))
data[data==0]=np.ma.masked
print(type(data))
print("数组形状:",data.shape)
print("查看数据信息:",f) ####其中'rainfall_end': datetime.datetime(2021, 12, 10, 2, 0)}为降水数据累计终止时间
print("查看扫描数据:",f.sym_block)
VV09 2021-1-25 11:29
您好,如何将ArcGI汉化?求10.2汉化包
croton 2019-1-29 09:38
您好,想咨询您一个问题。我已经把arcdesktop设为中文了,但是我的arcmap还是英文的,想问一下这个问题怎么解决?谢谢
深度狙杀 2017-7-20 14:22
忽忽儿: 我不是云南的呢,有个好朋友是云南的,你是云南哪里的呀?
大理呢
你朋友玉溪呢给事
aa22556688 2017-7-15 14:38
忽忽儿: 你好,能问个问题吗,那grads站点插值成格点后如何设置缺测值的数值呢?可以不按默认的走吗
前提要是知道 缺测值是多少,然后在ctl里面修改。
深度狙杀 2017-6-13 10:31
忽忽儿: 嘿嘿我当时也被那个问题困扰了很久,现在自己用.bat写gs,批量处理还挺快的
一看就是个大神
深度狙杀 2017-6-8 13:47
你好你好
哈哈
深度狙杀 2017-6-8 13:38
你也是云南呢噶
芊纤随风 2016-11-8 08:49
忽忽儿: 是呀~
我和她是老乡,都是玉溪人      ;也是她的学弟,比她小一届
芊纤随风 2016-10-21 08:27
你是周英姐的同学么?

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

返回顶部