Площадь поверхности четырехугольника долготы широты
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 сферу.