地球磁場の値の範囲

新関数のCoordinateBoundsとCoordinateBoundingBoxを使うと，多次元配列の各次元の極値を求めることができる．この例では，指定された地域における地球磁場の値の範囲を求める．

北極近くの地域を経緯度で定義する．

In[1]:=

lat = {65, 85}; lon = {-30, 30};

世界地図上にこの地域をプロットする．

In[2]:=

GeoGraphics[GeoBoundsRegion[{lat, lon}], GeoGridLines -> Automatic, GeoZoomLevel -> 5]

Out[2]=

この地域における地球磁場のデータを生成する．

In[3]:=

data = GeomagneticModelData[Transpose[{lat, lon}], GeoZoomLevel -> 1]

Out[3]=

このデータから地球磁場の成分を取り出す．

In[4]:=

components = {"NorthComponent", "EastComponent", "DownComponent"}; arr = Lookup[data, components]

Out[4]=

配列のリストを3次元ベクトルの配列に変換する．

In[5]:=

qa = Transpose[QuantityArray[arr], {3, 1, 2}]

Out[5]=

これが，指定地域における地球磁場の成分範囲である．

In[6]:=

Thread[components -> CoordinateBounds[qa]]

Out[6]=