[R-br] ?Problema para mudar projeção em um raster [Resolvido]
ASANTOS
alexandresantosbr em yahoo.com.br
Quinta Julho 25 08:38:06 BRT 2013
Obrigado Éder,
Problema resolvido, deu certo falatava um espaço mesmo no
comando CMS() e o topodata esta correto agora,
Redobrados agradecimentos,
Alexandre
Em 24/07/2013 22:56, Éder Comunello escreveu:
> Complementando a resposta anterior (abaixo), o arquivo topodata do CMR
> não corresponde à área de interesse definida. A área coincidente é
> dada pelo arquivo 19S435SN.zip.
>
> ###
> require(raster)
> url=("http://www.dsr.inpe.br/topodata/data/geotiff/19S435SN.zip")
> download.file(url, destfile = "19S435SN.zip")
> dem <- raster('19S435SN.tif')
> ###
>
> Observe as saídas transcritas abaixo:
> > extent(SP.geo)
> class : Extent
> xmin : -42.40624
> xmax : -42.37724
> ymin : -19.51573
> ymax : -19.48824
> > dem <- raster('19S435SN.tif')
> > extent(dem)
> class : Extent
> xmin : -43.50024
> xmax : -42
> ymin : -20
> ymax : -19
>
>
> 2013/7/24 Eder Comunello <ecomunel em gmail.com <mailto:ecomunel em gmail.com>>
>
> Prezado Alexandre,
>
> Quanto ao erro de spTransform(), faltou dar espaço antes do
> argumento '+south'.
>
> ### INICIO
>
> ### bnds fornecido no CMS
> bnds <- structure(c(772259.1, 775259.1, 775259.1, 772259.1, 772259.1,
> 7843046, 7843046, 7840046, 7840046, 7843046), .Dim =
> c(5L, 2L
> ), .Dimnames = list(NULL, c("x", "y")))
>
> ### CRS UTM
> SP <- SpatialPolygons(list(Polygons(list(Polygon(bnds)), "1")))
> proj4string(SP) = CRS("+proj=utm +zone=23 +south +datum=WGS84
> +units=m +no_defs") ## Projeção
>
> ### CRS GEO
> SPgeo<- spTransform(SP, CRS("+proj=longlat +datum=WGS84")) ###
> aqui não é necessário +south
> SPgeo em bbox
>
> ### FIM
>
> Com esse código você obterá:
>
> > SPgeo em bbox
> min max
> x -42.40624 -42.37724
> y -19.51573 -19.48824
>
> At.te,
>
>
> --
> Éder Comunello <comunello.eder em gmail.com
> <mailto:comunello.eder em gmail.com>>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
>
>
> --
> Éder Comunello <c
> <mailto:comunello.eder em gmail.com>omunello.eder em gmail.com
> <mailto:omunello.eder em gmail.com>>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
--
======================================================================
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/20130725/be19a4ab/attachment.html>
Mais detalhes sobre a lista de discussão R-br