Determining UTM zone (to convert) from longitude/latitude

by

Edit: For (non-R) code that works for all non-polar areas on earth, see here or here.

Unless you are dealing with data from a couple of exceptional areas (Svalbard and parts of Norway), this is a simple enough calculation that you might as well just do it yourself in R. Here is Wikipedia’s description of how longitude relates to UTM Zone number:

The UTM system divides the surface of Earth between 80°S and 84°N latitude into 60 zones, each 6° of longitude in width. Zone 1 covers longitude 180° to 174° W; zone numbering increases eastward to zone 60 that covers longitude 174 to 180 East.

So, assuming that in your data longitudes to the west of the Prime Meridian are encoded as running from -180 to 0 degrees, here’s an R-code version of the above:

long2UTM <- function(long) {
    (floor((long + 180)/6) %% 60) + 1
}

# Trying it out for San Francisco, clearly in UTM Zone 10 
# in the figure in the Wikipedia article linked above
SFlong <- -122.4192
long2UTM(SFlong)
# [1] 10

That expression could obviously be simplified a bit, but I think in this form the logic underlying its construction is most clear. The %% 60 bit is in there just in case some of your longitudes are greater than 180 or less than -180.

More Related Contents:

  1. Geographic / geospatial distance between 2 lists of lat/lon points (coordinates)
  2. R – Finding closest neighboring point and number of neighbors within a given radius, coordinates lat-long
  3. Latitude Longitude Coordinates to State Code in R
  4. Function to calculate geospatial distance between two points (lat,long) using R [duplicate]
  5. Cannot install R-forge package using install.packages
  6. Change the Blank Cells to “NA”
  7. Create a group number for each consecutive sequence
  8. Proper/fastest way to reshape a data.table
  9. Global variables in R
  10. Cumulative sum that resets when 0 is encountered
  11. How can I remove empty factors from ggplot2 facets?
  12. Efficiently sum across multiple columns in R
  13. Error in grid.Call(L_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : Polygon edge not found
  14. What can R do about a messy data format?
  15. Remove duplicates keeping entry with largest absolute value
  16. Plot data in descending order as appears in data frame [duplicate]
  17. R ifelse statement
  18. R knitr: Possible to programmatically modify chunk labels?
  19. Switch displayed traces via plotly dropdown menu
  20. R grep: is there an AND operator?
  21. Split up a dataframe by number of rows
  22. stat_contour not able to generate contour lines
  23. Make R exit with non-zero status code
  24. Add x and y axis to all facet_wrap
  25. display a matrix, including the values, as a heatmap
  26. How to automatically adjust the width of each facet for facet_wrap?
  27. Calculate mean across rows with NA values in R
  28. Is it possible to jitter two ggplot geoms in the same way?
  29. Convert numeric to date
  30. How to combine two lists in R

Leave a Comment