[R-br] extracao da Longitude correta para dados de SST

[Jônatan]

## usando rotate
require(raster)
b <- brick("~/Downloads/SST_81_14.nc", *stopIfNotEqualSpaced=FALSE*)
b
extent(b)
## fix long
extent(b) <- c(0, 360, -51, 11)
b
plot(b)
## rotate
br <- rotate(b)
extent(br)
br <- setValues(br, getValues(b))
plot(br)

2014-12-05 0:55 GMT-02:00 Humberto <hghazin em hotmail.com>:

> Boa noite pessoal,
>
>
>
> Estou extraindo de um arquivo ncdf os valores de latitude (lat) longitude(
> lon)  e temperatura da superfície da agua do mar( SST), entretanto a
> longitude (lon) ela está na projeção 0 a 360 graus. Minha pergunta é como
> faço para converter ela para -180 a 180 graus? Segue o script que estou
> utilizando juntamente com o link dos dados.
>
> Desde já agradeço a ajuda
>
> Humberto
>
>
>
>           lon
>
>              lat
>
>      data
>
>           sst
>
>                Lon como deveria ser
>
> 299.5
>
> 10.5
>
> 11/30/1981
>
> 27.85
>
>
>
> -60.000
>
> 300.5
>
> 10.5
>
> 11/30/1981
>
> 27.96
>
>
>
> -59.500
>
> 301.5
>
> 10.5
>
> 11/30/1981
>
> 28.02
>
>
>
> -58.500
>
> 302.5
>
> 10.5
>
> 11/30/1981
>
> 28.02
>
>
>
> -57.500
>
> 303.5
>
> 10.5
>
> 11/30/1981
>
> 27.97
>
>
>
> -56.500
>
> 304.5
>
> 10.5
>
> 11/30/1981
>
> 27.84
>
>
>
> -55.500
>
> 305.5
>
> 10.5
>
> 11/30/1981
>
> 27.81
>
>
>
> -54.500
>
> 306.5
>
> 10.5
>
> 11/30/1981
>
> 27.78
>
>
>
> -53.500
>
> 307.5
>
> 10.5
>
> 11/30/1981
>
> 27.75
>
>
>
> -52.500
>
>
>
>
>
> Link para acesso dos dados
> https://www.dropbox.com/s/xnab2g365mduykx/SST_81_14.nc?dl=0
>
>
>
> Script:
>
> require(ncdf) #'ncdf' versao 1.6.6
>
> require(reshape)
>
> require(chron)
>
>
>
> ###sst
>
> nc <- open.ncdf("SST_81_14.nc")
>
> nc
>
>
>
> lat <- get.var.ncdf(nc, "lat")
>
> lon <- get.var.ncdf(nc, "lon")
>
> tunits <- att.get.ncdf(nc,"time","units")
>
>
>
> # print tunits
>
> tunits$value
>
>
>
> # criar tempo "real" Times
>
> time<- get.var.ncdf(nc,"time")
>
> Time<-chron(time,origin=c(month=1, day=1, year=1900))
>
>
>
> ## array da variavel sst
>
> wm <- get.var.ncdf(nc, "sst")
>
> str(wm)
>
>
>
> ## atribuindo dimnames conforme lon, lat e tempo
>
> dimnames(wm)[[1]] <- as.character(lon)
>
> dimnames(wm)[[2]] <- as.character(lat)
>
> dimnames(wm)[[3]] <- as.character(Time)
>
> str(wm)
>
>
>
> ## melt para reestruturação dos dados
>
> mwm <- melt(wm)
>
> head(mwm)
>
> _______________________________________________
> 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.
>



-- 
###############################################################
##  Jônatan Dupont Tatsch
##  Professor do Departamento de Física
##  Centro de Ciências Exatas e Naturais (CCNE)
##  Universidade Federal de Santa Maria
##  Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil -
97105-900
##  Telefone: +55(55)33012083
##  www.ufsm.br/meteorologia
###############################################################
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141205/f65f3e4e/attachment.html>