如何在露天包裝R中添加我自己的形狀文件trajLevel（）

Question

我正在使用以下代碼繪制軌跡...........

 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

在這里，我想用自己的shape file更改基本world map 。

我的形狀文件如下。

 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？

Answer 1

為了避免在露天場所對“世界”進行硬編碼調用，您將不得不以與“地圖”包相同的基於文件的格式創建一個新的世界數據庫。

可能最簡單的方法是使用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")將繪制您的世界地圖版本。