Как я могу извлечь координаты x, y и z из географических данных Python? - PullRequest
Новая
       15

Как я могу извлечь координаты x, y и z из географических данных Python?

0 голосов
/ 29 января 2009

У меня есть географические данные, которые имеют 14 переменных. Данные представлены в следующем формате:

QUADNAME: rockport_colony_SD РАЗРЕШЕНИЕ: 10 ULLAT: 43,625
ULLON: -97,87527466 LRLAT: 43,5
LRLON: -97,75027466 HDATUM: 27
ZMIN: 361,58401489 ZMAX: 413.38400269 ZMEAN: 396.1293335 ZSIGMA: 12.36359215 PMETHOD: 5
QUADDATE: 20001001

Целые данные имеют много предыдущих переменных в последовательности.

Как извлечь из данных координаты ULLAT, ULLON и LRLAT из трех списков, чтобы каждая строка соответствовала одному местоположению?

Этот вопрос был задан проблемой в посте .

Ответы [ 2 ]

4 голосов
/ 29 января 2009

Примерно так может работать, если все данные находятся в большом плоском текстовом файле: 

import re

data = """
QUADNAME: rockport_colony_SD RESOLUTION: 10 ULLAT: 43.625
ULLON: -97.87527466 LRLAT: 43.5
LRLON: -97.75027466 HDATUM: 27
ZMIN: 361.58401489 ZMAX: 413.38400269 ZMEAN: 396.1293335 ZSIGMA: 12.36359215 PMETHOD: 5
QUADDATE: 20001001
"""

regex = re.compile(
    r"""ULLAT:\ (?P<ullat>-?[\d.]+).*?
    ULLON:\ (?P<ullon>-?[\d.]+).*?
    LRLAT:\ (?P<lrlat>-?[\d.]+)""", re.DOTALL|re.VERBOSE)

print regex.findall(data) # Yields: [('43.625', '-97.87527466', '43.5')]
2 голосов
/ 29 января 2009

С учетом StreamReader именованного читателя, это должно дать вам список (float, float, float). Я предлагаю список из 3-х кортежей, потому что, вероятно, будет удобнее и эффективнее проходить через него, если только по какой-то причине вы не хотите получать все очки индивидуально. 

coords = []
reader
while line=reader.readline():

  index_ullat = line.find("ULLAT")
  if index_ullat >= 0:
    ullat = float(line[ index_ULLAT+7 : ])

    line = reader.readline()

    index_ullon = line.find("ULLON")
    index_lrlat = line.find("LRLAT")
    if index_ullon >= 0 and index_lrlat >= 0:
      ullon = float(line[ index_ullon+7 : index_lrlat-1 ])
      lrlat = float(line[ index_lrlat+7 : ])
    else:
      raise InputError, "ULLON and LRLAT didn't follow ULLAT."

    coords.append(ullat, ullon, lrlat)

Это может сработать, но безобразно. Я не эксперт в разборе строк.

...