В этом примере показано, как вычислить географические области для векторных данных в формате многоугольника с помощью 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 возвращает три области: самую большую для многоугольника, представляющего Superior, Мичиган и Гурон вместе, две другие для Эри и Онтарио.
earthradius = almanac('earth','radius','miles'); area = areaint(gtlakelat,gtlakelon,earthradius)
area = 3×1
104 ×
8.0120
1.0382
0.7634