В этом примере показано, как вычислить географические области для векторных данных в формате многоугольника с помощью 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