[R-br] Função com vários if e else if que não funciona [RESOLVIDO]
ASANTOS
alexandresantosbr em yahoo.com.br
Quarta Fevereiro 3 22:45:07 BRST 2016
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 em 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 em yahoo.com.br
alexandre.santos em 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 em 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 em gmail.com <mailto:comunello.eder em 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 em 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160203/1ecd3f94/attachment-0001.html>
Mais detalhes sobre a lista de discussão R-br