R-скрипт для нахождения оптимального разрешения растра

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

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

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

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

С readShapePoly вроде бы cправился:
вместо

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

plots <- readShapePoly("plots.shp", proj4string=CRS(as.character(NA)))

сделал

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

plots <- readShapePoly(list.files(pattern="plots"), force_ring=TRUE, proj4string=CRS(as.character(NA)))

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

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

plots$SCI <- plots$PERIMETER / (2*sqrt( sapply(plots@polygons, getPolygonAreaSlot)/pi)*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.

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

pixel_gis-lab_version.R.7z

(6.34 КБ) 443 скачивания

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

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

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

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