MODIS. MOD02/03: связать координаты и время пролёта

[Эдуард Казаков]

Всем привет.
Может кто подскажет. Не могу разобраться с тем, как автоматически получить данные MOD021KM (и соответственно привязочные данные MOD03) на нужную территорию в нужный день. Они ещё не приведены в sinusoidal grid и на насовском FTP-сервере именуются по времени съемки. Время съемки одной и той же территории в разные дни не совпадает. Таким образом я не могу определить, имея координату и дату, какой файл мне нужен.
Вот на этом сайте можно найти гифки с нарисованной траекторией и подписанными отметками времени - когда и где снимал Terra. Отлично. Но можно ли найти эту же информацию, но в интерпретируемом виде, текстовом, например, "суточное время - координаты центра сцены"? Чтобы можно было автоматически её парсить.

Нужны именно сырые данные radiance по всем каналам, в продуктах с более высоким уровнем обработки такого нет.

[Максим Дубинин]

прямо помочь не могу, но расскажу как я это делал примерно 6 лет назад сидя в SDSU.

писал программу на VB которая использовала библиотеку SGP4 (orbital propagation model - орбитальное моделирование) и телеметрию спутников Terra и Aqua для генерации орбит (TLE, данные открыты).

Использовал для проверки именно эти данные с SSEC (благо они были под боком), получал кучу орбитальных треков в шейпах.

Накладывал на нужную мне территорию и смотрел когда именно там будет появляться спутник.

Качал нужные данные с LAADSWeb.

Осталась заготовка статьи (в основном начало про TLE) и можно откопать код. Жаль, что не дописал, похоже что никто так это толком и не решил с тех пор.

Вообще интересно было бы как-нибудь это восстановить, наверняка теперь можно легко наколбасить что-нибудь на питоне для того же самого.

[Эдуард Казаков]

Написал небольшой скрипт на эту тему
orbital_track_shapefile_creator (GitHub)

Мою подзадачу пока решает, можно что-нибудь добавлять.

Генерирует точечный шейп (WGS84) по набору исходных данных:
Год, месяц, день, шаг отслеживания положения спутника в минутах, две строки TLE, путь до результирующего файла.

Точки в атрибутах хранят широту, долготу и время UTC.

Можно использовать из командной строки, например:

Код: Выделить всё

python.exe orbital_track_shapefile_creator.py -y 2014 -m 6 -d 21 -s 5 -u "1 25994U 99068A 14172.15370853 .00000261 00000-0 68078-4 0 9990" -l "2 25994 098.2112 246.8346 0001762 090.7383 269.4028 14.57122372771643" -o "E:\track1.shp"

Можно обращаться к функции в коде:

Код: Выделить всё

import orbital_track_shapefile_creator

orbital_track_shapefile_creator.create_orbital_track_shapefile_for_day (2014,6,21,5,'1 25994U 99068A 14172.15370853 .00000261 00000-0 68078-4 0 9990', '2 25994 098.2112 246.8346 0001762 090.7383 269.4028 14.57122372771643','E:\track.shp')

Ну и дальше делать что вам нужно: искать ближайшую точку и т.д.

Использую pyorbital и shapefile, положил вместе со скриптом, чтобы не заморачиваться с установкой.

В итоге у этого скрипта результат это такой точечный шейп (даже совпадает с картинкой SSEC):