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