Boa noite pessoal,

Tenho um novo problema, tenho um raster (dem.crop.final.tif) que representa a elevação do terreno (parte de um objeto Topodata do INPE) e preciso georreferenciar o raster, sendo que tenho as coordenadas atuais (x.orig e y.orig) que estão erradas e preciso corrigir com as coordenadas certas (x.new e y.new) para todos os pixels do raster baseado nestes 6 pontos, alguém sabe se é possível realizar essa transformação sem utilizar o GRASS ou QGIS, segue CRM dos dados:

### <code r>

lapply(c('rgdal', 'maptools', 'shapefiles'), require, character.only=TRUE)

#Download do controno da área em shapefile

links <- c(

"https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif")

if (sum(sapply(basename(links), file.exists))<1) {

lapply(links, function(a) tryCatch(download.file(a, dest=basename(a), mode='wb'),

error=function(...) message("ERRO!")))} else message('OK!')

dem <- readGDAL('dem.crop.final.tif')

x.orig<-c(-42.33960, -42.33937, 42.34546, -42.34738, -42.21652, -42.40492)

y.orig<-c(-19.23349, -19.23514,-19.23217, -19.22818,-19.05137,-19.50393)

x.new<-c(-42.33929, -42.33921, -42.34515, -42.34714, -42.21620, -42.40471)

y.new<-c(-19.24486, -19.24655,-19.24350,-19.23955,-19.06276,-19.51528)

### </code>

Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal 
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial 
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr@yahoo.com.br 
        alexandre.santos@cas.ifmt.edu.br 
Lattes: http://lattes.cnpq.br/1360403201088680                   
======================================================================