Форумы GIS-Lab.info

Всем привет.
Может кто подскажет. Не могу разобраться с тем, как автоматически получить данные 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.

Можно использовать из командной строки, например:
Можно обращаться к функции в коде:

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

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

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