Re: [R-br] ?Problema para mudar projeção em um raster

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@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@bbox
### FIM
Com esse código você obterá:
SPgeo@bbox min max x -42.40624 -42.37724 y -19.51573 -19.48824
At.te,
-- Éder Comunello <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
-- Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

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@gmail.com <mailto:ecomunel@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@bbox
### FIM
Com esse código você obterá:
> SPgeo@bbox min max x -42.40624 -42.37724 y -19.51573 -19.48824
At.te,
-- Éder Comunello <comunello.eder@gmail.com <mailto:comunello.eder@gmail.com>> Dourados, MS - [22 16.5'S, 54 49'W]
-- Éder Comunello <c <mailto:comunello.eder@gmail.com>omunello.eder@gmail.com <mailto:omunello.eder@gmail.com>> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing list R-br@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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
participantes (2)
-
ASANTOS
-
Éder Comunello