Формат даних
Я сприйматиму це як невелику вправу програмування зчитувача даних. Погляньте на документацію :
Дані SRTM розподіляються у двох рівнях: SRTM1 (для США та її територій та володінь) з даними, відібраними з інтервалом дуги-секунди за широтою та довготою, та SRTM3 (для світу), відібраними у три дуги-секунди.
Дані поділяються на "по одному градусу широти та довготи" в "географічній" проекції, тобто растрове представлення з рівними інтервалами широти та довготи взагалі без проекції, але просто маніпулювати та мозаїчно.
Назви файлів відносяться до широти та довготи нижнього лівого кута плитки - наприклад, N37W105 має нижній лівий кут на 37 градусах північної широти та 105 градусів на західній довготі. Якщо точніше, ці координати відносяться до геометричного центру лівого нижнього пікселя, який у випадку даних SRTM3 становитиме близько 90 метрів.
Файли висоти мають розширення .HGT і підписані двома байтовими цілими числами. Байти перебувають у Motorola в "великому ендіанському" порядку з першим найважливішим байтом, який легко читається такими системами, як Sun SPARC, Silicon Graphics та Macintosh, що використовують процесори Power PC. DEC Alpha, більшість комп'ютерів та комп'ютерів Macintosh, побудовані після 2006 року, використовують замовлення Intel ("малопомітний"), тому може знадобитися деяка заміна байтів. Висоти в метрах посилаються на геоїд WGS84 / EGM96. Пустоті даних присвоюється значення -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 рядків, а в 701-му взяти зразок 1104). Два байти на зразок означає, що ми повинні пропустити перші 1683606 байтів у файл, а потім прочитати два байти, щоб отримати нашу комірку сітки. Дані є big-endian, що означає, що вам потрібно поміняти два байти, наприклад, на платформах 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 є деякі примітки щодо формату файлу , але конкретна покрокова відповідь залежить від наявного програмного забезпечення.