[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