Подложка нанесена растровым ASCII с границами страны - R - PullRequest
Новая
       27

Подложка нанесена растровым ASCII с границами страны - R

0 голосов
/ 18 июня 2019

Я хотел бы заложить карту под построенным ASCII-файлом (в виде матрицы) или хотя бы включить границы страны в график. Набор данных был получен из GAEZ (http://www.gaez.iiasa.ac.at/),, но для получения доступа к данным необходимо зарегистрироваться.

Это ASCII-файл, который я извлек: Ссылка

Я использовал пакет SDMTools для чтения файла и заполнил .contour для чтения и построения набора данных Финляндии. 


library(SDMTools)
library(viridis)

read.asc("...\\grid_temperature_2080\\data.asc", gz = FALSE)
x <- read.asc("...\\grid_temperature_2080\\data.asc", gz = FALSE)

filled.contour(x, col = viridis(24, option = "D"), 
              levels = seq(100, 220, 5),
              plot.title = {par(cex.main=1);title(main = "Estimated Thermal Regime 2080 (CCCma CGCM2 B2)",
                                                   xlab = "Longitude", ylab = "Latitude")},
              key.title = {par(cex.main=0.8);title(main="days > 5oC")},
              plot.axes = {contour(x, levels = seq(100,220,10), add = TRUE, lwd=2)}
)

Это дает мне эту карту (извините, у меня недостаточно репутации для публикации изображений
https://i.ibb.co/7vhDWNt/Rplot.jpg

И я бы хотел, чтобы фон был похож на этот

https://i.ibb.co/dB2DXzQ/map.png

Я пытался добавить map


library(maps)

map(database = "finland", regions = ".", exact = FALSE, boundary = TRUE,
    interior = TRUE, projection = "", parameters = NULL, orientation = NULL,
    fill = FALSE, col = 1, plot = TRUE, add = TRUE, namesonly = FALSE,
    xlim = NULL, ylim = NULL, wrap = FALSE, resolution = if (plot) 1 else 0,
    type = "l", bg = par("bg"), mar = c(4.1, 4.1, par("mar")[3], 0.1),
    myborder = 0.01, namefield="name", lforce="n")

И получил «Ошибка в as_mapper (.f, ...): отсутствует аргумент« .f », по умолчанию нет»

Я попробовал некоторые варианты, но больше ничего не получил.

Спасибо за ваши усилия!

...