[R-br] Mover coordenadas de um raster (Georreferenciamento)

ASANTOS alexandresantosbr em yahoo.com.br
Domingo Maio 4 23:51:18 BRT 2014


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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140504/fa4bfc8b/attachment.html>


Mais detalhes sobre a lista de discussão R-br