Обработка в GRASS лидарных данных

АлексЮстасу

Скачал себе GRASS. Что-то даже удалось загрузить, посмотреть картинку.
Загрузил один фрагмент воздушной лазернойной съемки из 135 000 точек и второй из 200 000 точек. Каждый всего на пару дециметров квадратных плана или до 200х400 м на местности. Т.е. ничтожно маленькие для планов.
Загрузил dxf-файл плана.
Но добиться от модулей для лазерных данных ничего пока не удалось. Надеюсь, разберусь.
Первые впечатления:
- очень медленное отображение. Экран при перемещениях, уменьшениях перерисовывается мучительно - забываешь, что хотел посмотреть.
- трехмерного отображения ни разу не смог получить. Не дождался?
- русские названия файлов, слоев необходимо переделывать латинскими буквами.
- как сделать разрез/посмотреть фрагмент в профиль, покрутить - не нашел.
Первые общие вопросы:
- сможет ли вообще работать при реальных объемах лазерных точек в миллион - несколько миллионов?
- может быть модули для обработки лазерной съемки рассчитаны на наземный метод с соответствующей плотностью точек, а для воздушного не годятся? (Хотя, точек при НЛС бывает и много больше.)

АлексЮстасу писал(а): - очень медленное отображение. Экран при перемещениях, уменьшениях перерисовывается мучительно - забываешь, что хотел посмотреть.

Не работал с лазернойной съемкой, но несколько знаком с GRASS. Первым делом вопрос, касающийся скорости отображения -- какая операционная система? какой интерфейс? Если не пробовали, то попробуйте TclTk -- на моей не очень мощной машине, linux, GRASS 6.4.2 миллион с небольшим точек отображается за 3-4 сек (если вид точки по умолчанию -- basic/circle, если упростить их отображение, то время падает до 1-2 сек.).

АлексЮстасу

GRASS 6.4.3RC2 - скачал просто последний.
Windows 7, x64, i7, 16 Гб, quadro 2000.
Вид точек - какой был по умолчанию. Т.е. крестик Basic - "x". А как еще упростить изображение?

АлексЮстасу писал(а):GRASS 6.4.3RC2 - скачал просто последний.
Windows 7, x64, i7, 16 Гб, quadro 2000.

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

АлексЮстасу писал(а): Вид точек - какой был по умолчанию. Т.е. крестик Basic - "x". А как еще упростить изображение?

Понятно. У меня по умолчанию кружок, который закрашивается да еще и обводка линии присутствует. Как раз при переходе на крестик Basic - "x" получается прирост скорости.

АлексЮстасу писал(а): Загрузил один фрагмент воздушной лазернойной съемки из 135 000 точек и второй из 200 000 точек

Как именно вы загрузили? Через модуль <v.in.ogr>? Приведите команду или скриншот с окном модуля.

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

А чего конкретно надо добиться?

АлексЮстасу писал(а): - очень медленное отображение. Экран при перемещениях, уменьшениях перерисовывается мучительно

Надо посмотреть на свойства импортированных точек (приведите вывод команды "v.info ваш_слой_точек").

Лидарные данные надо импортировать без построения топологии (она там ни к чему).

АлексЮстасу писал(а): - трехмерного отображения ни разу не смог получить. Не дождался?

Используйте NVIZ.

АлексЮстасу писал(а): - русские названия файлов, слоев необходимо переделывать латинскими буквами.

Да, GRASS приучает к правильному обращению с именами векторов — только SQL-совместимые

АлексЮстасу писал(а): - как сделать разрез/посмотреть фрагмент в профиль, покрутить - не нашел.

И не найдёте — спец. модуля для создания разрезов нет, надо программировать или смотреть в чём-то ещё.

АлексЮстасу писал(а): - сможет ли вообще работать при реальных объемах лазерных точек в миллион - несколько миллионов?

При правильном обращении сможет, в сети есть несколько статей по обработке данных лазерн. скан-я в GRASS.

АлексЮстасу писал(а): - может быть модули для обработки лазерной съемки рассчитаны на наземный метод с соответствующей плотностью точек, а для воздушного не годятся? (Хотя, точек при НЛС бывает и много больше.)

Скорее, с точностью до наоборот, спец. модули "v.lidar.*" были созданы для обработки данных воздушного сканирования.

АлексЮстасу

Что такое "интерфейс TclTk"? Как в него заходят?

Добиться хотелось демонстрации возможностей обработки лазерных точек в GRASS. Из-за различий в терминологии я в целом не очень понимаю, что делают все три модуля работы с лазерными точками: edgedetection, growing, correction. И других модулей работы с лазерными точками вообще нет?

Я пока не знаю, как привести вывод команды.

АлексЮстасу

Александр Мурый писал(а):
АлексЮстасу писал(а): Загрузил один фрагмент воздушной лазернойной съемки из 135 000 точек и второй из 200 000 точек
Как именно вы загрузили? Через модуль <v.in.ogr>? Приведите команду или скриншот с окном модуля.

Загрузил через v.in.ascii (через меню). Файл с точками xyz.

Александр Мурый писал(а):
АлексЮстасу писал(а):Но добиться от модулей для лазерных данных ничего пока не удалось. Надеюсь, разберусь.
А чего конкретно надо добиться?

Пока вообще понимаю очень-очень смутно, что эти модули вообще могут делать - терминология непривычная. До сих пор работал только в Terrasolid.

Александр Мурый писал(а):
АлексЮстасу писал(а): - очень медленное отображение. Экран при перемещениях, уменьшениях перерисовывается мучительно
Надо посмотреть на свойства импортированных точек (приведите вывод команды "v.info ваш_слой_точек").

Не знаю пока как сделать это "v.info ваш_слой_точек".

Александр Мурый писал(а):Лидарные данные надо импортировать без построения топологии (она там ни к чему).

Посмотрел загрузку - по умолчанию (а в первые разы грузил по умолчанию) стоит "Не строить топологию в точечном режиме".

Александр Мурый писал(а):
АлексЮстасу писал(а): - трехмерного отображения ни разу не смог получить. Не дождался?
Используйте NVIZ.

Как это NVIZ загружать?

Александр Мурый писал(а):
АлексЮстасу писал(а): - русские названия файлов, слоев необходимо переделывать латинскими буквами.
Да, GRASS приучает к правильному обращению с именами векторов — только SQL-совместимые

Но комментс.

Александр Мурый писал(а):
АлексЮстасу писал(а): - как сделать разрез/посмотреть фрагмент в профиль, покрутить - не нашел.
И не найдёте — спец. модуля для создания разрезов нет, надо программировать или смотреть в чём-то ещё.

А что тогда такое "3D view "?

Александр Мурый писал(а):
АлексЮстасу писал(а): - сможет ли вообще работать при реальных объемах лазерных точек в миллион - несколько миллионов?
При правильном обращении сможет, в сети есть несколько статей по обработке данных лазерн. скан-я в GRASS.

Александр Мурый писал(а):
АлексЮстасу писал(а): - может быть модули для обработки лазерной съемки рассчитаны на наземный метод с соответствующей плотностью точек, а для воздушного не годятся? (Хотя, точек при НЛС бывает и много больше.)
Скорее, с точностью до наоборот, спец. модули "v.lidar.*" были созданы для обработки данных воздушного сканирования.

Ок. Посмотрим.

Вот страница, на которой собрана информация по лидарным данным в GRASS.

АлексЮстасу писал(а): Из-за различий в терминологии я в целом не очень понимаю, что делают все три модуля работы с лазерными точками: edgedetection, growing, correction. И других модулей работы с лазерными точками вообще нет?

Модулей для анализа лидарных данных три, те, что вы перечислили. Они, как я понял, созданы были прежде всего для выделения зданий по облаку точек. Ещё есть модуль v.outlier для фильтрации "отлетающих" точек.
Кроме того, работать с облаком точек с помощью других "векторных" модулей никто не запрещает.

АлексЮстасу писал(а): Я пока не знаю, как привести вывод команды.

В меню "Vector --> Reports and statistics --> Basic vector metadata (v.info)", там выбрать векторный файл с точками.

АлексЮстасу писал(а):Что такое "интерфейс TclTk"? Как в него заходят?

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

В результате в стандартной поставке GRASS несколько таких оболочек, по умолчанию под Windows используется оболочка на python. Интерфейс на TclTk --- еще одна подобная оболочка, вызывать ее, если мне не изменяет память, можно из главного меню программ Windows, там же где вызывается питоновская оболочка, должен быть еще как минимум один пункт меню, с вызовом интерфейса TclTk .

АлексЮстасу

Пытался поставить GRASS (WinGRASS-6.4.3RC2-1-Setup) дома, но при запуске сначала появилась системная ошибка про отсутствие MSVCP100.dll. Переустановил. Потом про MSVCR100.dll. Переустановил.
Теперь появляется заставка запуска GRASS в виде архаичной карты, но GRASS не запускается.
На работе все встало на раз.

АлексЮстасу

Что-то судорожно поделал, и еще раз установил GRASS/
Удалось запустить через "...with Msys"
Попытался загрузить файл на 3 600 000 точек, но получил ответ:

GRASS_INFO_WARNING(2704,1): Файлы coor векторного слоя <prim_160_180@PERMANENT> больше, чем они должны быть (на 134287000 байтов)

с вылетом.
Это точек на 1 кв. км, т.е. обычный лист 1:2000.
Правда, пытался загрузить как 3D.

На попытку запустить через TclTk получил такой ответ:

: Ошибка_TclTk.png (22.06 КБ) 12530 просмотров

АлексЮстасу писал(а): Попытался загрузить файл на 3 600 000 точек, но получил ответ:
GRASS_INFO_WARNING(2704,1): Файлы coor векторного слоя <prim_160_180@PERMANENT> больше, чем они должны быть (на 134287000 байтов)
с вылетом.

Вы можете приложить файл с точками? Просто посмотреть, как он будет себя вести на заведомо рабочей машине?

АлексЮстасу

KolesovDmitry писал(а):Вы можете приложить файл с точками? Просто посмотреть, как он будет себя вести на заведомо рабочей машине?

На работе файл на 3 600 000 точек загрузился, точки отобразились.
Но результатов работы v.lidar.edgedetection добиться не могу.
Пишет:

WARNING: No data within this subregion. Consider changing the spline step.

Это относится к значению шагов интерполяции на восток/север? По умолчанию стоит 4 (чего?). Я ставил и 1, и 8 - результата никакого.

Как понимаю,v.lidar.edgedetection выполняет классификацию точек земли и сразу же возвышающихся над ней объектов в свои классы? И каким-то образом выделяет грани - линии смыкания объектов с землей? В виде чего - векторов или тоже специально классифицированных точек? Или вообще что?

АлексЮстасу писал(а):Пытался поставить GRASS (WinGRASS-6.4.3RC2-1-Setup) дома, но при запуске сначала появилась системная ошибка про отсутствие MSVCP100.dll. Переустановил. Потом про MSVCR100.dll. Переустановил.

Переустановка не решит проблемы отсутствующих в системе программных библиотек. Скачайте и установите пакет библиотек Microsoft Visual C++ 2010 (для 32-битной системы, для 64-битной системы). Правда, не понимаю, зачем они для GRASS...

АлексЮстасу писал(а): На попытку запустить через TclTk получил такой ответ:

Читаем то, что написано: "имя пользователя не должно содержать нелатинские символы". У вас юзер назван, очевидно, по-русски. Можно, как вариант, создать новую учётку с латинским именем специально для GRASS и подобных капризных программ

АлексЮстасу писал(а): Этот файл тоже на 3600000 точек, но на бОльшую площадь.http://yadi.sk/d/VMIYSvHp2s4QU

Можете ещё куда-то залить файл? На работе доступ к http://disk.yandex.ru закрыт.

Форумы GIS-Lab.info

Обработка в GRASS лидарных данных

Обработка в GRASS лидарных данных

Re: есть инструменты для лидарных данных или близких задач?

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Re: Обработка в GRASS лидарных данных

Кто сейчас на конференции