<div dir="ltr"><div><div><div><div><div><div>Alexandre,<br><br></div>Fiz isso para imagens de satélite<br><br></div>O procedimento que usei para criar o laço foi com um array de 3 dimensões, onde as duas primeiras eram o x e o y e coloquei um matriz "a traz" da outra como se fosse construindo um cubo então tinha um objeto com x,y e t . O registro do nome de cada figura, eu armazenei em um vetor<br>
<br></div>Pra ser sincero, não entendi muito bem sua rotina pois não domino algumas das funções que vc usou mesmo assim mando um exemplo do que fiz, que espero que te ajude.<br><br></div><div># carregar os arquivos de interesse no diretorio<br>
</div>my.dir = dir("/home/tito/tcc/tau",pattern=".nc",full.names=T)<br></div><div><br></div><div># criar a variável nula do array para não dar problema<br></div><div>TAU =NULL<br></div><div><br># loop para ler cada um dos arquivos<br>
</div>for(arquivo in my.dir){<br></div><br># carregar o arquivo pelo nome<br><div>nc = open.ncdf(arquivo)<br></div><div><br># extrair variavel de tau em x<br></div><div>taux = get.var.ncdf(arquivo,"x")<br><br></div>
<div># armazenar a variavel em um array segundo as dimenções e o numero da figura e de acordo com o nome do arquivo<br></div><div>TAU = array(c(TAU,taux),dim=c(nrow(taux),ncol(tauy),which(my.dir==arquivo))))<br></div><div>
<br><div><div>}<br><div><div><div><br></div><div>Espero ter ajudado<br><br></div><div><br><br><br><br><br><br><div><div><div><div><br><br><br></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_extra">
<br clear="all"><div>Tito Conte<br><br></div>
<br><br><div class="gmail_quote">2013/12/4 ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Caros Membros,<br>
<br>
Estou precisando de ajuda para criar um laço para extrair as coordenadas e o tempo (no nome do arquivo) em múltiplas figuras. No meu exemplo abaixo e usando apenas a figura 1, tenho uma regiões de pixels de interesse em branco e o nome da figura Fig1-0_00_00.png, preciso das coordenadas x e y das regiões em branco (Valor digital = 0), representado pelo objeto fig.data.pix e juntei ao tempo que é 0:00:00, pois o nome da Figura é 0_00_00. Conforme CRM abaixo, fazendo figura por figura esta OK, mas gostaria de fazer um laço para extrair essas informações das três figuras abaixo e armazená-las em um único objeto, alguém teria alguma sugestão?<br>
<br>
<br>
<br>
#-----------------------------<u></u>------------------------------<u></u>--------------------<br>
# Rotina para extração de coordenadas e tempo em imagens<br>
#-----------------------------<u></u>------------------------------<u></u>--------------------<br>
<br>
# Abertura das imagens<br>
<br>
require(RCurl)<br>
<br>
#Figura 1 ------------------------------<u></u>------------------------------<u></u>----------<br>
get.dropbox <- function(file, token) {<br>
require(RCurl, quietly=TRUE)<br>
bin <- getBinaryURL(paste0("<a href="https://www.dropbox.com/s/8b99gn9ve1yo7is/Fig1-0_00_00.png" target="_blank">https://<u></u>www.dropbox.com/s/<u></u>8b99gn9ve1yo7is/Fig1-0_00_00.<u></u>png</a>", token, "/", file),<br>
ssl.verifypeer = FALSE)<br>
con <- file(file, open = "wb")<br>
writeBin(bin, con)<br>
close(con)<br>
message(noquote(paste(file, "read into", getwd())))<br>
}<br>
#<br>
get.dropbox("Fig1-0_00_00.png"<u></u>, "Fig1-0_00_00.png")<br>
<br>
#Figura 2 ------------------------------<u></u>------------------------------<u></u>----------<br>
get.dropbox <- function(file, token) {<br>
require(RCurl, quietly=TRUE)<br>
bin <- getBinaryURL(paste0("<a href="https://www.dropbox.com/s/8b99gn9ve1yo7is/Fig2-0_20_00.png" target="_blank">https://<u></u>www.dropbox.com/s/<u></u>8b99gn9ve1yo7is/Fig2-0_20_00.<u></u>png</a>", token, "/", file),<br>
ssl.verifypeer = FALSE)<br>
con <- file(file, open = "wb")<br>
writeBin(bin, con)<br>
close(con)<br>
message(noquote(paste(file, "read into", getwd())))<br>
}<br>
#<br>
get.dropbox("Fig2-0_20_00.png"<u></u>, "Fig2-0_20_00.png")<br>
<br>
#Figura 3 ------------------------------<u></u>------------------------------<u></u>---------<br>
get.dropbox <- function(file, token) {<br>
require(RCurl, quietly=TRUE)<br>
bin <- getBinaryURL(paste0("<a href="https://www.dropbox.com/s/bo6gbt3ca6xwxm0/Fig3-0_40_00.png" target="_blank">https://<u></u>www.dropbox.com/s/<u></u>bo6gbt3ca6xwxm0/Fig3-0_40_00.<u></u>png</a>", token, "/", file),<br>
ssl.verifypeer = FALSE)<br>
con <- file(file, open = "wb")<br>
writeBin(bin, con)<br>
close(con)<br>
message(noquote(paste(file, "read into", getwd())))<br>
}<br>
#<br>
get.dropbox("Fig3-0_40_00.png"<u></u>, "Fig3-0_40_00.png")<br>
<br>
# Extraido as coordenadas dos pixels para figura 1 -----------------------------<br>
<br>
require(raster)<br>
r <- raster("Fig1-0_00_00.png")<br>
image(r)<br>
r2<-as.data.frame(r)<br>
<br>
# Cria as coordenadas ------------------------------<u></u>----------------------------<br>
grid0 <- expand.grid(x=seq(0,length.<u></u>out=r@nrows),<br>
y=seq(0,length.out=r@ncols))<br>
<br>
# Coordenadas + valor dos pixels<br>
<br>
fig.data<-cbind(grid0,r2)<br>
<br>
# Região de interesse - Número digital = 0<br>
<br>
fig.data.pix<-fig.data[fig.<u></u>data[,3]==0,]<br>
<br>
# Extrai o tempo do nome do arquivo<br>
tempo2<-r@data@names<br>
tempo1<-gsub("Fig1.", "", tempo2)<br>
tempo<-gsub("_", ":", tempo1)<br>
<br>
# Cria o arquivo final<br>
<br>
dados<-cbind(fig.data.pix[,1:<u></u>2],rep(tempo))<br>
<br>
head(dados)<br>
<br>
#-----------------------------<u></u>------------------------------<u></u>--------------------<br>
<br>
Obrigado,<br>
<br>
-- <br>
==============================<u></u>==============================<u></u>==========<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@<u></u>yahoo.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.<u></u>br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/<u></u>1360403201088680</a><br>
==============================<u></u>==============================<u></u>==========<br>
<br>
______________________________<u></u>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/<u></u>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br></div>