Кривая-укорачивание потока - Curve-shortening flow

Сходимость выпуклой кривой к окружности при потоке укорачивания кривой. Внутренние кривые (более светлый цвет) представляют собой плавные версии внешних кривых. Временные шаги между кривыми неодинаковы.

В математике поток сокращения кривой - это процесс, который изменяет гладкую кривую на евклидовой плоскости , перемещая ее точки перпендикулярно кривой со скоростью, пропорциональной кривизне . Поток с сокращением кривой является примером геометрического потока и одномерным случаем потока средней кривизны . Другие названия того же процесса включают евклидово сокращение потока , геометрический тепловой поток и эволюцию длины дуги .

Поскольку точки любой гладкой простой замкнутой кривой перемещаются таким образом, кривая остается простой и гладкой. Он теряет площадь с постоянной скоростью, а его периметр уменьшается настолько быстро, насколько это возможно для любой непрерывной эволюции кривой. Если кривая невыпуклая, ее полная абсолютная кривизна монотонно уменьшается, пока не станет выпуклой. После выпуклости изопериметрическое отношение кривой уменьшается по мере того, как кривая сходится к круглой форме, прежде чем схлопнуться в единственную точку сингулярности. Если две непересекающиеся простые гладкие замкнутые кривые эволюционируют, они остаются непересекающимися, пока одна из них не схлопнется в точку. Круг - единственная простая замкнутая кривая, которая сохраняет свою форму в потоке, укорачивающем кривую, но некоторые кривые, которые пересекаются друг с другом или имеют бесконечную длину, сохраняют свою форму, в том числе мрачная кривая жнеца, бесконечная кривая, которая перемещается вверх, и спирали, которые вращаются. при этом оставаясь того же размера и формы.

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

Изначально течение укорачивания кривой исследовалось как модель отжига металлических листов. Позже он был применен в анализе изображений, чтобы дать многомасштабное представление форм. Он также может моделировать системы реакции-диффузии и поведение клеточных автоматов . Поток сокращения кривой можно использовать для поиска замкнутых геодезических на римановых многообразиях и в качестве модели поведения многомерных потоков.

Определения

Поток представляет собой процесс , в котором точка пространства непрерывно изменяет свои места или свойство с течением времени. Более конкретно, в одномерном геометрическом потоке, таком как поток, укорачивающий кривую, точки, через которые проходит поток, принадлежат кривой , и что изменяется, так это форма кривой, ее вложение в евклидову плоскость, определяемую местоположением каждого из них. своих точек. В потоке сокращения кривой каждая точка кривой перемещается в направлении вектора нормали к кривой со скоростью, пропорциональной кривизне . Для развивающейся кривой, представленной двухпараметрической функцией C ( s , t ), где s параметризует длину дуги вдоль кривой, а t параметризует время в эволюции кривой, поток сокращения кривой может быть описан параболическим частным дифференциальное уравнение

форма уравнения теплопроводности , где κ - кривизна, а n - единичный вектор нормали.

Поскольку компоненты этого уравнения, длина дуги, кривизна и время, не зависят от сдвигов и поворотов евклидовой плоскости, из этого следует, что поток, определяемый этим уравнением, инвариантен относительно сдвигов и вращений (или, точнее, эквивариантен ). . Если масштабировать самолет с постоянным коэффициентом расширения, поток остается практически неизменным, но замедляется или ускоряется с тем же коэффициентом.

Негладкие кривые

Чтобы поток был хорошо определен, данная кривая должна быть достаточно гладкой, чтобы иметь непрерывную кривизну. Однако, как только поток начинается, кривая становится аналитической и остается таковой до тех пор, пока не достигнет сингулярности, в которой кривизна резко возрастет. Для гладкой кривой без пересечений единственная возможная сингулярность возникает, когда кривая схлопывается в точку, но погруженные кривые могут иметь особенности других типов. В таких случаях с некоторой осторожностью можно продолжить обтекание этих сингулярностей до тех пор, пока вся кривая не сузится до одной точки.

Для простой замкнутой кривой с использованием расширения потока на негладкие кривые на основе метода установки уровня есть только две возможности. Кривые с нулевой мерой Лебега (включая все многоугольники и кусочно-гладкие кривые) мгновенно превращаются в гладкие кривые, после чего они развиваются, как любая гладкая кривая. Однако кривые Осгуда с ненулевой мерой вместо этого немедленно превращаются в топологическое кольцо с ненулевой площадью и гладкими границами. В синусоида тополога в примере , который мгновенно становится гладким, несмотря даже не будучи локально связно ; Примеры, подобные этому, показывают, что обратная эволюция потока, укорачивающего кривую, может привести кривые с хорошим поведением к сложным сингулярностям за конечный промежуток времени.

Неевклидовы поверхности

Поток сокращения кривой и многие результаты о потоке сокращения кривой можно обобщить с евклидовой плоскости на любое двумерное риманово многообразие . Чтобы избежать дополнительных типов особенностей, важно, чтобы многообразие было выпуклым на бесконечности ; это означает, что каждый компакт имеет компактную выпуклую оболочку , как определено с помощью геодезической выпуклости . Поток, сокращающий кривую, не может вызвать отклонение кривой от ее выпуклой оболочки, поэтому это условие не позволяет частям кривой достигать границы многообразия.

Космические кривые

Поток сокращения кривой также изучался для кривых в трехмерном евклидовом пространстве . Вектор нормали в этом случае может быть определен (как на плоскости) как производная касательного вектора по длине дуги, нормированная как единичный вектор; это один из компонентов системы Френе – Серре . Он плохо определен в точках нулевой кривизны, но произведение кривизны и вектора нормали остается четко определенным в этих точках, что позволяет определить поток сокращения кривой. Кривые в пространстве могут пересекать друг друга или сами себя в соответствии с этим потоком, и течение может приводить к сингулярностям на кривых; каждая особенность асимптотична плоскости. Однако известно, что сферические кривые и кривые, которые могут быть ортогонально спроецированы в правильную выпуклую плоскую кривую, остаются простыми. Поток сокращения кривой для пространственных кривых использовался как способ определения обтекания сингулярностей на плоских кривых.

За пределами кривых

Можно расширить определение потока на более общие входные данные, чем кривые, например, используя выпрямляемые варифолды или метод установки уровня . Однако эти расширенные определения могут позволить частям кривых мгновенно исчезнуть или разрастаться до наборов ненулевой площади.

Для сетей кривых продолжение потока, укорачивающего кривую, за сингулярность может привести к неоднозначности или полноте.

Обычно изучаемый вариант проблемы включает сети непересекающихся внутри гладких кривых с соединениями, в которых пересекаются три или более кривых. Когда все стыки имеют ровно три кривые, пересекающиеся под углами 2 π / 3 (те же условия, наблюдаемые в оптимальном дереве Штейнера или двумерной пене из мыльных пузырей ), течение краткосрочного периода становится четко определенным. Однако он может в конечном итоге достичь сингулярного состояния с четырьмя или более кривыми, пересекающимися на стыке, и может быть более одного способа продолжить поток мимо такой сингулярности.

Поведение

Принцип уклонения, радиус и коэффициент растяжения

Если две непересекающиеся гладкие простые замкнутые кривые одновременно подвергаются потоку, сокращающему кривую, они остаются не пересекающимися по мере развития потока. Причина в том, что если две плавные кривые движутся таким образом, что создает пересечение, то во время первого пересечения кривые обязательно будут касаться друг друга, но не пересекаться. Но в такой ситуации кривизна двух кривых в точке касания неизбежно будет раздвигать их, а не сталкивать вместе, образуя пересечение. По той же причине простая замкнутая кривая никогда не может пересечь сама себя. Это явление известно как принцип избегания.

Принцип избегания подразумевает, что любая гладкая замкнутая кривая в конечном итоге должна достичь сингулярности, например точки бесконечной кривизны. Ведь, если данная гладкая кривая C окружена кругом, обе они останутся не пересекающимися, пока существуют обе. Но окружающий круг сжимается под действием потока кривизны, оставаясь круглым, пока не схлопнется, и по принципу избегания C должен оставаться внутри него. Итак, если бы C никогда не достигал сингулярности, он был бы захвачен в единственной точке в момент схлопывания круга, что невозможно для гладкой кривой. Это может быть определено количественно, наблюдая, что радиус наименьшего круга, который охватывает C, должен уменьшаться со скоростью, по крайней мере, так же быстро, как уменьшение радиуса круга, подвергающегося тому же потоку.

Huisken (1998) количественно определяет принцип избегания для одной кривой в терминах отношения между длиной дуги (меньшей из двух дуг) и евклидовым расстоянием между парами точек, иногда называемым коэффициентом растяжения . Он показывает, что коэффициент растяжения строго уменьшается в каждом из своих локальных максимумов, за исключением случая двух концов диаметра окружности, когда коэффициент растяжения постоянен на π . Это свойство монотонности подразумевает принцип избегания, поскольку, если кривая когда-либо коснется самой себя, коэффициент растяжения станет бесконечным в двух точках касания.

Длина

По мере того, как кривая подвергается потоку, сокращающему кривую, ее длина L уменьшается со скоростью, определяемой формулой

где интервал берется по кривой, κ - кривизна, а s - длина дуги вдоль кривой. Подынтегральная функция всегда неотрицательна, и для любой гладкой замкнутой кривой существуют дуги, внутри которых она строго положительна, поэтому длина монотонно уменьшается. В более общем смысле, для любой эволюции кривых, нормальная скорость которых равна f , скорость изменения длины равна

что можно интерпретировать как отрицательный внутренний продукт между данной эволюцией и потоком, укорачивающим кривую. Таким образом, поток, сокращающий кривую, можно описать как градиентный поток по длине, поток, который (локально) уменьшает длину кривой как можно быстрее относительно L 2 нормы потока. Это свойство и дало название потоку, сокращающему кривую.

Площадь

Для простой замкнутой кривой площадь, ограниченная кривой, сжимается с постоянной скоростью 2 π единиц площади в единицу времени, независимо от кривой. Следовательно, полное время, за которое кривая сузится до точки, пропорционально ее площади, независимо от ее исходной формы. Поскольку площадь кривой уменьшается с постоянной скоростью, и (в соответствии с изопериметрическим неравенством ) круг имеет максимально возможную площадь среди простых замкнутых кривых данной длины, из этого следует, что круги - самые медленные кривые, чтобы схлопнуться в точку ниже поток, укорачивающий кривую. Для свертывания всех остальных кривых требуется меньше времени, чем для окружности такой же длины.

Постоянная скорость уменьшения площади - единственный закон сохранения, которому удовлетворяет поток, сокращающий кривую. Это означает, что невозможно выразить «точку схода», в которой кривая в конечном итоге схлопывается, как интеграл по кривой любой функции от ее точек и их производных, потому что такое выражение привело бы к запрещенному второму закону сохранения. Однако, комбинируя постоянную скорость потери площади с принципом избегания, можно доказать, что точка схода всегда находится внутри круга, концентричного с минимальным охватывающим кругом, площадь которого является разницей в площадях между охватывающим кругом и окружностью. заданная кривая.

Полная абсолютная кривизна

Общее абсолютное кривизна гладкой кривой является интегралом от абсолютного значения кривизны по длине дуги кривой,

Его также можно выразить как сумму углов между векторами нормалей в последовательных парах точек перегиба . Это 2 π для выпуклых кривых и больше для невыпуклых кривых, служащее мерой невыпуклости кривой.

Новые точки перегиба не могут быть созданы потоком, укорачивающим кривую. Каждый из углов в представлении полной абсолютной кривизны в виде суммы монотонно уменьшается, за исключением моментов, когда две последовательные точки перегиба достигают того же угла или положения, что и друг друга, и обе исключаются. Следовательно, общая абсолютная кривизна никогда не может увеличиваться по мере развития кривой. Для выпуклых кривых он постоянен при 2 π, а для невыпуклых кривых монотонно убывает.

Теорема Гейджа – Гамильтона – Грейсона.

Если гладкая простая замкнутая кривая подвергается потоку сокращения кривой, она остается гладко вложенной без самопересечений. Со временем он станет выпуклым , и после этого останется выпуклым. По истечении этого времени все точки кривой переместятся внутрь, и форма кривой будет сходиться к кругу, так как вся кривая сузится до одной точки. Такое поведение иногда резюмируют, говоря, что каждая простая замкнутая кривая сжимается до «круглой точки».

Этот результат принадлежит Майклу Гейджу , Ричарду С. Гамильтону и Мэтью Грейсону. Гейдж ( 1983 , 1984 ) доказал сходимость к окружности для выпуклых кривых, стягивающихся в точку. В частности, Гейдж показал, что изопериметрическое отношение (отношение квадрата длины кривой к площади, число, равное 4 π для круга и больше для любой другой выпуклой кривой) уменьшается монотонно и быстро. Гейдж и Гамильтон (1986) доказали, что все гладкие выпуклые кривые в конечном итоге стягиваются в точку, не образуя никаких других особенностей, а Грейсон (1987) доказал, что каждая невыпуклая кривая в конечном итоге станет выпуклой. Эндрюс и Брайан (2011) предоставляют более простое доказательство результата Грейсона, основанное на монотонности фактора растяжения.

Предельная форма для всех сетей из двух коллинеарных лучей и двух кривых, соединяющих концы двух лучей. Центральная линза имеет форму vesica piscis .

Подобные результаты могут быть распространены с замкнутых кривых на неограниченные кривые, удовлетворяющие локальному условию Липшица . Для таких кривых, если обе стороны кривой имеют бесконечную площадь, то полученная кривая остается гладкой и без сингулярностей все время. Однако, если одна сторона неограниченной кривой имеет конечную площадь, а кривая имеет конечную общую абсолютную кривизну, то ее эволюция достигает сингулярности во времени, пропорциональной площади на стороне конечной площади кривой, с неограниченной кривизной вблизи сингулярности. . Для кривых, которые являются графиками достаточно хороших функций, асимптотических по лучу в каждом направлении, решение сходится по форме к уникальной форме, которая является асимптотической по отношению к тем же лучам. Для сетей, образованных двумя непересекающимися лучами на одной прямой, вместе с двумя гладкими кривыми, соединяющими концы двух лучей, имеет место аналог теоремы Гейджа – Гамильтона – Грейсона, согласно которой область между двумя кривыми становится выпуклой, а затем сходится к форме vesica piscis .

Особенности самопересечения кривых

Кривые, которые имеют самопересечение, могут достигать сингулярностей, прежде чем сузиться до точки. Например, если лемниската (любая плавная погруженная кривая с одним пересечением, напоминающая цифру 8 или символ бесконечности ) имеет неравные площади в двух долях, то в конечном итоге меньшая доля схлопнется в точку. Однако, если два лепестка имеют равные площади, то они будут оставаться равными на протяжении всей эволюции кривой, а изопериметрическое соотношение будет расходиться по мере того, как кривая сжимается до сингулярности.

Когда локально выпуклая кривая самопересечения приближается к сингулярности, когда одна из ее петель сжимается, она либо сжимается самоподобным образом, либо асимптотически приближается к кривой мрачного жатки (описанной ниже) по мере сжатия. Когда петля схлопывается до сингулярности, потеря полной абсолютной кривизны составляет не менее 2 π или ровно π .

О римановых многообразиях

На римановом многообразии любая гладкая простая замкнутая кривая останется гладкой и простой по мере развития, как и в евклидовом случае. Он либо схлопнется до точки за конечный промежуток времени, либо навсегда останется гладким и простым. В последнем случае кривая обязательно сходится к замкнутой геодезической поверхности.

Погруженные кривые на римановых многообразиях с конечным числом самопересечений становятся самокасающимися только в дискретный набор моментов времени, в каждый из которых они теряют пересечение. Как следствие, количество точек самопересечения не увеличивается.

Сокращение кривой на сфере можно использовать как часть доказательства теоремы о теннисном мячике . Эта теорема утверждает, что каждая гладкая простая замкнутая кривая на сфере, которая делит поверхность сферы на две равные области (как шов теннисного мяча ), должна иметь как минимум четыре точки перегиба . Доказательство основано на наблюдении, что сокращение кривой сохраняет свойства гладкости и деления площади пополам, а не увеличивает количество точек перегиба. Следовательно, это позволяет свести проблему к проблеме для кривых, близких к предельной форме сокращения кривой - большому кругу .

Формула монотонности Хьюскена

Согласно формуле монотонности Хьюскена , свертка развивающейся кривой с обращенным во времени тепловым ядром не возрастает. Этот результат можно использовать для анализа особенностей эволюции.

Конкретные кривые

Кривые с автомодельной эволюцией

Мрачная кривая жатки и ее переведенные копии, созданные потоком, сокращающим кривую.

Поскольку любая другая простая замкнутая кривая сходится к окружности, окружность - единственная простая замкнутая кривая, которая сохраняет свою форму в потоке сокращения кривой. Однако есть много других примеров кривых, которые либо непростые (они включают самопересечения), либо незамкнутые (они простираются до бесконечности) и сохраняют свою форму. Особенно,

который имеет положительные периодические решения (соответствующие кривым с автомодельной эволюцией) для любого периода, лежащего строго между π и .
  • Другие кривые, включая некоторые бесконечные спирали , остаются самоподобными с более сложными движениями, включая вращение или комбинации вращения, сжатия или расширения и перемещения.
  • Для сетей гладких кривых, встречающихся тройками на стыках с углами 2 π / 3, самоподобные сжимающиеся решения включают двойной пузырь, окружающий две равные области, форму линзы ( vesica piscis ), ограниченную двумя конгруэнтными дугами окружностей вместе с два коллинеарных луча, вершины которых находятся в углах линзы, и сеть «в форме рыбы», ограниченная отрезком прямой, двумя лучами и выпуклой кривой. Любые другие самоподобные сжимающиеся сети включают большее количество кривых. Другое семейство сетей растет гомотетически и остается самоподобным; это древовидные сети кривых, пересекающихся под углами 2 π / 3 в тройных стыках, асимптотические вееру из двух или более лучей, которые встречаются в общей конечной точке. Двухлучевой случай этих форм представляет собой неограниченную гладкую кривую; для трех или более лучей эволюция этих форм может быть определена с использованием обобщенных вариантов потока укорачивания кривой, таких как вариант для варифолдов. Данный веер из четырех или более лучей может быть асимптотическим по отношению к нескольким различным решениям этого типа, поэтому эти решения не обеспечивают однозначного определения потока, укорачивающего кривую, начинающегося с веера лучей.

Древние решения

Древнее решение к задаче потока является кривым, эволюция которой может быть экстраполирована назад на все время, без особенностей. Все самоподобные решения, которые сжимаются или остаются того же размера, а не растут, в этом смысле являются древними решениями; они могут быть экстраполированы назад, обращая преобразование самоподобия, которое они претерпели бы при прямом потоке, сокращающем кривую. Таким образом, например, круг, мрачный жнец и кривые Абреша – Лангера - все это древние решения.

Есть также примеры, которые не являются самоподобными. Ярким примером является решение овала Angenent по работе Angenent (1992) . Это семейство кривых можно параметризовать, задав кривизну как функцию касательного угла с помощью формулы

и имеют в качестве своей ограничивающей формы при обратной эволюции пару мрачных изгибов жнеца, приближающихся друг к другу с противоположных направлений. В декартовой системе координат они могут быть заданы неявным уравнением кривой

В физической литературе такие же формы известны как модель скрепки .

Решения Angenent в форме овала и стягивающейся окружности - единственные древние решения, временные интервалы которых ограничивают ограниченные выпуклые множества. Grim Reaper, стационарное полупространство и стационарное ленточное решение - единственные примеры, временные интервалы которых ограничивают неограниченные выпуклые множества. Существует много дополнительных (невыпуклых) локально выпуклых примеров, а также множество дополнительных (невыпуклых) вложенных примеров.

Численные приближения

Чтобы эффективно вычислить поток сокращения кривой, как непрерывную кривую, так и непрерывную эволюцию кривой необходимо заменить дискретным приближением.

Переднее отслеживание

Методы отслеживания фронта давно используются в гидродинамике для моделирования и отслеживания движения границ между различными материалами, крутых градиентов свойств материалов, таких как погодные фронты , или ударных волн внутри одного материала. Эти методы включают вывод уравнений движения границы и их использование для прямого моделирования движения границы, а не моделирование подстилающей жидкости и рассмотрение границы как возникающего свойства жидкости. Те же методы можно использовать для моделирования потока, сокращающего кривую, даже если кривая, на которой проходит поток, не является границей или скачком уплотнения.

В методах слежения вперед для сокращения кривой кривая, претерпевающая эволюцию, дискретизируется в виде многоугольника. Метод конечных разностей используется для вывода формул для приблизительного вектора нормали и кривизны в каждой вершине многоугольника, и эти значения используются для определения того, как перемещать каждую вершину на каждом временном шаге. Хотя поток сокращения кривой определяется движением кривой перпендикулярно самой себе, некоторые параметризации потока сокращения кривой могут позволить вершинам, которые аппроксимируют кривую, двигаться неперпендикулярно. Фактически, это позволяет вершинам перемещаться по кривой по мере ее развития. Выбор тщательной повторной параметризации может помочь более равномерно перераспределить вершины вдоль кривой в ситуациях, когда перпендикулярное движение может привести к их объединению. Мерриман, Бенс и Ошер (1992) пишут, что эти методы являются быстрыми и точными, но гораздо сложнее распространить их на версии потока сокращения кривой, которые применяются к более сложным входным данным, чем простые замкнутые кривые, где необходимо заниматься особенностями и изменениями топологии.

Для большинства таких методов Цао (2003) предупреждает, что «условия стабильности не могут быть легко определены, и временной шаг должен выбираться специально». Другой метод конечных разностей, разработанный Crandall & Lions (1996), изменяет формулу кривизны в каждой вершине, добавляя к ней небольшой член, основанный на операторе Лапласа . Эта модификация называется эллиптической регуляризацией , и ее можно использовать для доказательства существования обобщенных потоков, а также при их численном моделировании. Используя его, можно доказать, что метод Крэндалла и Лайонса сходится, и это единственный численный метод, указанный Цао, который снабжен границами скорости его сходимости. Для эмпирического сравнения прямого метода Эйлера , обратного Эйлера и более точных методов конечных разностей Кранка – Николсона см. Balažovjech & Mikula (2009) .

Свертка с повторной выборкой

Mokhtarian & Mackworth (1992) предлагают численный метод вычисления приближения к потоку сокращения кривой, который поддерживает дискретное приближение к кривой и чередует два шага:

  • Выполните повторную выборку текущей кривой, разместив новые точки выборки с одинаковым интервалом, измеренным по нормализованной длине дуги.
  • Сверните положения точек с помощью функции Гаусса с небольшим стандартным отклонением, фактически заменяя местоположение каждой точки средневзвешенным значением местоположения ближайших точек вдоль кривой с гауссовыми весами. Стандартное отклонение Гаусса должно быть выбрано достаточно малым, чтобы после этого шага точки выборки все еще имели почти равномерный интервал.

Как они показывают, этот метод сходится к распределению сокращения кривой в пределе, когда количество точек выборки растет, а нормализованная длина дуги радиуса свертки сокращается.

Медианная фильтрация

Мерриман, Бенс и Ошер (1992) описывают схему, работающую на двумерной квадратной сетке - фактически массив пикселей . Кривая, которая должна развиваться, представлена ​​путем присвоения значения 0 (черный) пикселям, находящимся вне кривой, и 1 (белый) пикселям, находящимся внутри кривой, что дает индикаторную функцию для внутренней части кривой. Это представление обновляется путем чередования двух шагов:

Чтобы эта схема была точной, временной шаг должен быть достаточно большим, чтобы кривая сместилась хотя бы на один пиксель даже в точках с низкой кривизной, но достаточно малым, чтобы радиус размытия был меньше минимального радиуса. кривизны. Следовательно, размер пикселя должен быть O (min κ / max κ 2 ) , достаточно малым, чтобы можно было выбрать подходящий промежуточный временной шаг.

Метод может быть обобщен на эволюцию сетей кривых, пересекающихся на стыках и разделяющих плоскость более чем на три области, путем одновременного применения одного и того же метода к каждой области. Вместо размытия и определения порога этот метод можно альтернативно описать как применение медианного фильтра с гауссовыми весами к каждому пикселю. Можно использовать ядра, отличные от теплового ядра, или адаптивно уточнять сетку, чтобы она имела высокое разрешение около кривой, но не тратила время и память на пиксели, далекие от кривой, которые не влияют на результат. Вместо использования только двух значений в пиксельном изображении версия этого метода, использующая изображение, значения пикселей которого представляют собой знаковое расстояние до кривой, может достичь субпиксельной точности и потребовать более низкого разрешения.

Приложения

Отжиг металлических листов

Раннее упоминание Уильяма В. Маллинза  ( 1956 ) об укорачивании кривой потока мотивирует его как модель физического процесса отжига , в котором термообработка вызывает смещение границ между зернами кристаллизованного металла. В отличие от мыльных пленок , которые из-за разницы в давлении воздуха становятся поверхностями постоянной средней кривизны , границы зерен при отжиге подвержены только локальным эффектам, которые заставляют их перемещаться согласно потоку средней кривизны. Одномерный случай этого потока, поток с сокращением кривой, соответствует листам отжига металла, которые достаточно тонкие, чтобы зерна стали фактически двумерными, а их границы - одномерными.

Анализ формы

В обработке изображений и компьютерное зрение , Mokhtarian & Макворт (1992) предлагает применять кривую укорачиванию поток к контуру формы , полученной из цифрового изображения, для того , чтобы удалить шум от формы и обеспечить масштабное пространство , что обеспечивает упрощенное описание формы на разных уровнях разрешения. Метод Мохтариана и Макворта включает в себя вычисление потока, сокращающего кривую, отслеживание точек перегиба кривой по мере их продвижения по потоку и построение графика, который отображает положения точек перегиба вокруг кривой в зависимости от параметра времени. Точки перегиба обычно удаляются с кривой попарно, когда кривая становится выпуклой (в соответствии с теоремой Гейджа – Гамильтона – Грейсона), а время жизни пары точек соответствует заметности особенности формы. Из-за метода свертки с повторной выборкой, который они описывают для вычисления численного приближения потока сокращения кривой, они называют свой метод масштабным пространством повторной выборки кривизны . Они отмечают, что это масштабное пространство инвариантно относительно евклидовых преобразований данной формы, и утверждают, что оно однозначно определяет форму и устойчиво к небольшим изменениям формы. Они сравнивают его экспериментально с несколькими связанными альтернативными определениями масштабного пространства для форм и обнаруживают, что масштабное пространство кривизны после повторной выборки менее требовательно к вычислениям, более устойчиво к неравномерному шуму и менее сильно зависит от мелкомасштабных различий форм.

Реакция – диффузия

В реакционно-диффузионных системах, моделируемых уравнением Аллена-Кана , предельное поведение для быстрой реакции, медленной диффузии и двух или более локальных минимумов энергии с одинаковым энергетическим уровнем друг у друга приводит к тому, что система оседает в областях с разными локальными значениями. минимумы, при этом фронты, ограничивающие границы между этими областями, эволюционируют согласно потоку, сокращающему кривую.

Клеточные автоматы

Клеточный автомат Anneal, 1600 шагов после случайного запуска

В клеточном автомате каждая ячейка в бесконечной сетке ячеек может иметь одно из конечного набора состояний, и все ячейки обновляют свои состояния одновременно, основываясь только на конфигурации небольшого набора соседних ячеек. Жизнь-как клеточный автомат правило является один , в котором сетка является бесконечной квадратной решеткой, существует ровно два состояния ячейки, множество соседей каждой ячейки восемь соседей по соседству Мура , и правило обновления зависит только от количество соседей с каждым из двух состояний, а не с какой-либо более сложной функцией этих состояний. В одном конкретном жизненном правиле, введенном Джерардом Вичниаком и называемом правилом искаженного большинства или правилом отжига, правило обновления устанавливает новое значение для каждой ячейки как большинство среди девяти ячеек, заданных ею и ее восемью соседями, за исключением случаев, когда эти ячейки делятся на четыре с одним состоянием и пять с другим состоянием, и в этом случае новое значение ячейки является меньшинством, а не большинством. Детальная динамика этого правила сложна, включая существование небольших стабильных структур. Однако в совокупности (при запуске со всеми ячейками в случайных состояниях) он имеет тенденцию формировать большие области ячеек, которые все находятся в том же состоянии, что и друг друга, причем границы между этими областями развиваются в соответствии с потоком сокращения кривой.

Строительство замкнутых геодезических

Поток сокращения кривой можно использовать для доказательства изопериметрического неравенства для поверхностей, гауссова кривизна которых является невозрастающей функцией расстояния от начала координат , например параболоида . На такой поверхности гладкий компакт, имеющий любую заданную площадь и минимальный периметр для этой области, обязательно представляет собой круг с центром в начале координат. Доказательство применяет поток сокращения кривой к двум кривым, метрической окружности и границе любого другого компакта и сравнивает изменение периметра двух кривых, поскольку они обе сводятся потоком к точке. Поток сокращения кривой также можно использовать для доказательства теоремы о трех геодезических , согласно которой каждое гладкое риманово многообразие, топологически эквивалентное сфере, имеет три геодезические, образующие простые замкнутые кривые .

Связанные потоки

К другим геометрическим потокам, связанным с потоком сокращения кривой, относятся следующие.

  • Для моделирования поведения кристаллов или других анизотропных материалов важно иметь варианты течения с сокращением кривой, для которых скорость потока зависит от ориентации кривой, а также от ее кривизны. Один из способов сделать это - определить энергию кривой как интеграл от гладкой функции γ ее нормальных векторов и сформировать градиентный поток этой энергии, согласно которому нормальная скорость, с которой течет кривая, пропорциональна анизотропный аналог кривизны. Этот поток можно смоделировать, дискретизируя кривую в виде многоугольника. В численных экспериментах кажется, что исходные кривые сходятся к форме Вульфа для γ, прежде чем сузиться до точки. В качестве альтернативы можно позволить кривой течь со скоростью a ( θ ) κ + b ( θ ), где κ - (обычная) кривизна, а a и b - гладкие функции ориентации θ . Когда a ( θ + π ) = a ( θ ) и b ( θ + π ) = - b ( θ ) (так что поток инвариантен относительно точечного отражения ), можно показать, что результирующий поток подчиняется принципу избегания и аналог теоремы Гейджа – Гамильтона – Грейсона.
  • Поток аффинного сокращения кривой был впервые исследован Alvarez et al. (1993) и Сапиро и Танненбаум (1993) . В этом потоке нормальная скорость кривой пропорциональна кубическому корню из кривизны. Результирующий поток инвариантен (с соответствующим масштабированием по времени) относительно аффинных преобразований евклидовой плоскости, большей группы симметрии, чем преобразования подобия, при которых инвариантен поток, сокращающий кривую. В рамках этого потока применяется аналог теоремы Гейджа – Гамильтона – Грейсона, согласно которому любая простая замкнутая кривая в конечном итоге становится выпуклой, а затем сходится к эллипсу, когда она схлопывается в точку.
  • Преобразование кривой с равными нормальными скоростями во всех точках было названо преобразованием травяного пожара . Кривые, образованные таким образом, обычно образуют острые углы, след которых образует среднюю ось кривой. Тесно связанная эволюция кривой, которая перемещает прямые сегменты многоугольной кривой с равной скоростью, но позволяет вогнутым углам двигаться быстрее, чем единичная скорость, вместо этого формирует другой тип топологического каркаса данной кривой, ее прямой каркас .
  • Для поверхностей в более высоких измерениях существует более одного определения кривизны, включая внешние (зависящие от внедрения) меры, такие как средняя кривизна, и внутренние меры, такие как скалярная кривизна и кривизна Риччи . Соответственно, существует несколько способов определения геометрических потоков на основе кривизны, включая поток средней кривизны (в котором нормальная скорость вложенной поверхности является ее средней кривизной), поток Риччи (собственный поток на метрике пространства, основанный на его кривизна Риччи), поток кривизны Гаусса и поток Уиллмора (градиентный поток для функционала энергии, сочетающего среднюю кривизну и гауссову кривизну). Поток сокращения кривой является частным случаем потока средней кривизны и потока кривизны Гаусса для одномерных кривых.
  • Вдохновленные потоком сокращения кривой на гладких кривых, исследователи изучили методы плавных многоугольников, чтобы они оставались многоугольными, с приложениями, включая формирование узоров и синхронизацию в распределенных системах роботов. Сохраняющие длину многоугольные потоки можно использовать для решения проблемы правила плотника .
  • В компьютерном зрении модель активного контура для обнаружения краев и сегментации изображения основана на сокращении кривой и развивает кривые на основе комбинации их кривизны и особенностей изображения.

Примечания

использованная литература