[Webinar] Streamline your web hosting managementRegister Today

x
  • Status: Solved
  • Priority: Medium
  • Security: Public
  • Views: 957
  • Last Modified:

Need to convert UTM to latitude / longitude

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
0
RUSS_EMI
Asked:
RUSS_EMI
1 Solution
 
Arthur_WoodCommented:
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

0
 
RUSS_EMIAuthor Commented:
Found the answer myself.  A square root term should have been squared.  Please close this question and refund the points.
Russ
0
Never miss a deadline with monday.com

The revolutionary project management tool is here!   Plan visually with a single glance and make sure your projects get done.

 
Arthur_WoodCommented:
Russ, you need to post a zero-point request to Community Support, asking them to delete this question, and refund your points.

AW
0
 
RUSS_EMIAuthor Commented:
Thank you, AW
0
 
RUSS_EMIAuthor Commented:
Thank you.
0
 
moduloCommented:
PAQed, with points refunded (500)

modulo
Community Support Moderator
0
 
LukeBCommented:
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 ?
0

Featured Post

Free Tool: Subnet Calculator

The subnet calculator helps you design networks by taking an IP address and network mask and returning information such as network, broadcast address, and host range.

One of a set of tools we're offering as a way of saying thank you for being a part of the community.

Tackle projects and never again get stuck behind a technical roadblock.
Join Now