- 积分
 - 7833
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 2011-8-30
 
	- 最后登录
 - 1970-1-1
 
 
 
 
 
 
 | 
	
 
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册 
 
 
 
x
 
def cinrad_plot_profile(file): 
 
    f = cinrad.io.read_auto(file) 
    rl = list(f.iter_tilt(230, "REF")) 
    cr = cinrad.calc.quick_cr(rl) 
    fig = PPI(cr, dpi=500, style="white", add_city_names=True) 
    fig.geoax.scatter( 
        x=122.34, y=30.3, s=200, c="r", marker="." 
    ) 
    fig.settings["is_inline"] = False  # notebook模式下需要这一行,正常不需要 
    vcs = cinrad.calc.VCS(rl) 
    sec = vcs.get_section(start_cart=(121., 30.3), end_cart=(123., 30.3))  # 传入经纬度坐标 
    # sec = vcs.get_section(start_polar=(113, 250), end_polar=(114, 28)) # 传入极坐标 
    # 附加操作 
    fig.plot_range_rings([50, 100, 150, 200, , color="white", linewidth=1)  # 用这个来画圈 
    fig.gridlines(draw_labels=True, linewidth=1, color="white")  # 用这个来画经纬度网格线 
    fig.plot_cross_section(sec, linecolor="red") 
    fig(file_path.split('\\')[5]+"_plot_profile.png") 
白线为30度,绘制30.3剖面为何跑北面去了 
 
 
 |   
 
 
 
 |