metdig

智能网格预报可视化诊断分析及权重推荐系统


剖面诊断分析绘图diag_crossection

<p>[TOC]</p> <h1>相当位温, 沿剖面风和风速</h1> <p>相当位温, 沿剖面风和风速产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_wsp(cross_wsp, cross_theta, cross_u, cross_v, cross_terrain, hgt,     st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),    h_pos=[0.125, 0.665, 0.25, 0.2],    wsp_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},     **pallete_kwargs)</code></pre> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">类型</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_wsp</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">风速</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度数据</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">wsp_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_3D_grid(...) psfc = get_model_3D_grid(...) wsp = mdgcal.other.wind_speed(u, v) cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_wsp = mdgcal.cross_section(wsp, st_point, ed_point) cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh) pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values) cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_time_rh_uv_theta(rh, u, v, theta, ...)</code></pre> <h1>相当位温, 锋生函数, 沿剖面风</h1> <p>相当位温, 锋生函数, 沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_fg( cross_u, cross_v, cross_theta, cross_fg, cross_terrain, hgt,     st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),     h_pos=[0.125, 0.665, 0.25, 0.2],     fg_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},     **pallete_kwargs)</code></pre> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">类型</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_fg</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color=&quot;red&quot;&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度数据</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">fg_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python"></code></pre> <h1>相当位温、垂直运动速度、沿剖面风</h1> <p>相当位温、垂直运动速度、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_w(cross_u, cross_v, cross_theta, cross_w, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], w_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_w</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">垂直运动速度</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度数据</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">w_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) w=read_w3d(...) cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_w = mdgcal.cross_section(w, st_point, ed_point) cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh) pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values) cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_theta_w(cross_u, cross_v, cross_theta, cross_w, cross_terrain, hgt, ...)</code></pre> <h1>相当位温、水平散度、沿剖面风</h1> <p>相当位温、水平散度、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_div(cross_div, cross_theta, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], div_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_div</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">水平散度</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度数据</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">div_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) div,u,v = read_div_uv_3d(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) # +form 3D psfc _, psfc_bdcst = xr.broadcast(tmp, psfc.squeeze()) psfc_bdcst = psfc_bdcst.where(psfc_bdcst &amp;gt; -10000, drop=True)  # 去除小于-10000 cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_div = mdgcal.cross_section(div, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh) pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values) cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td) cross_terrain = pressure - cross_psfc # plot draw_cross.draw_wind_theta_div(cross_div, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre> <h1>云水比、雪水比、水平风、温度时间剖面</h1> <p>云水比、雪水比、水平风、温度时间剖面产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_wind_qcld_qsn_tmp(qcld, qsn, tmp, u, v, terrain, mean_area=None, qcld_contour_kwargs={},qice_contour_kwargs={}, tmp_contourf_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">qcld</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">云水混合比</td> </tr> <tr> <td style="text-align: left;">qsn</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">雨水混合比</td> </tr> <tr> <td style="text-align: left;">tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">mean_area</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;"></td> </tr> <tr> <td style="text-align: left;">qcld_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">qice_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中brabs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data u = get_model_3D_grids(...) v = get_model_3D_grids(...) qcld = get_model_3D_grids(...) qsn = get_model_3D_grids(...) psfc = get_model_3D_grids(...) tmp = get_model_3D_grids(...) _, pressure = xr.broadcast(v, v['level']) terrain = pressure - psfc.values # plot draw_cross.draw_time_wind_qcld_qsn_tmp(qcld, qsn, tmp, u, v, terrain, ...)</code></pre> <h1>云水比、冰水比、水平风、温度时间剖面</h1> <p>云水比、冰水比、水平风、温度时间剖面产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_wind_qcld_qice_tmp(qcld, qice, tmp, u, v, terrain,mean_area=None, qcld_contour_kwargs={},qice_contour_kwargs={}, tmp_contourf_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">qcld</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">云水混合比</td> </tr> <tr> <td style="text-align: left;">qice</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">云冰混合比</td> </tr> <tr> <td style="text-align: left;">tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">mean_area</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;"></td> </tr> <tr> <td style="text-align: left;">qcld_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">qice_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中brabs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data u = get_model_3D_grids(...) v = get_model_3D_grids(...) qcld = get_model_3D_grids(...) qice = get_model_3D_grids(...) psfc = get_model_3D_grids(...) tmp = get_model_3D_grids(...) _, pressure = xr.broadcast(v, v['level']) terrain = pressure - psfc.values # plot draw_cross.draw_time_wind_qcld_qice_tmp(qcld, qice, tmp, u, v, terrain, ...)</code></pre> <h1>温度、温度平流、沿剖垂直环流</h1> <p>温度、温度平流、沿剖垂直环流产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_w_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_t, cross_w, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], tmpadv_contourf_kwargs={}, tmp_contour_kwargs={}, wind_quiver_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_tmpadv</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度平流</td> </tr> <tr> <td style="text-align: left;">cross_tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">cross_t</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">cross_w</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">风速</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">tmpadv_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">wind_quiver_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中quiver函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data u = get_model_3D_grid(...) v = get_model_3D_grid(...) w=read_w3d(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) tmpadv = mdgcal.var_advect(tmp, u, v) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_w = mdgcal.cross_section(w, st_point, ed_point) cross_t, cross_n = mdgcal.cross_section_components(cross_u, cross_v) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_tmpadv = mdgcal.cross_section(tmpadv, st_point, ed_point) pressure = mdgstda.gridstda_full_like_by_levels(cross_tmp, cross_tmp['level'].values) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_w_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_t, cross_w*ratio, cross_terrain, hgt, ...)</code></pre> <h1>温度、温度平流、沿剖面风</h1> <p>温度、温度平流、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], tmpadv_contourf_kwargs={}, tmp_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_tmpadv</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度平流</td> </tr> <tr> <td style="text-align: left;">cross_tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">tmpadv_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) tmpadv = mdgcal.var_advect(tmp, u, v) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_tmpadv = mdgcal.cross_section(tmpadv, st_point, ed_point) pressure = mdgstda.gridstda_full_like_by_levels(cross_tmp, cross_tmp['level'].values) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre> <h1>温度、垂直涡度平流、沿剖面风</h1> <p>温度、垂直涡度平流、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_vortadv_tmp(cross_vortadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], vortadv_contourf_kwargs={}, tmp_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_vortadv</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">垂直涡度平流</td> </tr> <tr> <td style="text-align: left;">cross_tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">vortadv_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) vort = mdgcal.vorticity(u, v) vortadv = mdgcal.var_advect(vort, u, v) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_vortadv = mdgcal.cross_section(vortadv, st_point, ed_point) pressure = mdgstda.gridstda_full_like_by_levels(cross_tmp, cross_tmp['level'].values) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_vortadv_tmp(cross_vortadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre> <h1>饱和相当位温, 准地转湿位涡, 沿剖面准地转风</h1> <p>饱和相当位温, 准地转湿位涡, 沿剖面准地转风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_thetaes_mpvg(cross_mpvg, cross_theta, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], mpv_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_mpvg</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">mpv_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python"></code></pre> <h1>相当位温、湿位涡、沿剖面风</h1> <p>相当位温、湿位涡、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_mpv(cross_mpv, cross_theta, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], mpv_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_mpv</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">湿位涡</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">mpv_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data theta = read_theta3d(...) mpv, _div, u, v = read_pv_div_uv(...) hgt = get_model_grid(...) psfc = get_model_grid(...) cross_theta = mdgcal.cross_section(theta, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_mpv = mdgcal.cross_section(mpv, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) pressure = mdgstda.gridstda_full_like_by_levels(cross_theta, cross_theta['level'].values) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_theta_mpv(cross_mpv, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre> <h1>相当位温、绝对涡度、沿剖面风</h1> <p>相当位温、绝对涡度、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_absv(cross_absv, cross_theta, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], absv_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_absv</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">绝对涡度</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">absv_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) absv = mdgcal.absolute_vorticity(u, v) cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_absv = mdgcal.cross_section(absv, st_point, ed_point) cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh) pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values) cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_theta_absv(cross_absv, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre> <h1>相当位温、相对湿度、沿剖面风</h1> <p>相当位温、相对湿度、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_rh(cross_rh, cross_theta, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], rh_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_rh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相对湿度</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">rh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh) pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values) cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_theta_rh(cross_rh, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre> <h1>相当位温、绝对湿度、沿剖面风</h1> <p>相当位温、绝对湿度、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_w_theta_spfh(cross_spfh, cross_theta, cross_t, cross_w, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], spfh_contourf_kwargs={}, theta_contour_kwargs={}, wind_quiver_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_spfh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">绝对湿度</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_t</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">cross_w</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">spfh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">wind_quiver_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中quiver函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) vvel = get_model_3D_grid(...) spfh = get_model_3D_grid(...) psfc = get_model_grid(...) #calcu w = mdgcal.vertical_velocity(vvel, tmp, spfh) cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_w = mdgcal.cross_section(w, st_point, ed_point) cross_t, cross_n = mdgcal.cross_section_components(cross_u, cross_v) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh) pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values) cross_spfh = mdgcal.specific_humidity_from_dewpoint(pressure, cross_td) cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_w_theta_spfh(cross_spfh, cross_theta, cross_t, cross_w*ratio, cross_terrain, hgt, ...)</code></pre> <h1>相当位温、绝对湿度、沿剖面风`</h1> <p>相当位温、绝对湿度、沿剖面风产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_theta_spfh(cross_spfh, cross_theta, cross_u, cross_v, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], spfh_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_spfh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">绝对湿度</td> </tr> <tr> <td style="text-align: left;">cross_theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">spfh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh) pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values) cross_spfh = mdgcal.specific_humidity_from_dewpoint(pressure, cross_td) cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_theta_spfh(cross_spfh, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre> <h1>温度、相对湿度、水平风场</h1> <p>温度、相对湿度、水平风场产品定制代码示例如下:</p> <pre><code class="language-python">draw_wind_tmp_rh_vvel(cross_rh, cross_tmp, cross_u, cross_v, cross_vvel, cross_terrain, hgt, st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65), h_pos=[0.125, 0.665, 0.25, 0.2], rh_contourf_kwargs={}, tmp_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">cross_rh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">水平散度</td> </tr> <tr> <td style="text-align: left;">cross_tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">cross_u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">cross_v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">cross_vvel</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> </tr> <tr> <td style="text-align: left;">cross_terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">hgt</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">高度</td> </tr> <tr> <td style="text-align: left;">st_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">起始点经纬度</td> </tr> <tr> <td style="text-align: left;">ed_point</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">结束点经纬度</td> </tr> <tr> <td style="text-align: left;">lon_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">经度剖面</td> </tr> <tr> <td style="text-align: left;">lat_cross</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">纬度剖面</td> </tr> <tr> <td style="text-align: left;">h_pos</td> <td style="text-align: left;">否</td> <td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td> <td style="text-align: left;">小图区域</td> </tr> <tr> <td style="text-align: left;">map_extent</td> <td style="text-align: left;">否</td> <td style="text-align: left;">(50, 150, 0, 65)</td> <td style="text-align: left;">绘图经纬度</td> </tr> <tr> <td style="text-align: left;">rh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data rh = get_model_3D_grid(...) u = get_model_3D_grid(...) v = get_model_3D_grid(...) tmp = get_model_3D_grid(...) hgt = get_model_grid(...) psfc = get_model_grid(...) cross_rh = mdgcal.cross_section(rh, st_point, ed_point) cross_u = mdgcal.cross_section(u, st_point, ed_point) cross_v = mdgcal.cross_section(v, st_point, ed_point) cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point) cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point) cross_u_t, cross_v_n = mdgcal.cross_section_components(cross_u, cross_v) _, pressure = xr.broadcast(cross_rh, cross_tmp['level']) cross_terrain = pressure - cross_psfc cross_terrain.attrs['var_units'] = '' # plot draw_cross.draw_wind_tmp_rh(cross_rh, cross_tmp, cross_u, cross_v, cross_u_t, cross_v_n, cross_terrain, hgt, ...)</code></pre> <h1>相当位温、相对湿度、水平风</h1> <p>相当位温、相对湿度、水平风产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_rh_uv_theta(rh, u, v, theta,terrain, rh_contourf_kwargs={}, uv_barbs_kwargs={}, theta_contour_kwargs={}, terrain_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">rh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相对湿度</td> </tr> <tr> <td style="text-align: left;">theta</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> </tr> <tr> <td style="text-align: left;">rh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">theta_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">terrain_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data tmp = get_model_3D_grids(...) u = get_model_3D_grids(...) v = get_model_3D_grids(...) rh = get_model_3D_grids(...) td = mdgcal.dewpoint_from_relative_humidity(tmp, rh) pressure = mdgstda.gridstda_full_like_by_levels(rh, rh['level'].values) theta = mdgcal.equivalent_potential_temperature(pressure, tmp, td) # plot draw_cross.draw_time_rh_uv_theta(rh, u, v, theta, ...)</code></pre> <h1>散度、垂直涡度、绝对湿度、水平风时间剖面</h1> <p>散度、垂直涡度、绝对湿度、水平风时间剖面产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_div_vort_spfh_uv(div, vort, spfh, u, v, terrain, spfh_contourf_kwargs={}, uv_barbs_kwargs={}, div_contour_kwargs={}, vort_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">div</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">水平散度</td> </tr> <tr> <td style="text-align: left;">vort</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">垂直涡度</td> </tr> <tr> <td style="text-align: left;">spfh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">绝对湿度</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">spfh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中brabs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">div_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">vort_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data div, u, v = read_div_uv_4d(...) vort, u, v = read_vort_uv_4d(...) spfh = read_spfh_4D(...) psfc = get_model_3D_grids(...) _, pressure = xr.broadcast(v, v['level']) terrain = pressure - psfc.values.repeat(pressure['level'].size, axis=1) # plot draw_cross.draw_time_div_vort_spfh_uv(div, vort, spfh, u, v, terrain, ...)</code></pre> <h1>温度、温度平流、水平风时间剖面</h1> <p>温度、温度平流、水平风时间剖面产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_wind_tmpadv_tmp(tmpadv, tmp, u, v, terrain,mean_area=None, tmpadv_contour_kwargs={'levels': np.arange(-15, 15, 1)}, tmp_contourf_kwargs={}, uv_barbs_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">tmpadv</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度平流</td> </tr> <tr> <td style="text-align: left;">tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">mean_area</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">tmpadv_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中brabs函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data u = get_model_3D_grids(...) v = get_model_3D_grids(...) psfc = get_model_3D_grids(...) tmp = get_model_3D_grids(...) tmpadv = mdgcal.var_advect(tmp, u, v) _, pressure = xr.broadcast(v, v['level']) terrain = pressure - psfc.values terrain.attrs['var_units'] = '' # plot draw_cross.draw_time_wind_tmpadv_tmp(tmpadv, tmp, u, v, terrain, ...)</code></pre> <h1>温度、垂直涡度平流、水平风时间剖面</h1> <p>温度、垂直涡度平流、水平风时间剖面产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_wind_vortadv_tmp(vortadv, tmp, u, v, terrain,mean_area=None, vortadv_contour_kwargs={}, tmp_contourf_kwargs={}, uv_barbs_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">vortadv</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">水平散度</td> </tr> <tr> <td style="text-align: left;">tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">mean_area</td> <td style="text-align: left;">否</td> <td style="text-align: left;">None</td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">vortadv_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中brabs函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data vort, u, v = read_vort_uv_4d(...) psfc = get_model_3D_grids(...) tmp = get_model_3D_grids(...) vortadv = mdgcal.var_advect(vort, u, v) _, pressure = xr.broadcast(v, v['level']) terrain = pressure - psfc.values terrain.attrs['var_units'] = '' # plot draw_cross.draw_time_wind_vortadv_tmp(vortadv, tmp, u, v, terrain, ...)</code></pre> <h1>散度、垂直涡度、相对湿度、水平风时间剖面</h1> <p>散度、垂直涡度、相对湿度、水平风时间剖面产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_div_vort_rh_uv(div, vort, rh, u, v, terrain, rh_contourf_kwargs={}, uv_barbs_kwargs={}, div_contour_kwargs={}, vort_contourf_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">div</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">水平散度</td> </tr> <tr> <td style="text-align: left;">vort</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相当位温</td> </tr> <tr> <td style="text-align: left;">rh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相对湿度</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">rh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中brabs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">div_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">vort_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data div, u, v = read_div_uv_4d(...) vort, u, v = read_vort_uv_4d(...) rh = get_model_3D_grids(...) psfc = get_model_3D_grids(...) _, pressure = xr.broadcast(v, v['level']) terrain = pressure - psfc.values.repeat(pressure['level'].size, axis=1) # plot draw_cross.draw_time_div_vort_rh_uv(div, vort, rh, u, v, terrain, ...)</code></pre> <h1>温度 相对湿度 水平风 气压垂直速度</h1> <p>温度 相对湿度 水平风 气压垂直速度产品定制代码示例如下:</p> <pre><code class="language-python">draw_time_rh_uv_tmp_vvel(rh, u, v, tmp, vvel, terrain, rh_contourf_kwargs={}, uv_barbs_kwargs={}, tmp_contour_kwargs={},vvel_contour_kwargs={}, **pallete_kwargs)</code></pre> <p>根据输入的数据,生成叠加产品并输出到指定目录。</p> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;参数说明&lt;/b&gt;&lt;/font&gt;</p> <table> <thead> <tr> <th style="text-align: left;">参数</th> <th style="text-align: left;">必选</th> <th style="text-align: left;">默认</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">rh</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">相对湿度</td> </tr> <tr> <td style="text-align: left;">u</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">U风</td> </tr> <tr> <td style="text-align: left;">v</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">V风</td> </tr> <tr> <td style="text-align: left;">vvel</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> </tr> <tr> <td style="text-align: left;">tmp</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">温度</td> </tr> <tr> <td style="text-align: left;">terrain</td> <td style="text-align: left;">&lt;font color='red'&gt;是&lt;/font&gt;</td> <td style="text-align: left;"></td> <td style="text-align: left;">地形</td> </tr> <tr> <td style="text-align: left;">rh_contourf_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contourf函数的可变参数</td> </tr> <tr> <td style="text-align: left;">tmp_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;">uv_barbs_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中barbs函数的可变参数</td> </tr> <tr> <td style="text-align: left;">vvel_contour_kwargs</td> <td style="text-align: left;">否</td> <td style="text-align: left;">{}</td> <td style="text-align: left;">matplotlib中contour函数的可变参数</td> </tr> <tr> <td style="text-align: left;"><code>**pallete_kwargs</code></td> <td style="text-align: left;"></td> <td style="text-align: left;"></td> <td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> &quot;graphics层画板设置draw_compose&quot;)</td> </tr> </tbody> </table> <table> <thead> <tr> <th style="text-align: left;">返回值</th> <th style="text-align: left;">说明</th> </tr> </thead> <tbody> <tr> <td style="text-align: left;">ret</td> <td style="text-align: left;">图像字典</td> </tr> </tbody> </table> <p>&lt;font size=&quot;4&quot;&gt;&lt;b&gt;调用示例&lt;/b&gt;&lt;/font&gt;</p> <pre><code class="language-python">from metdig.io import get_model_grid from metdig.io import get_model_3D_grid from metdig.io import get_model_3D_grids from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D from metdig.onestep.complexgrid_var.theta import read_theta3d from metdig.onestep.complexgrid_var.w import read_w3d from metdig.products import diag_crossection as draw_cross import metdig.cal as mdgcal import metdig.utl as mdgstda # get data tmp = get_model_3D_grids(...) u = get_model_3D_grids(...) v = get_model_3D_grids(...) rh = get_model_3D_grids(...) psfc = get_model_3D_grids(...) pressure = xr.broadcast(v, v['level']) terrain = pressure - psfc.values.repeat(pressure['level'].size, axis=1) terrain.attrs['var_units'] = '' # plot draw_cross.draw_time_rh_uv_tmp(rh, u, v, tmp, terrain, ...)</code></pre>

页面列表

ITEM_HTML