2012-06-01 14 views
7

Ciao Ho un db legacy con alcuni dati posizionali. I campi sono solo campi di testo con stringhe come questa 0°25'30"S, 91°7'W. C'è un modo in cui posso convertire questi in due numeri in virgola mobile per Decimal Latitude e Decimal Longitude?Python - Lotto converte le posizioni GPS in decimali Lat Lon

EDIT:

Quindi, un esempio potrebbe essere: 0°25'30"S, 91°7'W ->0.425, 91.116667 dove la posizione originale singolo campo produce due galleggianti.

Qualsiasi aiuto molto apprezzato.

+0

Quale sarebbe l'uscita corrispondente per la coordinata fornito. Inoltre, la latitudine ha 3 campi numerici, la longitudine 2 .. è quella tipica di come questi sono specificati (e costante nei tuoi dati)? Cosa hai provato fino ad ora? – Levon

+0

Quindi voglio generare su float: '0 ° 25'30" S, 91 ° 7'W' -> '0.425',' 91.116667'. Sembra che i dati possano o meno avere associati minuti. nessuno poi può assumere 0. –

+0

Longitudine 180 W = -180 180 E = 180 Latitudine 90 N = 90 90 S = -90 Esempio bisogna dare -0,425, -91,116667. Controllare questo cercando 0 ° Ad esempio, 25'30 "S, 91 ° 7'W su Google Maps. –

risposta

2

Converte la stringa di input nell'output previsto. Può gestire minuti e secondi in cui non è presente.

Attualmente non conta per Nord/Sud, Est/Ovest. Se mi dirai come vorresti che vengano gestiti, aggiornerò la risposta.

# -*- coding: latin-1 -*- 
import re 

PATTERN = re.compile(r"""(?P<lat_deg>\d+)°  # Latitude Degrees 
         (?:(?P<lat_min>\d+)')? # Latitude Minutes (Optional) 
         (?:(?P<lat_sec>\d+)")? # Latitude Seconds (Optional) 
         (?P<north_south>[NS]) # North or South 
         ,[ ] 
         (?P<lon_deg>\d+)°  # Longitude Degrees 
         (?:(?P<lon_min>\d+)')? # Longitude Minutes (Optional) 
         (?:(?P<lon_sec>\d+)")? # Longitude Seconds (Optional) 
         (?P<east_west>[EW]) # East or West 
         """, re.VERBOSE) 

LAT_FIELDS = ("lat_deg", "lat_min", "lat_sec") 
LON_FIELDS = ("lon_deg", "lon_min", "lon_sec") 

def parse_dms_string(s, out_type=float): 
    """ 
    Convert a string of the following form to a tuple of out_type latitude, longitude. 

    Example input: 
    0°25'30"S, 91°7'W 
    """ 
    values = PATTERN.match(s).groupdict() 

    return tuple(sum(out_type(values[field] or 0)/out_type(60 ** idx) for idx, field in enumerate(field_names)) for field_names in (LAT_FIELDS, LON_FIELDS)) 


INPUT = """0°25'30"S, 91°7'W""" 

print parse_dms_string(INPUT) # Prints: (0.42500000000000004, 91.11666666666666) 
+0

Grazie. Fammi vedere cosa posso fare con questo così com'è. –

16

Questo approccio può trattare con secondi e minuti assente, e credo che gestisce le direzioni della bussola in modo corretto:

# -*- coding: latin-1 -*- 

def conversion(old): 
    direction = {'N':1, 'S':-1, 'E': 1, 'W':-1} 
    new = old.replace(u'°',' ').replace('\'',' ').replace('"',' ') 
    new = new.split() 
    new_dir = new.pop() 
    new.extend([0,0,0]) 
    return (int(new[0])+int(new[1])/60.0+int(new[2])/3600.0) * direction[new_dir] 

lat, lon = u'''0°25'30"S, 91°7'W'''.split(', ') 
print conversion(lat), conversion(lon) 
#Output: 
0.425 91.1166666667 
+1

+1 per non farmi male alla testa .. "Semplice è meglio che complesso." – Levon

+1

Nessuna regex = +1. – Droogans

+0

E se i secondi sono in formato float? – ratata

1

Un approccio semplice (visto che ho imparato sulle espressioni regolari solo oggi a causa della questo problema). Si occupa di campi mancanti e indicazioni bussola.

# -*- coding: latin-1 -*- 
import re 
s = """0°25'30"S, 91°7'W""" 

def compLat_Long(degs, mins, secs, comp_dir): 
    return (degs + (mins/60) + (secs/3600)) * comp_dir 

def extract_DegMinSec(data): 
    m = re.search(r'(\d+°)*(\d+\')*(\d+")*', data.strip()) 
    deg, mins, secs = [0.0 if m.group(i) is None else float(m.group(i)[:-1]) for i in range(1, 4)] 
    comp_dir = 1 if data[-1] in ('W', 'S') else -1 
    return deg, mins, secs, comp_dir 

s1, s2 = s.split(',') 
dms1 = extract_DegMinSec(s1) 
dms2 = extract_DegMinSec(s2) 
print('{:7.4f} {:7.4f}'.format(compLat_Long(*dms1), compLat_Long(*dms2))) 

cede

0.4250 91.1167 
Problemi correlati