<div dir="ltr">Bom dia colegas programadores.<div><br></div><div>Estive trabalhando na obtenção de de dados climáticos a partir de um arquivo de extensão .nc, com um gride de 0,25 grau para todo o território nacional.</div><div>Minha necessidade é de extrair essas informações climáticas, dia a dia, em formato numérico para alguns pontos específicos do mapa, que irão coincidir com os pontos das estações meteorológicas de onde coletei os dados observados.</div><div>Após algumas tentativas e com algumas ajudas dos colegas, cheguei no código que segue ao final do texto junto com o link para baixar o arquivo analisado.</div><div>Gostaria de validar esse código, pois não estou familiarizado com esse tipo de arquivo e com o pacote do R que usei para abri-lo. Como não conheço o conteúdo do arquivo, também não sei se houve uma limitação na visualização da saída do comando View(e) pelo R.</div><div>Desde já, agradeço pela ajuda de todos!</div><div><br></div><div>Link para download do arquivo:</div><div><br></div><div><a href="https://drive.google.com/file/d/0B6tnaf2fmoCcOXhnRDIxWnM5UDA/view">https://drive.google.com/file/d/0B6tnaf2fmoCcOXhnRDIxWnM5UDA/view</a><br></div><div><br></div><div>Código do R:</div><div><br></div><div>#seleção de diretorio</div><div>setwd("C:/Users/user/Desktop/Novo_BD")</div><div>#seleção de pacote</div><div>library(raster)</div><div>#seleção do arquivo de interesse</div><div>r=raster("Rs_daily_UT_Brazil_v2_19800101_19891231.nc")</div><div>#seleção da variavel de interesse dentro do arquivo</div><div>s=stack("Rs_daily_UT_Brazil_v2_19800101_19891231.nc", varname="Rs")</div><div>#criar lista de pontos a serem extraidos os valores de Rs</div><div>longitude=c(-47.5,-47.5,-44.5,-40)</div><div>latitude=c(-18.5,-19.5,-20.5,-22)</div><div>xy=cbind(longitude,latitude)</div><div>sp=SpatialPoints(xy)</div><div>#extração dos valores</div><div>e=extract(s, sp, method= "bilinear", df=T) #ou method="simple"</div><div>#View(e)</div><div># Consertar data frame</div><div>df=data.frame(t(e)[,])</div><div>View(df) </div><div><br></div><div><br></div><div><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div>Yury Duarte<br></div>Engenheiro Agrônomo - ESALQ/USP<br></div></div></div>
</div></div>