Created
September 2, 2012 16:25
-
-
Save tomwhipple/3601130 to your computer and use it in GitHub Desktop.
convert DMS coordinates to decimal in Python
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/env/python | |
# coding=utf8 | |
""" | |
Converting Degrees, Minutes, Seconds formatted coordinate strings to decimal. | |
Formula: | |
DEC = (DEG + (MIN * 1/60) + (SEC * 1/60 * 1/60)) | |
Assumes S/W are negative. | |
""" | |
import re | |
def dms2dec(dms_str): | |
"""Return decimal representation of DMS | |
>>> dms2dec(utf8(48°53'10.18"N)) | |
48.8866111111F | |
>>> dms2dec(utf8(2°20'35.09"E)) | |
2.34330555556F | |
>>> dms2dec(utf8(48°53'10.18"S)) | |
-48.8866111111F | |
>>> dms2dec(utf8(2°20'35.09"W)) | |
-2.34330555556F | |
""" | |
dms_str = re.sub(r'\s', '', dms_str) | |
if re.match('[swSW]', dms_str): | |
sign = -1 | |
else: | |
sign = 1 | |
(degree, minute, second, frac_seconds, junk) = re.split('\D+', dms_str, maxsplit=4) | |
return sign * (int(degree) + float(minute) / 60 + float(second) / 3600 + float(frac_seconds) / 36000) |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Also assumes a very rigid format in line 39. This also chokes on something like
S26d15'10"
as thedegree
variable would then be an empty string. I changed this to:I put up a fork with this fix and the fixes mentioned by others here: https://gist.github.com/jeteon/89c41e4081d87b798d8006b16a52c695.