7-2 断裂させる

 

     グート図法に代表される断裂図法を試みました。

 

     7-2-1 サンソン図法を断裂

     あらかじめ断裂経度を決め、ブロック単位でデータをわけておく。ブロックの中央で

     サンソン図法の投影式で x , y を計算し、x軸のシフト値を加算して描画する。

 

                       ┌─ブロック─┐

                     左端   中央  右端

           ①   -20          20         60

           ②    60        130     -160    x のシフト値 = ①のブロック幅

           ③  -160        -90       -20    x のシフト値 = ①のブロック幅+②のブロック幅

 [7-2-1] 断裂サンソン図法 ( Interrupted Sinusoidal Projection )

【断裂<>中心】東経 20゚、緯度 0

                                        【断裂<>中心】東経130゚、緯度 0

                                                                                           【断裂<>中心】西経90゚、緯度 0

【経度間隔】10゚【緯度間隔】10

ダウンロード
7-2-1_Interrupted Sinusoidal.pdf
PDFファイル 286.6 KB

    7-2-2 モルワイデ図法を断裂

         データや描画については、前述『7-2-1 サンソン図法を断裂』を参照してください。

 

[図 7-2-2] 断裂モルワイデ図法 ( Interrupted Mollweide Projection )

【断裂<>中心】東経 20゚、緯度 0

                                        【断裂<>中心】東経130゚、緯度 0

                                                                                           【断裂<>中心】西経90゚、緯度 0

【経度間隔】10゚【緯度間隔】10

ダウンロード
7-2-2_Interrupted Mollweide.pdf
PDFファイル 321.7 KB

   7-2-3 エッケルト Ⅲ 図法を断裂

 

        データや描画については、前述『7-2-1 サンソン図法を断裂』を参照してください。

 

  [7-2-3] 断裂エッケルト Ⅲ 図法 ( Interrupted EckertProjection )

【断裂<>中心】東経 20゚、緯度 0

                                        【断裂<>中心】東経130゚、緯度 0

                                                                                           【断裂<>中心】西経90゚、緯度 0

【経度間隔】10゚【緯度間隔】10

ダウンロード
7-2-3_Interrupted Eckert Ⅲ.pdf
PDFファイル 304.5 KB

   7-2-4 McBryde S3 図法を断裂

         データや描画については、前述『7-2-1 サンソン図法を断裂』を参照してください。

 

  [7-2-4] 断裂McBryde S3 図法 ( Interrupted McBryde S3 Projection )

【断裂<>中心】東経 20゚、緯度 0

                                        【断裂<>中心】東経130゚、緯度 0

                                                                                           【断裂<>中心】西経90゚、緯度 0

【経度間隔】10゚【緯度間隔】10

ダウンロード
7-2-4_Interrupted McBryde S3.pdf
PDFファイル 291.0 KB

   7-2-5 Bonne Projection(ボンヌ図法)を断裂

         【 投影式 】赤道で北半球と南半球を断裂させる例です。

                                    北半球

                                     A = cotφ1 + φ1 - φ

                                   B = (λ - λ0)(cosφ) / A

                                    x = RA sin B

                                    y = R (cotφ1 - A cos B)

                南半球

                                    A' = cot( -φ2) + ( -φ2) - ( -φ)

                                   B' = (λ - λ0) [ cos ( -φ)] / A'

                                    x' = RA' sin B'

                                    y' = -R (cot(-φ2) - A' cos B') +

                                                                       2 x ( 緯度 中心経度の y の値 )

          【経度・緯度】経線:曲線 、緯線:曲線

          【 ポイント 】・データを北半球と南半球に分けて描画します。

  <基準経度・標準緯度・予備計算・シフト値 など>                
     記号  セルNo.         記述式                
      λ0      C2     ???     中心経度(希望経度を入力)                
      φ1     C3      ??     標準緯度北半球(希望緯度を入力)                
      φ2     C4  = -C3    標準緯度南半球                
      R     C5     1                
                                                                         
  <南半球へのシフト値計算>                 
     記号  セルNo.         記述式                
     経度      E3   ???  ← 実際の経度データ                
     緯度      F3   ???  ← 実際の緯度データ                
     補正
  λ-λ
0
    G3  =IF($C$2<0,IF(E3-$C$2>180,E3-$C$2-360,E3-$C$2),
                             IF(E3-$C$2<-180,E3+360-$C$2,E3-$C$2))
               
                 
      A     H3  =COT(RADIANS($C$3))+RADIANS($C$3)-RADIANS(F3)                
      B     I3  =IF(E3-$C$2<-180,RADIANS(E3-$C$2+360)*
      COS(RADIANS(F3))/H3,RADIANS(E3-$C$2)*COS(RADIANS(F3))/H3)
               
                 
      x     J3  =H3*SIN(I3)                
      y     K3  =COT(RADIANS($C$3))-H3*COS(I3)                
                                                                         
  <データ計算 など>  経度緯度のデータ毎に記述式を複写(行コピー)する。
   区分    記号  セルNo.         記述式
   

 

 
    No.      B9   ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      C9   ???  ← 実際の経度データ
     緯度      D9   ???  ← 実際の緯度データ
     補正
  λ-λ
0
    E9  =IF($C$2<0,IF(C9-$C$2>180,C9-$C$2-360,C9-$C$2),
                                             IF(C9-$C$2<-180,C9+360-$C$2,C9-$C$2))
 
      A     F9  =COT(RADIANS($C$3))+RADIANS($C$3)-RADIANS(D9)
      B     G9  =RADIANS(E9)*COS(RADIANS(D9))/F9
      x     H9  =F9*SIN(G9)
      y     I9  =COT(RADIANS($C$3))-F9*COS(G9)
   描画CH     J10  =IF(H10="",0,IF(AND(ABS(H10)>0.2,H9*H10<0),"●",
                              IF(AND(H9<>"",H10>PI()/2,ABS(H9-H10)>0.2),"●",0)))
 
   

 

 
    No. 

    K9

  ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      L9   ???  ← 実際の経度データ
     緯度      M9   ???  ← 実際の緯度データ
     補正
  λ-λ
0
    N9  =IF($C$2<0,IF(L9-$C$2>180,L9-$C$2-360,L9-$C$2),
                                                 IF(L9-$C$2<-180,L9+360-$C$2,L9-$C$2))
 
      A     O9  =COT(RADIANS(-$C$4))+RADIANS(-$C$4)-RADIANS(-M9)
      B     P9  =RADIANS(N9)*COS(RADIANS(-M9))/O9
      x     Q9  =O9*SIN(P9)
      y     R9  =-(COT(RADIANS(-$C$4))-O9*COS(P9))+2*$K$3
   描画CH     S10  =IF(Q10="",0,IF(AND(ABS(Q10)>0.2,Q9*Q10<0),"●",
                               IF(AND(Q9<>"",Q10>PI()/2,ABS(Q9-Q10)>0.2),"●",0)))
 

 

[7-2-5] 断裂ボンヌ図法 ( Interrupted Bonne Projection )

                               【地図主点(中心)】東経 135   【断裂位置】赤道

                                                 【標準緯度】北緯30゚、南緯30 【経度間隔】15  【緯度間隔】10

ダウンロード
7-2-5_Bonne_赤道断裂.pdf
PDFファイル 358.3 KB

   7-2-6 Collignon Projectionを断裂

              Collignon は外形が三角形で歪がおおきいため HEALPi などに利用される以外使い道

              が少ない図法です。そこで断裂変形させた 投影図を描画してみました。

    7-2-6-1 Collignon Projection

                    ----▲と▼を赤道で断裂(通常は緯度90-90゚まで描画するが、

                         赤道で分断し)そして合体させ形に

                  【 投影式 】 北半球

  x = 2 R (λ - λ0) (1 - sinφ)1/2 / π1/2

  y = π1/2R [l - (1- sinφ)1/2]

  南半球

  x' = 2 R (λ - λ0) (1 - sin(-φ)1/2 / π1/2

  y' = - π1/2R{l - [1- sin(-φ)]1/2 }

                 【経度・緯度】経線:直線 、緯線:直線

                  【 ポイント 】・データを北半球と南半球に分けて描画します。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セルNo.         記述式                
      λ0      C2    ???    中心経度(希望経度を入力)                
      R     C3      1                
                                                                         
  <データ計算 など>  経度緯度のデータ毎に記述式を複写(行コピー)する。
   区分    記号  セルNo.         記述式
   

 

 
    No.      B9   ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      C9   ???  ← 実際の経度データ
     緯度      D9   ???  ← 実際の緯度データ
     補正
  λ-λ
0
    E9  =IF($C$2<0,IF(C9-$C$2>180,C9-$C$2-360,C9-$C$2),
                                             IF(C9-$C$2<-180,C9+360-$C$2,C9-$C$2))
 
      x     F9  =COT(RADIANS($C$3))+RADIANS($C$3)-RADIANS(D9)
      y     G9  =RADIANS(E9)*COS(RADIANS(D9))/F9
   描画CH     H10  =IF(H10="",0,IF(AND(ABS(H10)>0.2,H9*H10<0),"●",
                              IF(AND(H9<>"",H10>PI()/2,ABS(H9-H10)>0.2),"●",0)))
 
   

 

 
    No.      I9    ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      J9    ???  ← 実際の経度データ
     緯度      K9    ???  ← 実際の緯度データ
     補正
  λ-λ
0
    L9   =IF($C$2<0,IF(J9-$C$2>180,J9-$C$2-360,J9-$C$2),
                                                        IF(J9-$C$2<-180,J9+360-$C$2,J9-$C$2))
 
      x'     M9   =2*RADIANS(L9)*SQRT(1-SIN(RADIANS(-K9)))/SQRT(PI())
      y'     N9   =-(SQRT(PI())*(1-SQRT(1-SIN(RADIANS(-K9)))))
   描画CH     O10   =IF(M10="",0,IF(AND(ABS(M10)>0.2,M9*M10<0),"●",
                                      IF(AND(M9<>"",M10>PI()/2,ABS(M9-M10)>0.2),"●",0))) 
 

 

[7-2-6-1] 断裂コリニョン図法 ( Interrupted Collignon Projection )

                 ▲と▼を赤道で断裂そして合体させひし形に

                              【断裂内容】北半球と南半球を別々に描画し、赤道にて合成 

                              【地図中心】東経135゚、緯度 0゚ 【経度間隔】15゚【緯度間隔】10

ダウンロード
7-2-6-1_Collignon_変形断裂①.pdf
PDFファイル 247.6 KB

    7-2-6-2 Collignon Projection

                     ----全球を経度で2分し、▲と▼を赤道で断裂そして合体させ◇◇形に

                【 投影式 】 北半球

 ①  x = 2R (λ - λ0)(1 - sinφ)1/21/2

                                      y = π1/2R [ l - (1- sinφ)1/2]

                                               A = 2π1/2 緯度0xの幅を求めた値

                                                                          ②と④に用いる

                           x = 2R (λ - λ1)(1 - sinφ)1/21/2 + A

                                       y = π1/2R [ l - (1- sinφ)1/2]

                     南半球

                                          x' = 2R (λ - λ0)(1 - sin()1/21/2

                                               y' = - π1/2R{l - [1- sin()]1/2}

 ④ x' = 2R (λ - λ1)(1 - sin()1/21/2 + A

                                               y' = - π1/2R{l - [1- sin()]1/2}

                【経度・緯度】経線:直線 、緯線:直線

                 【 ポイント 】・データを北半球と南半球に2分割さらにそれを2分割、つまり4つ

                                               に分割して描画します。

  <基準経度・標準緯度・予備計算・シフト値 など>                
     記号  セルNo.         記述式                
      λ0      C2    ???    左側の ◇ の中心経度(希望経度を入力)                
      λ1     C3    ???    右側の ◇ の中心経度(希望経度を入力)                
      A     C4  =2*SQRT(PI())                
      R     C5        1                
                                                                         
  <データ計算 など>  経度緯度のデータ毎に記述式を複写(行コピー)する。
   区分    記号  セルNo.         記述式
   
 
 

 
    No.      B9   ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      C9   ???  ← 実際の経度データ
     緯度      D9   ???  ← 実際の緯度データ
     補正
  λ-λ
0
    E9  =IF($C$2<0,IF(C9-$C$2>180,C9-$C$2-360,C9-$C$2),
                                             IF(C9-$C$2<-180,C9+360-$C$2,C9-$C$2))
 
      x     F9  =2*RADIANS(E9)*SQRT(1-SIN(RADIANS(D9)))/SQRT(PI())
      y     G9  =SQRT(PI())*(1-SQRT(1-SIN(RADIANS(D9))))
   
 
 

 
    No.      H9   ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      I9   ???  ← 実際の経度データ
     緯度      J9   ???  ← 実際の緯度データ
     補正
  λ-λ
0
    K9  =IF($C$3<0,IF(I9-$C$3>180,I9-$C$3-360,I9-$C$3),
                                                  IF(I9-$C$3<-180,I9+360-$C$3,I9-$C$3))
 
      x     L9  =2*RADIANS(K9)*SQRT(1-SIN(RADIANS(J9)))/SQRT(PI())+$C$4
      y     M9  =SQRT(PI())*(1-SQRT(1-SIN(RADIANS(J9))))
   
 
 

 
    No.      N9    ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      O9    ???  ← 実際の経度データ
     緯度      P9    ???  ← 実際の緯度データ
     補正
  λ-λ
0
    Q9   =IF($C$2<0,IF(O9-$C$2>180,O9-$C$2-360,O9-$C$2),
                                           IF(O9-$C$2<-180,O9+360-$C$2,O9-$C$2))
 
      x'     R9   =2*RADIANS(Q9)*SQRT(1-SIN(RADIANS(-P9)))/SQRT(PI())
      y'     S9   =-(SQRT(PI())*(1-SQRT(1-SIN(RADIANS(-P9)))))
   
 
 

 
    No.      T9    ???  ← 式の複写時、一旦ソートしたり,解析時に役立つので連番号を付与
     経度      U9    ???  ← 実際の経度データ
     緯度      V9    ???  ← 実際の緯度データ
     補正
  λ-λ
0
    W9   =IF($C$3<0,IF(U9-$C$3>180,U9-$C$3-360,U9-$C$3),
                                            IF(U9-$C$3<-180,U9+360-$C$3,U9-$C$3))
 
      x'     X9   =2*RADIANS(W9)*SQRT(1-SIN(RADIANS(-V9)))/SQRT(PI())+$C$4
      y'     Y9   =-(SQRT(PI())*(1-SQRT(1-SIN(RADIANS(-V9)))))

 

[7-2-6-2] 断裂コリニョン図法 ( Interrupted Collignon Projection )

             【断裂内容】西経135゚にて断裂、北半球と南半球を別々に描画し、赤道にて合成

                  【地図中心】東経135゚、緯度 0゚ 【経度間隔】15゚【緯度間隔】10

ダウンロード
7-2-6-2_Collignon_変形断裂②.pdf
PDFファイル 255.0 KB