Понижение частоты дискретизации GPS

Ваш код можно немного реорганизовать, чтобы сделать его более понятным. Вам нужно добавить дополнительную запись всякий раз, когда отображается EOL для случая, когда расстояние меньше 100 м:

from math import sin, cos, sqrt, atan2, radians    

def distance(coord1, coord2): #Haversin
    lat1,lon1=coord1
    lat2,lon2=coord2
    dlat = radians(lat2-lat1)
    dlon = radians(lon2-lon1)
    a = sin(dlat/2) * sin(dlat/2)
    + cos(radians(lat1))*cos(radians(lat2))*sin(dlon/2)*sin(dlon/2)
    c = 2 *atan2(sqrt(a),sqrt(1-a))
    s = (6367*c)*1000 #meter
    return s

def get_coordinates(data):
    # Importing only coordinates from asko input file (Row 2 and 4)
    # Converting the coordinates from DDMM.MMMM to DD.DDDDDD

    LAT = (data[2])
    LAT_D = LAT[0:2]               
    LATID = float(LAT_D)

    LAT_M = LAT[2:]
    LATM = float(LAT_M)
    LATIM = float(LATM) / 60.0

    latitude = (LATID + LATIM)                  

    LON = (data[4])
    LON_D = LON[1:3]
    LONGD = float(LON_D)

    LON_M = LON[3:]
    LONM = float(LON_M)
    LONGM = float(LONM) / 60.0

    longitude = (LONGD + LONGM)

    return (latitude, longitude)


coord1 = None

# with open as data will close itself after reading each line. so you don't need to close it yourself

with open('asko_nav_2013.nmea', 'r') as indata, open('asko_nav_out.txt', 'w') as outdata:
    for line in indata:
        if line.startswith('EOL'):  #if the line starts with EOL(end of line) it writes it in the output
            if dist < 100:
                outdata.write('%0.6f\t%0.6f\t%f\n' % (latitude, longitude, dist))
            outdata.write("\nEOL\n")
            coord1 = None   # Reset the first coordinate
        elif line.startswith('Line'): 
            outdata.write('\n%s' % line)
        elif line.startswith('$GPGGA'):  #when the fist line starts with $GPGGA it splits the columns
            data=line.split(",")        #the for loop reads the file line by line
            latitude, longitude = get_coordinates(data)

            if coord1:
                coord2 = (latitude, longitude)
                dist = distance(coord1, coord2)

                if dist >= 100:
                    outdata.write('%0.6f\t%0.6f\t%f\n' % (latitude, longitude, dist))
                    coord1 = coord2         
            else:
                # The first time through the loop "coord1" is None
                outdata.write('%0.6f\t%0.6f\t0.0 \n' % (latitude, longitude))
                coord1 = (latitude, longitude)

Для вашего ввода это создает следующий выходной файл:

Line 20130904_0848.nmea
58.822923   17.664710   0.0 
58.828186   17.664946   584.888514

EOL

Line 20130904_0926.nmea
58.827928   17.666404   0.0 
58.827936   17.666413   0.870480

EOL

Вам также необходимо сбросить coord1 всякий раз, когда обнаруживается EOL, чтобы убедиться, что 0 снова отображается для первой записи.

Немного сложно понять, решит ли это проблему полностью, поскольку ваши образцы данных, похоже, не соответствуют ожидаемому результату.

Решение второй проблемы, связанной с меньшим количеством строк результатов, чем ожидалось: вы предоставляете слишком мало информации о характере вашей проблемы и входных данных, которые вы обрабатываете. Выборка вашего ввода «каждые 100 м» может означать что-то другое, если ваши входные данные отбираются из траектории, по которой движется движущийся объект, особенно если движение не является чисто линейным.

Представьте, что ваши входные данные описывают координаты, полученные путем измерения GPS-координат через равные промежутки времени при движении по окружности с радиусом меньше, скажем, 15 м. Тогда независимо от того, сколько точек данных предоставляет ваш ввод, вывод для предлагаемого вами решения никогда не будет длиннее двух строк, потому что никакие две точки вдоль этой кривой не могут иметь абсолютное расстояние более 100 м. Это может объяснить, почему вы видите меньше строк в выводе, чем ожидалось.

Если вы хотите сэмплировать входные данные через каждые 100 м пройденного, вам придется суммировать все расстояния между входными семплами, начиная с последней точки, выбранной для выходных данных, и использовать это вместо dist. Модифицируя реорганизованный код Мартина, это можно было бы сделать так (некоторые строки для краткости опущены):

coord1 = None
coord_last = None  # holds coordinate of last input sample
dist = 0.0         # total distance travelled since coord1
# [...]
with open('asko_nav_2013.nmea', 'r') as indata, open('asko_nav_out.txt', 'w') as outdata:
    for line in indata:
    # [...]
            if coord1:
                coord2 = (latitude, longitude)
                delta = distance(coord_last, coord2)
                dist += delta
                coord_last = coord2

                if dist >= 100:
                    outdata.write('%0.6f\t%0.6f\t%f\n' % (latitude, longitude, dist))
                    coord1 = coord2
                    dist = 0.0
            else:
                # The first time through the loop "coord1" is None
                outdata.write('%0.6f\t%0.6f\t0.0 \n' % (latitude, longitude))
                coord1 = (latitude, longitude)
                coord_last = coord1
                dist = 0.0