Файл: Мауна-Лоа среднемесячная концентрация CO2.svg

# Изменение CO2 в атмосфере Мауна-Лоа # многоязычный график для Википедии # с использованием другого возможного источника (Скриппса, если NOAA недоступен)# Обязательные пакеты. используется для: библиотеки ( tidyverse )  # библиотеки обработки данных ( gridExtra )  # библиотеки позиционирования графика ( шкалы )  # красивых разрывов в библиотеке осей графика ( lubridate )  # библиотеки обработки даты ( httr )  # проверки URL# установить языковой стандарт заранее # sudo locale-gen ca_ES.UTF-8# Переводы ------------------------------------------------ ------------язык  <-  список (  ca_ES  =  list (  locale_lc_time  =  "ca_ES.UTF-8" ,  title  =  bquote ( "Mitjana mensual de la contratió de" ~ CO [ 2 ]),  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "Данные: RF Keeling, SJ Walker, SC Piper и AF Bollenbacher \ nScripps Программа CO2 (http://scrippsco2.ucsd.edu). Visitada " ,  Sys.Date ()),  caption_noaa  =  paste ("dades: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ n д-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). Visitada" ,  Sys.Date ()),  x  =  "любой" ,  y  =  bquote ( "fracció de" ~ CO [ 2 ] ~ "en aire sec (" * mu * "моль / моль)" ),  x2  =  "mes" ,  y2  =  bquote ( atop ( "fracció de" ~ CO [ 2 ] ~"en aire sec (" * mu * "mol / mol)" ,  "desviació de la mitjana anual" )),  title2  =  "Variació mensual"  ),  cs_CZ  =  list (  locale_lc_time  =  "cs_CZ.UTF-8" ,  title  =  "Průměrné měsíční koncentraceoxidu uhličitého" ,  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "данные: RF Keeling, SJ Walker, SC Пайпер Программа nScripps CO2 (http://scrippsco2.ucsd.edu). Přístupné " ,  Sys.Date ()),  caption_noaa  =  paste ( " данные: доктор Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd / ccgg / Trends /) a \ nД-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/)., Sys.Date ()),  x  =  "rok" ,  y  =  bquote ( "koncentrace" ~ CO [ 2 ] ~ "v suchém vzduchu (" * mu * "mol / mol)" ),  x2  =  "měsíc" ,  y2  =  bquote ( наверху ( "концентратор" ~ CO [ 2 ] ~ "в такой вздучу (" * mu * "моль / моль)" ,  "одчилка од рочного прмеру" )), title2  = "Měsíční změna (прůмерна рочная одчилка)"  ),  de_DE  =  list (  locale_lc_time  =  "de_DE.UTF-8" ,  title  =  bquote ( "Monatliche durchschnittliche" ~ CO [ 2 ] * "-Konzentration" ),  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "Datei: RF Keeling , SJ Walker, SC Piper und AF Bollenbacher \ nScripps Программа CO2 (http://scrippsco2.ucsd.edu). Zugänglich am " ,  Sys.Date ()),  caption_noaa  =  paste ("Дейте: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ nД-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). Zugänglich am" ,  Sys.Date ()),  x  =  "Яр" ,  y  =  bquote ( CO [ 2 ] * "-Anteil in trockener Luft (" * mu * "моль / моль)" ),  x2  =  "Monate" ,  y2  =  bquote ( поверх ( CO [ 2 ] * "-Anteil in trockener Luft (" * mu* "моль / моль)" ,  "Abweichung vom Jahresmittel" )),  title2  =  "Monatliche Variation"  ),   en_US  =  list (  locale_lc_time  =  "en_US.UTF-8" ,  title  =  bquote ( "Среднее значение за месяц" ~ CO [ 2 ] ~ "концентрация" ),  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "Данные: RF Keeling, SJ Walker, SC Piper и AF Bollenbacher \ nScripps Программа CO2 (http://scrippsco2.ucsd.edu). Доступ " ,  Sys.Date ()),  caption_noaa  =  paste («Данные: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ n д-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). Доступно» ,  Sys.Date ()),  x  =  "Год" ,  y  =  bquote ( CO [ 2 ] ~ "доля в сухом воздухе (" * mu * "моль / моль)" ),  x2  =  "Месяц" ,  y2  =  bquote ( наверху ( CO [ 2 ] ~ "фракция в сухом воздухе (" * mu *«моль / моль)» ,  «отклонение от среднегодового значения» )),  title2  =  «Сезонное изменение»  ),  es_ES  =  list (  locale_lc_time  =  "es_ES.UTF-8" ,  title  =  bquote ( "Media mensual de la contración de" ~ CO [ 2 ]),  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "dato: RF Keeling, SJ Walker, SC Piper и AF Bollenbacher \ nScripps Программа CO2 (http://scrippsco2.ucsd.edu). Visitada " ,  Sys.Date ()),  caption_noaa  =  paste ("дата: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ nД-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). Visitada" ,  Sys.Date ()),  x  =  "Año" ,  y  =  bquote ( "Fraccion de" ~ CO [ 2 ] ~ "en aire secco (" * mu * "моль / моль)" ),  x2  =  "Mes" ,  y2  =  bquote ( наверху ( "Fraccion de" ~ CO [ 2 ] ~"en aire secco (" * mu * "mol / mol)" ,  "Desviación de la media anual" )),  title2  =  "Variación mensual"  ),   fr_FR  =  list (  locale_lc_time  =  "fr_FR.UTF-8" ,  title  =  bquote ( "Моя менсуэла концентрации де" ~ CO [ 2 ]),  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "données: RF Keeling, SJ Walker, SC Piper et AF Bollenbacher \ nScripps CO2 Program (http://scrippsco2.ucsd.edu). Accédé le " ,  Sys.Date ()),  caption_noaa  =  paste ("données: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ n д-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). Accédé le" ,  Sys.Date ()),  x  =  "année" ,  y  =  bquote ( "фракция de" ~ CO [ 2 ] ~ "dans l'air sec (" * mu * "моль / моль)" ),  x2  =  " mois " ,  y2  =  bquote ( atop ( " дробь de "~ CO [ 2 ] ~"dans l'air sec (" * mu * "mol / mol)" ,  "en écart à la moyenne annuelle" )),  title2  =  "Saisonnière вариации"  ),  nn_NO  =  list (  locale_lc_time  =  "nn_NO.UTF-8" ,  title  =  bquote ( "Gjennomsnittlig månedlig" ~ CO [ 2 ] * "-konsentrasjon" ),  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "data: RF Keeling , SJ Walker, SC Piper и AF Bollenbacher \ nScripps Программа CO2 (http://scrippsco2.ucsd.edu). Vist " ,  Sys.Date ()),  caption_noaa  =  paste ("данные: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ nД-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). Vist" ,  Sys.Date ()),  x  =  "År" ,  y  =  bquote ( CO [ 2 ] * "-andel i tørr luft (" * mu * "моль / моль)" ),  x2  =  "Måned" ,  y2  =  bquote ( поверх ( CO [ 2 ] * "-andel i tørr luft (" * mu * "моль / моль) ",  "Avvik fra årlig gjennomsnitt" )),  title2  =  "Årlig variasjon"  ),  pl_PL  =  list (  locale_lc_time  =  "pl_PL.UTF-8" ,  title  =  bquote ( "rednie miesięczne stężenie" ~ CO [ 2 ]),  subtitle  =  "Mauna Loa" ,  caption_scripps  =  paste ( "Датчанин: RF Keeling, SJ Walker, SC Piper i AF Bollenbacher \ nScripps Программа CO2 (http://scrippsco2.ucsd.edu). Dostęp na " ,  Sys.Date ()),  caption_noaa  =  paste ("Датчанин: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ nД-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). Dostęp na" ,  Sys.Date ()),  x  =  "Rok" ,  y  =  bquote ( "Frakcja" ~ CO [ 2 ] ~ "w suchym powietrzu (" * mu * "mol / mol)" ),  x2  =  "Miesiąc" ,  y2  =  bquote ( наверху ( "Frakcja" ~ CO [ 2 ] ~"w suchym powietrzu (" * mu * "mol / mol)" ,  "Odejście od średniej rocznej" )),  title2  =  "Wahania sezonowe"  ),  ru_RU  =  list (  locale_lc_time  =  "ru_RU.UTF-8" ,  title  =  bquote ( "Среднемесячная и среднегодовая оценка" ~ CO [ 2 ]),  subtitle  =  "Мауна Лоа, Гаваи" ,  caption_scripps  =  paste ( "данные: RF Keeling, SJ Walker, SC Piper и AF Bollenbacher \ nScripps CO2 Program (http://scrippsco2.ucsd.edu). Доступ на " ,  Sys.Date ()),  caption_noaa  =  paste ("данные: д-р Питер Танс, NOAA / ESRL (www.esrl.noaa.gov/gmd/ccgg/trends/) и \ nД-р Ральф Килинг, Институт океанографии Скриппса (scrippsco2.ucsd.edu/). доступ на" ,  Sys.Date ()),  x  =  "Годы" ,  y  =  bquote ( "Доля" ~ CO [ 2 ] ~ "в сухом воздухе (" * mu * "моль / моль)" ),  x2  =  "Месяцы" ,  y2  =  bquote ( atop ( "Доля" ~ CO [ 2 ] ~"в сухом воздухе (" * mu * "моль / моль)" ,  "Отправление от среднегодовой" )),  title2  =  "Сезонные отклонения от среднегодовой"  ) )# Данные ------------------------------------------------ -------------------- # Сначала попробуйте данные NOAA, если их нет, используйте данные Scrippsif ( ! http_error ( "ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_mm_mlo.txt" ))  {  # https://www.esrl.noaa.gov/gmd/ccgg/trends/ data.html  source  <-  "noaa"  co2ml  <-  read_table2 (  "ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_mm_mlo.txt" ,  col_names  =  c ( "год" ,  "месяц" ,  «десятичный» ,  «средний» ,  «сезон» ,  «ndays» ,  «stddays» ,  «unc»),  col_types  = "iidddidd" ,  na  =  c ( "-99.99" ,  "-1" ),  comment  =  "#" )  %>%  group_by ( год )  %>%  mutate ( year_mean  =  mean ( average ,  na.rm  =  TRUE ),  delta  =  average  -  year_mean ,  vdate  =  ymd ( paste0 ( "2015-" ,  month ,  "-01" ))) %>%  переименовать( co2  =  средний )}  else  {  # http://scrippsco2.ucsd.edu/data/atmospheric_co2/primary_mlo_co2_record  # используется во время завершения работы правительства США.  Источник  <-  "scripps"  co2ml  <-  read_csv (  "http://scrippsco2.ucsd.edu/assets/data /atmospheric/stations/in_situ_co2/monthly/monthly_in_situ_co2_mlo.csv " ,  col_names  =  c ( " год " ,  " месяц " ,  " xls_date " ,  " десятичный " ,  " co2 " ,  " co2_seas_adj " ,  " fit " ,  " fit_seas_adj " ,  "co2_filled " ,  " co2_filled_seas_adj "),  col_types  =  "iiiddddddd" ,  skip  =  57 ,  na  =  "-99.99" ,  comment  =  "\" " )  %>%  group_by ( год )  %>%  mutate ( year_mean  =  mean ( co2_filled ,  na.rm  =  TRUE ) ,  delta  =  co2_filled  -  год_среднее ,  vdate  =  ymd ( paste0 ( "2015-" ,  месяц , "-01" )))  %>%  переименование ( ш2  =  co2_filled ) }# Сгенерировать график для каждого языка -------------------------------------for ( l  в  именах ( язык ))  {  сообщение ( l )  текущий  <-  язык [[ l ]]  # форматируем дату в локальных именах  Sys.setlocale ( "LC_TIME" ,  current $ locale_lc_time )  # основной график  p1  <-  ggplot ( co2ml ,  aes ( decimal ,  co2 ))  +  geom_line ( color  =  "pink" )  +  geom_point ( color  =  "red" ,  size  =  0.6 )  +  stat_smooth ( span  =  0.1 )  +  scale_x_continuous ( breaks  =  pretty_breaks ())  +  scale_y_continuous ( разрывы  = pretty_breaks ())  +  labs (  x  =  текущий $ x ,  y  =  текущий $ y ,  title  =  текущий $ title ,  subtitle  =  paste ( текущий $ subtitle ,  min ( co2ml $ year ),  "-" ,  max ( co2ml $ year ) ),  caption  =  current [[ paste0 ( "caption_" , источник )]])  +  theme_bw ()  +  theme ( plot.caption  =  element_text ( size  =  7 ))  # вставка графика  p2  <-  ggplot ( co2ml ,  aes ( vdate ,  delta ))  +  geom_hline ( yintercept  =  0 )  +  stat_smooth ( span  =  0.4 ,  se  =  FALSE )  +  stat_summary ( fun.data  =  "mean_cl_boot" ,  color  =  "red " ,  size  =  0,3 )  +  scale_x_date ( разрывы  = pretty_breaks ( 4 ),  minor_breaks  =  pretty_breaks ( 12 ),  labels  =  date_format ( "% b" ))  +  labs (  x  =  текущий $ x2 ,  y  =  текущий $ y2 ,  title  =  текущий $ title2 )  +  theme_bw ()  # объединить графики и экспортировать в SVG  p1  +  annotation_custom ( grob  =  ggplotGrob ( p2 ),  xmin  =  1957 ,  xmax  =  1991 ,  ymin  =  361 ,  ymax  =  412 )  ggsave ( file  =  paste ( "co2_mauna_loa" ,  l ,  Sys.Date (),  "wp.svg" ,  sep  =  "_" ),  ширина  =  20 ,  высота  = 20 ,  units  =  "cm" ,  device  =  svg ) }