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

프리웨어인 구글어스에는, 프로버전에서 있는 GPS연동 기능이 없다.
그래서 간단하게 살짝 만들어봤다..

GPS 신호를 수신하고, 구글어스를 조정해서 현재 위치를 보여주는것 까지다.
별로 신경을 안쓰고 발로 만들었으므로 자잘한 버그는 있지만 동작은 한다.
다운받아 테스트정도는 가능할 것이다.

구글어스를 실행시켜 놓은 상태에서 프로그램을 실행하고
GPS수신기가 연결된 포트 선택후에 연결 버튼만 누르면 현재 위치가 보여진다.

내 두번째(첫번째는 국산 리눅스PDA Yopy였다..) 리눅스 PDA Zaurus..
전에 구입한 GPS 수신기와 연결을 해봐야겠다는 생각이 들어서 꺼내봤더니..
어허...  USB 호스트 케이블이 어디있지?? -_-...

결국 찾다 찾다 못찾아서, 인터넷 쇼핑몰에서 포인트만으로(택배비 포함;;) 1,000원짜리 호스트 케이블을
구입하게 되었다.

팔때 아이리버용이라고 되있지만, USB 호스트케이블이 다 같은놈 아닌가?? -_-.. 아님 말고..
아무튼 USB호스트 케이블을 연결하고, GPS를 연결해 봤다.


자꾸 내용이 날아가네 -_-... IE 문제인가... 다시 쓰기 짜증나니 중간 생략하고..

아무튼 땜질하는 삽질 없이.. 적당히 성공.. 드라이버는 적절한 놈 깔아줘야 한다.
어떤거 깔았는지 궁금하면 질문하삼..



내가 설치한 드라이버는 좀 불안한 것 같다..
프로그램 짜는 도중에 자꾸 GPS와 연결이 끊기는 현상 발생..

나머지는 다음에 계속...

드라이버 다른거 찾아본 후에 네비게이션 프로그램 만들어볼 예정..(아주 천천히...)

갑자기 왠 방문수가 평소의 다섯배가 넘나 해서 유입 경로를 봤더니...

거의 자바, 서버, 파이썬, 패킷분석, 구글어스 등등 관련 자료를 찾으러 오셨구나...;;

어서와요.. :)

어제 약소하게 질렀다고 하는 놈이.. 이놈이다.
가격은 얼마 안되지만, 두둥… 무려 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로 내보내기 기능이 있네 -_-...
   나 왠 쓸데없는 삽질을 한거지?

전에, 뭐좀 만들어볼라다가 기능을 잘못 알아 헛 구입한 것을 깨닳은 후, 잡다구리 박스에 넣어 둔 것중에
Lantronix사의 Xport라는 놈이 있었다.

아주 작은 장치 안에, RTOS를 내장하여.. 전원만 공급하면 웹서버로 동작하고,
LAN을 통해 XPort에 접속 가능하도록 하며, 연결된 장비들을 시리얼 및 I/O포트를 통해 제어할 수 있도록 한다.

이거 구입할 당시에 원하던 기능은 -_-... 시리얼로 인터넷의 특정 페이지를 GET하는거였다.
아무튼, 고놈의 메뉴얼을 오랜만에 봐보니 안보이던 문장이 보여, 한번 가지고 놀아 보았다.

TCP/UDP를 통해 GPIO(General Purpose Input/Output)제어 가능한 포트가 3개(CP1~CP3) 있는데,
그놈들을 제어하는 것을 해보았다. 이것만 가지고도 최소 3개의 장비를 원격으로 ON/OFF할 수 있다.

프로그램은 이렇게 작성해보았다.(파이썬 스크립트)