APRS Lati/Long Decoding: Satellites information update by JE9PEL

Duchifat-1 の管制局から送られてきた compressed packet は、!"<X?%qh$$y!!_
です。｢個人｣ を特定するものですが 悪用のしようがないと判断し、ここに公開
して具体的な計算を行ってみます。このコードをアスキー変換して計算した結果
は、自局の緯度経度(Lati,Long)と小数点以下第5位まで正確に一致しました。

APRS compressed packet is: !"<X?%qh$$y!!_

Lati = <X?%
                y1  y2  y3  y4
     < X ? % =  60  88  63  37

     Lati = 90-((y1-33)*91^3+(y2-33)*91^2+(y3-33)*91+(y4-33))/380926
          = 90-((60-33)*91^3+(88-33)*91^2+(63-33)*91+(37-33))/380926
          = 90-(27*753571+55*8281+30*91+4)/380926
          = 90-(20346417+455455+2730+4)/380926
          = 90-20804606/380926
          = 90-54.61587
          = 35.38413 (北緯35.38413度)

Long = qh$$
                x1  x2  x3  x4
     q h $ $ = 113 104  36  36

     Long = -180+((x1-33)*91^3+(x2-33)*91^2+(x3-33)*91+(x4-33))/190463
          = -180+((113-33)*91^3+(104-33)*91^2+(36-33)*91+(36-33))/190463
          = -180+(80*753571+71*8281+3*91+3)/190463
          = -180+(60285680+587951+273+3)/190463
          = -180+60873907/190463
          = -180+319.61014
          = 139.61014 (東経139.61014度)

[参考]
http://www.aprs.org/doc/APRS101.PDF
http://www.qsl.net/py4zbz/satelite/aprscompr.gif
http://web-apps.nbookmark.com/ascii-converter/
http://www.asahi-net.or.jp/~ei7m-wkt/numbr526.htm
http://www.asahi-net.or.jp/~ei7m-wkt/numbr879.htm
http://blog.goo.ne.jp/je1cvl/d/20151005