Как добавить свой собственный файл формы trajLevel () в openair пакет R - PullRequest
Новая
       46

Как добавить свой собственный файл формы trajLevel () в openair пакет R

0 голосов
/ 12 декабря 2018

Я использую этот следующий код для построения траекторий ........... 

library(openair)

load("GDASNDL1000m.Rdata")

trajLevel(traj,method="hexbin",col="jet",xbin=40,parameters=NULL,
orientation=c(90,0,0),projection="mercator")

результат https://github.com/adeckmyn/maps/files/2667752/GDASNDL1000m.zip

Здесь я хотел бы изменить базу world map на свою собственную shape file.

мой файл формы выглядит следующим образом ... 

z1=maptools::readShapePoly("/home/sateeshm/shapefiles/ncmrwf/india_map")
library(maps)
map(z1)

https://github.com/adeckmyn/maps/files/2667336/World-India.zip

#

Теперь актуальный вопрос как связать z1 с trajLevel?

1 Ответ

0 голосов
/ 12 декабря 2018

Чтобы избежать жестко закодированного вызова «world» в openair, вам нужно будет создать новую базу данных world в том же файловом формате, что и пакет «maps».

Вероятно, самый простойспособ сделать это - использовать пакет mapMaker .Этот пакет не на CRAN, но его можно найти на github.Это пакет, который я использовал для создания стандартной карты мира.Документация минимальна, но если вам не нужны названия полигонов и т. Д., Вы можете создать «быструю и грязную» карту мира следующим образом: 

# get your new map as a simple list of polygons (or lines)
z1=maps::map(maptools::readShapePoly("india_map"), plot=FALSE)
# create internal representation
z2=mapMaker::map.make(z1)
# write binary files:
mapMaker::map.export.bin(z2, "/my/path/to/world")

# To make map() call this new database:
library(maps)
worldMapEnv="MYMAP"
Sys.setenv("MYMAP"="/my/path/to/") # don't add the "world" !

, теперь map("world") нарисует вашу версиюкарта мира.

...