python 畫雷達回波PPI圖

秋刀鱼CCC發表於2024-08-14
Python

雷達回波的程式碼挺簡單的，關鍵是執行環境的問題，之前用python3.8，根本畫不出來圖，尤其是地圖類的，即使安裝了cartopy也不行。

後來在一番折騰之後，嘗試用了python3.9，一切都變得簡單了起來，首先弄個虛擬環境python3.9，然後按裝一系列的包：cartopy,cinrad等包

具體詳情參考：PyCINRAD保姆級教程 - Heywhale.com

中間還包括一系列的github註冊等，github註冊和git配置環境變數參考：

具體詳情參考：Git的安裝和環境變數的配置_git環境變數-CSDN部落格

#!usr/bin/env python
# -*- coding:utf-8 -*-
"""
@author: Suyue
@file: leida.py
@time: 2024/07/30
@desc:
"""
import warnings
import cinrad
import matplotlib
matplotlib.use('TkAgg')
from cinrad.visualize import PPI
warnings.filterwarnings("ignore")
import numpy as np


basePath = "F:/"
cinrad.__version__

nFiles = basePath + "/Z_RADR_I_Z9471_20230804000007_O_DOR_CD_CAP_FMT.bin.bz2"
f = cinrad.io.read_auto(nFiles)
data = f.get_data(0,230,"REF")
print(data)

f.available_product(0) #第0個仰角有哪些產品可以讀取
f.available_tilt('REF')
f.available_tilt('VEL')
vel0 = f.get_data(1,230,"VEL")
print(vel0)

f = cinrad.io.read_auto(nFiles)
print(type(f).__name__)
data = f.get_data(0, 230, "REF")
f.available_product(0)
data["REF"].values = np.ma.masked_less(data["REF"].values, 0)
fig = cinrad.visualize.PPI(data, style="black")
print(123)
fig('F:/abc.png')

相關文章