Оценка объема земляных выемок в QGIS

[slon]

Добрый день!
По работе столкнулся с такой задачей:
Есть карьер, необходимо его засыпать, для этого необходимо оценить объем завозимого грунта??
основой исходных данных послужил топоплан местности 1:500, по отметкам существующих и проектных высот были построены гриды двумя методами Kriging и Radial Basis function. Для построения поверхностей использовался Surfer и в нем же был рассчитан объем.
Вопрос заключается в следующем - можно ли средствами Qgis произвести такой же расчет??

[ericsson]

Средствами SAGA из QGIS.
http://www.saga-gis.org/saga_tool_doc/2 ... lus_2.html
https://docs.qgis.org/2.6/en/docs/user_ ... olume.html

[trir]

В Civil 3D очень удобно это делать. А главное там можно в 3D рассмотреть TIN и проверить правильность её построения. А ещё там есть структурные линии и т. п.
И в таких масштабах надо использовать TIN, а не GRID.

[slon]

Спасибо, будем пробывать

[slon]

Попробовал средства SAGA, не совсем то что нужно, верхний слой выбирается значение, а у меня две поверхности. Попробовал еще средства Grass - r.volume , объем существенно отличается от рассчитанного в Surfer, и при "загрузке результирующего слоя" в логе ошибка:

The following layers were not correctly generated.
Centroids
You can check the log messages to find more information about the execution of the algorithm

Пробовал как TIN так и GRID.

[slon]

ребят помогите, разъясните порядок расчета объема в Qgis-grass.
в архиве исходники - точечные шэйпы высот нижнего и верхнего слоя, и маска слоя в котором оценивается объем
в surfer получены вот- такие вот цифры:

metod V,м3
Radial bas.func 176114
Kriging 173994
TIN 176762

[Игорь Белов]

Давайте начнём с Surfer'а. На результат влияет геометрия грида. Поскольку точки на нижней поверхности расположены очень плотно, нужно строить грид с шагом не более одного метра, чтобы можно было говорить о сходимости. Построим гриды в рамках X = [498050..498650], Y = [5102800..5103400] с шагом метр на метр методом TIN, отбланкуем маской, вычислим объёмы. Полученные насыпи/выемки:
Positive Volume [Cut]: 174758.89023836
Negative Volume [Fill]: 294.70535759341
Net Volume [Cut-Fill]: 174464.18488077

В других программах, использующих гриды, геометрия должна быть идентична. Это значит, в QGIS следует строить гриды в рамках X = [498049.5..498650.5], Y = [5102799.5..5103400.5] с шагом метр на метр. Используем тот же метод TIN. Один из возможных путей:
Интерполяция:
[Растр] → [Интерполяция] → [Интерполяция]
Обрезать растр по маске:
[Растр] → [Извлечение] → [Обрезка]
Получить разности объёмов по ячейкам:
[Растр] → [Калькулятор растров]
Вычислить объём:
В инструментах анализа [Геоалгоритмы QGIS] → [Инструменты обработки растра] → [Статистика растрового слоя]
Получаем объём:
Sum: 174472.996208

С учётом того, что разные программы по-разному выполняют триангуляцию Делоне, можно говорить о совпадении результатов.

[slon]

Игорь большое спасибо за развернутый ответ, все понятно.
Есть вопрос в продолжение темы о шаге грида, не знаю создать отдельную тему или здесь продолжить на имеющемся примере?!
Не силен в аналитической геометрии, Вы написали о шаге не более метра при такой плотности точек. Каким способом можно определить наилучший (оптимальный) шаг грида, чтобы обеспечить наилучший(достоверный) результат??
Попробовал данный пример с шагом 0,25 м, разница получается всего в 4 м3)) конечно для данных объемов в 170 тыс. это смешная погрешность)

[Игорь Белов]

Минимальный шаг грида должен быть порядка расстояния между близкими точками.

Добавим к таблице данных колонку с порядковым номером id. Вычислим расстояния между ближайшими точками:
[Вектор] → [Анализ] → [Матрица расстояний]

dist_martix.png (13.66 КБ) 6834 просмотра

Откроем полученную электронную таблицу CSV, отсортируем колонку MEAN и полюбуемся на наименьшие значения, отличные от нуля:

stats.png (13.88 КБ) 6834 просмотра

Похоже, шаг грида для сходимости результатов должен быть не более 0,5 метра.

[slon]

Спасибо большое))!

[trir]

Civil 3D
Объем выемки (отрегулировано) 306.71 Куб. м
Объем насыпи (отрегулировано) 174767.32 Куб. м
Чистый объем (отрегулировано) 174460.61 Куб. м<Насыпь>