Prezado Elias, Fornecendo dados junto com seu código, permitindo torná-lo reproduzível (CMR), você obtém ajuda mais facilmente.
Suponho que o 17 seja em função da variável que você quer representar no grid que está exportando (predição). Estou supondo que 17 seja o índice da última coluna adicionada (cbind). Um objeto SpatialPixelsDataFrame permite multivariáveis, mas pode ser que o formato escolhido exija a seleção de uma única variável pra compor um grid. ================================================ Éder Comunello Em 27 de setembro de 2013 07:43, ELIAS MENDES COSTA MENDES escreveu: > > Obrigado Éder, era isso mesmo, agora que rodei um comando para tirar todos > os NAs deu certo. > Já aproveitando o email talvez pudesse me ajudar mais uma vez. No script > que estou usando o comando que uso para salvar o mapa predito em um formato > que > abra no SagaGis é: writeGDAL(grid_V[,17], "V_rna37.sdat", > drivername="SAGA", mvFlag=-99999) Você sabe explicar o porque do 17? > > Os comandos anteriores a esse foram. > > # fazer o mapa para V%, ou aplicar o modelo para a area de estudo > grid_V <- cbind(grid1, predito) # usar "predito" ou "test1" > names(grid_V) > gridded(grid_V) <- ~X+Y # transformar "grid_V" de DataFrame para > SpatialPixelsDataFrame > proj4string(grid_V) <- CRS(utm) # definir projecao > names(grid_V) > spplot(grid_V, "predito") # plotar no R apenas a coluna "predito" > >
_______________________________________________ R-br mailing list [email protected] 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.
