GPSの研究 その6
概要
GPSを理解したかった。
RINEXから、読み出した衛星の軌道をEARTHで、眺めてみた。
写真
サンプルコード
from lxml import etree
from pykml.factory import KML_ElementMaker as KML
import pyproj
ecef = pyproj.Proj(proj = 'geocent', ellps = 'WGS84', datum = 'WGS84')
lla = pyproj.Proj(proj = 'latlong', ellps = 'WGS84', datum = 'WGS84')
sats = [
"G02 6262234.031370984 -14613110.981922923 21254935.229655903",
"G04 19040661.998708855 -8311756.524602333 16532658.634675292",
"G07 26078664.69836943 5588070.077162541 -239561.27794465935",
"G08 24470813.305349953 -95731.17177840695 -10965308.78942738",
"G10 10946425.106999923 -13217422.981333993 19960118.04486387",
"G13 17400028.443231095 1461399.2922087098 19887536.270703305",
"G17 21078475.93615281 -13351300.07070003 -8639709.183587706",
"G20 14385911.569008492 20514420.345839 8765413.838559393",
"G23 9336918.80677975 12048291.269695656 21723935.800126303"
]
doc = KML.kml(
KML.Folder(
)
)
for sat in sats:
l = sat.split(' ')
x = float(l[1])
y = float(l[2])
z = float(l[3])
lon, lat, alt = pyproj.transform(ecef, lla, x, y, z, radians = False)
pm = KML.Placemark(
KML.name(l[0]),
KML.Point(
KML.altitudeMode("absolute"),
KML.coordinates("{},{},{}".format(lon, lat, alt)),
),
KML.Style(
KML.IconStyle(
KML.Icon(
KML.href('http://maps.google.com/mapfiles/kml/paddle/A.png')
),
KML.scale(40)
)
)
)
doc.Folder.append(pm)
print (etree.tostring(doc, pretty_print = True).decode())
from lxml import etree
from pykml.factory import KML_ElementMaker as KML
import pyproj
ecef = pyproj.Proj(proj = 'geocent', ellps = 'WGS84', datum = 'WGS84')
lla = pyproj.Proj(proj = 'latlong', ellps = 'WGS84', datum = 'WGS84')
sats = [
"G02 6262234.031370984 -14613110.981922923 21254935.229655903",
"G04 19040661.998708855 -8311756.524602333 16532658.634675292",
"G07 26078664.69836943 5588070.077162541 -239561.27794465935",
"G08 24470813.305349953 -95731.17177840695 -10965308.78942738",
"G10 10946425.106999923 -13217422.981333993 19960118.04486387",
"G13 17400028.443231095 1461399.2922087098 19887536.270703305",
"G17 21078475.93615281 -13351300.07070003 -8639709.183587706",
"G20 14385911.569008492 20514420.345839 8765413.838559393",
"G23 9336918.80677975 12048291.269695656 21723935.800126303"
]
doc = KML.kml(
KML.Folder(
)
)
for sat in sats:
l = sat.split(' ')
x = float(l[1])
y = float(l[2])
z = float(l[3])
lon, lat, alt = pyproj.transform(ecef, lla, x, y, z, radians = False)
pm = KML.Placemark(
KML.name(l[0]),
KML.Point(
KML.altitudeMode("absolute"),
KML.coordinates("{},{},{}".format(lon, lat, alt)),
),
KML.Style(
KML.IconStyle(
KML.Icon(
KML.href('http://maps.google.com/mapfiles/kml/paddle/A.png')
),
KML.scale(40)
)
)
)
doc.Folder.append(pm)
print (etree.tostring(doc, pretty_print = True).decode())
以上。
Author And Source
この問題について(GPSの研究 その6), 我々は、より多くの情報をここで見つけました https://qiita.com/ohisama@github/items/c779354f72d2eb895cc9著者帰属:元の著者の情報は、元の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 .