Форумы GIS-Lab.info

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

Нужна помощь знатоков в "осовременивании" R-скрипта (был написан в 2005 г.) для вычисления оптимального разрешения растра: pixel.R.
Я честно пытался (не зная R), но увяз в ошибках. Например, "Ошибка в Polygon(coords = crds) : ring not closed"
Подозреваю, что это связано с readShapePoly, поиск в сети не особо помог..

С readShapePoly вроде бы cправился:
вместо сделал Далее ------------------------------------------->
Теперь загвоздка со строкой, которая считает т.н. "shape complexity index" (SCI = P/2*r *pi , где r = sqrt(A/pi)) Ошибка:

Ошибка в `[[<-.data.frame`(`*tmp*`, name, value = numeric(0)) :
replacement has 0 rows, data has 147
Вдобавок: Было 50 или более предупреждений (введите warnings() чтобы просмотреть первые 50)

Рад приветствовать
Я свой вопрос присоединю, он в каком-то смысле одного порядка. Пытаюсь прочесть слой постгис
pon.ogr<-readOGR("PG: dbname=crab_sdb","ponet_wgs84")

Получаю ошибку

Ошибка в ogrInfo(dsn = dsn, layer = layer, input_field_name_encoding = input_field_name_encoding) :
Cannot open file

Промежуточные итоги: при поддержке gamm'a (за что ему большое человеческое спасибо и вообще респект) подготовлен вариант скрипта, который должен работать в современных версиях R (>=12.x), мной испытывался только на R 2.14.1.
Для работы скрипта необходимы библиотеки: maptools, spatstat, gstat, RSAGA, rgdal, akima.

Любые отзывы/предложения/идеи приветствуются.

может сразу оформить с описанием что оно делает и как в вики?

Постараюсь в ближайшие дни, но не как обычно, а реально..