Scribd
Upload
62 views3 pages

Geodetik2UTM Transformation

This document contains a script to transform geographic coordinates to UTM coordinates. It provides the name, student number, subject, and home point coordinates in degrees. The script defines ellipsoid parameters, converts to radians, determines UTM zone and applies formulas to calculate UTM coordinates of X = 680189.9655 m and Y = 9192541.1147 m.

Uploaded by

rizka
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
62 views3 pages

Geodetik2UTM Transformation

This document contains a script to transform geographic coordinates to UTM coordinates. It provides the name, student number, subject, and home point coordinates in degrees. The script defines ellipsoid parameters, converts to radians, determines UTM zone and applies formulas to calculate UTM coordinates of X = 680189.9655 m and Y = 9192541.1147 m.

Uploaded by

rizka
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

GEOGRAPHIC TO UTM COORDINATE TRANSFORMATION

Name : Rizka Amelia Dwi Safira


NRP : 03311940000044
Subject : Coordinate System and Map Projection

Home Point Coordinate:

(Lakarsantri Gg. II No.8A, Lakarsantri, Surabaya)


7°18' 7.15'' LU;112°37' 56.13'' LS

Script:
clear all
clc
format longG

disp("=====Geographic to UTM Coordinate Transformation=====")


disp("Name : Rizka Amelia Dwi Safira")
disp("NRP : 03311940000044")
disp("Subject : Coordinat System and Map Projection")
disp("=====================================================")

% Reference ellipsoid parameters: WGS1984


a = 6378137; %semi-major axis
f = 1/298.257223563; %flattening
b = a-a*f; %semi-minor axis
e1 = f*(2-f); %square eccentricity-1 (e^2)
e2 = (a^2-b^2)/b^2; %square eccentricity-2 (e'^2)

% Geographic coordinate
Lat = -(7+(18/60)+(7.15/3600)); %latitude(deg) of home point
Latrad = Lat*pi/180; %latitude(rad) of home point
Lon = 112+(37/60)+(56.13/3600); %longitude(deg) of home point
disp("Known:")
disp(['Latitude of Home Point = ' , num2str(Lat), ' degree'])
disp(['Longitude of Home Point = ' num2str(Lon), ' degree'])

% Determining l and b
% Central meridian (Long0)
Long0 = (floor(Lon/6)*6) + 3;
dLong0 = Lon-Long0;
% Central paralel (Lat0)
Lat0 = Lat;

% Scale for UTM transformation


k01 = 0.9996;

% Deviation parameter (East-West)


miu_k = e2*cosd(Lat)*cosd(Lat); %square deviation parameter
p = 180/pi; %constanta

% Ellipsoid normal radius (N)


X = sqrt(1-e1*(sind(Lat))^2);
N = a/X; %ellipsoid normal radius

% Meridian Length (G)


E0 = a*(1-e1/4-3*e1^2/64-5*e1^3/256);
E2 = 3*a*(e1+e1^2/4+15*e1^3/128)/8;
E4 = a*15*(e1^2+3*e1^3/4)/256;
E6 = a*35*e1^3/3072;
G = E0*Latrad-E2*sind(2*Lat)+E4*sind(4*Lat) -E6*sind(6*Lat);

% Geograpic coordinate to UTM


a0 = k01*G;
a1 = k01*N*cosd(Lat)/p;
a2 = k01*N*sind(Lat)*cosd(Lat)/(2*p^2);
a3 = k01*N*cosd(Lat)*cosd(Lat)*cosd(Lat)*(1 -(tand(Lat))^2+miu_k)/(6*p^3);
a4 = k01*N*sind(Lat)*cosd(Lat)*cosd(Lat)*cosd(Lat)*(5-
(tand(Lat))^2+9*miu_k+4*miu_k^2)/(24*p^4);
a5 = k01*N*(cosd(Lat))^5*(5-(18*(tand(Lat))^2)+(tand(Lat))^4+(14*miu_k)-
(58*e2*(sind(Lat))^2))/(120*p^5);
a6 = k01*N*sind(Lat)*(cosd(Lat))^5*(61-58*(tand(Lat))^2+(tand(Lat))^4+270*miu_k -
330*e2*(sind(Lat))^2)/(270*p^6);

% UTM Coordinate
x = (a1*dLong0 + a3*dLong0^3 + a5*dLong0^5) + 500000;
y = (a0+a2*dLong0^2+a4*dLong0^4+a6*dLong0^6) + 10000000;

disp(' ')
disp('Calculation Result:')
disp(['X UTM = ', num2str(x), ' m'])
disp(['Y UTM = ', num2str(y), ' m'])

disp(' ')
disp('Conclusion:')
disp(['Home point in geographic coordinate (degree): ', num2str(Lat), ' and ',
num2str(Lon), ' has UTM coordinate: X = ',...
num2str(x), ' m and Y = ' num2str(y), ' m'])

disp(' ')
disp('Thank you')
Result:
=====Geographic to UTM Coordinate Transformation=====
Name : Rizka Amelia Dwi Safira
NRP : 03311940000044
Subject : Coordinat System and Map Projection
=====================================================
Known:
Latitude of Home Point = -7.302 degree
Longitude of Home Point = 11 2.6323 degree

Calculation Result:
X UTM = 680189.9655 m
Y UTM = 9192541.1147 m

Conclusion:
Home point in geographic coordinate (degree): -7.302 and 112.6323 has UTM
coordinate: X = 680189.9655 m and Y = 9192541.1147 m

Thank you

Computation from SRGI Converter:

You might also like