Это команда snaphu, которую можно запустить в бесплатном хостинг-провайдере OnWorks, используя одну из наших многочисленных бесплатных онлайн-рабочих станций, таких как Ubuntu Online, Fedora Online, онлайн-эмулятор Windows или онлайн-эмулятор MAC OS.
ПРОГРАММА:
ИМЯ
snaphu - алгоритм развёртки фазы для SAR-интерферометрии
СИНТАКСИС
снэпху [параметры] [infile] [длина строки] [параметры]
ОПИСАНИЕ
снэпху - это sстатистическая стоимость nсетевой поток aалгоритм для phаза unwrapping. Учитывая ввод
интерферограмма и другие наблюдаемые данные, снэпху пытается вычислить конгруэнтную фазу-
развернутые решения, которые максимально вероятны в приближенном a задний смысл.
Процедура решателя алгоритма основана на оптимизации сети. По умолчанию, снэпху предполагает,
что его входом является поверхность измерения интерферограммы радара с синтезированной апертурой (SAR)
топография. Измерения деформации предполагаются, если -d предоставляется вариант. Гладкий; плавный,
общие данные предполагаются, если -s предоставляется вариант.
Только документы на этой странице руководства снэпхусинтаксис и использование. Его теоретические основы:
обсуждается в ссылках, цитируемых ниже.
Наиболее распространенные входные параметры могут быть указаны в командной строке, в то время как многие другие
параметры twiddle обрабатываются через -f файлы параметров и конфигурации. В самом
по крайней мере, необходимо указать имя входного файла с фазовой оболочкой и длину его строки.
Дальность должна увеличиваться вправо на интерферограмме, и фаза плоской земли
рампа должна быть удалена с входной интерферограммы перед снэпху запущен. Для деформации
На интерферограммах также должны быть удалены фазовые вариации из-за топографии.
За исключением имени входного файла и длины строки, все входные параметры принимают значения по умолчанию.
значения, если не указаны. Однако по возможности эти параметры следует настраивать.
поскольку точность решения зависит от того, насколько хорошо статистика оценки
проблемы смоделированы. Чтобы избежать некачественных решений, пользователям настоятельно рекомендуется
предоставить свои наилучшие оценки соответствующих параметров проблемы. Параметры устанавливаются в
порядок, в котором они указаны в командной строке, поэтому несколько файлов конфигурации или
могут быть заданы параметры, при этом более поздние значения имеют приоритет перед более ранними.
Допустимые форматы файлов подробно описаны ниже. Формат по умолчанию для входного файла:
COMPLEX_DATA, но можно использовать любой из описанных форматов. Если любой из
Используются форматы ALT_LINE_DATA или ALT_SAMPLE_DATA, величина и фаза (в радианах)
интерферограмма должна быть в первом и втором каналах файла соответственно.
Если используется формат FLOAT_DATA, входной файл должен содержать только фазу
интерферограмма (в радианах); величина может быть передана с -m опцию.
ДОПОЛНИТЕЛЬНЫЕ ОПЦИИ
-a amp-файл
Считать данные яркости из файла amp-файл. Файл должен содержать амплитуды
(не мощности) двух отдельных изображений SAR, образующих интерферограмму, если
используются форматы ALT_SAMPLE_DATA (по умолчанию) или ALT_LINE_DATA. Он должен содержать
среднее значение этих двух изображений, если используется формат FLOAT_DATA. Если (1)
доступны амплитуды обоих изображений, (2) величина интерферограммы также
доступны, и (3) -c не используется, то оценка согласованности
автоматически формируется из имеющихся данных. Количество использованных образов для этого
смету можно задать в файле конфигурации. Если нет данных об амплитуде или мощности
задано, то величина входной интерферограммы используется как среднее
амплитуды, и оценка когерентности не формируется. Обратите внимание, что величина
интерферограмма не равна средней амплитуде изображений РСА. В
данные об амплитуде должны быть в той же системе единиц, что и для входных данных.
интерферограмму, а также зарегистрированную к ней.
-A pwrфайл
Как и в случае -a вариант, за исключением того, что данные в указанном файле предполагается
представляют мощность двух отдельных изображений SAR.
-b Бперп
Для режима топографии используйте Бперп (десятичное значение в метрах) как значение
перпендикулярная составляющая интерферометрической базовой линии. Знак определяется таким
которые Бперп отрицателен, если развернутая фаза увеличивается с высотой. К
по умолчанию предполагается режим повторного прохода или пинг-понг; для передачи с одной антенной
данные, значение Бперп следует уменьшить вдвое или установить режим передачи
соответственно в файле конфигурации (см. -f вариант). Базовое значение
используется только в режиме топографии.
-c Corrfile
Прочитать данные корреляции из файла Corrfile. Данные корреляции должны быть
того же размера, что и входная интерферограмма, и зарегистрированы в ней. Следовательно, сырой
оценка корреляции может потребовать повышения дискретизации, если она включает в себя больше взглядов, чем
интерферограмма. Если -c вариант не дан, формируется оценка согласованности
если возможно, из имеющихся данных. В противном случае единообразная согласованность по умолчанию
предполагается для всей интерферограммы. Если ALT_LINE_DATA (по умолчанию) или
Используются форматы ALT_SAMPLE_DATA, данные корреляции должны быть во вторых данных
канал файла; первый канал игнорируется. Формат FLOAT_DATA также может
использоваться. Значения корреляции должны быть от нуля до единицы включительно.
-d Запуск в режиме деформации. Статистика проблемы и результирующие функции затрат:
основанный на предположении, что истинная развернутая фаза представляет собой поверхность
смещение, а не высота.
-e файл оценки
Сглаживание с использованием развернутой оценки фазы в файле файл оценки. Оценка
вычитается из входной интерферограммы перед разворачиванием и вставляется обратно
в решение непосредственно перед записью вывода. Оценка также влияет на
используемые функции затрат, поскольку вычитание константы из случайной величины сдвигает
функция плотности вероятности случайной величины. Если форматы ALT_LINE_DATA
(по умолчанию) или ALT_SAMPLE_DATA, развернутая оценка (в радианах) должна
находиться во втором канале данных файла; первый канал игнорируется. В
Также можно использовать формат FLOAT_DATA.
-f файл конфигурации
Прочитать параметры конфигурации из файла файл конфигурации. Файл анализируется строкой
строка для пар "ключ-значение". Файлы конфигурации шаблона включены в
снэпху исходный код: snaphu.conf.полный содержит все допустимые пары ключ-значение;
snaphu.conf.brief содержит наиболее важные параметры. Строки, не начинающиеся с
буквенно-цифровые символы рассматриваются как строки комментариев. Параметры командной строки
указано после -f переопределит параметры, указанные в файл конфигурации и порок
наоборот. В -f опция может быть указана несколько раз с разной конфигурацией
файлы, с параметрами в указанных позже файлах, переопределяющими параметры в более ранних.
-g файл-маска
Вырастите маску связного компонента для развернутого решения и запишите маску на
файла файл-маска. Связанный компонент - это область пикселей в решении.
что, как полагают, было развернуто в относительном, внутренне самосогласованном
способом в соответствии с используемыми статистическими затратами. Регионы меньше
предварительно выбранный порог маскируется. Параметры для этой опции могут быть установлены в
конфигурационный файл. Файл связанного компонента состоит из беззнаковых
символы, причем все пиксели одного значения принадлежат одному и тому же подключенному
компонент и ноль, соответствующие замаскированным пикселям.
-G файл-маска
Увеличьте маску связного компонента (см. -g option) для массива входных данных,
предполагая, что он уже развернут, и запишите маску в файл файл-маска.
Статистические функции затрат вычисляются для формирования маски, но новый развернутый
решение не вычислено.
-h Распечатайте справочное сообщение с описанием параметров командной строки и выйдите.
-i Запуск в режиме только инициализации. Как обычно, снэпху использует приблизительный минимум
алгоритм связующего дерева (MST) или алгоритм потока минимальной стоимости (MCF) для генерации
инициализация его итеративного модифицированного решателя сетевой симплексной системы. Если -i is
задано, инициализация записывается в выход, и программа завершается без
запускает итеративный решатель.
-l журнальный файл
Записывать все параметры времени выполнения и некоторую другую информацию о среде в
указанный файл. Файл журнала - это текстовый файл того же формата, что и конфигурация.
.
-m магнитный файл
Считать данные магнитуды интерферограммы из указанного файла. Эта опция полезна
в основном, если входной файл обернутой фазы задан как набор значений реальной фазы
а не комплексные значения интерферограммы. Величина интерферограммы используется для
сформировать оценку когерентности, если также будут даны соответствующие данные об амплитуде. В
формат файла по умолчанию - FLOAT_DATA. Если форматы ALT_LINE_DATA или ALT_SAMPLE_DATA
используются, величина должна быть в первом канале данных файла; секунда
канал игнорируется. Если используется формат COMPLEX_DATA, информация о фазе
игнорируются.
-n Запускайте в режиме без статистических затрат. Если -i or -p даны варианты, снэпху предусматривает
не использовать статистические затраты. Информация из файла веса (-w вариант) все равно будет
использоваться, если дано.
-o Outfile
Запишите распакованный вывод в файл с именем Outfile. Если форматы файлов
Используются ALT_LINE_DATA (по умолчанию) или ALT_SAMPLE_DATA, записывается развернутая фаза
во второй канал данных, а величина интерферограммы записывается в
первый канал. Также можно использовать формат FLOAT_DATA.
-p ценностное
Запустите в режиме Lp-norm с p =ценностное , Где ценностное является неотрицательной десятичной дробью. Вместо того
статистические функции затрат, программа использует функции стоимости Lp со статистически
на основе весов (если -n также дано). Решения по-прежнему всегда совпадают.
Более того, конгруэнтность обеспечивается в рамках процедуры решателя, а не после
шаг обработки оптимизации. Следовательно, если p = 2, например, стоимость наименьших квадратов
функции, но решение, вероятно, будет более точным, чем одно
генерируется на основе алгоритма наименьших квадратов на основе преобразования.
-q Запускать только в режиме количественной оценки. Предполагается, что входные данные уже развернуты, и
рассчитывается и распечатывается общая стоимость этого решения. Развернутая фаза
завернутый в предположении соответствия для расчета стоимости. Ошибки округления могут ограничивать
точность количественной оценки стоимости. Увидеть -u вариант для допустимого файла
форматов.
-s Запускайте в режиме плавного решения. Статистика проблемы и результирующие функции затрат
основаны на предположении, что истинная развернутая фаза представляет собой общий
поверхность без разрывов. Это то же самое, что и режим деформации с
Параметр DEFOMAX равен нулю.
-t Запуск в режиме топографии. Статистика проблемы и результирующие функции затрат:
на основе предположения, что истинная развернутая фаза представляет собой отметку поверхности.
Это значение по умолчанию.
-u Предположим, что входной файл развернут, а не упакован. Алгоритм делает
итеративные улучшения этого решения вместо использования процедуры инициализации.
Входной файл может иметь формат ALT_LINE_DATA (по умолчанию) или ALT_SAMPLE_DATA;
величина интерферограммы должна быть в первом канале данных, а развернутый
фаза должна быть во втором канале данных. Формат FLOAT_DATA также может быть
используемый.
-v Запустить в подробном режиме. Дополнительная информация о прогрессе алгоритма выводится на
стандартный вывод.
-w файл веса
Чтение внешних скалярных весов из файла файл веса. Веса, который должен быть
положительные короткие целые числа применяются к любым используемым функциям стоимости. Там
- одно значение веса для каждой дуги в сети, поэтому файл веса должен быть
объединение растровых весов дуги горизонтального и вертикального потока. Таким образом, для
интерферограмма из N строк по M столбцов, файл веса будет состоять из растеризованного (N-1)
на массив M, за которым следует растеризованный массив N на (M-1) коротких целочисленных данных. Этот
вариант не очень хорошо протестирован.
--аа ampfile1 ampfile2
Данные амплитуды считываются из указанных файлов. Данные двух человек
Предполагается, что изображения РСА, образующие интерферограмму, хранятся отдельно в файлах.
ampfile1 и ampfile2. Эти файлы должны быть в формате FLOAT_DATA. Этот
вариант аналогичен -a опцию.
--АА pwrfile1 pwrfile2
Как и в случае --аа вариант, но данные о питании читаются из указанных файлов.
- собрать имя_директории
Соберите временные файлы тайлового режима в указанном каталоге. Самый
параметры конфигурации (из командной строки и любых файлов конфигурации) должны быть
указано. Эта опция полезна, если пользователь желает изменить сборку плитки.
параметры без повторного разворачивания отдельных плиток.
--авторское право, --Информация
Распечатайте уведомление об авторских правах на программное обеспечение и информацию об отчете об ошибке, затем выйдите.
--costinfile файл стоимости
Чтение массивов статистических затрат из файла файл стоимости. Этот файл должен быть в формате
написано --costoutfile вариант. Файл стоимости не контролирует, снэпху
работает в режиме топографии, деформации или плавного раствора; последние два должны быть
указано явно, даже если файл стоимости был сгенерирован во время работы в этих режимах.
--costoutfile файл стоимости
Записать массивы статистических затрат в файл файл стоимости. Эту опцию можно использовать с
--costinfile возможность сэкономить время на генерацию статистических затрат, если то же самое
затраты используются многократно.
--отлаживать, --dumpall
Сбрасывать в файлы всевозможные промежуточные массивы.
--мст Используйте алгоритм минимального остовного дерева (MST) для инициализации. Это
по умолчанию.
--mcf Используйте алгоритм потока минимальной стоимости (MCF) для инициализации. Решатель cs2 от
Гольдберга и Черкасского. Модифицированный сетевой симплексный решатель в режиме L1
может дать другие результаты, чем решатель cs2, хотя в принципе оба должны быть
L1 оптимальный.
--nproc n
Используйте n параллельные процессы в режиме плитки. Программа создает новый процесс для
каждую плитку, чтобы плитки можно было развернуть параллельно; в большинстве n процессы будут работать
одновременно. Разветвление выполняется до чтения данных. Стандартные выходные потоки
дочерние процессы направляются в файлы журнала во временном каталоге плиток.
--кусок Первая строка первый столбец Nrow нкол
Прочтите и разверните только часть или часть входной интерферограммы. Читаемый фрагмент
это Nrow by нкол прямоугольник, верхний левый угол которого - пиксель в строке Первая строка
и колонка первый столбец (проиндексировано с 1). Все входные файлы (например, амплитуда,
согласованность и т. д.), как предполагается, имеют тот же размер, что и входной файл фазы. Все
выходные файлы Nrow by нкол.
--плитка Ntilerow Нтилкол рядовврлп Colovrlp
Разверните интерферограмму в тайловом режиме. Интерферограмма разбита на
Ntilerow by Нтилкол плитки, каждая из которых разворачивается независимо. Плитки
перекрывать рядовврлп и Colovrlp пикселей в строках и столбцах. В
плитки затем сегментируются на надежные области на основе функций стоимости, а
регионы собраны заново. Программа создает подкаталог для временных файлов в
каталог конечного выходного файла. Эта опция в настоящее время включена только
для статистических функций затрат.
ФАЙЛОВ ФОРМАТЫ
Форматы входных файлов могут быть указаны в файле конфигурации. Все эти форматы
состоят из растровых данных с плавающей запятой одинарной точности (с плавающей запятой, вещественное * 4 или комплексное * 8)
типы в собственном порядке байтов платформы. Данные читаются построчно (поперек, затем вниз).
Независимо от формата файла, все массивы входных данных должны иметь одинаковое количество
образцы по ширине и глубине и должны быть совмещены друг с другом. Обратите внимание, что вес
файлы и файлы стоимости имеют собственные форматы. Допустимые форматы для других файлов данных
описаны ниже.
КОМПЛЕКС_ДАННЫХ
Чередующиеся поплавки соответствуют реальной (синфазной) и мнимой (квадратурной)
компоненты сложных выборок данных. Указанная длина строки должна быть числом
комплексных отсчетов (пары реальных и мнимых отсчетов) на строку.
ALT_LINE_DATA
Чередующиеся строки (строки) данных соответствуют строкам чисто реальных данных из двух
отдельные массивы. Первый массив часто представляет собой величину интерферограммы, а
вторым может быть развернутая фаза, когерентность и т. д. Это также иногда называют
хгт или формат с чередованием строк.
ALT_SAMPLE_DATA
Чередующиеся выборки соответствуют чисто реальным выборкам из двух отдельных массивов.
Этот формат иногда используется для амплитуд двух изображений SAR.
FLOAT_DATA
Файл содержит данные только для одного канала или массива, и эти данные являются чисто реальными.
ПРИМЕРЫ
Разверните завернутую топографическую интерферограмму под названием `` завернутый файл '', длина строки которого равна
1024 сложных семпла (вывод будет записан в файл, имя которого скомпилировано в
программа):
snaphu обернутый файл 1024
Разверните тот же файл, что и выше, но используйте информацию о яркости из файла ampfile,
установите перпендикулярную базовую линию на -165 м на средней ширине полосы и поместите результат в файл
называется `` развернутый файл '' (данные согласованности генерируются автоматически, если `` обернутый файл ''
содержит сложные данные, а файл ampfile содержит данные об амплитуде обоих изображений РСА):
snaphu wrapfile 1024 -a ampfile \
-b -165 -o распакованный файл
Разверните интерферограмму, как указано выше, но прочтите информацию о корреляции из файла.
`` corrfile '' вместо генерации его из данных интерферограммы и амплитуды:
snaphu wrapfile 1024 -a ampfile -c corrfile \
-b -165 -o распакованный файл
Следующее эквивалентно предыдущему примеру, но входные параметры считываются из
файл конфигурации, и отображается подробный вывод:
cat> configfile
# Это строка комментария, которая будет проигнорирована
AMPFILE файл усилителя
КОРРФАЙЛ
БПЕРП -165
OUTFILE развернутый файл
snaphu -v -f файл конфигурации обернутый файл 1024
Разверните ту же интерферограмму, но используйте только инициализацию MST (со скалярной
статистические веса) и запишите результат в `` mstfile '':
snaphu -f файл конфигурации -i файл-оболочка 1024 -o mstfile
Прочтите развернутые данные в `` mstfile '' и используйте их как инициализацию измененного
сетевой симплексный решатель:
snaphu -f configfile -u mstfile 1024 -o распакованный файл
Обратите внимание, что в двух предыдущих примерах имя выходного файла в файле конфигурации:
заменяется тем, что указано в командной строке. Две предыдущие команды вместе находятся в
принцип эквивалентен предыдущему, хотя ошибки округления в межфазном потоке
преобразования могут вызывать незначительные различия
Разверните интерферограмму, как указано выше, но используйте алгоритм MCF для инициализации:
snaphu -f файл конфигурации wrappedfile 1024 --mcf
Разверните интерферограмму еще раз, но сначала сгладьте ее с развернутыми данными в
`` estfile '', затем повторно вставьте вычтенную фазу после разворачивания:
snaphu -f файл конфигурации wrappedfile 1024 -e файл конфигурации
Далее предполагается, что на обернутой входной интерферограмме измеряется деформация, а не
топография. Разверните интерферограмму с данными корреляции:
snaphu -d обернутый файл 1024 -c corrfile
Разверните входную интерферограмму, минимизируя невзвешенную конгруэнтную норму L2:
snaphu -p 2 -n обернутый файл 1024
Разверните интерферограмму в виде набора плиток размером три на четыре, перекрывающихся на 30 пикселей, с
указанный файл конфигурации с использованием двух процессоров:
snaphu wrappedfile 1024 -f файл конфигурации \
--tile 3 4 30 30 --nproc 2
HINTS И СОВЕТЫ
Программа может напечатать предупреждающее сообщение о сокращении затрат, чтобы избежать переполнения. Если
слишком много затрат сокращено, значение РАСХОДНОГО МАТЕРИАЛА может потребоваться уменьшить в
файл конфигурации (через -f вариант). Если программа выводит предупреждающее сообщение о
неожиданное увеличение общей стоимости решения, это показатель того, что слишком много затрат
обрезаны. Обычно это нормально, если сокращаются лишь некоторые расходы.
В режиме топографии, если развернутый результат содержит слишком много разрывов, попробуйте
увеличение значения LAYMINEI или уменьшение значения LAYCONST. Первый
определяет нормированный порог интенсивности для стоянки, и последний является относительным
вероятность остановки. Если по азимуту слишком много разрывов, попробуйте
уменьшение значения AZDZFACTOR, которое влияет на соотношение стоимости азимута и дальности. Если
базовая линия не известна, попробуйте угадать ее и убедитесь, что ее знак правильный. Указать
геометрические параметры изображения SAR, насколько это возможно. По умолчанию используются данные ERS.
с пятью взглядами, снятыми по азимуту.
В режиме деформации, если развернутый результат содержит слишком много разрывов, попробуйте
увеличение значения DEFOTHRESHFACTOR или уменьшение значения DEFOCONST. Если
смещение поверхности меняется медленно, и истинных разрывов не ожидается вообще,
DEFOMAX_CYCLE можно установить равным нулю. Это поведение также вызывается с помощью -s вариант.
результирующие функции затрат будут аналогичны функциям затрат L2, взвешенным по корреляции, хотя
первые не обязательно центрируются на обернутых градиентах. Конгруэнтность все еще
применяется во время, а не после оптимизации.
Программа может быть запущена только в режиме инициализации (-i) режим для быстрой очистки MST или MCF
решений.
СИГНАЛЫ
После запуска итеративного решателя снэпху улавливает прерывание (INT) и зависание (HUP)
сигналы. При получении прерывания, например, если пользователь нажимает Ctrl-C, программа
завершает малую итерацию, выгружает свое текущее решение на выходе и выходит. Если
второе прерывание выдается после первого (пойманного) прерывания, программа завершается
немедленно. Если получен сигнал зависания, программа сбрасывает свое текущее решение, а затем
продолжает нормально работать.
ВЫХОД статус
При успешном завершении программа выходит с кодом 0. Ошибки приводят к коду выхода 1.
Используйте snaphu в Интернете с помощью сервисов onworks.net