[R-br] exportar para SAGAGIS classificação via writeGDAL

[Jônatan]

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 em 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 em 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 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.
>



-- 
###############################################################
##  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
###############################################################
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130716/211bd716/attachment.html>