GEODETIC2NED convert from geodetic to NED coordinates

Inputs
outputs

Inputs

lat,lon, alt: ellipsoid geodetic coordinates of point under test (degrees, degrees, meters)
lat0, lon0, alt0: ellipsoid geodetic coordinates of observer/reference (degrees, degrees, meters)
spheroid: referenceEllipsoid
angleUnit: string for angular units. Default 'd': degrees

outputs

n,e,d: North, East, Down coordinates of test points (meters)

function [north, east, down] = geodetic2ned(lat, lon, alt, lat0, lon0, alt0, spheroid, angleUnit)
arguments
  lat {mustBeReal}
  lon {mustBeReal}
  alt {mustBeReal}
  lat0 {mustBeReal}
  lon0 {mustBeReal}
  alt0 {mustBeReal}
  spheroid (1,1) matmap3d.referenceEllipsoid = matmap3d.wgs84Ellipsoid()
  angleUnit {mustBeTextScalar} = 'd'
end

[x1,y1,z1] = matmap3d.geodetic2ecef(spheroid, lat,lon,alt,angleUnit);
[x2,y2,z2] = matmap3d.geodetic2ecef(spheroid, lat0,lon0,alt0,angleUnit);

dx = x1-x2;
dy = y1-y2;
dz = z1-z2;

[east, north, up] = matmap3d.ecef2enuv(dx, dy, dz, lat0, lon0, angleUnit);

down = -up;

end

GEODETIC2NED convert from geodetic to NED coordinates

Contents

Inputs

outputs