รูปแบบข้อมูล
ฉันจะใช้เป็นแบบฝึกหัดเล็กน้อยในการตั้งโปรแกรมตัวอ่านข้อมูล ดูที่เอกสารประกอบ :
ข้อมูล SRTM มีการกระจายในสองระดับ: SRTM1 (สำหรับสหรัฐอเมริกาและดินแดนและการครอบครอง) โดยมีการสุ่มตัวอย่างข้อมูลในหนึ่งช่วงอาร์ควินาทีในละติจูดและลองจิจูดและ SRTM3 (สำหรับโลก) สุ่มตัวอย่างที่สามอาร์ควินาที
ข้อมูลถูกแบ่งออกเป็นหนึ่งโดยหนึ่งองศาละติจูดและลองจิจูดแผ่นในประมาณการ "ทางภูมิศาสตร์" ซึ่งจะกล่าวว่างานนำเสนอแรสเตอร์ที่มีช่วงเวลาที่เท่าเทียมกันของละติจูดและลองจิจูดโดยไม่ต้องฉายภาพเลย แต่ง่ายต่อการจัดการและโมเสก
ชื่อไฟล์อ้างถึงละติจูดและลองจิจูดของมุมซ้ายล่างของแผ่นกระเบื้อง - เช่น N37W105 มีมุมซ้ายล่างที่ละติจูด 37 องศาเหนือและลองจิจูดละติจูด 105 องศาทางตะวันตก เพื่อความแม่นยำยิ่งขึ้นพิกัดเหล่านี้อ้างถึงศูนย์กลางทางเรขาคณิตของพิกเซลซ้ายล่างซึ่งในกรณีของข้อมูล SRTM3 จะอยู่ที่ประมาณ 90 เมตร
ไฟล์ Height มีนามสกุล. HGT และมีการเซ็นชื่อเป็นจำนวนเต็มสองไบต์ ไบต์อยู่ในลำดับโมโตโรล่า "big-endian" ที่มีไบต์ที่สำคัญที่สุดก่อนสามารถอ่านได้โดยตรงโดยระบบเช่น Sun SPARC, Silicon Graphics และคอมพิวเตอร์ Macintosh โดยใช้โปรเซสเซอร์ Power PC DEC Alpha คอมพิวเตอร์ส่วนใหญ่และคอมพิวเตอร์ Macintosh ที่สร้างขึ้นหลังจากปี 2006 ใช้คำสั่ง Intel ("little-endian") ดังนั้นอาจจำเป็นต้องมีการแลกเปลี่ยนไบต์ ความสูงมีหน่วยเป็นเมตรอ้างอิงถึงพิกัด WGS84 / EGM96 Data voids ถูกกำหนดค่า -32768
วิธีการดำเนินการ
สำหรับตำแหน่งของคุณ 50 ° 24'58.888 "N 14 ° 55'11.377" E คุณพบไพ่ที่ถูกต้องแล้ว N50E14.hgt เรามาดูกันว่าพิกเซลตัวไหนที่คุณสนใจละติจูดแรก 50 ° 24'58.888 "N:
24'58.888" = (24 * 60)" + 58.888" = 1498.888"
ส่วนโค้งวินาที หารด้วยสามและปัดเศษเป็นจำนวนเต็มที่ใกล้เคียงที่สุดให้แถวกริด 500 การคำนวณเดียวกันสำหรับผลลัพธ์ลองจิจูดในคอลัมน์กริด 1104
เอกสารฉบับย่อเริ่มขาดข้อมูลเกี่ยวกับวิธีการจัดเรียงแถวและคอลัมน์ในไฟล์ แต่ในเอกสารฉบับเต็มมีการระบุไว้ว่า
ข้อมูลจะถูกจัดเก็บในลำดับหลักของแถว (ข้อมูลทั้งหมดสำหรับแถว 1 ตามด้วยข้อมูลทั้งหมดสำหรับแถว 2 เป็นต้น)
แถวแรกในไฟล์น่าจะอยู่เหนือสุดมากเช่นถ้าเราสนใจแถว 500 จากขอบล่างเราต้องดูแถว
1201 - 500 = 701
จากจุดเริ่มต้นถ้าไฟล์ เซลล์กริดของเราคือจำนวน
(1201 * 700) + 1104 = 841804
จากจุดเริ่มต้นของไฟล์ (เช่นข้าม 700 แถวและใน 701st หนึ่งใช้ตัวอย่าง 1104) สองไบต์ต่อตัวอย่างหมายความว่าเราต้องข้าม 1683606 ไบต์แรกในไฟล์แล้วอ่านสองไบต์เพื่อรับกริดเซลล์ของเรา ข้อมูลมีขนาดใหญ่ซึ่งหมายความว่าคุณต้องสลับทั้งสองไบต์เช่นแพลตฟอร์มของ Intel
โปรแกรมตัวอย่าง
โปรแกรม Python ที่ง่ายต่อการดึงข้อมูลที่ถูกต้องจะมีลักษณะดังนี้ (ดูเอกสารสำหรับการใช้โมดูล struct):
import struct
def get_sample(filename, n, e):
i = 1201 - int(round(n / 3, 0))
j = int(round(e / 3, 0))
with open(filename, "rb") as f:
f.seek(((i - 1) * 1201 + (j - 1)) * 2) # go to the right spot,
buf = f.read(2) # read two bytes and convert them:
val = struct.unpack('>h', buf) # ">h" is a signed two byte integer
if not val == -32768: # the not-a-valid-sample value
return val
else:
return None
if __name__ == "__main__":
n = 24 * 60 + 58.888
e = 55 * 60 + 11.377
tile = "N50E14.hgt" # Or some magic to figure it out from position
print get_sample(tile, n, e)
โปรดทราบว่าการดึงข้อมูลที่มีประสิทธิภาพจะต้องดูซับซ้อนกว่าเล็กน้อย (เช่นไม่เปิดไฟล์สำหรับตัวอย่างแต่ละตัวอย่าง)
ทางเลือก
คุณยังสามารถใช้โปรแกรมที่สามารถอ่านไฟล์. hgt ได้จากกล่อง แต่นั่นน่าเบื่อ
.hgtรูปแบบไฟล์ในเอกสารคู่มือSRTMแต่คำตอบทีละขั้นตอนที่เฉพาะเจาะจงขึ้นอยู่กับซอฟต์แวร์ที่คุณมี