Точность сгенерированной модели рельефа

[gamm]

Кстати gamm хочется услышать Ваше мнение по поводу выбора размера ячейки грида?!

легко!

Берем масштаб оригинала, потом знаменитые "0.1мм на бумаге", потом умножаем на эмпирический коэффициент (примерно 3-5), и получаем шаг сетки. Для масштаба 1:200 000 получаем 20м*(3..5)=60-100, для масштаба 1:10 000 получаем 3-5м. Это если мы в средней полосе России, и забыли про овраги, обрывы, карстовые воронки, и пр. прелести.

[BES59]

потом умножаем на эмпирический коэффициент

нееее с такими коэф. не интересно работать, как-то не научно... Вообщем судя по всему Все подбирают размер ячейки ЦМР методом тыка. Признаю, зачастую это проще.

[gamm]

потом умножаем на эмпирический коэффициент

нееее с такими коэф. не интересно работать, как-то не научно... Вообщем судя по всему Все подбирают размер ячейки ЦМР методом тыка. Признаю, зачастую это проще.

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

[eucrariano]

Поскольку выводов по предыдущим постам еще нет, задам вопрос по ходу:

... для масштаба 1:10 000 получаем 3-5м. Это если мы в средней полосе России, и забыли про овраги, обрывы, карстовые воронки, и пр. прелести.

А если не забыли, и местность реально насыщена оврагами, обрывами и карстовыми об`ектами, то какой шаг Вы выбираете? А в случае предгорий и гор? При масштабе 1:10000 не приближаемся ли мы в таком случае к "знаменитым "0.1мм на бумаге" (2 м для предгорий и 1 м - горы)? И если местность сильно рассечена в тех же горах, то не стоит ли выбирать шаг сетки меньше 1 м (0.5 м например), для более адекватного построения той же модели и более полного учитывания "структурных линий".

[gamm]

А если не забыли, и местность реально насыщена оврагами, обрывами и карстовыми об`ектами, то какой шаг Вы выбираете?

никакой - ничто уже не поможет, нужно строить векторную модель, например TIN, и в нее встраивать "структурные элементы", что и делается в геологическом моделировании

[Александр Фадеев]

Дамы и господа меня убивают Ваши 0,1 мм на бумаге. Задайте себе вопрос: 1) технически возможно в 1 мм напечатать на бумаге 10 линий, или хотя бы их различить или отличить. 2) и совсем убивает, чем выше разрешение сканирования, тем точность выше. Вопрос к чему? К ошибке заложенной в п. 1, т.е. её описания (прорисовки).

[gamm]

Дамы и господа меня убивают Ваши 0,1 мм на бумаге.

вы просто не в теме, речь не о линейном масштабе

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

И с каким разрешением ни сканируй, детали, отсутствующие в данном масштабе, не появятся.

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

* eucrariano

Указанный выше по теме R-скрипт (для расчёта оптимального разрешения растра) в обновлённом виде выложил здесь.

[eucrariano]

Спасибо amuriu и gamm, вы молодцы. Завтра попробую запустить скрипт, докачав нужные либы к R, поскольку с моим домашним интернетом (declared 115 kbit, а реально меньше) только форум со скрипом и читать

[eucrariano]

Ну вот докачал все что надо.
Имеем:
R version 2.14.1 (2011-12-22)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-pc-linux-gnu (64-bit)
и оболочка - R Commander Version 1.8-1
Скрипт загрузил и дальше ...
Нужен например example 4, для него нужны contours.shp и contours2.shp.
как их подгрузить?
Как запустить сам скрипт?
Короче понятно, что R-statistic в статистике - это аналог GRASS в ГИС, со всеми вытекающими последствиями...
Очень поможет вот это: http://r-analytics.blogspot.com/p/blog-page_09.html - много чего есть по теме.
Не оставляя за бортом постепенное изучение R (а это действительно нужно, необходимо и познавательно), хотелось бы у уважаемых (без кавычек) модификаторов скрипта Hengl’я попросить выложить более детальную инструкцию к использованию их труда. Не все обладают, к сожалению, необходимыми познаниями/временем для быстрого освоения/вникания в суть, и такая "инструкция" очень бы пригодилась многим. Может даже лучше перекинуть ее в "Рецепты"
Спасибо.

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

От лица НЕглавного модификатора скрипта (я просто взялся за это дело, а gamm сделал почти все изменения) скажу, что в R абсолютный новичок, а в математике сильно плаваю, но постараюсь-таки описать работу со скриптом в статье (сначала на вики).
С R Commander дела не имел.
Из командной строки скрипт запускается, например, так:
-- запускается R
-- выполняем source "полный/путь/к/pixel.R" source("полный/путь/к/pixel.R")
-- скрипт должен заработать, по пути данные должны сами скачаться с http://spatial-analyst.net/DATA/ и разархивироваться в текущую директорию
-- после окончания работы скрипта то, что он насчитал, можно посмотреть командой ls(), каждый из выведенных объектов можно посмотреть, введя его имя + Enter, например:

Код: Выделить всё

> plots.pixsize3
       5% 
0.8732125

-- в тексте скрипта есть пояснения, что же это всё значит
-- .....
-- надеюсь, что PROFIT

[eucrariano]

Спасибо за пинок в правильную сторону. Профит получился 50%
Из запущенного R выполняем согласно инструкции:
> source "/home/eucrariano/R/pixel_gis-lab_version.R"
получаем ошибку:
Error: unexpected string constant in "source "/home/eucrariano/R/pixel_gis-lab_version.R""

Читаем документацию и, поставив скобки, получаем:
source("/home/eucrariano/R/pixel_gis-lab_version.R") - скрипт запускается, рисует один график и останавливается со множественными ошибками. Вины gamma и amuriu здесь нет, просто не все пакеты у меня установлены, по ходу их доустановки, ошибок вылазит меньше.
А вот запуск в R-Studio и построчное выполнение скрипта показывает уже другое. Графиков намного больше. Хотя ругань еще есть, несмотря на полную дебиановскую установку. Не хватает rgdal, rsaga и всякого прочего. Вот еще попробую скомпилить R из исходников и отпишусь.

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

получаем ошибку:
Error: unexpected string constant in "source "/home/eucrariano/R/pixel_gis-lab_version.R""

Пардоньте, забыл скобки.

Не хватает rgdal, rsaga и всякого прочего. Вот еще попробую скомпилить R из исходников и отпишусь.

Зачем же ж такие поспешные решения? Может, просто доустановить нужные пакеты из самого R? Про установку пакетов есть прямо чуть ли не на главной странице оф.сайта.

[gamm]

Может, просто доустановить нужные пакеты из самого R?

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

[Максим Дубинин]

Код: Выделить всё

?require
?install.packages