Линии экстремумов

[nirhoh]

Подскажите пожалуйста. В какой-нибудь ГИС реализована функция поиска условного и безусловного минимума-максимума (построение линий экстремумов по полю; переход от площадных характеристик к линейным)? Пример на рис. 1

рис.1.jpg (90.05 КБ) 14589 просмотров

[Александр Мурый]

Не скажу насчёт "линий экстремумов по полю", но на картинке лично я вижу горизонтали и профиль, проведённый по водоразделу. Если нужно выделить линии водоразделов, то этого можно добиться несколькими похожими способами.
1) В растровом калькуляторе в любой ГИС умножаем ЦМР на (-1), получаем "обращённый" рельеф. Далее с помощью соответствующих инструментов строим речную сеть — это будет сеть водоразделов.
2) Строим бассейны рек по исходной ЦМР. С помощью растровой алгебры выделяем границы бассейнов, векторизуем их —> получаем водоразделы.

[nirhoh]

Спасибо. Но первый способ однозначно не подходит. На водоразделах бывают седловины (разнонаправленное изменение гипсометрического уровня), а водоток - односторонне направленное изменение гипсометрического уровня.

[nirhoh]

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

[gamm]

пунктирная линия - линия экстремума

это ни разу не экстремум (либо у вас вся линия имеет один уровень). Если уровень действительно один, то линия трассируется по равенству этому значению. Если нет, то у вас либо тальвег, либо водораздел (говоря геоморфологически). Если теперь весь "рельеф" сильно наклонить (чтобы гарантировать "сток" в одну сторону), то на вогнутых участках вы получите водотоки, на выпуклых - водоразделы. Что и строится в какой-нибудь SAGA

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

[nirhoh]

Попробую сформулировать задачу. Работаю в ArcGis 9.3 c ЦМР. Выделение «водотоков» и определение их порядка инструментами Гидрология (Hydrology) мне не подходит (сейчас не буду объяснять почему). Наиболее приемлемо для меня выделение «водотоков» через горизонтальную кривизну. Области аномальной горизонтальной кривизны (рис.2) нужно представить в виде линий.

[gamm]

Работаю в ArcGis 9.3

причем здесь это чудо? если вам нужна инструкция, какие кнопки нажимать, это не ко мне. Задача не сформулирована - горизонтальную кривизну посчитать проблем нет, что мешает? и потом строить "водоток" на этой кривизне

[nirhoh]

Инструкция, какие кнопки нажимать в ArcGis мне не нужна. Я сформулировал задачу - подскажите какими средствами ее можно решить.

[nirhoh]

строить "водоток" на этой кривизне... Подскажите каким образом... В том то и дело, что "аномальная" кривизна носит "фрагментарный" характер...Примерно как на рис.3

[gamm]

Я сформулировал задачу

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

Успехов в достижении цели!

[nirhoh]

Спасибо. Видимо у нас разные понятия экстремума. Посмотрите еще раз мое сообщение Сегодня, 13:34 и Ваш ответ на него.

[gamm]

Спасибо. Видимо у нас разные понятия экстремума. Посмотрите еще раз мое сообщение Сегодня, 13:34 и Ваш ответ на него.

экстремум - это точка. У нее кривизны нет.

судя по всему, вам нужно что-то подобное этому, там и вторая статья есть. И с терминами нужно аккуратнее

[nirhoh]

Вначале стандартная процедура расчета горизонтальной кривизны по ЦМР. Горизонтальная кривизна имеет числовую характеристику. Далее стандартная процедура классической статистики: среднее и стандартное отклонение. Аномальные значения - среднее +_ стандартное отклонение. По числовым показателям аномальной кривизны рисуются изолинии (рис.2) Вопрос как в "автомате" нарисовать линию экстремума?

[gamm]

Вначале стандартная процедура расчета горизонтальной кривизны по ЦМР. Горизонтальная кривизна имеет числовую характеристику. Далее стандартная процедура классической статистики: среднее и стандартное отклонение. Аномальные значения - среднее +_ стандартное отклонение. По числовым показателям аномальной кривизны рисуются изолинии (рис.2) Вопрос как в "автомате" нарисовать линию экстремума?

еще раз - это не экстремум (по крайней мере с точки зрения математики). Ссылку я дал, там есть код для Матлаба, потом нужно трассировать, например в ET. В Арке вы вряд ли что-то найдете, тем более готовую кнопку.

с точки зрения классификации форм это валы (если они идут вверх - параболоиды при локальной квадратичной аппроксимации), можно попробовать что-то типа GIS ECO им Шарого, или самому прописать по статье Вуда.

[nirhoh]

экстремум - это точка. У нее кривизны нет.

Видимо мы точно путаемся в понятиях.