How would i convert latitude and longitude to meters from a point?

I have lots of data points (over 50000) of measured positions from a gps antenna where i have the latitude and longitude of every point. What I would like to do is to plot the points relative to a reference points where the scale of the axis are in meters. I would assume that since the area of the points is very small (approx. 10-20 meters) it is possible to assume that the surface is not a part of a sphere but is a flat surface. Please help me on how to plot this.

Preferably I would like the plot to look something like this:

http://www.google.co.uk/imgres?start=127&um=1&hl=sv&sa=N&authuser=0&biw=1366&bih=631&tbm=isch&tbnid=MTkiipN3vyE-5M:&imgrefurl=http://gps.ece.cornell.edu/scintplot_aopr.html&docid=_IpDd14QuOtYBM&imgurl=http://gps.ece.cornell.edu/aopr_elazA.png&w=256&h=256&ei=pjMIUJPxMZSKhQfm8O3mAw&zoom=1&iact=rc&dur=220&sig=118366609262092269785&page=6&tbnh=143&tbnw=143&ndsp=24&ved=1t:429,r:21,s:127,i:241&tx=88&ty=63

with meters as the scales on the axis. Atm the only thing I have achieved is scatter(ChipLatLong(:,1),ChipLatLong(:,2))

Which is just a plot of the latitude and longitude on a regular graph.

I would be very thankful for any response

Viktor

Answer 1

[Meysam Mahooti]

%--------------------------------------------------------------------------

%

% Position: Position vector (r [m]) from geodetic coordinates

% (Longitude [rad], latitude [rad], altitude [m])

%

% Last modified: 2018/01/27 M. Mahooti

%

%--------------------------------------------------------------------------

function r = Position(lon, lat, h)

R_equ = 6378.137e3; % Earth's radius [m]; WGS-84

f = 1/298.257223563; % Flattening; WGS-84

e2 = f*(2-f); % Square of eccentricity

CosLat = cos(lat); % (Co)sine of geodetic latitude

SinLat = sin(lat);

% Position vector

N = R_equ/sqrt(1-e2*SinLat*SinLat);

r(1) = (N+h)*CosLat*cos(lon);

r(2) = (N+h)*CosLat*sin(lon);

r(3) = ((1-e2)*N+h)*SinLat;