其他
GSOD全球逐日气象站点数据介绍与下载
前面给大家分享过用R语言获取NOAA ISD气象数据,那个是逐小时的,想要做逐日或者逐月的研究需要的合成比较复杂。在这里给大家推荐一套逐日的气象数据。
ISD气象数据获取可以看以前的推文:
Global Surface Summary of the Day - GSOD数据,本套数据衍生自The Integrated Surface Hourly (ISH)数据集,从1929年至今,包括全球超9000个站点的气象观测数据。
GSOD官网:https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C00516#Documentation
GSOD包含的主要内容
从这套数据的内容来看,和ISD数据差异还是比较大的,这个字段比较少,字段单位比较奇葩,是英制单位,不是常见的国际公制单位,使用时需要进行单位换算。
字段说明
STATION 气象站代码,WMO/DATSAV3合并w/WBAN码 DATE 日期,mm/dd/yyyy格式 LATITUDE 小数点格式的纬度,南半球负值 LONGITUDE 小数点格式的经度,西半球负值 ELEVATION 高程,单位米 NAME, 气象站名称 TEMP, Mean temperature平均温度,单位Fahrenheit华氏度,精确到一位小数,缺失值9999.9 TEMP_ATTRIBUTES 参与平均温度计算的观测值个数 DEWP,Mean dew point平均露点,单位Fahrenheit华氏度,精确到一位小数,缺失值9999.9 DEWP_ATTRIBUTES 参与平均露点计算的观测值个数 SLP, Mean sea level pressure平均海平面压强,单位mb毫巴mbar,1mb=1hPa百帕斯卡,缺失值9999.9 SLP_ATTRIBUTES,参与平均海平面压强计算的观测值个数 STP,Mean station pressure平均测站压强,单位mb 毫巴mbar,缺失值9999.9 STP_ATTRIBUTES,参与平均测站压强计算的观测值个数 VISIB,Mean visibility平均能见度,单位miles英里,缺失值999.9 VISIB_ATTRIBUTES,参与平均能见度计算的观测值个数 WDSP,Mean wind speed平均风速,单位knots节(海里/小时),缺失值999.9 WDSP_ATTRIBUTES,参与平均风速计算的观测值个数 MXSPD,Maximum sustained wind speed最大持续风速,单位knots节,缺失值999 GUST,Maximum wind gust最大阵风,单位knots节,缺失值999.9 MAX,Maximum temperature最高温度,单位Fahrenheit华氏度,缺失值9999.9 MAX_ATTRIBUTES 空值,自动观测到的最高温度 *
星号,从逐小时数据中提取得到的最高温度MIN, Minimum temperature最低温度,单位Fahrenheit华氏度 MIN_ATTRIBUTES 空值,自动观测到的最低温度 *
星号,从逐小时数据中提取得到的最低温度PRCP,Precipitation amount降水量(雨或融化的雪),单位inches英寸,缺失值99.99,“0”值表示当日没有测量到降水 PRCP_ATTRIBUTES A,1 report of 6-hour precipitation amount. B,Summation of 2 reports of 6-hour precipitation amount. C,3次报告6小时降水量的和 D,4次报告6小时降水量的和 E,1次报告12小时降水量 F,2次报告12小时降水量的和 G,1次报告24小时降水量 H,气象站报告当日降水量为0,但是在逐小时观测中报告了至少一次降水,可能是当日数据不完整 I,气象站当日没有报告任何降水,逐日观测中也没有报告任何降水,但是仍然可能存在降水只是没有被报告 SNDP,Snow depth积雪深度,单位inches英寸,取当日最后一次观测值(如果观测了多次),缺失值999.9 FRSHTT,指示器,当日报告的天气情况,是否有雾Fog ('F' - 1st digit)、雨Rain or Drizzle ('R' - 2nd digit)、雪Snow or Ice Pellets ('S' - 3rd digit)、冰雹Hail ('H' - 4th digit)、雷Thunder ('T' - 5th digit)、龙卷风或漏斗云Tornado or Funnel Cloud ('T' - 6th digit)。 1,有 0,没有
单位换算
压强单位:1hPa = 1mb 温度单位: ,常见温度对照: 距离单位: 英里转公里 英寸转毫米 速度单位: 1节=1海里/小时=1.852km/h
数据下载
全球数据下载地址:https://www.ncei.noaa.gov/data/global-summary-of-the-day/access/
按条件筛选数据下载:https://www.ncei.noaa.gov/access/search/data-search/global-summary-of-the-day
R语言下载GSOD数据
既然有处理ISD数据的R包,我猜肯定也有处理GSOD数据的R包,还真让我找到了,目前CRAN上面公开的就是这个GSODR。
加载GSODR包,包中内置isd_history.rda数据,内含气象站基本情况。
library(GSODR)
library(tidyverse)
load(system.file("extdata", "isd_history.rda", package = "GSODR"))
as.tibble(isd_history)
关于GSODR更多内容可以查看GSODR网站:
https://ropensci.github.io/GSODR/articles/GSODR.html
或者等下期内容我再介绍如何使用GSODR和tidyverse下载和整理气象数据。
多多点赞和再看,或者加入培训班会员据说可以让下期推文更快出来。