[R-br] ajuda normatização de grid
Eder Comunello
comunello.eder em gmail.com
Sexta Setembro 27 11:44:29 BRT 2013
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"
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130927/285028d4/attachment.html>
Mais detalhes sobre a lista de discussão R-br