## 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@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@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 ###############################################################