【Mosaic-HAT】Raspberry Pi+RTKLIBでRTK移動局を作ってみよう2
前回、移動局(Rover)を作りましたが、補正処理をRTKLIBに任せました。
Mosaic-HAT自身も補正処理機能を備えているので、そちらで動かしてみましょう。
※移動局その1はこちら。https://qiita.com/KIT-tokunaga/items/f9a7249bdba8b1aceb3b
1. 構成
構成は下図の通り。基準局の測位情報を、RTK2goなどのNtripサーバを経由し、さらにラズパイを経由し、mosaicHATに入れてあげるだけ。基準局はRTCMv3などのフォーマットで情報を配信しているので、このフォーマットのままシリアルポートに流し込んであげればMosaicHATがフォーマットを理解し、誤差補正してくれるのです。近くに基準局があれば、わざわざRTK2goを通さずとも、シリアルからシリアルへと伝達してあげれば、同じことができます。
2. RTKLIB/str2strの設定
RTK2goから基準局の情報を読み出し、Mosaic-HATに流し込むツールは、RTKLIB/str2strを使います。これは基準局の作り方でも出てきました。
※基準局はこちら。https://qiita.com/KIT-tokunaga/items/5df5854f1471528f31ee
基準局では、Mosaic-HATの情報をRTK2goに配信するために使いましたが、移動局はその逆になります。つまり、inとoutを変えます。
cd /home/pi/myRTKLIB/RTKLIB/app/str2str/gcc
./str2str \
-in ntrip://rtk2go.com:2101/MOUNT_POINT#rtcm3 \
-out serial://ttyS0:9600:8:n:1:off#rtcm3
cd先はstr2strの実態のあるパスにしてください。
また、-inにあるMOUNT_POINTには、使用する基準局のマウントポイントにしてください。-outは、シリアルポートを示します。Mosaic-HATと繋がるコネクタにアサインされているポートはttyS0になります。スピードを9,600bpsにしていますが、インターネット側が遅れた場合を想定しています。なお、注意点として、このシリアルポートのMosaic-HAT側(COM1)は、デフォルト119,200bpsとなっていますので、上記を実行する前に、速度を変更しないといけません。次節にMosaic-HATの設定コマンドをまとめて示します。
3. Mosaic-HATの設定
Mosaic-HATに設定すべき内容は、以下の3点です。
1) 出力データのフォーマット指定
2) 送るメッセージ種別の設定
3) 速度を9600に変更
1)は、setDataInOutコマンドでCOM1の入力をコマンド(CMD)、出力をNMEAメッセージに指定します。2)は、NMEAメッセージのうち、どれを出力させるかの設定となり、移動局は補正済の位置情報が最低限必要となるので、setNMEAOutputコマンドを使ってGGAメッセージを1秒間隔で出力する設定をします。3)は、setCOMSettingsコマンドで設定します。以下にスクリプト例を示します。この例では、最初にリセット(GPIO5をLOW)して、コマンドモードのおまじない(S連打)の後に、3つの設定を行っています。なお、コマンドを確実に認識したかを見るために、Mosaic-HATからの反復メッセージ($Rで始まるコメンドのオウム返し)を見るようにしました。
なお、2のstr2strよりも先にこのスクリプトを実行してください。2のシェルの最初に、"python3 ./setNMEA.py"を入れておけば、まとめて実行できます。
import serial
import RPi.GPIO as GPIO
from time import sleep
#Module Reset
print("Module Reset")
GPIO.setwarnings(False)
GPIO.setmode(GPIO.BCM)
GPIO.setup(5, GPIO.OUT, initial=GPIO.LOW)
sleep(1)
GPIO.output(5, GPIO.HIGH)
sleep(10)
#COM Initialize
print("COM Initialize")
ser = serial.Serial('/dev/ttyS0', 115200)
sleep(1) # wait for connection
#Module Initialize
def PutCmd(command) :
ser.write(command)
while 1 :
msg = ser.readline().decode()
print(msg)
if '$R' in msg :
print(msg)
break
print("Module Initialize")
ser.write(b'SSSSSSSSSSSSS\n') # push MOSAIC to run in command mode
sleep(0.1)
PutCmd(b'sdio ,COM1,CMD,NMEA\n') # SetDataInOut command
sleep(0.1)
PutCmd(b'sno ,Stream1,COM1,GGA,sec1\n') # SetNMEAOutput command
sleep(0.1)
ser.write(b'scs ,COM1,baud9600\n') # SetCOMSettings command
print("Finish")
ser.close()
3. 測位情報の読み出し
最後に、測位情報を取り出します。Mosaic-HATは、NMEAフォーマットでGGAメッセージを出力しつづけているので、これをラズパイで採取します。以下のスクリプト例は、取得した位置情報を表示するとともに、sol.posというファイル名でログを出力します。RTKLIBのrtkplotで読み込むことができます。
注意点ですが、GGAメッセージの緯度経度は、dddmm.mmmmというフォーマットです。例えば、3612.3456という値は、36°12.3456分を指します。分秒ではありません。なので、10進に変換するには、ddd+mm.mmmm/60という計算をします。min2decでその計算をしています。
import serial
import time
from datetime import datetime, timezone, timedelta
import chardet
def min2dec(value) :
val_deg = int(value)
val_min = value * 100 - val_deg * 100
val_dec = float(val_deg + val_min / 60)
return val_dec
#Initialize log file
with open('sol.pos', mode='w') as f :
f.write("% JST latitude(deg) longitude(deg) height(m) Q \n")
while True:
with serial.Serial('/dev/ttyS0', 9600) as ser: # Mosaic's COM1
msg_bytes = ser.readline()
if chardet.detect(msg_bytes)['encoding'] != "ascii" :
continue
msg_string = str(msg_bytes.decode())
msg_string = msg_string.rstrip()
if (msg_string.startswith('$GPGGA')):
nmea_array = [element.strip() for element in msg_string.split(',')]
Quality_Indicator = int(nmea_array[6])
if Quality_Indicator == 0:
print("No GPS Fix Available!")
else:
# parse NMEA GGA message
UTC_Time = datetime.fromtimestamp(float(nmea_array[1]), timezone.utc)
Latitude = float(nmea_array[2])*0.01
Latitude_direction = nmea_array[3]
Longitude = float(nmea_array[4])*0.01
Longitude_direction = nmea_array[5]
Height = float(nmea_array[9])
Height_unit = nmea_array[10]
Lat_deg = min2dec(Latitude)
Long_deg = min2dec(Longitude)
# print coordicates
print('Quality: ' + str(Quality_Indicator))
print('UTC Time: ' + UTC_Time.strftime("%H:%M:%S.%f"))
print(' Latitude: ' + str(Latitude) + Latitude_direction)
print(' Longitude: ' + str(Longitude) + Longitude_direction)
print(' Height: ' + str(Height) + Height_unit)
print(' Lat_deg: ' + str(Lat_deg))
print(' Long_deg: ' + str(Long_deg))
# pos.sol output
with open('sol.pos', mode='a') as f :
f.write("0 {0} {1} {2} {3} {4}\n".format(UTC_Time.strftime("%H:%M:%S.%f"), Lat_deg, Long_deg, Height, Quality_Indicator))
time.sleep(0.1)
else:
continue
Author And Source
この問題について(【Mosaic-HAT】Raspberry Pi+RTKLIBでRTK移動局を作ってみよう2), 我々は、より多くの情報をここで見つけました https://qiita.com/KIT-tokunaga/items/3a9b46843c847fb4157b著者帰属:元の著者の情報は、元のURLに含まれています。著作権は原作者に属する。
Content is automatically searched and collected through network algorithms . If there is a violation . Please contact us . We will adjust (correct author information ,or delete content ) as soon as possible .