2 вопроса про карты и gstat

[Игорь Черниенко]

Рад приветствовать!
Имеются интерполированные в ячейки регулярной сетки (при помощи gstat) значения определенного параметра. Необходимо отобразить их в виде схемы распределения с привязкой к схеме побережья. Функция contour отрисовывать мою сетку категорически не хочет, зато она прекрасно рисуется contourplot и levelplot из lattice.
Карта побережья загружена из шейпа при помощи maptools, может быть представлена в виде полигонов или преобразована в точки - узлы полигонов, которые можно нарисовать при помощи plot(oj,type='l') или lines. Хотелось бы эти карты совместить. При использовании функции lines схема побережья на карту накладывается, но очертания побережья и схема распределения не имеют между собой ничего общего.
Пользуясь случаем, еще вопрос, для тех кто использует gstat.
Имеется код:
создаем объект:
cmm.gs<-gstat(id=names(cmm)[4],formula=cmm[,4]~1,locations=lon+lat,data=cmm)
(cmm - таблица данных со значениями некоторых параметров на станциях)
далее, добавляем данные в объект:
for(i in c(5:12)) {
cmm.gs<-gstat(cmm.gs,id=names(cmm),formula=cmm[,i]~1,locations=lon+lat,data=cmm)
}
имена переменных в объект записываются как положено, а вот данные - только из одного столбца. Может, у меня в цикле что-то не так?
Заранее спасибо.