Площадь поверхности четырехугольника широты-долготы
area = areaquad(lat1,lon1,lat2,lon2)
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid)
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid,units)
area = areaquad(lat1,lon1,lat2,lon2) возвращает площадь поверхности, ограниченную параллелями lat1 и lat2 и меридианы lon1 и lon2. Продукция area - доля площади единичной сферы, равная 4δ, поэтому результат находится в диапазоне от 0 до 1.
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid) позволяет специфицировать эллипсоидную модель с ellipsoid. ellipsoid является referenceSphere, referenceEllipsoid, или oblateSpheroid объект или вектор формы [semimajor_axis eccentricity]. Когда ellipsoid является входным, результирующий area задаётся в терминах (квадратичных) единиц эллипсоида. Например, если эллипсоид referenceEllipsoid('grs80','kilometers') используется, результирующая площадь в км2.
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid,units) где units определяет единицы ввода. Значение по умолчанию: 'degrees'.
Найдите долю поверхности Земли, лежащую между 30ºN и 45ºN, а также между 25ºW и 60ºE:
area = areaquad(30,-25,45,60)
area =
0.0245Предполагая сферический эллипсоид, найдите площадь поверхности Земли в квадратных километрах.
earthellipsoid = referenceSphere('earth','km');
area = areaquad(-90,-180,90,180,earthellipsoid)
area =
5.1006e+08Для сравнения:
earthellipsoid.SurfaceArea ans = 5.1006e+08
areaquad точный расчет, основанный на простой сферической геометрии. Для несферических эллипсоидов данные преобразуются во вспомогательную авторитарную сферу.