Работа с растровыми данными в R

Обсуждение материалов сайта: вопросы, замечания, предложения
xen87
Активный участник
Сообщения: 168
Зарегистрирован: 08 апр 2009, 14:44
Проекты: 1
Репутация: 18
Откуда: Кострома

Re: Работа с растровыми данными в R

Сообщение xen87 » 06 авг 2015, 11:24

Максим Дубинин писал(а):вот такое решение есть, более простое, создать world-файл, но разумеется без прописывания всех тэгов геотиф:
после чтения растра

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

test <- "D:\\Programming\\R\\raster-open-save-georef\\s_2000.tif"
pixsize = as.numeric(GDALinfo(test)[6])
originx = as.numeric(GDALinfo(test)[4])
originy = as.numeric(GDALinfo(test)[5])
tfwinfo = rbind(pixsize,0,0,-1*pixsize,originx,originy)
Как я понял, при таком способе в ворд-файле сохраняются координаты левого нижнего угла растра, но зачастую гис-программы (в частности КуГис) требует левый верхний угол, поэтому следует чуть подкорректировать:

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

heigth = as.numeric(GDALinfo(test)[1])
tfwinfo = rbind(pixsize,0,0,-1*pixsize,originx,originy+pixsize*heigth)
где heigth - "высота" сохраняемого растра в пикселях

П.С. большое спасибо за статью!

Ответить

Вернуться в «Материалы сайта»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость