图像质量指数
<p>[TOC]</p>
<pre><code class="language-python">%matplotlib inline
%load_ext autoreload
%autoreload 2
import meteva.base as meb
import meteva.method as mem
import meteva.perspact as mps
import numpy as np
import pandas as pd
import datetime
import meteva</code></pre>
<h1>UQI</h1>
<p><font face="黑体" color=Blue size=3>universal image quality index (通用图像质量指数)</font></p>
<p><img src="https://www.showdoc.com.cn/server/api/attachment/visitFile?sign=a8ce03f3f6e27dcdf1a04d4c70bdb417&amp;file=file.png" alt="" /></p>
<p><font face="黑体" color=Blue size=3><strong>uqi(grd_ob, grd_fo)</strong></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;"><strong><font face="黑体" color=Blue size=5>grd_ob</font></strong></td>
<td style="text-align: left;">网格数据形式的观测数据,只支持包含单一平面场的网格数据</td>
</tr>
<tr>
<td style="text-align: left;"><strong><font face="黑体" color=Blue size=5>grd_fo</font></strong></td>
<td style="text-align: left;">网格数据形式的预报数据,只支持包含单一平面场的网格数据</td>
</tr>
<tr>
<td style="text-align: left;"><font face="黑体" color=blue size=5>return</font></td>
<td style="text-align: left;">字典,内容包括uqi和3个分量的取值</td>
</tr>
</tbody>
</table>
<p>参考示例</p>
<pre><code class="language-python">grid1 = meb.grid([100, 120, 0.05], [24, 40, 0.05])
path_ob = r'H:\test_data\input\mem\mode\ob\rain03\20072611.000.nc'
path_fo_03 = r'H:\test_data\input\mem\mode\ec\rain03\20072608.003.nc'
grd_ob = meb.read_griddata_from_nc(path_ob, grid=grid1, time=&quot;2020072611&quot;, dtime=0, data_name=&quot;OBS&quot;)
grd_fo_03 = meb.read_griddata_from_nc(path_fo_03, grid=grid1, time=&quot;2020072608&quot;, dtime=3, data_name=&quot;ECMWF&quot;)
result = mem.space.uqi(grd_ob,grd_fo_03)
print(result)</code></pre>
<pre><code>{'corrcoef_score': 0.4353316462756795, 'brightness_score': 0.9998681593415543, 'variability_score': 0.9745083524098113, 'uqi': 0.4241783940311719}</code></pre>