[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