Azimuth¶

Calculates the great circle azimuth from a point 1 to a point 2. The latitude and longitude of these two points can either be given independently or as columns of the matrices point 1 and point 2 in the form [latitude longitude].

>> azimuth([10, 10], [10, 40])
ans = 87.336
>> azimuth([0, 10], [0, 40])
ans = 90
>> azimuth(pi / 4, 0, pi / 4, -pi / 2, "radians")
ans = 5.3279

Reckon¶

Compute the coordinates of the end-point of a displacement on a sphere. lat e lon are the coordinates of the starting point, range is the covered distance of the displacements along a great circle and azimuth is the direction of the displacement relative to the North.

This function can also be used to define a spherical coordinate system with rotated poles.

>> [lato, lono] = reckon(0, 10, 30, 90)
lato = 0
lono =  40.000
>> [lato, lono] = reckon(0, 0, pi, pi / 2, "radians")
lato = 0
lono = -3.1416

Distance¶

Calculates the distance (in degrees) between point 1 and point 2.

>> distance([37,-76], [37,-9])
ans = 52.309
>> distance([37,-76], [67,-76])
ans = 30.000

Departure¶

Calculates the distance in degrees between longitudes lon1 and lon2 at latitude defined by lat.

>> departure(0, 10, 0)
ans = 10
>> departure(0, 10, 60)
ans = 5

Vincenty¶

Calculates the distance (in kilometers) between point 1 and point 2 using the formula devised by Thaddeus Vincenty, with an accurate ellipsoidal model of the earth. The default ellipsoidal model is WGS-84, which is the most globally accurate model.

>> vincenty([37, -76], [37, -9])
ans = 5830.081
>> vincenty([37, -76], [67, -76], referenceEllipsoid(7019))
ans = 3337.843

Vincenty (Direct Form)¶

Compute the coordinates of the end-point of a displacement on a geodesic. The parameters lat and lon are the coordinates of the starting point, The parameter range is the covered distance of the displacements along a specified geodesic and azi is the direction of the displacement relative to the North.

>> [lat, lon] = vincentyDirect(0, 0, pi, pi / 2)
lat = 0
lon = 3.1311
>> [lat, lon, az] = vincentyDirect(0, 0, pi, pi / 2)
lat = 0
lon = 3.1311
az = 1.5708

Great-Circle to Small-Circle¶

Converts a great circle to small circle notation.

Degree to Kilometer¶

Convert a distance along a great circle of the Earth from degrees to kilometer. Great-circle distance uses a spherical model of the earth, using the average great-circle radius of 6372.795 kilometers, resulting in an error of up to about 0.5%.

Kilometer to Degree¶

Convert a distance along a great circle of the Earth from kilometer to degrees. A radius Great-circle distance uses a spherical model of the earth, using the average great-circle radius of 6372.795 kilometers, resulting in an error of up to about 0.5%.

Degree to Radians¶

Converts angles input in degrees to the equivalent in radians.

Radians to Degree¶

Converts angles input in radians to the equivalent in degrees.

Geocode¶

Geolocation service from Google Maps API. Define a geocode method for resolving a location from a string.

Geocoding is the process of converting addresses (like “1600 Amphitheatre Parkway, Mountain View, CA”) into geographic coordinates (like latitude 37.423021 and longitude -122.083739), which you can use to place markers or position the map.

Reverse geocoding is the process of converting geographic coordinates into a human-readable address.

>> [lat, lon] = geocode('new york,ny')
lat = 40.714
lon = -74.006
>> [lat, lon] = geocode('brasilia,df')
lat = -15.827
lon = -47.922

Reverse¶

Geolocation service from Google Maps API. Define a geocode method for resolving a address from a location.

The term geocoding generally refers to translating a human-readable address into a location on a map. The process of doing the converse, translating a location on the map into a human-readable address, is known as reverse geocoding.

>> reverse([40.714 -74.006])
ans = '58-68 Chambers Street, Nova York, NY 10007, EUA'
>> reverse([-15.827 -47.922])
ans = 'Asa Sul, Brasília, DF, 70383-070, Brasil'

nPI to PI¶

This brings the angin into the \([-pi pi]\) interval.

ZERO to 2PI¶

This brings the angin into the \([0 2*pi[\) interval.

Reference Ellipsoid¶

This function returns a reference ellipsoid object corresponding to the specified code (numerical EPSG). The values of the SemimajorAxis and SemiminorAxis properties are in kilometers. The reference ellipsoid has five properties: Code, Name, SemimajorAxis, SemiminorAxis and Flattening.

The form code can receive a valid EPSG code. 46 codes are currently implemented between 7001 and 7053 (except for 7017, 7023, 7026 and 7037-7040).

The valid values for name form are as follows: sphere, unitsphere, earth, moon, mercury, venus, mars, jupiter, saturn, uranus, neptune and pluto.