인간 남자케릭 아이스의 블로그

어제 약소하게 질렀다고 하는 놈이.. 이놈이다.
가격은 얼마 안되지만, 두둥… 무려 4개월 할부(무이자로;;)로 질렀다.
왜냐면 할부로 지른놈들이 4개월 후에 모두 할부가 끝나기 때문..
한꺼번에 할부가 다 떨어져 버리면 기분도 더 좋겠지...

이놈 이름이 큐비랩에서 만든, Location Plus A1이라고 하더라….

 
요렇게 생겼는데, 한손에 잡히는정도의 크기다. 사진을 찍다 보니.. 남는것은 사진밖에 없다는 생각이 들어..
'사진에 더 많은 의미를 부여하기위해' 라는 변명으로 구입했다. 이놈으로 사진에 찍은 위치 정보를 추가할 수 있다.

이놈을 PC에 USB로 연결해 외장 GPS수신기로 사용할 수도 있지만 주로 GPS데이터를 그대로 로깅하는 기능을 한다.(저 회사에서 내세우는, 사진에 위치 추가하는 기능은 포함된 SW의 역할이다.)

로깅되는 GPS데이터에는 시간정보도 포함되어 있으므로.. 사진을 찍은 시간과 GPS데이터를 이용하면 그 사진이 어디서 찍혔는지 파악할 수 있다. 그 점을 저 회사에서 착안하여 만든 것 같다.

저놈 설명서를 보니, 외장 GPS수신기로 사용 가능하다고 해서 VisualGPS라는 GPS테스트 프로그램으로 체크해봤더니.. 방안 창문가에서 GPS위성이 8개가 잡힐 정도로 수신률이 괜찮았다.

이놈의 주요 사용법은.. 전원을 켠 후 그냥 가지고 다니면 된다.. 그러면 GPS데이터가 파일로 저장되고..
PC에 연결하면 그 파일을 보고 조작할 수 있다.(대용량저장장치)

어떤 사람들을 보니, 이런 데이터로 구글어스에 이동궤적을 표시했던데.. 어떤 프로그램을 쓰는지는 모르겠고..
그냥 간단하게 구글어스용 KML생성기를 만들어서 적용해 보았다.

GPS데이터를 이용해, 궤적을 그려주는 KML파일로 생성해 준다.

언어는 외계어(-_-..)가 아니고 python이니, 알아서 이해하기 바란다.

# NMEA.py import string import math class NMEA: def __init__(self, line=None): self.prn = range(12) self.elevation = range(12) self.azimuth = range(12) self.ss = range(12) self.zs = range(12) self.zv = range(12) self.time = '?' self.mode = 0 self.lat = 0.0 self.lon = 0.0 self.altitude = 0.0 self.track = 0.0 self.speed = 0.0 self.in_view = 0 self.SAT = 0 self.ZCH = 0 self.ZCHseen = 0 self.LATLON = 0 # added by heisice(08-05-14) if line is not None: self.handle_line(line) def add_checksum(self,sentence): csum = 0 for c in sentence: csum = csum ^ ord(c) return sentence + "%02X" % csum + "\r\n" def checksum(self,sentence, cksum): csum = 0 for c in sentence: csum = csum ^ ord(c) return "%02X" % csum == cksum def update(self, lval, value, category): if lval != value: return (value, 1) else: return (lval, category) #$GPGGA,000032.997,0000.0000,N,00000.0000,E,0,00,50.0,0.0,M,,.j...???.??.?|.V??jj.h..00032.997,V,0000.0000,N,00000.0006 #Lat: 0.000000 Lon: 0.000000 Alt: 0.000000 Sat: 0 Mod: 1 Time: 11/26/2000 00:00:31 #$GPGGA,000033.997,0000.0000,N,00000.0000,E,0,00,50.0,0.0,M,,,,0000*3C def do_lat_lon(self, words): if words[0][-1] == 'N': words[0] = words[0][:-1] words[1] = 'N' if words[0][-1] == 'S': words[0] = words[0][:-1] words[1] = 'S' if words[2][-1] == 'E': words[2] = words[2][:-1] words[3] = 'E' if words[2][-1] == 'W': words[2] = words[2][:-1] words[3] = 'W' if len(words[0]): lat = string.atof(words[0]) frac, intpart = math.modf(lat / 100.0) lat = intpart + frac * 100.0 / 60.0 if words[1] == 'S': lat = -lat (self.lat, self.LATLON) = self.update(self.lat, lat, self.LATLON) if len(words[2]): lon = string.atof(words[2]) frac, intpart = math.modf(lon / 100.0) lon = intpart + frac * 100.0 / 60.0 if words[3] == 'W': lon = -lon (self.lon, self.LATLON) = self.update(self.lon, lon, self.LATLON) #$GPRMC,024932.992,V,4443.7944,N,07456.7103,W,,,270402,,*05 #$GPGGA,024933.992,4443.7944,N,07456.7103,W,0,00,50.0,192.5,M,,,,0000*27 #$GPGSA,A,1,,,,,,,,,,,,,50.0,50.0,50.0*05 #$GPGSV,3,1,09,14,77,023,,21,67,178,,29,64,307,,30,42,095,*7E #$GPGSV,3,2,09,05,29,057,,11,15,292,,18,08,150,,23,08,143,*7A #$GPGSV,3,3,09,09,05,052,*4B #$GPRMC,024933.992,V,4443.7944,N,07456.7103,W,,,270402,,*04 #$GPGGA,024934.991,4443.7944,N,07456.7103,W,0,00,50.0,192.5,M,,,,0000*23 #$GPGSA,A,1,,,,,,,,,,,,,50.0,50.0,50.0*05 # RMC - Recommended minimum specific GPS/Transit data # RMC,225446,A,4916.45,N,12311.12,W,000.5,054.7,191194,020.3,E*68 # 225446 Time of fix 22:54:46 UTC # A Navigation receiver warning A = OK, V = warning # 4916.45,N Latitude 49 deg. 16.45 min North # 12311.12,W Longitude 123 deg. 11.12 min West # 000.5 Speed over ground, Knots # 054.7 Course Made Good, True # 191194 Date of fix 19 November 1994 # 020.3,E Magnetic variation 20.3 deg East # *68 mandatory checksum def processGPRMC(self, words): global seconds # the Navman sleeve's GPS firmware sometimes puts the direction in the wrongw ord. day = string.atoi(words[8][0:2]) month = string.atoi(words[8][2:4]) year = 2000 + string.atoi(words[8][4:6]) hours = string.atoi(words[0][0:2]) minutes = string.atoi(words[0][2:4]) seconds = string.atoi(words[0][4:6]) if words[1] == "V" or words[1] == "A": self.time = ("%02d/%02d/%04d %02d:%02d:%02d" % (day, month, year, hours, minutes, seconds)) if words[6]: self.speed = string.atof(words[6]) if words[7]: self.track = string.atof(words[7]) self.do_lat_lon(words[2:]) # GGA - Global Positioning System Fix Data # GGA,123519,4807.038,N,01131.324,E,1,08,0.9,545.4,M,46.9,M, , *42 # 123519 Fix taken at 12:35:19 UTC # 4807.038,N Latitude 48 deg 07.038' N # 01131.324,E Longitude 11 deg 31.324' E # 1 Fix quality: 0 = invalid # 1 = GPS fix # 2 = DGPS fix # 08 Number of satellites being tracked # 0.9 Horizontal dilution of position # 545.4,M Altitude, Metres, above mean sea level # 46.9,M Height of geoid (mean sea level) above WGS84 # ellipsoid # (empty field) time in seconds since last DGPS update # (empty field) DGPS station ID number def processGPGGA(self,words): self.do_lat_lon(words[1:]) # sometimes I get the following, which of course produces an error when parsing the status: #$GPGGA,051122.091,0000.0000,N,00000.0000,E0,,00,50.0,0.0,M,,,,0000*36 self.status = string.atoi(words[5]) self.satellites = string.atoi(words[6]) self.altitude = string.atof(words[8]) # GSA - GPS DOP and active satellites # GSA,A,3,04,05,,09,12,,,24,,,,,2.5,1.3,2.1*39 # A Auto selection of 2D or 3D fix (M = manual) # 3 3D fix # 04,05... PRNs of satellites used for fix (space for 12) # 2.5 PDOP (dilution of precision) # 1.3 Horizontal dilution of precision (HDOP) # 2.1 Vertical dilution of precision (VDOP) # DOP is an indication of the effect of satellite geometry on # the accuracy of the fix. def processGPGSA(self,words): (self.mode, self.LATLON) = self.update(self.mode, string.atof(words[1]), self.LATLON) self.pdop = string.atof(words[14]) self.hdop = string.atof(words[15]) self.vdop = string.atof(words[16]) # GSV - Satellites in view # GSV,2,1,08,01,40,083,46,02,17,308,41,12,07,344,39,14,22,228,45*75 # 2 Number of sentences for full data # 1 sentence 1 of 2 # 08 Number of satellites in view # 01 Satellite PRN number # 40 Elevation, degrees # 083 Azimuth, degrees # 46 Signal strength - higher is better # <repeat for up to 4 satellites per sentence> # There my be up to three GSV sentences in a data packet def processGPGSV(self,words): n = string.atoi(words[1]) in_view = string.atoi(words[2]) (self.in_view, self.SAT) = self.update(self.in_view, in_view, self.SAT) f = 3 n = (n - 1) * 4; m = n + 4; while n < in_view and n < m: if words[f+0]: (self.prn[n], self.SAT) = self.update(self.prn[n], string.atoi(words[f+0]), self.SAT) if words[f+1]: (self.elevation[n], self.SAT) = self.update(self.elevation[n], string.atoi(words[f+1]), self.SAT) (self.azimuth[n], self.SAT) = self.update(self.azimuth[n], string.atoi(words[f+2]), self.SAT) if f + 3 < len(words) and words[f+3]: (self.ss[n], self.SAT) = self.update(self.ss[n], string.atoi(words[f+3]), self.SAT) f = f + 4 n = n + 1 #PRWIZCH - Rockwell Zodiac Proprietary #Channel Information #$PRWIZCH ,00,0,03,7,31,7,15,7,19,7,01,7,22,2,27,2,13,0,11,7,08,0,02,0*4C #SATELLITE #IDENTIFICATION NUMBER - 0-31 #SIGNAL QUALITY - 0 low quality - 7 high quality Repeats 12 tims. def processPRWIZCH(self,words): for i in range(12): (self.zs[i], self.ZCH) = self.update(self.zs[i], string.atoi(words[2*i+0]), self.ZCH) (self.zv[i], self.ZCH) = self.update(self.zv[i], string.atoi(words[2*i+1]), self.ZCH) self.ZCHseen = 1; def handle_line(self, line): if line[0] == '$': line = string.split(line[1:-1], '*') if len(line) != 2: return if not self.checksum(line[0], line[1]): return "Bad checksum" words = string.split(line[0], ',') if NMEA.__dict__.has_key('process'+words[0]): NMEA.__dict__['process'+words[0]](self, words[1:]) else: return "Unknown sentence" else: return "Not NMEA" def get_status(self,satellite): if self.ZCHseen: for i in range(12): if satellite == self.zs[i]: return (self.zv[i] & 7) | 8 return 0 else: for i in range(12): if satellite == self.prn[i]: s = self.ss[i] / 6 return min(s, 7) return 0 if __name__ == '__main__': # Self-testing code goes here. nmea = NMEA() lines = [ "$GPGGA,000033.997,0000.0000,N,00000.0000,E,0,00,50.0,0.0,M,,,,0000*3C\n", "$GPRMC,024932.992,V,4443.7944,N,07456.7103,W,,,270402,,*05\n", "$GPGSA,A,1,,,,,,,,,,,,,50.0,50.0,50.0*05\n", "$GPGSV,3,1,09,14,77,023,,21,67,178,,29,64,307,,30,42,095,*7E\n", "$GPGSV,3,2,09,05,29,057,,11,15,292,,18,08,150,,23,08,143,*7A\n", "$GPGSV,3,3,09,09,05,052,*4B\n", ] for line in lines: nmea.handle_line(line) print nmea.__dict__
# coding: utf-8 # log2kml # 2008-05-14 by HEISICE(heisice@gmail.com) import NMEA kml_name = "오늘의 이동 경로" kml_desc = "오늘 이동 경로를 선으로 표시합니다." input_filename = "LP20080514_091223.log" output_filename = "output.kml" output = open(output_filename,"w") output.write("""<?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://earth.google.com/kml/2.2"> <Document> <name>%s</name> <description>%s</description> <Style id="lineStyle"> <LineStyle> <color>9fffff00</color> <width>3</width> </LineStyle> </Style> <Placemark> <name>%s</name> <description>%s</description> <styleUrl>#lineStyle</styleUrl> <LineString> <altitudeMode>relativeToGround</altitudeMode> <coordinates>"""%(kml_name, kml_desc, kml_name, kml_desc)) for line in open(input_filename): try: nmea = NMEA.NMEA(line.replace("\r","")) if nmea.lat <> 0: output.write("%s,%s,0\n"%(nmea.lon, nmea.lat)) except: continue output.write("""</coordinates> </LineString> </Placemark> </Document> </kml>""") output.close()

- 추가: 다른 사람들은 GPSBabel이라는 프로그램 쓰나보네.. NMEA 파일 -> 구글어스(keyhole)로 변환 되네..
- 또 추가: 기본제공하는 (이미지에 위치정보 넣는) 프로그램에도 kml로 내보내기 기능이 있네 -_-...
   나 왠 쓸데없는 삽질을 한거지?