Skip to content
Advertisement

degree decimal to utm conversion, python

I got surprised when I tried to convert degree decimal coordinates into utm by using available library in anacondautm.

my coordinates looks like:

lon = array([83.71666667, 83.7       , 83.6       , 83.65      , 83.6       ,
       83.88333333, 83.56666667, 83.88333333, 83.96666667, 83.75      ,
       83.7       , 83.1       , 83.21666667, 83.73333333, 83.65      ,
       83.4       , 83.56666667, 84.        , 83.78333333, 83.68333333,
       83.6       , 83.48333333, 83.3       , 83.38333333, 83.61666667,
       83.98333333, 83.43333333, 83.53333333, 84.05      , 84.21666667,
       83.15      , 83.06666667, 83.26666667, 83.25      , 83.8       ,
       84.9       , 84.36666667, 84.        , 83.88333333, 84.61666667,
       84.35      , 84.41666667, 84.61666667, 83.81666667, 84.11666667,
       83.81666667, 83.8       , 84.1       , 84.23333333, 84.28333333,
       83.96666667, 84.01666667, 83.8       , 84.61666667, 84.1       ,
       83.76666667, 84.13333333, 83.75      , 83.78333333, 83.91666667,
       85.        , 84.28333333, 84.41666667, 84.53333333, 85.13333333,
       85.05      , 85.01666667, 84.81666667, 84.98333333, 84.43333333,
       85.38333333, 84.81666667, 85.01666667, 84.93333333, 85.25      ,
       85.18333333, 85.31666667, 85.3       , 85.11666667, 85.41666667,
       85.25      , 85.54722   ])

and lat values are:

lat = array([28.78333333, 28.75      , 28.26666667, 28.48333333, 28.63333333,
       28.81666667, 28.35      , 29.05      , 29.18333333, 28.18333333,
       28.21666667, 28.4       , 28.6       , 28.4       , 28.03333333,
       28.38333333, 28.15      , 29.1       , 28.96666667, 28.9       ,
       28.46666667, 28.38333333, 28.05      , 28.56666667, 28.13333333,
       29.18333333, 27.95      , 27.86666667, 27.68333333, 27.68333333,
       27.93333333, 27.55      , 28.01666667, 28.06666667, 27.86666667,
       28.36666667, 28.28333333, 28.21666667, 28.1       , 28.06666667,
       28.13333333, 27.93333333, 28.        , 27.88333333, 28.11666667,
       28.26666667, 28.3       , 28.03333333, 28.55      , 27.96666667,
       28.26666667, 28.06666667, 28.38333333, 28.2       , 28.36666667,
       27.98333333, 27.86666667, 28.26666667, 28.26666667, 28.08333333,
       28.48333333, 28.76666667, 27.61666667, 27.58333333, 27.55      ,
       27.41666667, 27.41666667, 27.55      , 27.43333333, 27.06666667,
       28.28333333, 28.05      , 27.91666667, 27.86666667, 27.8       ,
       27.71666667, 28.01666667, 28.1       , 28.01666667, 27.75      ,
       27.75      , 28.20946   ])

when I made scatter plot it looks like normal:

enter image description here

After converting coordinates in to utm distribution of the points change completely and looks strange like in this figure:

enter image description here

the code that i used to conver is:

import utm
X = []
Y = []
for i in np.arange(len(lat)):
    LAT,LON,Z,S = utm.from_latlon(lat[i],lon[i])
    X = np.append(X,LON)
    Y = np.append (Y,LAT)

as suggested in comments: I switch lon lat positn i.e, plt.scatter (Y,X) it gives:

enter image description here

Am I doing something wrong, or is this result correct?

Advertisement

Answer

The answer might be helpful for others so I have posted my solution. I found pyproj works better then utm. In pyproj we can specify utm zone.

import pyproj

from pyproj import Proj

myProj = Proj("+proj=utm +zone=45U, +north +ellps=WGS84 +datum=WGS84 +units=m +no_defs")
lon_,lat_ = myProj(lon, lat) 
plt.scatter(lon_,lat_)

and the figure looks good!

enter image description here

User contributions licensed under: CC BY-SA
9 People found this is helpful
Advertisement