贡献一个读swan产品（bin.bz2格式）的函数

schliezer

站里大神太多，也贡献了很多源代码，雷达方面不是很多，我就不揣冒昧贡献一个自己写的读SWAN产品CR组合反射率的小程序，需要事先安装好numpy，struct，bz2三个包：

import numpy as np
from bz2 import BZ2File as bz
import struct as sr
def ReadCR(filename):
    f=bz(filename)
    # Read parameters from file header
    f.seek(66)        #skip to date tag
    year=sr.unpack('H',f.read(2))
    month=sr.unpack('H',f.read(2))
    day=sr.unpack('H',f.read(2))
    hour=sr.unpack('H',f.read(2))
    minute=sr.unpack('H',f.read(2))
    time=year+month+day+hour+minute
    f.seek(2,1)
    XNumGrids=sr.unpack('H',f.read(2)) #grid number in X direction
    YNumGrids=sr.unpack('H',f.read(2)) #grid number in Y direction
    ZNumGrids=sr.unpack('H',f.read(2)) #grid number in Z direction
    XYZdim=XNumGrids+YNumGrids+ZNumGrids
    NumGrids=1
    for n in XYZdim:
        NumGrids*=n
    f.seek(4,1)
    StartLon=sr.unpack('f',f.read(4))
    StartLat=sr.unpack('f',f.read(4))
    CentreLon=sr.unpack('f',f.read(4))
    CentreLat=sr.unpack('f',f.read(4))
    XReso=sr.unpack('f',f.read(4))
    YReso=sr.unpack('f',f.read(4))
    EndLon=np.array(StartLon)+np.array(XReso)*(np.array(XNumGrids)-1)
    EndLat=np.array(StartLat)-np.array(YReso)*(np.array(YNumGrids)-1)
    f.seek(1024,0)
    fmt='%sB'%NumGrids
    CR=sr.unpack(fmt,f.read(NumGrids))
    CR=(np.array(CR).reshape((-1,XNumGrids[0]))-66)/2 # quality control on CR
    CR[CR==-33]=0                                     # quality control on CR
    Lat=np.linspace(np.array(StartLat),EndLat,np.array(YNumGrids))
    Lon=np.linspace(np.array(StartLon),EndLon,np.array(XNumGrids))
    return CR,Lat,Lon,XReso,YReso

ququ

楼主厉害！

你管呢

楼主想问您一下，我现在有了雷达基数据，用SWAN应该怎么拼图呢？可以加我QQ879762500

edwardli

你管呢 发表于 2021-8-3 12:00
楼主想问您一下，我现在有了雷达基数据，用SWAN应该怎么拼图呢？可以加我QQ879762500

pycinrad就可以。如果有能力pyart或者wralib也可以。

你管呢

edwardli 发表于 2021-8-3 12:53
pycinrad就可以。如果有能力pyart或者wralib也可以。

谢谢您，那应该怎么做呢？我用cinrad.easycalc.GridMapper就可以吗？pyart要怎么做呢？麻烦您了！

HUYAN

楼主请问有没有做过把某一高度比如3km的反射率提取出来的呢

schliezer

HUYAN 发表于 2022-10-8 17:23
楼主请问有没有做过把某一高度比如3km的反射率提取出来的呢

抱歉还没有，正在尝试，所以翻出来了以前自己发的贴子