missInfo
package is used to detect daily meteorological station data
missing information and interpolate missing values.
Tasklist:
You can install the released version of missInfo from github with:
devtools::install_github("kongdd/missInfo")
This is a basic example which shows you how to solve a common problem:
library(missInfo)
## basic example code
data("prcp")
data("st840")
r <- interp_main(prcp, st840, smax = 200, verbose = FALSE)
#> [0 ] Detect missing information of original data ...
#> [1.1 ] runing na.approx interpolate ...
#> [1.2 ] Detect missing info after na.approx interpolate ...
#> [2.1.1] runing linear lm interpolate ...
#> [2.1.2] Detect missing info after first lm interpolate ...
#> [2.2.1] runing second linear lm interpolate ...
#> [2.2.2] Detect missing info after second lm interpolate ...
#> [3.1 ] History average interpolate...
#> [3.2 ] Detect missing info after hisavg interpolate ...
台站海拔高度 +100000 当台站海拔高度为估测值时,在估测数据基础上加100000
各要素项 32766 数据缺测或无观测任务
气压日极值 +20000 气压极值取自定时值,在原值上加20000
日最小相对湿度 +300 最小相对湿度取自定时值,在原值上加300
风速 +1000 当风速超过仪器测量上限时,在上限数据基础上加1000
风向 1-17 用数字表示风向方位,17表示静风 +100 当表示风向为八风向时,在原值上加100 90X 风向出现X个时,风向数据用个数X表示 95X 风向至少出现X个时,风向数据用个数X表示
降水量
| flag | Description | | ----- | -------------------------- | | 32700 | 表示降水“微量” (0.1mm) | | 32XXX | XXX为纯雾露霜 | | 31XXX | XXX为雨和雪的总量 | | 30XXX | XXX为雪量(仅包括雨夹雪,雪暴) |
蒸发量 32700 表示蒸发器结冰 +1000 蒸发器中注入的水全部蒸发,在注入的水量数据基础上加1000
0cm地温 +10000 实际温度(零上)超仪器上限刻度,在上限数据基础上加10000 -10000 实际温度(零下)超仪器下限刻度,在下限数据基础上减10000
Quality Control (QC) flags
| flag | Description | | ---- | ----------- | | 0 | 正确 | | 1 | 可疑 | | 2 | 错误 | | 4 | 定制数据 | | 7 | 无观测任务 | | 8 | 缺测 | | 9 | 未进行质量控制 |
[1] Dongdong Kong, R package: missing information detect and interpolate,
missInfo
version 0.1.1, https://github.com/kongdd/missInfo
Keep in mind that this repository is released under a GPL-3 license, which permits commercial use but requires that the source code (of derivatives) is always open even if hosted as a web service.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.