2016-03-08 25 views
1

這應該是一個容易的。Lambert conformal conic projection不能在R中使用mapproj?

library(maps) 
map("world", projection = "mercator", wrap = TRUE) 
map("world", projection = "cylequalarea", param = 45, wrap = TRUE) 
map("world", projection = "albers", parameters = c(lat0 = 0, lat1 = 0), wrap = TRUE) 
map("world", projection = "lambert", parameters = c(lat0 = 0, lat1 = 0), wrap = TRUE) 

前三個投影有效,但最後一個沒有 - 它只是顯示墨卡託投影。我究竟做錯了什麼?

我知道我可以使用spTransform來實際進行投影,但我希望能做些簡單的事情。

+0

什麼是你期待LCC投影看起來像在lat0 = 0,LAT1 = 0? – MikeJewski

+0

@MikeJewski與墨卡託不同嗎?我試圖做一個這樣的比較:(https://multimedia.journalism.berkeley.edu/media/upload/tutorials/qgis-basics/albers.jpg) – beanumber

+1

我不是非常精通GIS的東西,但就我所知,在lat0 = 0和lat1 = 0處進行LCC投影對於實際圖像沒有任何作用。這更多是你在找什麼? [Image](http://imgur.com/ADXcUS0) – MikeJewski

回答

0

你只需要調整您的LCC標準緯線讓它看起來錐形

library(maps) 
map("usa", projection = "mercator", wrap = TRUE) 
map("usa", projection = "cylequalarea", param = 45, wrap = TRUE) 
map("usa", projection = "albers", parameters = c(lat0 = 0, lat1 = 0), wrap = TRUE) 
map("usa", projection = "lambert", parameters = c(lat0 = 20, lat1 = 50), wrap = TRUE)