地球磁场范围值

全新函数 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]=