Obtain Ranges of Values of the Earth's Magnetic Field

Define a region near the North Pole from its latitudes and longitudes:

In[1]:=

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

Plot this region on a world map:

In[2]:=

plot=GeoGraphicsGeoBoundsRegion[{lat,lon}],



Out[2]=

Get data about the Earth's magnetic field in this region:

In[3]:=

data=GeomagneticModelDataTranspose[{lat,lon}],

...



Out[3]=

Extract the components of the Earth's magnetic field:

In[4]:=

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

Out[4]=

QuantityArray

Dimensions: {40,94}

Unit: Nanoteslas

,QuantityArray

Dimensions: {40,94}

Unit: Nanoteslas

,QuantityArray

Dimensions: {40,94}

Unit: Nanoteslas



Convert the list of arrays into an array of 3-vectors:

In[5]:=

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

Out[5]=

QuantityArray

Dimensions: {40,94,3}

Unit: Nanoteslas



Use

CoordinateBounds

to obtain the ranges of the components of the Earth's magnetic field in the specified region:

In[6]:=

Thread[componentsCoordinateBounds[qa]]//TableForm

Out[6]//TableForm=

NorthComponent{

1823.47

nT

,

14720.3

nT

}

EastComponent{

-3963.88

nT

,

3266.14

nT

}

DownComponent{

49264.4

nT

,

56798.8

nT

}

Related Symbols

Publisher Information