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