Страница 2 из 2

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

Добавлено: 06 авг 2015, 11:24
xen87
Максим Дубинин писал(а):вот такое решение есть, более простое, создать 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 - "высота" сохраняемого растра в пикселях

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