ホーム  •  ニュース  •  フォーラム  •  アカウント情報  •  サイト内検索  •  新規登録
 ログイン
ユーザ名:

パスワード:


パスワード紛失

新規登録
 カウントダウンブロック
カウントダウンイベントはありません
 メニュー
メイン
   地球科学のための R
     maps & mapproj パッケージ		投稿するにはまず登録を

スレッド表示 | 新しいものから 前のトピック | 次のトピック | 下へ
投稿者 スレッド
森
投稿日時: 2008-1-22 16:34
登録日: 2004-7-29
居住地: 地球
投稿: 303
maps & mapproj パッケージ
R の地図を描くためのパッケージ maps と、地図投影のためのパッケージ mapproj を入れた。
描いた地図を図法ごとに描くと次のようになる。
(描画のためのコードはこちら)

・ aitoff



・ albers(30,60)



・ azequalarea



・ azequidistant



・ bicentric(30)



・ bicentric(60)



・ bonne(60)



・ conic(30)



・ conic(60)



・ cylequalarea(0)



・ cylequalarea(60)



・ cylindrical



・ eisenlohr



・ elliptic(30)



・ elliptic(60)



・ fisheye(1.5)



・ fisheye(2)



・ gall(0)



・ gall(60)



・ gilbert



・ globular



・ gnomonic



・ guyou



・ harrison(2,10)



・ hex



・ homing(60)



・ lagrange



・ lambert(30,60)



・ laue



・ lune(10,60)



・ mecca(60)



・ mercator



・ mollweide



・ newyorker(10)



・ orthographic



・ perspective(10)



・ perspective(2)



・ polyconic



・ simpleconic(30,60)



・ sinusoidal



・ square



・ stereographic



・ tetra



・ trapezoidal(30,60)



・ vandergrinten
AM
投稿日時: 2008-1-22 17:21
管理人
登録日: 2004-7-27
居住地:
投稿: 16
Re: maps & mapproj パッケージ
描画に用いたコードは次の通り
library(mapproj)

# par(ask=T)

plot(0,0,type='n',xlim=c(-180,180),ylim=c(-90,90)); map(add=T, interior=F, projection="")                                ; map.grid()
dev.off(); png("mercator.png");
  plot(0,0,type='n',xlim=c(-pi,pi),ylim=c(-3,3)); map(add=T, interior=F, projection="mercator")                            ; map.grid()
dev.off(); png("sinusoidal.png");
  plot(0,0,type='n',xlim=c(-pi,pi),ylim=c(-3,3)); map(add=T, interior=F, projection="sinusoidal")                          ; map.grid()
dev.off(); png("cylequalarea_p=0.png");
  plot(0,0,type='n',xlim=c(-pi,pi),ylim=c(-1.5,1.5)); map(add=T, interior=F, projection="cylequalarea", parameter=0)       ; map.grid()
dev.off(); png("cylequalarea_p=60.png");
  plot(0,0,type='n',xlim=c(-pi/2,pi/2),ylim=c(-1.5,1.5)); map(add=T, interior=F, projection="cylequalarea", parameter=60)  ; map.grid()
dev.off(); png("cylindrical.png");
  plot(0,0,type='n',xlim=c(-pi,pi),ylim=c(-3,3)); map(add=T, interior=F, projection="cylindrical")                         ; map.grid()
dev.off(); png("gall_p=0.png");
  plot(0,0,type='n',xlim=c(-pi,pi),ylim=c(-1,1)); map(add=T, interior=F, projection="gall", parameter=0)                   ; map.grid()
dev.off(); png("gall_p=60.png");
  plot(0,0,type='n',xlim=c(-pi,pi),ylim=c(-1,1)); map(add=T, interior=F, projection="gall", parameter=60)                  ; map.grid()
dev.off(); png("mollweide.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-1,1)); map(add=T, interior=F, projection="mollweide")                             ; map.grid()
dev.off(); png("gilbert.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="gilbert")                               ; map.grid()
dev.off(); png("azequidistant.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="azequidistant")                         ; map.grid()
dev.off(); png("azequalarea.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="azequalarea")                           ; map.grid()
dev.off(); png("gnomonic.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="gnomonic")                              ; map.grid()
dev.off(); png("perspective_p=10.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="perspective",parameter=10)              ; map.grid()
dev.off(); png("perspective_p=2.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="perspective",parameter=2)               ; map.grid()
dev.off(); png("orthographic.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="orthographic")                          ; map.grid()
dev.off(); png("stereographic.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-2,2)); map(add=T, interior=F, projection="stereographic")                         ; map.grid()
dev.off(); png("laue.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-2,2)); map(add=T, interior=F, projection="laue")                                  ; map.grid()
dev.off(); png("fisheye_p=1.5.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="fisheye", parameter=1.5)                ; map.grid()
dev.off(); png("fisheye_p=2.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="fisheye", parameter=2)                  ; map.grid()
dev.off(); png("newyorker_p=10.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="newyorker", parameter=10)               ; map.grid()
dev.off(); png("conic_p=30.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="conic", parameter=30)                   ; map.grid()
dev.off(); png("conic_p=60.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="conic", parameter=60)                   ; map.grid()
dev.off(); png("simpleconic_p=30_60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="simpleconic", parameter=c(30,60))       ; map.grid()
dev.off(); png("lambert_p=30_60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="lambert", parameter=c(30,60))           ; map.grid()
dev.off(); png("albers_p=30_60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="albers", parameter=c(30,60))            ; map.grid()
dev.off(); png("bonne_p=60.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-3,1)); map(add=T, interior=F, projection="bonne", parameter=60)                   ; map.grid()
dev.off(); png("polyconic.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-2,2)); map(add=T, interior=F, projection="polyconic")                             ; map.grid()
dev.off(); png("aitoff.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-2,2)); map(add=T, interior=F, projection="aitoff")                                ; map.grid()
dev.off(); png("lagrange.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="lagrange")                              ; map.grid()
dev.off(); png("bicentric_p=60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="bicentric",parameter=60)                ; map.grid()
dev.off(); png("bicentric_p=30.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="bicentric",parameter=30)                ; map.grid()
dev.off(); png("elliptic_p=30.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="elliptic",parameter=30)                 ; map.grid()
dev.off(); png("elliptic_p=60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="elliptic",parameter=60)                 ; map.grid()
dev.off(); png("globular.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-2,2)); map(add=T, interior=F, projection="globular")                              ; map.grid()
dev.off(); png("vandergrinten.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="vandergrinten")                         ; map.grid()
dev.off(); png("eisenlohr.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1,1)); map(add=T, interior=F, projection="eisenlohr")                             ; map.grid()
dev.off(); png("guyou.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="guyou")                                 ; map.grid()
dev.off(); png("square.png");
  plot(0,0,type='n',xlim=c(0,3),ylim=c(-3,0)); map(add=T, interior=F, projection="square")                                 ; map.grid()
dev.off(); png("tetra.png");
  plot(0,0,type='n',xlim=c(-6,6),ylim=c(-6,6)); map(add=T, interior=F, projection="tetra")                                 ; map.grid()
dev.off(); png("hex.png");
  plot(0,0,type='n',xlim=c(0,3),ylim=c(-2,2)); map(add=T, interior=F, projection="hex")                                    ; map.grid()
dev.off(); png("harrison_p=2_10.png");
  plot(0,0,type='n',xlim=c(-2,2),ylim=c(-2,2)); map(add=T, interior=F, projection="harrison",parameter=c(2,10))            ; map.grid()
dev.off(); png("trapezoidal_p=30_60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-2,0)); map(add=T, interior=F, projection="trapezoidal",parameter=c(30,60))        ; map.grid()
dev.off(); png("lune_p=10_60.png");
  plot(0,0,type='n',xlim=c(-1,1),ylim=c(-1/2,1/2)); map(add=T, interior=F, projection="lune",parameter=c(10,60))           ; map.grid()
dev.off(); png("mecca_p=60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="mecca",parameter=60)                    ; map.grid()
dev.off(); png("homing_p=60.png");
  plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3)); map(add=T, interior=F, projection="homing",parameter=60)                   ; map.grid()
dev.off()
スレッド表示 | 新しいものから 前のトピック | 次のトピック | トップ

投稿するにはまず登録を
 

WWW を検索 meteorology.jp を検索

Powered by XOOPS 2.0 © 2001-2006 The XOOPS Project, Maitained by A. Mori
FI Theme :: XOOPS 2 Theme by ImageSquare :: Costomized by matchan and A.Mori