6-7-1  描画形状が〇、一部が〇、∞のようなもの

 

 6-7-1-1 Bacon

        【 投影式 】F = [ ( π / 2) 2 / | λ - λ0 | + | λ - λ0 | ] / 2

y = R ( π / 2 ) sinφ

 x = ±R [ | λ - λ0 | - F + (F2 - y2 / R2)1/2]

 ± の符号は ( λ - λ0 ) を用いる

 もし  λ = λ0 , x = 0

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

        【 ポイント 】・極は点にならず、平となる。

                                          ・投影中心から±90°の経線は、円弧となります。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      F     F9  =((PI()/2)^2/RADIANS(ABS(E9))+RADIANS(ABS(E9)))/2
      x     G9  =IF(E9=0,0,SIGN(E9)*(RADIANS(ABS(E9))-F9+SQRT(F9^2-H9^2)))
      y     H9  =PI()/2*SIN(RADIANS(D9))
     端CH     I10  =IF(G10="",0,IF(AND(ABS(G10)>0.2,G9*G10<0),"●",
                              IF(AND(G9<>"",G10>PI()/2,ABS(G9-G10)>0.2),"●",0))) 
 
        左右の端をチェックしている。"●"の付いた行の前に一行挿入すると、勝手に左端と右端が線で結ばれない。

 

[6-7-1-1] Bacon 図法 ( Bacon Projection )

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

     【備考】中心経度の±90゚で描画すると半球となる。それを Bacon Globular 図法という。

ダウンロード
6-7-1-1_Bacon_描画_A.pdf
PDFファイル 276.6 KB

 

  6-7-1-2 Bacon Globular

      Bacon の一部を円(半球)として切り出し、さらに並列に並べ全球とした。

 

 

[6-7-1-2]  Bacon Globular 図法 ( Bacon Globular Projection )

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

ダウンロード
6-7-1-2_Bacon Globular.pdf
PDFファイル 147.9 KB

 

[6-7-1-2]  Bacon Globular 図法 - 半球並列 ( Bacon Globular Projection )

【地図主点(中心)】左半球:東経 135゚、緯度 0゚ 右半球:西経45゚、緯度 0

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

ダウンロード
6-7-1-2_Bacon Globular_半球並列.pdf
PDFファイル 271.2 KB

  6-7-1-3 Apian I

         【 投影式 】y = Rφ

 F = [ ( π / 2 )2 / | λ - λ0 | + | λ - λ0 | ] / 2

 x = ± R [ | λ - λ0 | - F + ( F2 - y2 / R2 )1/2 ]

 ± の符号は ( λ - λ0 ) を用いる

 もし  λ = λ0  ,  x = 0

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

        【 ポイント 】・極は点にならず、平となる。

                                         ・投影中心から±90°の経線は、円弧となります。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      F     F9  =((PI()/2)^2/ABS(RADIANS(E9))+ABS(RADIANS(E9)))/2
      x     G9  =IF(C9=$C$2,0,IF(E9<0,-(ABS(RADIANS(E9))-F9+SQRT(F9^2-H9^2)),                                                         (ABS(RADIANS(E9))-F9+SQRT(F9^2-H9^2))))
      y     H9  =RADIANS(D9)
     端CH     I10  =IF(G10="",0,IF(AND(ABS(G10)>0.2,G9*G10<0),"●",
                              IF(AND(G9<>"",G10>PI()/2,ABS(G9-G10)>0.2),"●",0))) 
 
        左右の端をチェックしている。"●"の付いた行の前に一行挿入すると、勝手に左端と右端が線で結ばれない。

 

[6-7-1-3]  Apian I 図法 ( Apian I Projection )

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

     【備考】中心経度の±90゚で描画すると半球となる。それを Apian Globular I 図法という。

ダウンロード
6-7-1-3_Apian I.pdf
PDFファイル 284.9 KB

  6-7-1-4 Apian Ⅰ Globular

     Apian I の一部を円(半球)として切り出し、さらに並列に並べ全球とした

 

[6-7-1-4]  Apian ⅠGlobular 図法 ( Apian Globular I Projection )

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

ダウンロード
6-7-1-4_Apian Globular I_半球.pdf
PDFファイル 153.8 KB

[6-7-1-4]  Apian ⅠGlobular 図法 - 半球並列 ( Apian Globular I Projection )

【地図主点(中心)】左半球:東経 135゚、緯度 0゚ 右半球:西経45゚、緯度 0

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

ダウンロード
6-7-1-4_Apian Globular I_半球並列.pdf
PDFファイル 278.9 KB

 

  6-7-1-5 Fournier Ⅰ

         【 投影式 】もし (λ - λ0) = 0 , または |φ| = π /2 , x = 0 , y = Rφ

 もしφ = 0,  x = R( λ - λ0) , y = 0

 もし | λ - λ0| = π / 2,

 x = R ( λ - λ0) cosφ

 y = R (π / 2) sinφ

 C = π2 / 4

 P = | π sin φ |

 S = (C - φ2) / (P - 2| φ | )

 A = ( λ - λ0)2 / C - 1

 y = ±R ({ S2 - A [ C- PS- ( λ - λ0)2]}1/2 - S ) / A

 符号は φ

 x = ±R ( λ - λ0)[ 1 - y2 / R2 C )]1/2

  参考文献では『符号は (λ - λ0)』とあるが、

 投影式内部に (λ - λ0) があるため記入式には

 符号を付与しなかった。

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

        【 ポイント 】・投影中心から±90°の経線は、円弧となります。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      C     F9  =PI()^2/4
      P     G9  =ABS(PI()*SIN(RADIANS(D9)))
      S     H9  =(F9-(RADIANS(D9))^2)/(G9-2*ABS(RADIANS(D9)))
      A     I9  =(RADIANS(E9))^2/F9-1
      x     J9  =IF(OR(E9=0,ABS(RADIANS(D9))=PI()/2),0,
       IF(D9=0,RADIANS(E9),IF(ABS(E9)=PI()/2,RADIANS(E9)*
                 COS(RADIANS(D9)),RADIANS(E9)*SQRT(1-K9^2/F9))))
 
 
      y     K9  =IF(OR(E9=0,ABS(RADIANS(D9))=PI()/2),RADIANS(D9),IF(D9=0,0,
             IF(ABS(RADIANS(E9))=PI()/2,(PI()/2)*SIN(RADIANS(D9)),
                    IF(D9<0,-(SQRT(H9^2-I9*(F9-G9*H9-(RADIANS(E9))^2))-H9)/I9,
                                      (SQRT(H9^2-I9*(F9-G9*H9-(RADIANS(E9))^2))-H9)/I9))))
 
 
 
     端CH     L9  =IF(OR(C9=$C$2-180,AND(E9<>"",C9<>"",E9>$C$2-180,C9<$C$2-180),
                                       AND(E9<>"",C9<>"",E9<$C$2-180,C9>$C$2-180)),"●",0)
 

 

[6-7-1-5]  Fournier I 図法 ( Fournier I Projection )

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

【備考】中心経度の±90゚で描画すると半球となる。それを Fournier I Globular 図法という。

 

ダウンロード
6-7-1-5_Fournier I_全球.pdf
PDFファイル 332.7 KB

  6-7-1-6 Fournier Ⅰ  Globular

 

[6-7-1-6]  Fournier I Globular 図法

                         ( Fournier Globular I Projection ) in Two Hemispheres

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

ダウンロード
6-7-1-6_Fournier I Globular_半球.pdf
PDFファイル 189.6 KB

[6-7-1-6]  Fournier I Globular 図法 半球並列

                         ( Fournier Globular I Projection ) in Two Hemispheres

【地図主点(中心)】左半球:東経 135゚、緯度 0゚ 右半球:西経45゚、緯度 0

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

ダウンロード
6-7-1-6_Fournier I Globular_半球並列.pdf
PDFファイル 349.9 KB

  6-7-1-7 Nicolosi Globular

        【 投影式 】もし  (λ - λ0) = 0  または  | φ | = π / 2

                                      x = 0 , y = Rφ

 もし  φ = 0 

 x = R (λ - λ0) , y = 0

 もし | λ - λ0| = π/2

 x = R (λ - λ0) cosφ, y = R (π / 2) sinφ

 それ以外の場合は、

 b = π / [ 2 (λ - λ0)] – 2 (λ - λ0) / π

 c = 2 φ / π

 d = (1- c2) / (sinφ- c )

 M = [b (sinφ) / d - b/ 2] / (1 + b2/ d2)

 N = [d2 (sin φ) / b2 + d/ 2] / (1 + d2 / b2)

 x = (π R/2){M ± [M2 + cos2 φ/ (1 + b2 /d2)]1/2}

                                                                                符号は (λ - λ0)

 y = (πR/2){N ± [N2 - ( d2(sin2φ) / b2 +

 d sin φ- 1) / (1+ d2 / b2 )]1/2}

                                                              符号はφの逆

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

         【 ポイント 】・±90°の外付近が特に歪みます。

                      ・投影中心から±90°の経線は、円弧となります。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      b     F9  =PI()/(2*RADIANS(E9))-2*RADIANS(E9)/PI()
      c     G9  =2*RADIANS(D9)/PI()
      d     H9  =(1-G9^2)/(SIN(RADIANS(D9))-G9)
      M     I9  =(F9*SIN(RADIANS(D9))/H9-F9/2)/(1+F9^2/H9^2)
      N     J9  =(H9^2*SIN(RADIANS(D9))/F9^2+H9/2)/(1+H9^2/F9^2)
      x     K9  =IF(ABS(RADIANS(E9))=PI()/2,RADIANS(E9)*COS(RADIANS(D9)),
      IF(OR(E9=0,ABS(RADIANS(D9))=PI()/2),0,IF(D9=0,RADIANS(E9),
                          IF(E9<0,(PI()/2)*(I9-SQRT(I9^2+COS(RADIANS(D9))^2/
                                                       (1+F9^2/H9^2))),(PI()/2)*(I9+SQRT(I9^2+
                                                            COS(RADIANS(D9))^2/(1+F9^2/H9^2)))))))
 
 
 
 
      y     L9  =IF(OR(E9=0,ABS(RADIANS(D9))=PI()/2),RADIANS(D9),IF(D9=0,0,
      IF(ABS(RADIANS(E9))=PI()/2,(PI()/2)*SIN(RADIANS(D9)),IF(D9<0,
           (PI()/2)*(J9+SQRT(J9^2-(H9^2*(SIN(RADIANS(D9)))^2/F9^2+H9*
               SIN(RADIANS(D9))-1)/(1+H9^2/F9^2))),(PI()/2)*(J9-SQRT(J9^2-(H9^2*
                   (SIN(RADIANS(D9)))^2/F9^2+H9*SIN(RADIANS(D9))-1)/(1+H9^2/F9^2)))))))
 
 
 
 
     端CH     M10  =IF(K10="",0,IF(AND(ABS(K10)>0.2,K9*K10<0),"●",IF(AND(K9<>"",
                                                                  K10>PI()/2,ABS(K9-K10)>0.2),"●",0)))
 

 

[6-7-1-7]  Nicolosi Globular 図法  ( Nicolosi Globular Projection )

                        通常は半球描画ですが、ここでは歪部分を含む全球描画

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

 

ダウンロード
6-7-1-7_Nicolosi Globular_全球.pdf
PDFファイル 319.3 KB

[6-7-1-7]  Nicolosi Globular 図法  ( Nicolosi Globular Projection )

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

ダウンロード
6-7-1-7_Nicolosi Globular_半球.pdf
PDFファイル 191.3 KB

[6-7-1-7]  Nicolosi Globular 図法 - 半球並列 ( Nicolosi Globular Projection )

【地図主点(中心)】左半球:東経 135゚、緯度 0゚ 右半球:西経45゚、緯度 0

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

ダウンロード
6-7-1-7_Nicolosi Globular_半球並列.pdf
PDFファイル 354.1 KB

  6-7-1-8 Lagrange

         【 投影式 】もし  |φ| = π / 2

                                             x = 0 , y = ±2R , 符号はφ

 A1 = [ ( 1 + sin φ1) / ( 1 - sin φ1) ]1/(2 W)

 A = [ (1 + sin φ) / ( 1- sin φ) ]1/(2W)

 V = A / A1

 C = ( V + 1/V ) / 2 + cos [ (π - π0 ) / W]

      x = 2R sin [ ( λ - λ0 ) / W] / C

      y = R (V - 1/V) / C

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

        【 ポイント 】・W =1,φ1 =0 にすると Stereographic Projection なる

 

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

【 ポイント 】

        ・W =1,φ1 =0 にすると

              Stereographic Projection なる

                 ・外周は、円となる。

         ・経線の高緯度(±8990°)

           は緯度ピッチを細かく

           しないと円弧になりません。

           左図が参考例です。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セルNo.         記述式                
      λ0     C2    ???  ← 中心経度(希望経度を入力)                
      W     C3        2    ← W=1,φ1=0 にすると Stereographic Projection なる                
      φ1     C4         0    ← W=1,φ1=0 にすると Stereographic Projection なる                
      R      C5         1                
      A1      F2 

 =((1+SIN(RADIANS($C$4)))/

                        (1-SIN(RADIANS($C$4))))^(1/(2*$C$3))

               
                                                                         
  <データ計算 など>  経度緯度のデータ毎に記述式を複写(行コピー)する。
     記号  セル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  =((1+SIN(RADIANS(D9)))/(1-SIN(RADIANS(D9))))^(1/(2*$C$3))
      V     G9  =F9/$F$2
      C     H9  =(G9+1/G9)/2+COS(RADIANS(E9/$C$3))
      x     I9  =IF(ABS(D9)=PI()/2,0,2*SIN(RADIANS(E9)/$C$3)/H9)
      y     J9  =IF(ABS(D9)=PI()/2,SIGN(D9)*2,(G9-1/G9)/H9)
     端CH     K10  =IF(I10="",0,IF(AND(ABS(I10)>0.2,I9*I10<0),"●",IF(AND(I9<>"",
                                                                   I10>PI()/2,ABS(I9-I10)>0.2),"●",0)))
 
        左右の端をチェックしている。"●"の付いた行の前に一行挿入すると、勝手に左端と右端が線で結ばれない。

 

[6-7-1-8]  Lagrange 図法 ( Lagrange Projection )

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

ダウンロード
6-7-1-8_Lagrange.pdf
PDFファイル 334.8 KB

  6-7-1-9 Vander Grinten

         【 投影式 】もし  φ = 0   ,  x = R (λ - λ0 )  ,  y = 0

 もし  |φ| = π/2 , x = 0  , y = ± π R

                      符号はφに同じ

 B = | 2φ / π |  , C = (1- B2 )1/2

 もし  ( λ - λ0 ) = 0  , x = 0  ,  y = ±π RB / ( 1 + C )

 符号はφに同じ

 A = | π / ( λ - λ0) - ( λ - λ0) / π | / 2

 G = C / ( B + C - 1)

 P = G (2 / B - 1)

 Q = A2 + G

 S = P2 + A2

 T = G - P2

 x = ± π R{AT + [A2T2- S(G2- P2)]1/2} / S

 符号は(π - π0)

 y = ± π R{PQ -A [(A2 + l )S- Q2 ]1/2} / S

 符号はφ

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

       【 ポイント 】・外周は、円となる。

                      ・経線の高緯度(±8990°)は緯度ピッチを細かくしないと

                       円弧になりません。

                       (6-7-8 Lagrange の図を参照してください)

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      A     F9  =ABS(PI()/RADIANS(E9)-RADIANS(E9)/PI())/2
      B     G9  =ABS(2*RADIANS(D9)/PI())
      C     H9  =SQRT(1-G9^2)
      G     I9  =H9/(G9+H9-1)
      P     J9  =I9*(2/G9-1)
      Q     K9  =F9^2+I9
      S     L9  =J9^2+F9^2
      T     M9  =I9-J9^2
      x     N9  =IF(D9=0,RADIANS(E9),IF(E9=0,0,IF(RADIANS(ABS(D9))=PI()/2,0,
                            SIGN(E9)*(PI()*(F9*M9+SQRT(F9^2*M9^2-L9*(I9^2-J9^2)))/L9))))
 
      y     O9   =IF(D9=0,0,IF(E9=0,SIGN(D9)*PI()*G9/(1+H9),IF(RADIANS(ABS(D9))=PI()/2,
           SIGN(D9)*PI(),SIGN(D9)*(PI()*(J9*K9-F9*SQRT((F9^2+1)*L9-K9^2))/L9))))
 
     端CH     P10   =IF(N10="",0,IF(AND(ABS(N10)>0.2,N9*N10<0),"●",IF(AND(N9<>"",
                                                               N10>PI()/2,ABS(N9-N10)>0.2),"●",0)))
 
        左右の端をチェックしている。"●"の付いた行の前に一行挿入すると、勝手に左端と右端が線で結ばれない。

 

[6-7-1-9]  Van der Grinten I  図法 ( Van der Grinten I Projection )

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

ダウンロード
6-7-1-9_Van der Grinten I.pdf
PDFファイル 324.3 KB

  6-7-1-10 Vander Grinten

        【 投影式 】A = | π / ( λ - λ0) - ( λ - λ0) / π | / 2

  B = | 2φ / π |

   C = (1- B2 )1/2

 もし (λ - λ0) = 0

 x = 0

 y = ± π R B / ( 1 + C )  符号はφに同じ

 それ以外は

 x1 = [ C (1+A2)1/2-A C2] / (1+A2B2)

 x = ± π R x1     符号は(π - π0)

 y = ± π R (1 - x12-2 A x1)1/2    符号はφに同じ

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

         【 ポイント 】・外周は、円となる。

                       ・経線の高緯度(±8990°)は緯度ピッチを細かくしないと

円弧になりません。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      A     F9  =ABS(PI()/RADIANS(E9)-RADIANS(E9)/PI())/2
      B     G9  =ABS(2*RADIANS(D9)/PI())
      C     H9  =SQRT(1-G9^2)
      x1     I9  =(H9*SQRT(1+F9^2)-F9*H9^2)/(1+F9^2*G9^2)
      x     J9  =IF(E9=0,0,SIGN(E9)*PI()*I9)
      y     K9  =IF(E9=0,SIGN(D9)*PI()*G9/(1+H9),SIGN(D9)*PI()*SQRT(1-I9^2-2*F9*I9))
     端CH     L10  =IF(J10="",0,IF(AND(ABS(J10)>0.2,J9*J10<0),"●",
                                          IF(AND(J9<>"",J10>PI()/2,ABS(J9-J10)>0.2),"●",0)))
 
        左右の端をチェックしている。"●"の付いた行の前に一行挿入すると、勝手に左端と右端が線で結ばれない。

 

[6-7-1-10]  Van der Grinten Ⅱ 図法 ( Van der GrintenProjection )

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

ダウンロード
6-7-1-10_Van der Grinten Il.pdf
PDFファイル 346.0 KB

  6-7-1-11 Van der Grinten

         【 投影式 】 A = | π / ( λ - λ0) - ( λ - λ0) / π | / 2

 B = | 2φ / π |

 C = (1- B2 )1/2

 もし  φ = 0

  x = R (λ - λ0 )

 y = 0

 もし  ( λ - λ0 ) = 0

 x = 0

 y = ± π R B / ( 1 + C )  符号はφに同じ

 それ以外なら

 y1 = B / ( 1 + C )

 x = ± π R [(A2 + l -y12)1/2 - A ]   符号は(π - π0)

                                             y = ± π R y1  符号はφに同じ

 

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

 

【 ポイント 】・外周は、円となる。

              ・経線の高緯度(±8990°)

                                緯度ピッチを細かくしないと

                円弧になりません。

                左図が参考例です。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      A     F9  =ABS(PI()/RADIANS(E9)-RADIANS(E9)/PI())/2
      B     G9  =ABS(2*RADIANS(D9)/PI())
      C     H9  =SQRT(1-G9^2)
      x     I9  =IF(E9=0,0,IF(D9=0,RADIANS(E9),SIGN(E9)*PI()*(SQRT(F9^2+1-J9^2)-F9)))
      y1     J9  =G9/(1+H9)
      y     K9  =SIGN(D9)*PI()*J9
   描画CH     L10  =IF(OR(E10=$C$2-180,AND(E9<>"",E10<>"",E9>$C$2-180,E10<$C$2-180),
                                        AND(E9<>"",E10<>"",E9<$C$2-180,E10>$C$2-180)),"●",0)
 

 

[6-7-1-11]  Van der Grinten  図法 ( Van der GrintenProjection )

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

ダウンロード
6-7-1-11_Van der Grinten Ill.pdf
PDFファイル 285.5 KB

  6-7-1-12 Van der Grinten IV

         【 投影式 】もし  φ = 0

 x = R ( λ - λ0)  ,  y = 0

 もし  ( λ - λ0) = 0 または | φ | = π / 2,

 x = 0   ,  y = Rφ

 その他の場合

 B = | 2φ/π |

 C = (8B- B4- 2B2 -5) / (2B3 -2B2)

 D = ±{ [ 2 ( λ - λ0) / π + π / (2( λ - λ0))]2- 4}1/2

                         符号は、 ( | λ - λ0 | - π/2) に同じ

 F = (B + C)2(B2 + C2D2- 1)+ (1- B2){B2[(B + 3C)2+ 4C2] +

  12BC3 + 4C4}

 x1 = {D[(B + C)2 + C2- 1] + 2F1/2} / [4(B+ C)2 + D2]

 x = ±π R x1 / 2   符号は、( λ - λ0) に同じ

 y = ± (πR/2)(1+D|x1| - x12)1/2  符号はφに同じ

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

        【 ポイント 】・外周は、∞のような形となる。

            ・中心経度から±90°の経線は円弧となるため 半球並列の

              図が投影できます。

 

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      B     F9  =ABS(2*RADIANS(D9)/PI())
      C     G9  =(8*F9-F9^4-2*F9^2-5)/(2*F9^3-2*F9^2)
      D     H9  =IF((ABS(RADIANS(E9))-PI()/2)<0,-SQRT((2*RADIANS(E9)/PI()+PI()/
       (2*RADIANS(E9)))^2-4),SQRT((2*RADIANS(E9)/PI()+PI()/(2*RADIANS(E9)))^2-4))
 
      F     I9  =(F9+G9)^2*(F9^2+G9^2*H9^2-1)+(1-F9^2)*(F9^2*
                                                      ((F9+3*G9)^2+4*G9^2)+12*F9*G9^3+4*G9^4)
 
      x1     J9  =(H9*((F9+G9)^2+G9^2-1)+2*SQRT(I9))/(4*(F9+G9)^2+H9^2)
      x     K9  =IF(D9=0,RADIANS(E9),IF(OR(E9=0,ABS(RADIANS(D9))=PI()/2),0,SIGN(E9)*PI()*J9/2))
      y     L9  =IF(D9=0,0,IF(OR(E9=0,ABS(RADIANS(D9))=PI()/2),RADIANS(D9),
                                                   SIGN(D9)*(PI()/2)*SQRT(1+H9*ABS(J9)-J9^2)))
 
   描画CH     M10  =IF(K10="",0,IF(AND(ABS(K10)>0.2,K9*K10<0),"●",IF(AND(K9<>"",
                                                              K10>PI()/2,ABS(K9-K10)>0.2),"●",0)))
 

 

[6-7-1-12]  Van der Grinten  図法 ( Van der GrintenProjection )

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

ダウンロード
6-7-1-12_Van der Grinten IV.pdf
PDFファイル 329.8 KB

[6-7-1-12]  Van der Grinten  図法 半球並列  ( Van der GrintenProjection )

【左半球主点(中心)】東経 135゚、緯度 0゚  【右半球主点(中心)】西経 45゚、緯度 0

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

ダウンロード
6-7-1-12_Van der Grinten IV_半球並列.pdf
PDFファイル 322.9 KB

  6-7-1-13 Eisenlohr

         【 投影式 】S1 = sin [ ( λ - λ0 )/2 ]

 C1 = cos [ (λ - λ0 )/2 ]

 T = sin (φ/2) / [ cos (φ/2) + ( 2 cosφ)1/2C1 ]

 C = [ 2 / ( 1+T2) ]1/2

 V = { [ cos (φ/2) + ((cos φ)/2)1/2(C1+S1) ] /

  [ cos (φ/2) + ((cosφ)/2)1/2(C1-S1) ] }1/2

 x = ( 3 + 81/2) R [ - 2 ln V + C (V - 1/V) ]

 y = ( 3 + 81/2) R [ - 2 arctan T + C T (V + 1/V) ]

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

         【 ポイント 】・外周は、∞のような形となる。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セルNo.         記述式                
      λ0     C2     ???  ← 中心経度(希望経度を入力)                
       R      C3    1                
     係数     C4  =3+SQRT(8)   ←3 + 81/2                
                                                                         
  <データ計算 など>  経度緯度のデータ毎に記述式を複写(行コピー)する。
     記号  セル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))
 
      S1     F9  =SIN(RADIANS(E9)/2)
      C1     G9  =COS(RADIANS(E9)/2)
      T     H9  =SIN(RADIANS(D9/2))/(COS(RADIANS(D9/2))+SQRT(2*COS(RADIANS(D9)))*G9)
      C     I9  =SQRT(2/(1+H9*H9))
      V     J9  =SQRT((COS(RADIANS(D9/2))+SQRT(COS(RADIANS(D9))/2)*
             (G9+F9))/(COS(RADIANS(D9/2))+SQRT(COS(RADIANS(D9))/2)*(G9-F9)))
 
      x     K9  =$C$4*(-2*LN(J9)+I9*(J9-1/J9))
      y     L9  =$C$4*(-2*ATAN(H9)+I9*H9*(J9+1/J9))
   描画CH     M10  =IF(K10="",0,IF(AND(ABS(K10)>0.2,K9*K10<0),"●",
                                     IF(AND(K9<>"",K10>PI()/2,ABS(K9-K10)>0.2),"●",0)))
 

 

[6-7-1-13]  アイゼンロール図法( Eisenlohr projection

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

ダウンロード
6-7-1-13_Eisenlohr projection.pdf
PDFファイル 349.9 KB

  6-7-1-14 August Epicycloidal

        【 投影式 】Cl = [ 1 - tan2(φ/2) ]1/2

C = 1 + Cl cos [ (λ – λ0) / 2 ]

x1 = sin [ (λ - λ0)/2 ] C1/C

 y1 = tan ( φ/2) / C

 x = 4 R x1 (3 + x12 - 3y12) / 3

 y = 4 R y1 (3 + 3x12 - y12) / 3

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

        【 ポイント 】・外周は、∞のような形となる。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      C1     F9  =SQRT(1-TAN(RADIANS(D9/2))*TAN(RADIANS(D9/2)))
      C     G9  =1+F9*COS(RADIANS(E9)/2)
      x1     H9  =SIN(RADIANS(E9)/2)*F9/G9
      y1     I9  =TAN(RADIANS(D9/2))/G9
      x     J9  =4*H9*(3+H9*H9-3*I9*I9)/3
      y     K9  =4*I9*(3+3*H9*H9-I9*I9)/3
   描画CH     L10  =IF(J10="",0,IF(AND(ABS(J10)>0.2,J9*J10<0),"●",
                               IF(AND(J9<>"",J10>PI()/2,ABS(J9-J10)>0.2),"●",0)))
 

 

[6-7-1-14]  August Epicycloidal Projection

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

ダウンロード
6-7-1-14_August Epicycloidal.pdf
PDFファイル 349.8 KB

  6-7-1-15 Ortelius Oval

        【 投影式 】F = [ (π/2)2 / | λ - λ0 | + | λ - λ0 | ] / 2

もし  | λ - λ0 | <= π / 2

 x = ± R [ | λ - λ0 | - F + ( F2 - y2/R2)1/2]

                                                         ± の符号は ( λ - λ0 ) を用いる

                                           y = Rφ

 もし  | λ - λ0 | > π / 2

  x = ± R {[ (π/2)2 - φ2]1/2 + | λ - λ0 | - π/2}

 ± の符号は ( λ - λ0 ) を用いる

                                            y = Rφ

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

         【 ポイント 】・極は点ではなく直線となります。

            ・中心経度から±90°の経線は円弧となるため半球並列の

            図が投影できます。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セル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))
 
      F      F9   =((PI()/2)^2/RADIANS(ABS(E9))+RADIANS(ABS(E9)))/2
      x      G9  =IF(ABS(RADIANS(E9))<=PI()/2,SIGN(E9)*(RADIANS(ABS(E9))-F9+
          SQRT(F9^2-H9^2)),SIGN(E9)*(SQRT((PI()/2)^2-
                  RADIANS(D9)^2)+ABS(RADIANS(E9))-PI()/2))
 
 
      y      H9   =RADIANS(D9)
    端CH     I10   =IF(G10="",0,IF(AND(ABS(G10)>0.2,G9*G10<0),"●",
                                       IF(AND(G9<>"",G10>PI()/2,ABS(G9-G10)>0.2),"●",0))) 
 
        左右の端をチェックしている。"●"の付いた行の前に一行挿入すると、勝手に左端と右端が線で結ばれない。

 

[6-7-1-15]  Ortelius Oval 図法 ( Ortelius Oval Projection )

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

                                              【備考】中心経度の±90゚で描画すると半球となる。

 

ダウンロード
6-7-1-15_Ortelius Oval.pdf
PDFファイル 285.0 KB

[6-7-1-15]  Ortelius Oval 図法 - 半球並列

                                  ( Ortelius Oval Projection,in Two Hemispheres )

【地図主点(中心)】左半球:東経 135゚、緯度 0゚ 右半球:西経45゚、緯度 0

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

ダウンロード
6-7-1-15_Ortelius Oval_半球並列.pdf
PDFファイル 284.1 KB

  6-7-1-16 Conformal World in an Ellipse 図法

        【 投影式 】もし ks = sin θ  and  kc = cos θ なら

                   K' = ∫0 π/2 ( 1 - kc2sin2λ')-1/2 dλ'

                                      θ = 23.8958° のとき K' = 2.34767

                                      u1 = 2 (1- kc) cos φ

                   v1 = (1 + kc) [ l + cos φ cos (λ - λ0)]

                                      A = ks2 [ 1 - cos φ cos (λ - λ0) ]

                                      B = u1 - v1

                                      C = u1 + v1

                                      R1 = A - (A2- 4kcBC)1/2

                                      λ' = -arcsin [R1 / (2Ckc)]

                                     φ'=arcsin ({ 1- [R1 / (2B)]2}1/2 / ks)

                                      u2 = ∫0 φ' (1 - ks2sin2φ' )-1/2dφ'

                                      v2 = K'- ∫0 λ' ( l - kc2sin2λ' )-1/2 dλ'

                                      u3 = exp [πu2 / (4K')]

                                      v3 = πv2 / (4K')

                                           x = ± R (u3 + 1 / u3)(sin v3) / 2

                                           y = ± R (u3 - l / u3)(cos v3) / 2

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

              【 ポイント 】・全球を楕円形に投影したものです。緯線の間隔は低緯度において

                          細かくした方が、楕円形の横左右端の経線がスムーズに描けます。

                         ・投影式の中に積分が含まれていますが、シンプソンの公式を利用

                                 して簡易積分をして見ました。

                         ・簡易積分のせいか、補助線などでエラーなるポイントが発生しま

                           す。端っこの経線などで起こる可能性があります。

                           そのポイントは、経度か緯度を微調整してください。

                                 <例>  緯度  135°       ←エラー発生

                                         135.000003° ←エラー回避

                         さらに、簡易積分の影響か 経線の xx 軸との隙間ができるので、

                         x の値を補正する必要があります。

                          左右端の経線の緯度 ±90°の x 計算値と x 軸との

                           差が隙間の値となっています。

 【入力イメージ】

     x y を求めるための係数が多いので入力セルの配置などを画像にしました。

  <基準経度・標準緯度・予備計算・共通数値 など>                
     記号  セルNo.         記述式                
      λ0     C2     ???  ← 中心経度(希望経度を入力)                
       R      C3    1                
       θ     E2   23.8958                
       K'     E3    2.34767                
       ks     G2  =SIN(RADIANS($E$2))                
       ks2     G3  =G2^2                
       kc      I2  =COS(RADIANS($E$2))                
       kc2     I3  =I2^2                
                                                                         
  <係数計算・共通数値 など>
        記号  セルNo.         記述式                                                
        経度     C6   -45
        緯度     D6    90
        補正
     λ-λ
0
    E6  =IF($C$2<0,IF(C6-$C$2>=180,C6-$C$2-360,
                                     C6-$C$2),IF(C6-$C$2<=-180,C6+360-$C$2,C6-$C$2))
 
         u1     F6  =2*(1-$I$2)*COS(RADIANS(D6))
         v1     G6  =(1+$I$2)*(1+COS(RADIANS(D6))*COS(RADIANS(E6)))
         A     H6  =$G$3*(1-COS(RADIANS(D6))*COS(RADIANS(E6)))
         B     I6  =F6-G6
         C     J6  =F6+G6
         R1     K6  =H6-SQRT(H6^2-4*$I$2*I6*J6)
         λ'     L6  =-ASIN(K6/(2*J6*$I$2))
         φ'     M6  =ASIN(SQRT(1-(K6/(2*I6))^2)/$G$2)
         u2     N6  =(M6/6)*(1/SQRT(1-$G$3*(SIN(RADIANS(0)))^2)+4*
                            (1/SQRT((1-$G$3*(SIN(M6/2))^2)))+1/SQRT((1-$G$3*(SIN(M6))^2)))
 
         v2     O6  =$E$3-(L6/6)*(1/SQRT((1-$I$3*(SIN(0))^2))+4*
                             1/SQRT((1-$I$3*(SIN(L6/2))^2))+1/SQRT((1-$I$3*(SIN(L6))^2)))
 
         u3     P6  =EXP(PI()*N6/(4*$E$3))                                          
         v3     Q6  =PI()*O6/(4*$E$3)                                              
    x補正係数     R6  =SIGN(E6)*(P6+1/P6)*SIN(Q6)/2                                      
      R7  =-R6                                                      
                                                                         
  <データ計算 など>  経度緯度のデータ毎に記述式を複写(行コピー)する。
     記号  セルNo.         記述式
      No.     B11   ???  ← 式を複写する時に、一旦ソートしたり、解析時に役立つので連番号を付与
     経度     C11   ???  ← 実際の経度データ
     緯度     D11   ???  ← 実際の緯度データ
     補正
  λ-λ
0
    E11  =IF($C$2<0,IF(C11-$C$2>180,C11-$C$2-360,
                                     C11-$C$2),IF(C11-$C$2<-180,C11+360-$C$2,C11-$C$2))
 
       u1     F11  =2*(1-$I$2)*COS(RADIANS(D11))
       v1     G11  =(1+$I$2)*(1+COS(RADIANS(D11))*COS(RADIANS(E11)))
       A     H11  =$G$3*(1-COS(RADIANS(D11))*COS(RADIANS(E11)))
       B     I11  =F11-G11
       C     J11  =F11+G11
       R1     K11  =H11-SQRT(H11^2-4*$I$2*I11*J11)
       λ'     L11  =-ASIN(K11/(2*J11*$I$2))
       φ'     M11  =ASIN(SQRT(1-(K11/(2*I11))^2)/$G$2)
       u2     N11  =(M11/6)*(1/SQRT(1-$G$3*(SIN(RADIANS(0)))^2)+4*
                         (1/SQRT((1-$G$3*(SIN(M11/2))^2)))+1/SQRT((1-$G$3*(SIN(M11))^2)))
 
       v2     O11  =$E$3-(L11/6)*(1/SQRT((1-$I$3*(SIN(0))^2))+4*
                              1/SQRT((1-$I$3*(SIN(L11/2))^2))+1/SQRT((1-$I$3*(SIN(L11))^2)))
 
       u3     P11  =EXP(PI()*N11/(4*$E$3))
       v3     Q11  =PI()*O11/(4*$E$3)
        x     R11  =SIGN(E11)*(P11+1/P11)*SIN(Q11)/2
     x補正     S11  =IF(R11=0,0,IF(R11<0,R11-$R$7,R11-$R$6))
        y     T11  =SIGN(D11)*(P11-1/P11)*COS(Q11)/2

 

 [6-7-1-16]  Conformal World in an Ellipse 図法

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

ダウンロード
6-7-1-16_正積楕円.pdf
PDFファイル 313.2 KB