規定のブラウザでURLを指定して開く場合のサンプルプログラム (google mapで緯度、経度を指定) import webbrowser #pythonにデフォルトで入っている webbrowser.open('https://www.google.com/maps?q='+str(blh[0])+','+str(blh[1])) #↑でURL指定すると規定…

pythonのpandasでプロットする場合とnumpy→matplotlibでプロットする場合 import pandas as pd import numpy as np import matplotlib.pyplot as plt df1 = pd.read_csv('****.csv') #csvを読み込み df1.plot(x='**x',y='**y') #pandasの機能でプロットする…

ちょっとマニアックだが…ecef座標系をENU(EAST-NORTH-UP)座標に変換するpythonプログラムを使って 移動経路などを3次元プロットするプログラム。例えばGPSから得たECEF座標から基準点を原点として飛行経路を描くときに使えます。 5行目に出てくるecef_enuと…

ECEF座標からENU座標に変換するプログラムです。 下記を参考にしてpythonでプログラムを作成しました。 www.enri.go.jppythonを使うと行列の扱いがとても楽です。 別プログラムで使用しやすいようecef2enuという関数にしています。 #WGS84 -> ENU import mat…

下記を参考にしてECEF座標系から緯度経度、その逆を変換するpythonプログラムを 作成したので、備忘録を兼ねて記載しておきます。 www.enri.go.jp汎用的に使えるよう、それぞれ関数にしておきました。 #緯度経度からECEF座標系を算出する import math import…

GPS時刻をJSTに変換するプログラム print("Please input gps-time") gps_time = int(input())+(9*60*60)-19 #うるう秒は19秒で計算 roll_over = gps_time // (1024*7*24*60*60) week_time = gps_time - (roll_over*1024*7*24*60*60) week_number = week_time…

緯度経度→ECEF座標系に変換する。 ECEF座標系→緯度経度に変換するプログラム #緯度経度からECEF座標系を算出する import math #東京タワー p_deg = 35.658581 #[deg] l_deg = 139.745433 #[deg] h = 333.0 #[m] p_rad = math.radians(p_deg) l_rad = math.ra…