Alexandre, Fiz isso para imagens de satélite 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 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. # carregar os arquivos de interesse no diretorio my.dir = dir("/home/tito/tcc/tau",pattern=".nc",full.names=T) # criar a variável nula do array para não dar problema TAU =NULL # loop para ler cada um dos arquivos for(arquivo in my.dir){ # carregar o arquivo pelo nome nc = open.ncdf(arquivo) # extrair variavel de tau em x taux = get.var.ncdf(arquivo,"x") # armazenar a variavel em um array segundo as dimenções e o numero da figura e de acordo com o nome do arquivo TAU = array(c(TAU,taux),dim=c(nrow(taux),ncol(tauy),which(my.dir==arquivo)))) } Espero ter ajudado Tito Conte 2013/12/4 ASANTOS <alexandresantosbr@yahoo.com.br>
Caros Membros,
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?
#----------------------------------------------------------- -------------------- # Rotina para extração de coordenadas e tempo em imagens #----------------------------------------------------------- --------------------
# Abertura das imagens
require(RCurl)
#Figura 1 ------------------------------------------------------------ ---------- get.dropbox <- function(file, token) { require(RCurl, quietly=TRUE) bin <- getBinaryURL(paste0("https://www.dropbox.com/s/ 8b99gn9ve1yo7is/Fig1-0_00_00.png", token, "/", file), ssl.verifypeer = FALSE) con <- file(file, open = "wb") writeBin(bin, con) close(con) message(noquote(paste(file, "read into", getwd()))) } # get.dropbox("Fig1-0_00_00.png", "Fig1-0_00_00.png")
#Figura 2 ------------------------------------------------------------ ---------- get.dropbox <- function(file, token) { require(RCurl, quietly=TRUE) bin <- getBinaryURL(paste0("https://www.dropbox.com/s/ 8b99gn9ve1yo7is/Fig2-0_20_00.png", token, "/", file), ssl.verifypeer = FALSE) con <- file(file, open = "wb") writeBin(bin, con) close(con) message(noquote(paste(file, "read into", getwd()))) } # get.dropbox("Fig2-0_20_00.png", "Fig2-0_20_00.png")
#Figura 3 ------------------------------------------------------------ --------- get.dropbox <- function(file, token) { require(RCurl, quietly=TRUE) bin <- getBinaryURL(paste0("https://www.dropbox.com/s/ bo6gbt3ca6xwxm0/Fig3-0_40_00.png", token, "/", file), ssl.verifypeer = FALSE) con <- file(file, open = "wb") writeBin(bin, con) close(con) message(noquote(paste(file, "read into", getwd()))) } # get.dropbox("Fig3-0_40_00.png", "Fig3-0_40_00.png")
# Extraido as coordenadas dos pixels para figura 1 -----------------------------
require(raster) r <- raster("Fig1-0_00_00.png") image(r) r2<-as.data.frame(r)
# Cria as coordenadas ------------------------------ ---------------------------- grid0 <- expand.grid(x=seq(0,length.out=r@nrows), y=seq(0,length.out=r@ncols))
# Coordenadas + valor dos pixels
fig.data<-cbind(grid0,r2)
# Região de interesse - Número digital = 0
fig.data.pix<-fig.data[fig.data[,3]==0,]
# Extrai o tempo do nome do arquivo tempo2<-r@data@names tempo1<-gsub("Fig1.", "", tempo2) tempo<-gsub("_", ":", tempo1)
# Cria o arquivo final
dados<-cbind(fig.data.pix[,1:2],rep(tempo))
head(dados)
#----------------------------------------------------------- --------------------
Obrigado,
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
_______________________________________________ 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.