Диапазон значений магнитного поля Земли
Новые функции 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]=
Переведите список массивов в массив трёх векторов.
In[5]:=
qa = Transpose[QuantityArray[arr], {3, 1, 2}]Out[5]=
Далее представлены диапазоны компонентов магнитного поля Земли в заданной области.
In[6]:=
Thread[components -> CoordinateBounds[qa]]Out[6]=
