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 ======================================================================