MeteoInfoLab脚本示例：模式检验

呼小喵喵

太强大了，要好好学习Meteoinfo啊~~发现编程还是困难户

静言_GRMC

好强大，刚看到就第一时间测试了
ddate,OBS_AQI,OBS_AQIL,F_AQI,F_AQIL
11,88,2,95,2
12,104,3,103,3
13,153,4,132,3
14,99,2,180,4
15,70,2,68,2
16,48,1,48,1
17,63,2,55,2
18,62,2,75,2
19,78,2,70,2
110,82,2,60,2
111,53,2,49,1
112,48,1,46,1
113,54,2,57,2
114,64,2,55,2
*******365天的数据

按照王老师的方法，运行结果如下：
Continuous verification...
Score        Value
MAE        15.01
R        0.72
MeanErr        -1.80
Bias        0.97
RMSE        21.81
MSE        475.83

Dichotomous verification...
Hit: 366
Miss: 0
False alarm: 0
Correct negative: 0
---------------------------
Score        Value
Accuracy        1.00
HK        NaN
OR        NaN
POD        1.00
ETS        NaN
Bias        1.00
FAR        0.00
ORSS        NaN
POFD        NaN
HSS        NaN
SR        1.00
TS        1.00



这部分结果貌似不对，请问王老师是输入文件中缺少什么吗？

MeteoInfo

静言_GRMC 发表于 2017-5-21 16:01
好强大，刚看到就第一时间测试了
ddate,OBS_AQI,OBS_AQIL,F_AQI,F_AQIL
11,88,2,95,2

有无检验需要设置阈值。

静言_GRMC

MeteoInfo 发表于 2017-5-21 22:00
有无检验需要设置阈值。

嗯嗯   懂了   谢谢

wuyu

好强大，学习了

liuxiunuli

学习了！！！

小小气象人@ouc

楼主你好，请问wrfda里有没有自带画这种图的脚本啊

ymy_00000

程序：
fn='H:/ymy/Precipitaion6h/2017/6/SURF_CHN_PRE_H6-2017060100.txt'
#Read station data
table=readtable(fn,delimiter=' ',format='%5i%5i%6f%f%f%f%i%i%i%i%i%i%i%f%i')
lon=table['Lon']
lat=table['Lat']
rain=table['PRE_6h']
#griddata function - interpolate
x = arange(75.09375, 135, 0.28125)
y = arange(10.125, 55.125, 0.28125)
gdata,gx,gy = griddata((lon, lat), rain, xi=(x, y), method='idw', radius=0.5)
#save data to file
gdata.savegrid(x,y,'H:/ymy/Precipitaion6h/IDW/2017/06/2017060100.txt')

运行结果：
>>> run script...
Traceback (most recent call last):
  File "E:\ymy\MetoeInfo\program\IDW.py", line 5, in <module>
    table=readtable(fn,delimiter=' ',format='%5i%5i%6f%f%f%f%i%i%i%i%i%i%i%f%i')
  File "E:\ymy\MetoeInfo\MeteoInfo_1.6.2\MeteoInfo\pylib\mipylib\dataset\midata.py", line 470, in readtable
    tdata = TableUtil.readASCIIFile(filename, delimiter, headerlines, format, encoding, readvarnames)
        at java.lang.NumberFormatException.forInputString(Unknown Source)

        at java.lang.Integer.parseInt(Unknown Source)

        at java.lang.Integer.valueOf(Unknown Source)

        at org.meteoinfo.table.DataColumn.convertTo(DataColumn.java:273)

        at org.meteoinfo.table.DataRow.setValue(DataRow.java:136)

        at org.meteoinfo.table.DataRow.setValue(DataRow.java:111)

        at org.meteoinfo.table.DataTable.setValue(DataTable.java:245)

        at org.meteoinfo.data.TableUtil.readASCIIFile(TableUtil.java:168)

        at sun.reflect.GeneratedMethodAccessor18.invoke(Unknown Source)

        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)

        at java.lang.reflect.Method.invoke(Unknown Source)

java.lang.NumberFormatException: java.lang.NumberFormatException: For input string: "52.9744"

原始数据：
Station_Id_C Station_Id_d Lat Lon Alti PRS_Sensor_Alti WIN_S_Sensor_Heigh Station_Type Station_levl Year Mon Day Hour PRE_6h  Q_PRE_6h
50136 50136 52.9744 122.5108 438.5 439.7 999999 4 11 2017 6 1 0 0 0
50137 50137 53.4692 122.3783 296 297.3 999999 4 13 2017 6 1 0 0 0
50246 50246 52.3481 124.7183 361.9 363 999999 4 12 2017 6 1 0 0 0。。。。

请问王老师，我这个错误是因为readtable的format编写不正确导致的吗？如果是的话，请问该如何修改。

MeteoInfo

ymy_00000 发表于 2018-12-23 16:36
程序：
fn='H:/ymy/Precipitaion6h/2017/6/SURF_CHN_PRE_H6-2017060100.txt'
#Read station data

format='%2s%5f%6i%2f'   试试， s指字符类型，i是整形，f是浮点型，如果几个同样类型的在一起可以用数字简写，比如%5f表示5个浮点类型的字段。

ymy_00000

MeteoInfo 发表于 2018-12-24 09:14
format='%2s%5f%6i%2f'   试试， s指字符类型，i是整形，f是浮点型，如果几个同样类型的在一起可以用数字 ...

谢谢王老师，程序可以运行了