组合反射率的数据如何变成格点数据 #95

liyuan3970 · 2023-08-29T13:33:11Z

您好，我用您的算法出来的组合反射率的数据自己用contour画图，回波地理位置存在偏移的情况，
path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/"
f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2')
tilt_number = 0
data_radius = 230
data_dtype = 'REF' # stands for reflectivity

数据加载核心

ra = f.get_data(tilt_number, data_radius, data_dtype)
rl = list(f.iter_tilt(230, 'REF'))
cr = cinrad.calc.quick_cr(rl)
data = cr.CR

颜色 ------

colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5']
levels1 = [20,25,30,35,40,45,50,55,60,65,70,75]
cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11)
contour = plt.contourf(data.longitude, data.latitude, data,cmap=cmaps,levels = levels1)

后来我验证了用单仰角绘制最低仰角的反射率是不存在位置偏移的情况的
path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/"
f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2')
v = f.get_data(0,230, 'REF')
fg =cinrad.calc.GridMapper([v])
output = fg(0.01)
colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5']
levels1 = [20,25,30,35,40,45,50,55,60,65,70,75]
cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11)
contour = plt.contourf(output.longitude, output.latitude, output.REF,cmap=cmaps,levels = levels1)

请问您这个包有组合反射率数据变成格点数据的方法么？或者调用您底层函数实现这一点您有什么建议?

CyanideCN · 2023-08-30T07:22:09Z

组合反射率有一点偏移很正常啊

…

On Tue, Aug 29, 2023 at 21:33 小飞侠 ***@***.***> wrote: 您好，我用您的算法出来的组合反射率的数据自己用contour画图，回波地理位置存在偏移的情况， path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/" f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2') tilt_number = 0 data_radius = 230 data_dtype = 'REF' # stands for reflectivity 数据加载核心 ra = f.get_data(tilt_number, data_radius, data_dtype) rl = list(f.iter_tilt(230, 'REF')) cr = cinrad.calc.quick_cr(rl) data = cr.CR 颜色 ------ colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5'] levels1 = [20,25,30,35,40,45,50,55,60,65,70,75] cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11) contour = plt.contourf(data.longitude, data.latitude, data,cmap=cmaps,levels = levels1) 后来我验证了用单仰角绘制最低仰角的反射率是不存在位置偏移的情况的 path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/" f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2') v = f.get_data(0,230, 'REF') fg =cinrad.calc.GridMapper([v]) output = fg(0.01) colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5'] levels1 = [20,25,30,35,40,45,50,55,60,65,70,75] cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11) contour = plt.contourf(output.longitude, output.latitude, output.REF,cmap=cmaps,levels = levels1) 请问您这个包有组合反射率数据变成格点数据的方法么？或者调用您底层函数实现这一点您有什么建议? — Reply to this email directly, view it on GitHub <#95>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AJWF65D7NH3MB4DZL5SN4UTXXXVSLANCNFSM6AAAAAA4C7AI2M> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

hehuan710 · 2023-09-28T05:46:00Z

组合反射率有一点偏移很正常啊
…
On Tue, Aug 29, 2023 at 21:33 小飞侠 @.> wrote: 您好，我用您的算法出来的组合反射率的数据自己用contour画图，回波地理位置存在偏移的情况， path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/" f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2') tilt_number = 0 data_radius = 230 data_dtype = 'REF' # stands for reflectivity 数据加载核心 ra = f.get_data(tilt_number, data_radius, data_dtype) rl = list(f.iter_tilt(230, 'REF')) cr = cinrad.calc.quick_cr(rl) data = cr.CR 颜色 ------ colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5'] levels1 = [20,25,30,35,40,45,50,55,60,65,70,75] cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11) contour = plt.contourf(data.longitude, data.latitude, data,cmap=cmaps,levels = levels1) 后来我验证了用单仰角绘制最低仰角的反射率是不存在位置偏移的情况的 path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/" f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2') v = f.get_data(0,230, 'REF') fg =cinrad.calc.GridMapper([v]) output = fg(0.01) colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5'] levels1 = [20,25,30,35,40,45,50,55,60,65,70,75] cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11) contour = plt.contourf(output.longitude, output.latitude, output.REF,cmap=cmaps,levels = levels1) 请问您这个包有组合反射率数据变成格点数据的方法么？或者调用您底层函数实现这一点您有什么建议? — Reply to this email directly, view it on GitHub <#95>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJWF65D7NH3MB4DZL5SN4UTXXXVSLANCNFSM6AAAAAA4C7AI2M . You are receiving this because you are subscribed to this thread.Message ID: @.>

从地图边界来看，组合反射率确实较0.5°有明显偏差，PS：这样比较是因为0.5°以上仰角基本没回波了，而且只用0.5°一个仰角的数据计算出的组合反射率没有偏移

hehuan710 · 2023-09-29T02:07:01Z

组合反射率有一点偏移很正常啊
…
On Tue, Aug 29, 2023 at 21:33 小飞侠 @.> wrote: 您好，我用您的算法出来的组合反射率的数据自己用contour画图，回波地理位置存在偏移的情况， path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/" f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2') tilt_number = 0 data_radius = 230 data_dtype = 'REF' # stands for reflectivity 数据加载核心 ra = f.get_data(tilt_number, data_radius, data_dtype) rl = list(f.iter_tilt(230, 'REF')) cr = cinrad.calc.quick_cr(rl) data = cr.CR 颜色 ------ colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5'] levels1 = [20,25,30,35,40,45,50,55,60,65,70,75] cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11) contour = plt.contourf(data.longitude, data.latitude, data,cmap=cmaps,levels = levels1) 后来我验证了用单仰角绘制最低仰角的反射率是不存在位置偏移的情况的 path ="/home/workspace/Data/My_Git/web_met/static/data/downfile/" f = cinrad.io.StandardData(path+'Z_RADR_I_Z9576_20230822032019_O_DOR_SAD_CAP_FMT.bin.bz2') v = f.get_data(0,230, 'REF') fg =cinrad.calc.GridMapper([v]) output = fg(0.01) colorslist = ['#00C800','#019000','#FFFF00','#E7C000','#FF9000','#D60000','#C00000','#FF00F0','#780084','#AD90F0','#AE0AF5'] levels1 = [20,25,30,35,40,45,50,55,60,65,70,75] cmaps = LinearSegmentedColormap.from_list('mylist',colorslist,N=11) contour = plt.contourf(output.longitude, output.latitude, output.REF,cmap=cmaps,levels = levels1) 请问您这个包有组合反射率数据变成格点数据的方法么？或者调用您底层函数实现这一点您有什么建议? — Reply to this email directly, view it on GitHub <#95>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJWF65D7NH3MB4DZL5SN4UTXXXVSLANCNFSM6AAAAAA4C7AI2M . You are receiving this because you are subscribed to this thread.Message ID: @.>

原因可能在此处，因为每个仰角的经纬度范围有差异，导致CR偏移

liyuan3970 · 2023-09-29T11:57:57Z

我回头看看，这个0.5仰角的位置偏移不明显，如果这里有问题，那经纬度按照0.5仰角的来应该来解决

CyanideCN · 2023-10-08T03:03:11Z

可以看看这个commit有没有修复，刚刚改了一下

CyanideCN added a commit that referenced this issue Oct 8, 2023

Possibly fix CR coords problem, #95

abc7d58

CyanideCN added the component:calc label Oct 8, 2023

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

组合反射率的数据如何变成格点数据 #95

组合反射率的数据如何变成格点数据 #95

liyuan3970 commented Aug 29, 2023

CyanideCN commented Aug 30, 2023 via email

hehuan710 commented Sep 28, 2023 •

edited

Loading

hehuan710 commented Sep 29, 2023

liyuan3970 commented Sep 29, 2023

CyanideCN commented Oct 8, 2023

组合反射率的数据如何变成格点数据 #95

组合反射率的数据如何变成格点数据 #95

Comments

liyuan3970 commented Aug 29, 2023

数据加载核心

颜色 ------

CyanideCN commented Aug 30, 2023 via email

hehuan710 commented Sep 28, 2023 • edited Loading

hehuan710 commented Sep 29, 2023

liyuan3970 commented Sep 29, 2023

CyanideCN commented Oct 8, 2023

hehuan710 commented Sep 28, 2023 •

edited

Loading