Função com vários if e else if que não funciona
Caros Listeiros, Estou tentando criar uma função abaixo para converter um SpatialGridDataFrame originado de uma imagem de satélite para reflectância, mas a função que eu criei (radconvL()) pelo excesso de else if não esta funcionando. Meu problema é que não dá erro, o que normalmente ajuda a detectar o problema, mas no meu caso a função nem funciona, segue meu CRM abaixo, alguém poderia me ajudar? require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para reflectância radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) Lmin = 0 Qmax = 127 if (is.na(band)) stop("invalid band") if (band == 2) { Lmax = 120.64 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 3) { Lmax = 151.31 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 4) { Lmax = 157.57 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if(band == 5) { Lmax = 69.03 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } } ###### teste<-radconvL(band2, band=2) ## Aqui não funciona e não sei porque teste Muito 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 ======================================================================
Alexandre, A melhor forma de você atacar esse problema é usar as facilidades de depuração do R:
?debug
para ver a documentação do depurador padrão. Um tutorialzinho sobre seu uso pode ser encontrado aqui: http://www.biostat.jhsph.edu/~rpeng/docs/R-debug-tools.pdf HTH -- Cesar Rabak 2015-12-15 18:47 GMT-02:00 ASANTOS <alexandresantosbr@yahoo.com.br>:
Caros Listeiros,
Estou tentando criar uma função abaixo para converter um SpatialGridDataFrame originado de uma imagem de satélite para reflectância, mas a função que eu criei (radconvL()) pelo excesso de else if não esta funcionando. Meu problema é que não dá erro, o que normalmente ajuda a detectar o problema, mas no meu caso a função nem funciona, segue meu CRM abaixo, alguém poderia me ajudar?
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para reflectância radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) Lmin = 0 Qmax = 127
if (is.na(band)) stop("invalid band")
if (band == 2) { Lmax = 120.64 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 3) { Lmax = 151.31 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 4) { Lmax = 157.57 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if(band == 5) { Lmax = 69.03 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } } ###### teste<-radconvL(band2, band=2) ## Aqui não funciona e não sei porque teste
Muito 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 fornea cdigo mnimo reproduzvel.
Uma opção mais eficiente (e provavelmente mais legível) é aproveitar os recursos de vetorização do R. O que eu tenho feito é criar um raster stack com todas as bandas do landsat. Depois, eu crio uma lista com os parametros de calibração de cada banda (que você vai ler nos arquivos de calibração). Funciona melhor se você deixar todos os arquivos - imagens e arquivos de calibração - em uma pasta só. Assumindo que o raster stack tem o mesmo tamanho da lista, você pode aplicar uma multiplicação direta para fazer a calibração. Greetings, -- Thiago V. dos Santos PhD studentLand and Atmospheric ScienceUniversity of Minnesota On Tuesday, December 15, 2015 1:46 PM, Cesar Rabak <cesar.rabak@gmail.com> wrote: Alexandre, A melhor forma de você atacar esse problema é usar as facilidades de depuração do R:
?debug para ver a documentação do depurador padrão. Um tutorialzinho sobre seu uso pode ser encontrado aqui: http://www.biostat.jhsph.edu/~rpeng/docs/R-debug-tools.pdf
HTH--Cesar Rabak 2015-12-15 18:47 GMT-02:00 ASANTOS <alexandresantosbr@yahoo.com.br>: Caros Listeiros, Estou tentando criar uma função abaixo para converter um SpatialGridDataFrame originado de uma imagem de satélite para reflectância, mas a função que eu criei (radconvL()) pelo excesso de else if não esta funcionando. Meu problema é que não dá erro, o que normalmente ajuda a detectar o problema, mas no meu caso a função nem funciona, segue meu CRM abaixo, alguém poderia me ajudar? require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para reflectância radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) Lmin = 0 Qmax = 127 if (is.na(band)) stop("invalid band") if (band == 2) { Lmax = 120.64 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 3) { Lmax = 151.31 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 4) { Lmax = 157.57 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if(band == 5) { Lmax = 69.03 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } } ###### teste<-radconvL(band2, band=2) ## Aqui não funciona e não sei porque teste Muito 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 fornea cdigo mnimo reproduzvel. _______________________________________________ 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.
Muito bem pesado Thiago, vou atentar a isso. E agradeço também ao Cesar, nunca havia usado a ferramenta de depuração e resolveu o meu problema. Segue minha função final, se alguém for trabalhar com imagens do Resourcesat e sensor LISS IIII: ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) if (is.na(band)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if (band == 2) { Lmax = 120.64 } else if (band == 3) { Lmax = 151.31 } else if (band == 4) { Lmax = 157.57 } else if(band == 5) { Lmax = 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######-------------------------------------------------------------------------- -- ====================================================================== 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 ====================================================================== Em 16/12/2015 21:11, Thiago V. dos Santos escreveu:
Uma opção mais eficiente (e provavelmente mais legível) é aproveitar os recursos de vetorização do R.
O que eu tenho feito é criar um raster stack com todas as bandas do landsat. Depois, eu crio uma lista com os parametros de calibração de cada banda (que você vai ler nos arquivos de calibração). Funciona melhor se você deixar todos os arquivos - imagens e arquivos de calibração - em uma pasta só.
Assumindo que o raster stack tem o mesmo tamanho da lista, você pode aplicar uma multiplicação direta para fazer a calibração. Greetings, -- Thiago V. dos Santos
PhD student Land and Atmospheric Science University of Minnesota
On Tuesday, December 15, 2015 1:46 PM, Cesar Rabak <cesar.rabak@gmail.com> wrote:
Alexandre,
A melhor forma de você atacar esse problema é usar as facilidades de depuração do R:
?debug
para ver a documentação do depurador padrão.
Um tutorialzinho sobre seu uso pode ser encontrado aqui: http://www.biostat.jhsph.edu/~rpeng/docs/R-debug-tools.pdf <http://www.biostat.jhsph.edu/%7Erpeng/docs/R-debug-tools.pdf>
HTH -- Cesar Rabak
2015-12-15 18:47 GMT-02:00 ASANTOS <alexandresantosbr@yahoo.com.br <mailto:alexandresantosbr@yahoo.com.br>>:
Caros Listeiros,
Estou tentando criar uma função abaixo para converter um SpatialGridDataFrame originado de uma imagem de satélite para reflectância, mas a função que eu criei (radconvL()) pelo excesso de else if não esta funcionando. Meu problema é que não dá erro, o que normalmente ajuda a detectar o problema, mas no meu caso a função nem funciona, segue meu CRM abaixo, alguém poderia me ajudar?
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para reflectância radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) Lmin = 0 Qmax = 127
if (is.na <http://is.na/>(band)) stop("invalid band")
if (band == 2) { Lmax = 120.64 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 3) { Lmax = 151.31 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 4) { Lmax = 157.57 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if(band == 5) { Lmax = 69.03 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } } ###### teste<-radconvL(band2, band=2) ## Aqui não funciona e não sei porque teste
Muito 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 <mailto:e-mails%3Aalexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.br <mailto:alexandre.santos@cas.ifmt.edu.br> Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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 fornea cdigo mnimo reproduzvel.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
_______________________________________________ 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.
Como é que eu faço para rodar esta função? Muito bem pesado Thiago, vou atentar a isso. E agradeço também ao Cesar, nunca havia usado a ferramenta de depuração e resolveu o meu problema. Segue minha função final, se alguém for trabalhar com imagens do Resourcesat e sensor LISS IIII: ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) if (is.na(band)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if (band == 2) { Lmax = 120.64 } else if (band == 3) { Lmax = 151.31 } else if (band == 4) { Lmax = 157.57 } else if(band == 5) { Lmax = 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######-------------------------------------------------------------------------- -- ====================================================================== 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 ====================================================================== Em 16/12/2015 21:11, Thiago V. dos Santos escreveu: Uma opção mais eficiente (e provavelmente mais legível) é aproveitar os recursos de vetorização do R. O que eu tenho feito é criar um raster stack com todas as bandas do landsat. Depois, eu crio uma lista com os parametros de calibração de cada banda (que você vai ler nos arquivos de calibração). Funciona melhor se você deixar todos os arquivos - imagens e arquivos de calibração - em uma pasta só. Assumindo que o raster stack tem o mesmo tamanho da lista, você pode aplicar uma multiplicação direta para fazer a calibração. Greetings, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota On Tuesday, December 15, 2015 1:46 PM, Cesar Rabak <mailto:cesar.rabak@gmail.com> <cesar.rabak@gmail.com> wrote: Alexandre, A melhor forma de você atacar esse problema é usar as facilidades de depuração do R:
?debug
para ver a documentação do depurador padrão. Um tutorialzinho sobre seu uso pode ser encontrado aqui: http://www.biostat.jhsph.edu/~rpeng/docs/R-debug-tools.pdf <http://www.biostat.jhsph.edu/%7Erpeng/docs/R-debug-tools.pdf> HTH -- Cesar Rabak 2015-12-15 18:47 GMT-02:00 ASANTOS <alexandresantosbr@yahoo.com.br>: Caros Listeiros, Estou tentando criar uma função abaixo para converter um SpatialGridDataFrame originado de uma imagem de satélite para reflectância, mas a função que eu criei (radconvL()) pelo excesso de else if não esta funcionando. Meu problema é que não dá erro, o que normalmente ajuda a detectar o problema, mas no meu caso a função nem funciona, segue meu CRM abaixo, alguém poderia me ajudar? require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para reflectância radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) Lmin = 0 Qmax = 127 if (is.na <http://is.na/> (band)) stop("invalid band") if (band == 2) { Lmax = 120.64 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 3) { Lmax = 151.31 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 4) { Lmax = 157.57 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if(band == 5) { Lmax = 69.03 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } } ###### teste<-radconvL(band2, band=2) ## Aqui não funciona e não sei porque teste Muito 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 <mailto:e-mails%3Aalexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus
Boa tarde Mauro, Segue a função abaixo a uma simulação de um raster, como se fosse uma imagem de satélite e o emprego da função radconvL(): require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) if (is.na(band)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if (band == 2) { Lmax = 120.64 } else if (band == 3) { Lmax = 151.31 } else if (band == 4) { Lmax = 157.57 } else if (band == 5) { Lmax = 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######-- teste<-radconvL(band2, band = 2) teste -- ====================================================================== 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 ====================================================================== Em 20/12/2015 22:04, Mauro Sznelwar escreveu:
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) if (is.na(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if (band == 2) { Lmax = 120.64 } else if (band == 3) { Lmax = 151.31 } else if (band == 4) { Lmax = 157.57 } else if(band == 5) { Lmax = 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######
Ainda não deu certo, veja:
require(raster)
Loading required package: raster
Loading required package: sp
> require(sp)
>
> # RasterLayer inventado
> r <- raster(nrows=10, ncols=10)
> r <- setValues(r, 1:ncell(r))
> plot(r)
>
> band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
>
>
> ##Função para conversão DN para radiância - LISS III
> ---------------------------
+ radconvL<-function(x, band = 2)
+ {
+ BAND <- c("2", "3", "4", "5")
+ band <- pmatch(band, BAND)
+
+ if (is.na(band))
+ stop("invalid band")
+
+ Lmin = 0
+ Qmax = 127
+ x <- as.vector(as.matrix(x))
+ results <- x
+
+ if (band == 2) {
+ Lmax = 120.64
+ }
+ else if (band == 3) {
+ Lmax = 151.31
+ }
+ else if (band == 4) {
+ Lmax = 157.57
+ }
+ else if (band == 5) {
+ Lmax = 69.03
+ }
+
+ x <- Lmin + ((Lmax-Lmin)*x)/Qmax
+ if (class(results) == "SpatialGridDataFrame")
+ results@data[, 1] <- x
+ else if (is.data.frame(x))
+ results <- data.frame(matrix(x, nrow = nrow(results),
+ ncol = ncol(results)))
+ else results <- x
+ results
+ }
Error in ---------------------------radconvL <- function(x, band = 2) { :
object 'radconvL' not found
> ######--
>
> teste<-radconvL(band2, band = 2)
Error: could not find function "radconvL"
> teste
Error: object 'teste' not found
Boa tarde Mauro,
Segue a função abaixo a uma simulação de um raster, como se fosse
uma imagem de satélite e o emprego da função radconvL():
require(raster)
require(sp)
# RasterLayer inventado
r <- raster(nrows=10, ncols=10)
r <- setValues(r, 1:ncell(r))
plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em
SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III
---------------------------
radconvL<-function(x, band = 2)
{
BAND <- c("2", "3", "4", "5")
band <- pmatch(band, BAND)
if (is.na(band))
stop("invalid band")
Lmin = 0
Qmax = 127
x <- as.vector(as.matrix(x))
results <- x
if (band == 2) {
Lmax = 120.64
}
else if (band == 3) {
Lmax = 151.31
}
else if (band == 4) {
Lmax = 157.57
}
else if (band == 5) {
Lmax = 69.03
}
x <- Lmin + ((Lmax-Lmin)*x)/Qmax
if (class(results) == "SpatialGridDataFrame")
results@data[, 1] <- x
else if (is.data.frame(x))
results <- data.frame(matrix(x, nrow = nrow(results),
ncol = ncol(results)))
else results <- x
results
}
######--
teste<-radconvL(band2, band = 2)
teste
--
======================================================================
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
======================================================================
---
Este email foi escaneado pelo Avast antivírus.
https://www.avast.com/antivirus
Mauro, Desculpe a demora, achei o erro, segue a função corrigida: require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) if (is.na(band)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######-- teste<-radconvL(band2, band = 4) teste -- ====================================================================== 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 22/12/2015 06:58, Mauro Sznelwar escreveu:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if (band == 2) { Lmax = 120.64 } else if (band == 3) { Lmax = 151.31 } else if (band == 4) { Lmax = 157.57 } else if (band == 5) { Lmax = 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste<-radconvL(band2, band = 2) teste
Caros Listeiros, A um mês atrás eu consegui criar a função abaixo com a ajuda de vocês, no entanto, descobri que ela funciona quando seleciono band=3, 4 ou 5, mas para band=2 que é a primeira condição da função não funciona (Error in radconvL(band2, band = 2) : object 'Lmax' not found ), alguém poderia me ajudar? Segue CRM: require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) if (is.na(band)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######-- teste2<-radconvL(band2, band = 2) teste2 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate:https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn:https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 28/12/2015 16:22, ASANTOS escreveu:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste<-radconvL(band2, band = 4) teste
Lendo a mensagem de erro, vemos que Lmax não foi definido. Olhando seu programa, Lmax só é definido quando band é 2, ou 3, ou 4 ou 5. Conclui-se que, na sua chamada de função, band não é igual a nenhum desses valores. Quando a linha x <- Lmin + ((Lmax-Lmin)*x)/Qmax é executada, acontece o erro. 2016-02-03 17:44 GMT-02:00 ASANTOS <alexandresantosbr@yahoo.com.br>:
Caros Listeiros,
A um mês atrás eu consegui criar a função abaixo com a ajuda de vocês, no entanto, descobri que ela funciona quando seleciono band=3, 4 ou 5, mas para band=2 que é a primeira condição da função não funciona (Error in radconvL(band2, band = 2) : object 'Lmax' not found ), alguém poderia me ajudar?
Segue CRM:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste2<-radconvL(band2, band = 2) teste2
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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================
Em 28/12/2015 16:22, ASANTOS escreveu:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste<-radconvL(band2, band = 4) teste
_______________________________________________ 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.
Obrigado Thiago, mas isso eu havia recebido, mas não entendi porque ele não lê band=2, mas lê band=3, band=4 e band=5!!!! Caros Listeiros, A um mês atrás eu consegui criar a função abaixo com a ajuda de vocês, no entanto, descobri que ela funciona quando seleciono band=3, 4 ou 5, mas para band=2 que é a primeira condição da função não funciona (Error in radconvL(band2, band = 2) : object 'Lmax' not found ), alguém poderia me ajudar? Segue CRM: require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) if (is.na(band)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######-- teste2<-radconvL(band2, band = 2) teste2 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 03/02/2016 17:10, Thiago Serafim escreveu:
Lendo a mensagem de erro, vemos que Lmax não foi definido. Olhando seu programa, Lmax só é definido quando band é 2, ou 3, ou 4 ou 5. Conclui-se que, na sua chamada de função, band não é igual a nenhum desses valores.
Quando a linha x <- Lmin + ((Lmax-Lmin)*x)/Qmax é executada, acontece o erro.
2016-02-03 17:44 GMT-02:00 ASANTOS <alexandresantosbr@yahoo.com.br <mailto:alexandresantosbr@yahoo.com.br>>:
Caros Listeiros,
A um mês atrás eu consegui criar a função abaixo com a ajuda de vocês, no entanto, descobri que ela funciona quando seleciono band=3, 4 ou 5, mas para band=2 que é a primeira condição da função não funciona (Error in radconvL(band2, band = 2) : object 'Lmax' not found ), alguém poderia me ajudar?
Segue CRM:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na <http://is.na>(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste2<-radconvL(band2, band = 2) teste2
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 <tel:%28%2B55%29%2065%208132-8112> (TIM)(+55) 65 9686-6970 <tel:%28%2B55%29%2065%209686-6970> (VIVO) e-mails:alexandresantosbr@yahoo.com.br <mailto:e-mails:alexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.br <mailto:alexandre.santos@cas.ifmt.edu.br> Lattes:http://lattes.cnpq.br/1360403201088680 OrcID:orcid.org/0000-0001-8232-6722 <http://orcid.org/0000-0001-8232-6722> Researchgate:https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn:https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================
Em 28/12/2015 16:22, ASANTOS escreveu:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na <http://is.na>(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste<-radconvL(band2, band = 4) teste
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
_______________________________________________ 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.
BAND <- c("2", "3", "4", "5"); band = 2; pmatch(band, BAND) ?pmatch Value An integer vector (possibly including NA if nomatch = NA) of the same length asx, *giving the indices of the elements* in table which matched, or nomatch. alternativa: band <- as.integer(BAND[pmatch(band, BAND)]) 2016-02-03 17:44 GMT-02:00 ASANTOS <alexandresantosbr@yahoo.com.br>:
Caros Listeiros,
A um mês atrás eu consegui criar a função abaixo com a ajuda de vocês, no entanto, descobri que ela funciona quando seleciono band=3, 4 ou 5, mas para band=2 que é a primeira condição da função não funciona (Error in radconvL(band2, band = 2) : object 'Lmax' not found ), alguém poderia me ajudar?
Segue CRM:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste2<-radconvL(band2, band = 2) teste2
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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================
Em 28/12/2015 16:22, ASANTOS escreveu:
require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND)
if (is.na(band)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } ######--
teste<-radconvL(band2, band = 4) teste
_______________________________________________ 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.
-- ############################################################### ## Jônatan Dupont Tatsch ## Professor do Departamento de Física ## Coordenador Substituto do Programa de Pós-Graduação em Meteorologia (PPGMET) ## Centro de Ciências Exatas e Naturais (CCNE) ## Universidade Federal de Santa Maria - UFSM ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900 ## Telefone: +55(55)33012083 ## www.ufsm.br/meteorologia ###############################################################
Senhores, boa noite! O problema é que você sobrescreve "band" com o resultado de pmatch (1:4). Pode observar que mesmo rodando com band 3:5 o resultado tá errado. Inclua um print(band) antes do if pra verificar. ================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
Senhores, segue o código rodando. O objeto "band" redundante foi alterado pra "ver". ######################################################################## require(raster); require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") ver <- pmatch(band, BAND) if (is.na(ver)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) return(results) } ######-- teste2<-radconvL(band2, band = 2) teste2<-radconvL(band2, band = 3) ================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] Em 3 de fevereiro de 2016 19:48, Éder Comunello <comunello.eder@gmail.com> escreveu:
Senhores, boa noite!
O problema é que você sobrescreve "band" com o resultado de pmatch (1:4). Pode observar que mesmo rodando com band 3:5 o resultado tá errado.
Inclua um print(band) antes do if pra verificar.
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
Obrigado Jonathan e Éder, Não ia conseguir jamais encontrar o erro ligado pmach() e Éder ficou uma importante dica em sempre fazer um print() para ver os elementos e resultados envolvidos para não incorrer em erros, Redobrados agradecimentos pessoal, Segue o problema resolvido: require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") Band <- pmatch(band, BAND) if (is.na(Band)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 } x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(results) results } ######-- teste2<-radconvL(band2, band = 2) teste3<-radconvL(band2, band = 3) -- ====================================================================== 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 03/02/2016 19:56, Éder Comunello escreveu:
Senhores, segue o código rodando. O objeto "band" redundante foi alterado pra "ver".
########################################################################
require(raster); require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") ver <- pmatch(band, BAND)
if (is.na <http://is.na>(ver)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
if(band == 2) { Lmax <- 120.64 } else if(band == 3) { Lmax <- 151.31 } else if(band == 4) { Lmax <- 157.57 } else if(band == 5) { Lmax <- 69.03 }
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) return(results)
} ######--
teste2<-radconvL(band2, band = 2) teste2<-radconvL(band2, band = 3)
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
Em 3 de fevereiro de 2016 19:48, Éder Comunello <comunello.eder@gmail.com <mailto:comunello.eder@gmail.com>> escreveu:
Senhores, boa noite!
O problema é que você sobrescreve "band" com o resultado de pmatch (1:4). Pode observar que mesmo rodando com band 3:5 o resultado tá errado.
Inclua um print(band) antes do if pra verificar.
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
_______________________________________________ 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.
Alexandre, bom dia! Embora você já tenha resolvido seu problema, lembrei de uma possibilidade de melhoria para o código: o uso de switch() no lugar do bloco if else. ### <code r> require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { Lmax <- switch(as.character(band), "2" = 120.64, "3" = 151.31, "4" = 157.57, "5" = 69.03, NA) if (is.na(Lmax)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) results } ######-- teste2<-radconvL(band2, band = 2) teste3<-radconvL(band2, band = 3) teste1<-radconvL(band2, band = 1) ### </code> ================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
Muito mais "limpo" fica o código final Éder, Muito obrigado, Abraços, -- ====================================================================== 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 05/02/2016 09:58, Éder Comunello escreveu:
Alexandre, bom dia!
Embora você já tenha resolvido seu problema, lembrei de uma possibilidade de melhoria para o código: o uso de switch() no lugar do bloco if else.
### <code r> require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { Lmax <- switch(as.character(band), "2" = 120.64, "3" = 151.31, "4" = 157.57, "5" = 69.03, NA)
if (is.na <http://is.na>(Lmax)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) results } ######-- teste2<-radconvL(band2, band = 2) teste3<-radconvL(band2, band = 3) teste1<-radconvL(band2, band = 1) ### </code>
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
_______________________________________________ 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.
Prezados Listeiros, Mais um problema com a minha função, apesar de eu especificar ao final que class(results) == "SpatialGridDataFrame", o resultado final não é um SpatialGridDataFrame mas um objeto numeric e não sei porque, alguém pode me ajudar? Segue a função nova desenvolvida com a ajuda do Éder: ### <code r> require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { Lmax <- switch(as.character(band), "2" = 120.64, "3" = 151.31, "4" = 157.57, "5" = 69.03, NA) if (is.na <http://is.na>(Lmax)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) results } ######-- teste2<-radconvL(band2, band = 2) teste3<-radconvL(band2, band = 3) teste1<-radconvL(band2, band = 1) ### </code> -- ====================================================================== 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 05/02/2016 09:58, Éder Comunello escreveu:
Alexandre, bom dia!
Embora você já tenha resolvido seu problema, lembrei de uma possibilidade de melhoria para o código: o uso de switch() no lugar do bloco if else.
### <code r> require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { Lmax <- switch(as.character(band), "2" = 120.64, "3" = 151.31, "4" = 157.57, "5" = 69.03, NA)
if (is.na <http://is.na>(Lmax)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- x
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) results } ######-- teste2<-radconvL(band2, band = 2) teste3<-radconvL(band2, band = 3) teste1<-radconvL(band2, band = 1) ### </code>
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
_______________________________________________ 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.
Alexandre, bom dia! Fiz apenas uma pequena alteração indicada no código. ### <code r> require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { Lmax <- switch(as.character(band), "2" = 120.64, "3" = 151.31, "4" = 157.57, "5" = 69.03, NA) if (is.na(Lmax)) stop("invalid band") Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- band2 ### <==== alterado aqui!!! x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) results } ######-- teste2<-radconvL(band2, band = 2) spplot(teste2) teste3<-radconvL(band2, band = 3) teste1<-radconvL(band2, band = 1) ### </code> ================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
Muito obrigado Éder, funcionou perfeitamente!!!!!! Também tem outra alteração que resolve o meu problema e aumenta a velocidade de processamento, uma vez que tenho rasters muito grandes, sendo: radconvL <- function(x, band = 2) { Lmax <- switch(band, "2" = 120.64, "3" = 151.31, "4" = 157.57, "5" = 69.03, NA) if (is.na(Lmax)) stop("invalid band") Lmin = 0 Qmax = 127 setValues(x, Lmin + ((Lmax-Lmin)*values(x))/Qmax) } library(raster) b <- brick(system.file("external/rlogo.grd", package="raster")) test <- radconvL(b[[2]], band = 2) Abraços, -- ====================================================================== 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 08/02/2016 12:14, Éder Comunello escreveu:
Alexandre, bom dia!
Fiz apenas uma pequena alteração indicada no código.
### <code r> require(raster) require(sp)
# RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r)
band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame
##Função para conversão DN para radiância - LISS III --------------------------- radconvL<-function(x, band = 2) { Lmax <- switch(as.character(band), "2" = 120.64, "3" = 151.31, "4" = 157.57, "5" = 69.03, NA)
if (is.na <http://is.na>(Lmax)) stop("invalid band")
Lmin = 0 Qmax = 127 x <- as.vector(as.matrix(x)) results <- band2 ### <==== alterado aqui!!!
x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x print(paste(band, Lmax)) print(results) results } ######--
teste2<-radconvL(band2, band = 2) spplot(teste2)
teste3<-radconvL(band2, band = 3) teste1<-radconvL(band2, band = 1)
### </code>
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
participantes (7)
-
ASANTOS -
Cesar Rabak -
Jônatan -
Mauro Sznelwar -
Thiago Serafim -
Thiago V. dos Santos -
Éder Comunello