Gerar arquivo de saída netCDF no R
Olá, eu preciso gerar um arquivo de saída .nc (netcdf) dos resultados da média da coluna de ozônio que carregue as informações de lat, lon e valor em cada ponto de grade da região analisada. Já consegui gerar um arquivo .txt, porém o programa que utilizo não lê arquivos modo texto. abaixo o arquivo de entrada também .nc Tentei utilizar a função create da biblioteca net https://www.dropbox.com/s/ntvmtx7wonwo839/teste.nc?dl=0 library(maps) library(ncdf4) dados <- nc_open('teste.nc') # lendo coordenadas espaço-temporal lat <- ncvar_get( dados, 'lat' ) lon <- ncvar_get( dados, 'lon' ) time <- ncvar_get( dados, 'time' ) # lendo dados coluna total de Ozônio ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' ) # dimensoes da variavel ColumnAmountO3 dims_ColumnAmountO3 <- dim(ColumnAmountO3) # tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de grade X tempo dim(ColumnAmountO3) <- c(dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] ) # calculando a média e retornado-a em 2D media_ColumnAmountO3 <- rowMeans( ColumnAmountO3) dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1],dims_ColumnAmountO3[2] ) # longitude varia de 0 a 360, convertendo para -180 a 180, essaconversão é feita para plotagem sobre o mapa for (i in 1:dim(lon)) { if (lon[i] > 180) { lon[i] <- lon[i]-360 } } desde já agradeço ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 81125154 ____________________________________________________________________________
A resposta já foi postada antes - veja a sua própria pergunta intitulada "Criar arquvo de saída a partir de dados NetCDF". Nenhuma das respostas funcionou? Em que ponto do seu código está a função "create" para nós tentarmos te ajudar com o erro? Greetings, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota On Thursday, February 25, 2016 1:00 PM, Mateus Dias Nunes <nunes.mateusdias@gmail.com> wrote: Olá, eu preciso gerar um arquivo de saída .nc (netcdf) dos resultados da média da coluna de ozônio que carregue as informações de lat, lon e valor em cada ponto de grade da região analisada. Já consegui gerar um arquivo .txt, porém o programa que utilizo não lê arquivos modo texto. abaixo o arquivo de entrada também .nc Tentei utilizar a função create da biblioteca net https://www.dropbox.com/s/ntvmtx7wonwo839/teste.nc?dl=0 library(maps) library(ncdf4) dados <- nc_open('teste.nc') # lendo coordenadas espaço-temporal lat <- ncvar_get( dados, 'lat' ) lon <- ncvar_get( dados, 'lon' ) time <- ncvar_get( dados, 'time' ) # lendo dados coluna total de Ozônio ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' ) # dimensoes da variavel ColumnAmountO3 dims_ColumnAmountO3 <- dim(ColumnAmountO3) # tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de grade X tempo dim(ColumnAmountO3) <- c(dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] ) # calculando a média e retornado-a em 2D media_ColumnAmountO3 <- rowMeans( ColumnAmountO3) dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1],dims_ColumnAmountO3[2] ) # longitude varia de 0 a 360, convertendo para -180 a 180, essaconversão é feita para plotagem sobre o mapa for (i in 1:dim(lon)) { if (lon[i] > 180) { lon[i] <- lon[i]-360 } } desde já agradeço ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 81125154 ____________________________________________________________________________ _______________________________________________ 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� c�igo m�imo reproduz�el.
participantes (2)
-
Mateus Dias Nunes -
Thiago V. dos Santos