We help IT Professionals succeed at work.

Check out our new AWS podcast with Certified Expert, Phil Phillips! Listen to "How to Execute a Seamless AWS Migration" on EE or on your favorite podcast platform. Listen Now

x

Need to convert UTM to latitude / longitude

RUSS_EMI
RUSS_EMI asked
on
Medium Priority
1,164 Views
Last Modified: 2007-12-19
Need to convert UTM to Lat / Long. It's close, but not right. This is needed VERY ugently.
Thanks,
Russ


'********************************************************************
a = 6378137
b = 6356752.3142
e12 = 6.69437999014138E-03
e = 8.18191908426219E-02
e1 = 1.6792203863838E-03

M = Northing / 0.9996
mu = M / (a * (1 - e12 / 4 - 3 * e ^ 4 / 64 - 5 * e ^ 6 / 256))



J1 = 3 * e1 / 2 - 27 * e1 ^ 3 / 32 '
J2 = 21 * Sqr(e1) / 16 - (55 * e1) ^ 4 / 32 '
J3 = (151 * e1) ^ 3 / 96 '
J4 = (1097 * e1) ^ 4 / 512 '
fp = mu + J1 * Sin(2 * mu) + J2 * Sin(4 * mu) + J3 * Sin(6 * mu) + J4 * Sin(8 * mu) '
C1 = e12 * Sqr(Cos(fp)) '
T1 = Sqr(Tan(fp)) '
R1 = a * (1 - Sqr(e)) / (1 - Sqr(e) * Sqr(Sin(fp))) ^ 1.5 '
N1 = a / (1 - Sqr(e) * Sqr(Sin(fp))) ^ 0.5 '
D = Easting / (N1 * 0.9996) '
Q1 = (N1 * Tan(fp)) / R1 '
Q2 = (Sqr(D) / 2) '
Q3 = (5 + 3 * T1 + 10 * C1 - 4 * Sqr(C1) - 9 * e12) * D ^ 4 / 24
Q4 = (61 + 90 * T1 + 298 * C1 + 45 * Sqr(T1) - 3 * Sqr(C1) - 252 * e12) * D ^ 6 / 720
Q5 = D
Q6 = (1 + 2 * T1 + C1) * D ^ 3 / 6
Q7 = (5 - 2 * C1 + 28 * T1 - 3 * Sqr(C1) + 8 * e12 + 24 * Sqr(T1)) * D ^ 5 / 120

Latitude = fp - Q1 * (Q2 - Q3 + Q4) '
Longitude = Central_Meridian + (Q5 - Q6 + Q7) / Cos(fp) '

Longitude = Longitude * 180 / PI
Latitude = Latitude * 180 / PI
Easting = Easting
Northing = Northing
Comment
Watch Question

g0arth, unfortunately, that code does EXACTLY the inverse of what RUSS EMI was asking for.  He has UTM and needs Latitude/Longitude.  The provided link returns UTM from Latitude/Longitude.

AW

Author

Commented:
Found the answer myself.  A square root term should have been squared.  Please close this question and refund the points.
Russ
Russ, you need to post a zero-point request to Community Support, asking them to delete this question, and refund your points.

AW

Author

Commented:
Thank you, AW

Author

Commented:
Thank you.
Commented:
Unlock this solution and get a sample of our free trial.
(No credit card required)
UNLOCK SOLUTION

Commented:
Hi Russ

could you post your final code? Also, I assume that your UTM paramters of a,b,e are specific to your 'zone'. Can you elaborate how I can use for my 'zone', here in Australia ?
Unlock the solution to this question.
Thanks for using Experts Exchange.

Please provide your email to receive a sample view!

*This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.

OR

Please enter a first name

Please enter a last name

8+ characters (letters, numbers, and a symbol)

By clicking, you agree to the Terms of Use and Privacy Policy.