[andrik75]

Доброго времени суток всем!
Есть у меня задача над решением которой я бьюсь уже год с лишним и все никак.
Выходные данные: растровая карта месности или же выкопировка с горизонталями.
Задача: По этой карте построить цифровую модель рельефа. "Оцифровать" можно и вручную, тоесть ввести в программу нужные горизонтали и отметки высот. В результате программа на области карты должна построить некую гладкую функцию вида z=f(x,y), имея которую я бы мог уже анализировать рельеф местности. Визуализация не нужна. Нужна идея и алгоритм самой реализации на любом языке программирования, так как я сам програмер и разобраться в этом смогу.
Интересует: алгоритм триангуляции, применимый для геодезии, или модернизированная триангуляция Делонэ, а также алгоритм быстрого поиска попадания точки P(x,y) в ккакой то из треугольников триангуляции.
А вопщем интересует все, что с этим связано!

[andrik75]

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

...вот вот и я о том же!

Вообще хранить в дереве точки неправильно ИМХО, нужно хотябы треугольники.

))) и я как не странно пока заснул до этого тоже додумался

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

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

[Владимир.]

Сообщение от Supermax Триангуляция говорите? - А зачем?

Зачем вдаваться в детали: что, почем и почему...
Это уже давно придумали и я не знаю другого, чтобы использовали в современных программах типа Геоникса, кроме триангуляции Делонэ

Сообщение от andrik75 Триангуляция Делонэ, где-то я читал, в случае с горизонталями дает не очень удовлетворительные результаты.

В смысле? Построенная по горизонталям? Или наоборот? Горизонтали, построенные по труиангуляции?
Для расчетов отметок/объемов горизонтали не используют. Они только для визуальной "проверки" рельефа. Считают по поверхности, т.е. по тем самым треугольникам

[andrik75]

в смысле построенная по горизонталям!
ну в принципе, то что нужно я уже нашел, спасибо Александру achawk.narod.ru, буду разбираться, так что тему можно считать решенной

[zamtmn]

этот топик был рекламой?

[andrik75]

Сообщение от zamtmn этот топик был рекламой?

а Вы видите в этом смысл? я нет!
это уже офтоп, но все же, это мои поиски решения, которые на данном этапе отобразились именно в вышеизложенное. и не более того!

[slogos]

Если тема еще актуальна, возможно смогу оказать конкретную помощь. У меня есть собственные наработки по этой теме. По набору пикетов (точек с xyz) cтроится триангуляция Делоне. Есть функция укладки проектной линии на поверхность (с математической точки зрения - проектирование плоской полилинии на поверхность, представленную треугольниками, а точнее 3DFace). а в ней, конечно z = f(x,y), только определение в какой треугольник попадает точка решено просто перебором.Приложения составлены на VLISPe. Поэкспериментировал со сглаживанием поверхности - пришел к выводу, что это не приводит к ощутимому улучшению результата. Проверял на стереомоделях в Photomode. Точнее в каких то случаях это происходит, но нет способа заранее сказать когда это случится. А вообще на эту тему много публикаций, в одной диссертации, например, доказывается, что точнее всего сглаживать поверхность рельефа методом крикинга. Замена сети треугольников регулярной сеткой целесообразна только для облегчения дальнейших вычислений, но это искажает результат. Если работаете по тополанам, где горизонтали строились на основе пикетов, то их не нужно использовать - горизонтали это производная от пикетов это раз, а второе - Вы ну уверены что они построены верно.

[ShaggyDoc]

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

Да, идеально было бы использовать "пикеты", т.е. точки, в которых непосредственно бралась отметка. К сожалению, отметок всех точек, которые брались в поле, на топопланах нет. Это уж внутренние хитрости топографов.

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

[Сергей Дубина]

А как быть с негладкими функциями: дороги с бордюрами, обрывы, откосы, подпорные стены и т.д.?

[slogos]

Это т.н. структурные линии рельефа. Есть два варианта решения: 1)Строить триангуляцию с учетом этого, а это значит в алгоритме д.б. предусмотрено, что треугольники не должны пересекать эти динии 2)строить триангуляцию "так как есть " а потом удалаять ребра пересекающие линии т.е. преобразовать триангуляцию. Вообще говоря, через перечисленные объекты не должны "проходить" горизонтали - только и всего, а ЦМР может быть построена правильно и без учета их, но в этом случае пикеты должны выбираться "грамотно".

[slogos]

Какова конечная цель? Для построения цифровой модели рельефа и дальнейшего ее использования я использую Lisp-программу. Поскольку это делаю в Acade, то не создаю никакого дерева и не храню информацию о соответствии точек и треугольников. Что называется построил 3D-Face и забыл. ЦМР в виде треугольников (3Dfase) наиболее точна, т.к. она строится по пикетам, выбираемым специальным образом на месиности. Всякая сглаживание не улучшает модель. А если все-таки это необходимо делать, то можно создать регулярную сеть, а можно сгладить и сплайнами. Все это легко делается на Lispe/

[slogos]

"каждая точка знает в каких треугольниках учавствует", "бинарное дерево", "матрица инцидентности" и пр - это хорошо для классических языков, для Vlispa (да и для самой задачи) все должно быть проще. И есть - например у Елпанова. По массиву точек строится поверхность в виде 3Dface. Всякое сглаживание поверхности вносит погрешность. А для определения z (т.е. z = f(x,y)) нужной точки достаточно прогу которая находит точку пересечени отрезка с плоскостью 3Dграни. Если время поиска попадания точки на ту или иную грань критично, можно придумать всякие хитрости. Для сведения - при количестве граней 500-700 (точек 400) ждать не долго.

[lema]

Здравствуйте все. Сразу скажу - я не работаю в автокаде. Меня интересует, есть ли возможность преобразовать топографическую карту автокада, в растровую карту высот?
Типа вот в такую:
<a href='http://hostingkartinok.com/show-image.php?id=07148561a947f824d251e4c428457e4d' title='imagehost'><img src='http://s7.hostingkartinok.com/uploads/thumbs/2015/09/07148561a947f824d251e4c428457e4d.png' /></a>
Если такое можно сделать, то на сколько это сложно и кто может такое осуществить? Могут ли геодезисты такую услугу оказать?

[Boxa]

Хм, сначала нужно посмотреть что Вы понимаете под "топографическую карту автокада", потом будет видно можно ли преобразовать.

ЗЫ.
ИМХО, если Вы не собираетесь писать код, то ошиблись разделом.

[lema]

Сообщение от Boxa Хм, сначала нужно посмотреть что Вы понимаете под "топографическую карту автокада", потом будет видно можно ли преобразовать.

Во вложении.

ИМХО, если Вы не собираетесь писать код, то ошиблись разделом.

Что за код? Я не знаю о чём речь. Увидел обсуждение о преобразовании из одного в другое, на сколько я понял из растрового в векторное, вот и задал вопрос здесь, ибо в этой теме скорее всего должны находиться люди которые в "теме"(простите за каламбур) Мне же вроде тоже самое но только наоборот нужо
Я же не знаю какова данная реализация, может там (грубо говоря) пару кнопок нажать для ковертирования, просто я об этом не знаю. Для начала хотел узнать существует ли такая возможность, если да, то на сколько сложно это реализовать.

[Boxa]

Сообщение от lema Во вложении.

Прикольно, а что должно на выходе получится, применительно к этому исходнику?

Сообщение от lema Что за код? Я не знаю о чём речь. Увидел обсуждение о преобразовании из одного в другое, на сколько я понял из растрового в векторное, вот и задал вопрос здесь, ибо в этой теме скорее всего должны находиться люди которые в "теме"(простите за каламбур) Мне же вроде тоже самое но только наоборот нужо Я же не знаю какова данная реализация, может там (грубо говоря) пару кнопок нажать для ковертирования, просто я об этом не знаю. Для начала хотел узнать существует ли такая возможность, если да, то на сколько сложно это реализовать.

Код- программный код на любом языке программирования.
Штатно, в автокаде, такого конвертера нет, на счет наличия такого в Civil 3D не знаю.
Так что или смотреть, что там есть в Civil 3D или писать код.

тоже самое но только наоборот -

[lema]

Сообщение от Boxa Прикольно, а что должно на выходе получится, применительно к этому исходнику?

Примерно такое:

Штатно, в автокаде, такого конвертера нет, на счет наличия такого в Civil 3D не знаю.

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

Так что или смотреть, что там есть в Civil 3D или писать код.

Civil 3D - я вообще не знаю, что это. Я работаю в скетчапе и в специальной программе по ландшафту, которая понимает только растровую карту высот. Вот такая боль.
тоже самое но только наоборот - [/quote]
А что не так?

Код- программный код на любом языке программирования.

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

[lema]

Сообщение от Boxa Прикольно, а что должно на выходе получится, применительно к этому исходнику?

То что выше я по ссылке указал
hostingkartinok.com/show-image.php?id=07148561a947f824d251e4c428457e4d
Мне не важно как это делается, меня интересует возможно ли это сделать и чего за это хотят?

[trir]

есть в Civil 3D - "Экспорт поверхности в DEM"

[Кулик Алексей aka kpblc]

lema, картинку не видно.

[trir]

в специальной программе по ландшафту

?

Откуда, кстати, такой план? там точки - анонимные блоки...