Surface area of latitude-longitude quadrangle
area = areaquad(lat1,lon1,lat2,lon2)
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid)
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid,
area = areaquad(lat1,lon1,lat2,lon2) returns
the surface area bounded by the parallels
area is a fraction of the unit sphere's
area of 4π, so the result ranges from 0 to 1.
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid) allows
the specification of the ellipsoid model with
ellipsoid is a
oblateSpheroid object, or a vector of the form
ellipsoid is input, the resulting
area is given in terms of the (squared) units of the ellipsoid. For
example, if the ellipsoid
used, the resulting area is in km2.
area = areaquad(lat1,lon1,lat2,lon2,ellipsoid, specifies
the units of the inputs. The default is
Find the fraction of the Earth's surface that lies between 30ºN and 45ºN, and also between 25ºW and 60ºE:
area = areaquad(30,-25,45,60) area = 0.0245
Assuming a spherical ellipsoid, find the surface area of the Earth in square kilometers.
earthellipsoid = referenceSphere('earth','km'); area = areaquad(-90,-180,90,180,earthellipsoid) area = 5.1006e+08
earthellipsoid.SurfaceArea ans = 5.1006e+08
A latitude-longitude quadrangle is a region bounded by two meridians and two parallels. In spherical geometry, it is the intersection of a lune (a section bounded by two meridians) and a zone (a section bounded by two parallels).
areaquad calculation is exact, being
based on simple spherical geometry. For nonspherical ellipsoids, the
data is converted to the auxiliary authalic sphere.