Theos Where are you

______สืบเนื่องมาจากมีคำถามเรื่อง ephemeris ของดาวเทียมทาง blog ผมเองก็ไม่ได้รู้อะไรมากมายหรอกครับ ออกตัวไว้ก่อน เพียงแต่สนใจและมีโจทย์จากอาจารย์มาให้ลองเล่นเลยมีโอกาสได้ศึกษา วันนี้เลยขอมาแชร์แนวทางการเรียนรู้ เพื่อให้ท่านได้นำไปพัฒนาต่อ จริงๆแล้วการที่เราจะเข้าใจวงโคจรของดาวเทียม ก็จะทำให้เราสามารถทราบถึงตำแหน่งของดาวเทียมนะเวลาต่างๆที่อยู่บนหัวเราได้ ผมคงไม่อธิบาย Keplerian Elements เพราะดูแล้วเยอะ ถ้าต้องการทำโปรแกรมที่ precise มากก็ไปหาหนังสือต่างประเทศอ่านน่าจะดี ทุกอย่างที่อธิบายวงโคจรของดาวเทียมอยู่ใน Keplerien Elements เช่น Inclination (มุมระนาบวงโคจรกับ equator) , Right Ascension of Ascending Node (อันนี้ลืม), Eccentricity (ค่าเยื้องศูนย์) , descending เป็นต้น ค่าเหล่านี้ใช้อธิบายการเคลื่อนที่ของดาวเทียมตามแบบ Kepler Model

_____แน่นอนว่าถ้าเราต้องการการคำนวณขั้นสูงการใช้สมการเพื่อคำนวณหาตำแหน่งบนโมเดลย่อมเป็นเรื่องที่ยุ่งยากและซับซ้อน แต่สำหรับการคำนวณพื้นฐาน ใน python มีโมดูลที่ชื่อ pyephem สามารถคำนวณหาพารามิเตอร์ต่างๆที่กล่าวมาได้ จากค่า TLE (Two line Element) ซึ่งมีแบบเปิดไม่ถูกต้องมากนักจาก NORAD ผมทดลองให้คำนวณค่าตำแหน่ง(x,y) และมุมสูง โดยผมใช้ TLE ของ THEOS จาก NORAD ซึ่งดาวน์โหลดได้จาก  http://celestrak.com/NORAD/elements/resource.txt

_____การคำนวณคราวๆนะครับ ตัวอย่างการเขียน code ก็ไม่ยากครับ ดังนี้

import ephem
Session = ephem.Observer()
Session.date = ‘2009/12/1 00:00:00’
LOCAL_THAI = 7
# NORAD two-line element
l1 = “Thai THEOS”
l2 = “1 33396U 08049A   08348.12285235 -.00000044  00000-0  51000-9 0  1071”
l3 = “2 33396 098.7734 051.3202 0001046 089.3375 270.7921 14.20009526 10371”
sat = ephem.readtle(l1, l2, l3)
print sat.name
sat.compute(Session)
print “Right ascension : “, sat.g_ra
print “Declination : “, sat.g_dec
print “Speed/Range: “, sat.range_velocity
print “Azimuth : “, sat.az
print “eclipse :”, sat.eclipsed
print “Magnitude :”, sat.mag
print “*********************************************”
print “position (lon,lat,attitude)”
print sat.sublong, sat.sublat, (sat.elevation)/1000
ผลจากการประมวลผล
_____ผลที่ได้ก็สามารถนำมาใช้วางแผนคราวๆในการหาตำแหน่งของดาวเทียม เพื่อถ่ายภาพหรือหาความเป็นไปได้ของการมีภาพถ่ายดาวเทียม ณ ช่วงเวลานั้นๆที่ตำแหน่งนั้น แน่นอนว่าสามารถแสดงเป็น animation เพื่อดูการเคลื่อนผ่านของดาวเทียมได้ดังภาพ

ศึกษาเพิ่มเติมที่

http://www.amsat.org/amsat/keps/kepmodel.html

http://celestrak.com/NORAD/elements/

Advertisements

ใส่ความเห็น

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / เปลี่ยนแปลง )

Twitter picture

You are commenting using your Twitter account. Log Out / เปลี่ยนแปลง )

Facebook photo

You are commenting using your Facebook account. Log Out / เปลี่ยนแปลง )

Google+ photo

You are commenting using your Google+ account. Log Out / เปลี่ยนแปลง )

Connecting to %s

สร้างเว็บไซต์หรือบล็อกฟรีที่ WordPress.com.

Up ↑

%d bloggers like this: