Этот пример показывает, как вычислить географические области для векторных данных в формате полигона с помощью функции areaint. areaint выполняет численное интегрирование с помощью Теоремы Зеленого для области на поверхности, заключенной полигоном. Поскольку это - дискретное интегрирование на дискретных данных, результаты не точны. Тем не менее, метод обеспечивает лучшие средние значения вычисления областей областей произвольной формы. Лучшие меры следуют из лучших данных. Для получения дополнительной информации смотрите areaint.
Загрузите MAT-файл континентальной части США, conus.mat, и вычислите радиус Земли.
load conus earthradius = almanac('earth','radius');
Вычислите область континентальной части США, наряду с областью Лонг-Айленда и Виноградника Марты. areaint как другие функции области, areaquad и areamat, возвращает область как часть поверхности целой планеты, если вы не обеспечиваете радиус. Поскольку Наземный радиус по умолчанию находится в километрах, область находится в квадратных километрах.
area = areaint(uslat,uslon,earthradius)
area = 3×1
106 ×
7.9256
0.0035
0.0004
Вычислите области Великих озер с помощью тех же переменных, на этот раз в квадратных милях. areaint возвращает три области: самое большое для полигона, представляющего Выше, Мичиган, и Гурон вместе, другие два для Эри и Онтарио.
earthradius = almanac('earth','radius','miles'); area = areaint(gtlakelat,gtlakelon,earthradius)
area = 3×1
104 ×
8.0120
1.0382
0.7634