Waldir, experimente usar o pacote raster: require(raster) ## dados para exemplo com fator r <- raster(ncol=10,nrow=5) r[]<- as.factor(sample(c("a","b","c"),50,replace=T)) plot(r) writeRaster(r,"r.sdat") r1 <- raster("r.sdat") plot(r1) 2013/7/16 Waldir De Carvalho Junior <waldir.carvalho@embrapa.br>
boa tarde a Todos Estou precisando exportar o resultado da classificação via randomforest de um SpatialPixelsDataFrame, onde a predição da classe é a 17ª coluna do objeto. Quero exportar para o SAGAGIS e estou tentando pelo "writeGDAL", mas a resposta é a seguinte:
class(grid_cls)[1] "SpatialPixelsDataFrame" attr(,"package") [1] "sp"> names(grid_cls) [1] "ID" "ascii.mde30" "Flow.Direct" "Watershed.B" "Analytical" "Slope" "Aspect" [8] "Plan.Curvat" "Profile.Cur" "Convergence" "Catchment.A" "Topographic" "LS.Factor" "Channel.Net" [15] "Vertical.Di" "Relative.Sl" "prediction" > grid_cls$prediction[1:15] [1] GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ Levels: CX GM GSJ GX LA LVA PA PVA RL RY> writeGDAL(grid_cls[,17], "Soil_Class.sdat", drivername="SAGA", mvFlag=-99999)Error in create2GDAL(dataset = dataset, drivername = drivername, type = type, : Numeric bands required
O erro indica que o grid (grid_cls[,17]) é uma categoria, mas é isso mesmo, não é numérico, é classe de solo. Qualquer ajuda será bem-vinda. obrigado
__________ Dr. Waldir de Carvalho Junior Pesquisador PhD - Mapeamento Digital Embrapa Solos Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Rio de Janeiro/RJ
Waldir.Carvalho@embrapa.br Telefone: +55 (21) 2179-4614 | Fax: +55 (21) 2274-5291 www.cnps.embrapa.br | twitter.com/embrapa Confira também: www.facebook.com/agrosustentavel
____________________________________________________________________________
Aviso de confidencialidade.
Esta mensagem da Empresa Brasileira de Pesquisa Agropecuária (Embrapa), empresa pùblica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, é enviada exclusivamente a seu destinatário e pode conter informações confidenciais, protegidas por sigilo profissional. Sua utilização desautorizada é ilegal e sujeita o infrator às penas da lei. Se você a recebeu indevidamente, queira, por gentileza, reenviá-la ao emitente, esclarecendo o equívoco.
Confidentiality note
This message from Empresa Brasileira de Pesquisa Agropecuária (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.
_______________________________________________ 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 ## Centro de Ciências Exatas e Naturais (CCNE) ## Universidade Federal de Santa Maria ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900 ## Telefone: +55(55)33012084 ## www.ufsm.br/meteorologia ###############################################################