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"