Take the 2-minute tour ×
Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

How to use the UTM coordinate system (northing, easting) in Mathematica using the built-in functions (GeoPositionENU, GeoGridPosition, etc.). There are UTM Zones implemented in Mathematica, e.g.:

GeoProjectionData["UTMZone33"]

gives

{"TransverseMercator", {"Centering" -> {0, 15},   "CentralScaleFactor" -> 1, "GridOrigin" -> {0, 0},   "ReferenceModel" -> "WGS84"}}

From N45, E15 one should get Easting 500000 and Northing 4982950.4 (zone 33), but I can't reproduce it.

GeoGridPosition[GeoPosition[{45, 15, 0}, "WGS84"], "UTMZone33"]

gives

GeoGridPosition[{0., 4.98494*10^6, 0}, "UTMZone33"]

(Easting offset 500000 is understandable - default, but Northing is way off.)

Any help would be appreciated.

share|improve this question
    
I've wanted to know this for a while too. Right now, I just use the function included in this blog post to go from LatLon to UTM: blog.wolfram.com/2009/04/17/mapping-gps-data –  kale Jan 8 '13 at 14:56
2  
The output of GeoProjectionData is not correct for the UTM system: the value of GridOrigin is wrong (it should be {500000,0} and so is the value of CentralScaleFactor (it should be 0.9996). –  whuber Jan 8 '13 at 17:04
    
@kale: the short version of the function would therefore be: ll2utm[coord_] := GeoGridPosition[ GeoPosition[coord, "WGS84Original"], {"UTMZone33", "CentralScaleFactor" -> 0.9996, "GridOrigin" -> {500000, 0}}] [[1, {1, 2}]] and then ll2utm[{45, 15}] // Round gives {500000, 4982950}. –  Boocko Jan 9 '13 at 11:08
    
@Boocko, Except, you would need to include the logic to figure out what zone you were in and build a string "UTMZone"<>ToString@zone. –  kale Jan 9 '13 at 15:19
add comment

1 Answer

up vote 4 down vote accepted

As in comments by WReach, the correct answer is

GeoGridPosition[ 
  GeoPosition[{45, 15, 0}, "WGS84"], 
  {"UTMZone33", "CentralScaleFactor" -> 0.9996, "GridOrigin" -> {500000, 0}}
]

The problem lies with the scale factor: 0.9996 and the grid origin

If you are familiar with projection systems used in GIS systems, you can check out http://spatialreference.org/ref/epsg/32633/

Details are here:

PROJCS["WGS 84 / UTM zone 33N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",15],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    AUTHORITY["EPSG","32633"],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]

More explanations of the UTM system can be found here: "The central meridian in each UTM zone has a scale factor of 0.9996"

share|improve this answer
1  
I feel this falls a bit short as an answer. Could you elaborate somewhat? –  Sjoerd C. de Vries Jan 8 '13 at 16:29
1  
I agree with your diagnosis, although the easting should also be 500000 as the OP observes. You might want to mention that the desired result can be obtained using the expression GeoGridPosition[ GeoPosition[{45, 15, 0}, "WGS84"], {"UTMZone33", "CentralScaleFactor" -> 0.9996, "GridOrigin" -> {500000, 0}}]. –  WReach Jan 8 '13 at 19:18
    
Thank you, this works. But to convert back, the simple GeoPosition on previous output should work, but I get errors. Looks like the GeoPosition is not happy with {"UTMZone33",...}. –  Boocko Jan 9 '13 at 10:51
1  
works in this way: GeoProjectionData["UTMZone33"] gives {"TransverseMercator", {"Centering" -> {0, 15}, "CentralScaleFactor" -> 1, "GridOrigin" -> {0, 0}, "ReferenceModel" -> "WGS84"}}, then use this: GeoPosition@ GeoGridPosition[ GeoPosition[{45, 15, 0}, "WGS84"], {"TransverseMercator", "Centering" -> {0, 15}, "CentralScaleFactor" -> 0.9996, "GridOrigin" -> {500000, 0}, "ReferenceModel" -> "WGS84"}] –  Tuku Jan 9 '13 at 15:30
    
@Tuku that works flawlessly. Thanks. –  Boocko Jan 10 '13 at 10:27
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.