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