<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Senhores, bom dia!</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Fiz mais uns testes, com 3 alternativas {raster, ncdf4, RNetCDF} para escrever o arquivo NetCDF a partir de uma matriz. Acredito que seja possível melhorar o script do RNetCDF. Se alguém puder colaborar, será bem-vindo.</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style=""><div class="gmail_default"><font face="monospace, monospace">### <code r></font></div><div class="gmail_default"><font face="monospace, monospace">### Dados</font></div><div class="gmail_default"><font face="monospace, monospace">require(raster)</font></div><div class="gmail_default"><font face="monospace, monospace">data(volcano); volcano <- 10*volcano</font></div><div class="gmail_default"><font face="monospace, monospace">image(volcano)</font></div><div class="gmail_default"><font face="monospace, monospace">x <- 100*(1:nrow(volcano))</font></div><div class="gmail_default"><font face="monospace, monospace">y <- 100*(1:ncol(volcano))</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">########## {ncdf4} ###########################################################</font></div><div class="gmail_default"><font face="monospace, monospace"># browseURL("<a href="https://www.getdatajoy.com/learn/Read_and_Write_NetCDF_from_R">https://www.getdatajoy.com/learn/Read_and_Write_NetCDF_from_R</a>")</font></div><div class="gmail_default"><font face="monospace, monospace">library(ncdf4)</font></div><div class="gmail_default"><font face="monospace, monospace">nc1_dim1  <- ncdim_def('EW', 'm', as.double(x))</font></div><div class="gmail_default"><font face="monospace, monospace">nc1_dim2  <- ncdim_def('SN', 'm', as.double(y))</font></div><div class="gmail_default"><font face="monospace, monospace">nc1_varz  <- ncvar_def('Elevation', 'm', list(nc1_dim1, nc1_dim2), -1, longname = 'Data about a Volcano')</font></div><div class="gmail_default"><font face="monospace, monospace">nc1_out   <- nc_create('<a href="http://volcano.nc">volcano.nc</a>', nc1_varz, force_v4 = TRUE)</font></div><div class="gmail_default"><font face="monospace, monospace">ncvar_put(nc1_out, nc1_varz, volcano)</font></div><div class="gmail_default"><font face="monospace, monospace">nc_close(nc1_out)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">print(nc_open("<a href="http://volcano.nc">volcano.nc</a>"))</font></div><div class="gmail_default"><font face="monospace, monospace">plot(raster("<a href="http://volcano.nc">volcano.nc</a>")) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">########## {RNetCDF} #########################################################</font></div><div class="gmail_default"><font face="monospace, monospace"># browseURL("<a href="https://journal.r-project.org/archive/2013-2/michna-woods.pdf">https://journal.r-project.org/archive/2013-2/michna-woods.pdf</a>")</font></div><div class="gmail_default"><font face="monospace, monospace">require(RNetCDF)</font></div><div class="gmail_default"><font face="monospace, monospace">nc2_out <- <a href="http://create.nc">create.nc</a>("<a href="http://volcano2.nc">volcano2.nc</a>")</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://dim.def.nc">dim.def.nc</a>(nc2_out, "EW", length(x))</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://dim.def.nc">dim.def.nc</a>(nc2_out, "SN", length(y))</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.def.nc">var.def.nc</a>(nc2_out, "EW", "NC_DOUBLE", "EW")</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.def.nc">var.def.nc</a>(nc2_out, "SN", "NC_DOUBLE", "SN")</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.def.nc">var.def.nc</a>(nc2_out, "Elevation", "NC_DOUBLE", c("EW", "SN"))</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://att.put.nc">att.put.nc</a>(nc2_out, "Elevation", "_FillValue", "NC_DOUBLE", -999.9)</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.put.nc">var.put.nc</a>(nc2_out, "EW", as.double(x))</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.put.nc">var.put.nc</a>(nc2_out, "SN", as.double(y))</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.put.nc">var.put.nc</a>(nc2_out, "Elevation", volcano)</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://close.nc">close.nc</a>(nc2_out)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">print(nc_open("<a href="http://volcano2.nc">volcano2.nc</a>"))    ### {ncdf4}</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://print.nc">print.nc</a>(<a href="http://open.nc">open.nc</a>("<a href="http://volcano2.nc">volcano2.nc</a>")) ### {RNetCDF}</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.get.nc">var.get.nc</a>(nc2, "EW")</font></div><div class="gmail_default"><font face="monospace, monospace"><a href="http://var.get.nc">var.get.nc</a>(nc2, "Elevation")</font></div><div class="gmail_default"><font face="monospace, monospace">plot(raster("<a href="http://volcano2.nc">volcano2.nc</a>"))</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">########## {raster} ##########################################################</font></div><div class="gmail_default"><font face="monospace, monospace">require(raster)</font></div><div class="gmail_default"><font face="monospace, monospace">tmp <- raster(volcano)</font></div><div class="gmail_default"><font face="monospace, monospace">r <- t(flip(tmp, 1))</font></div><div class="gmail_default"><font face="monospace, monospace">plot(r)</font></div><div class="gmail_default"><font face="monospace, monospace">extent(r) <- c(range(x), range(y))</font></div><div class="gmail_default"><font face="monospace, monospace">plot(r)</font></div><div class="gmail_default"><font face="monospace, monospace">writeFormats()</font></div><div class="gmail_default"><font face="monospace, monospace">writeRaster(r, "<a href="http://volcano3.nc">volcano3.nc</a>", "CDF")</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">print(nc_open("<a href="http://volcano3.nc">volcano3.nc</a>")) ### {ncdf4}</font></div><div class="gmail_default"><font face="monospace, monospace">plot(raster("<a href="http://volcano3.nc">volcano3.nc</a>"))</font></div><div class="gmail_default"><font face="monospace, monospace">### </code></font></div><div style="font-family:verdana,sans-serif"><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="arial, helvetica, sans-serif" style="font-size:small"><div style="font-family:'trebuchet ms',sans-serif;display:inline">​</div></font><font face="arial, helvetica, sans-serif" style="font-size:12.8px"><span style="font-size:small">================================================</span><br>Éder Comunello</font><div style="font-size:12.8px"><span style="font-family:arial,helvetica,sans-serif">Agronomist (UEM), MSc in Environ. Sciences (UEM)</span><font face="arial, helvetica, sans-serif"><br></font></div><div style="font-size:12.8px"><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">DSc in Agricultural Systems Engineering (USP/Esalq)</span><br></div><div><span style="font-size:small">Brazilian Agricultural Research Corporation (</span><font face="arial, helvetica, sans-serif" style="font-size:small">Embrapa)</font></div><div><div><font face="arial, helvetica, sans-serif" style="font-size:small">Dourados, MS, Brazil </font><font face="arial, helvetica, sans-serif" size="2">|<O>|</font></div><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">================================================</span><br></div></div></div><div style="font-size:small"><font face="arial, helvetica, sans-serif">GEO, </font>-22.2752, -54.8182, 408m<span style="font-family:arial,helvetica,sans-serif"><br></span></div><div style="font-size:12.8px"><span style="font-size:small;font-family:arial,helvetica,sans-serif">UTC-04:00 / DST: UTC-03:00</span></div><div><div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div>