Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

Ambos lados da revisão anteriorRevisão anterior
Próxima revisão
Revisão anterior
software:art:curso:chuvapr [2007/12/12 00:43] paulojussoftware:art:curso:chuvapr [2007/12/12 00:56] (atual) paulojus
Linha 5: Linha 5:
 (banco geográfico no formato TerraLib) (banco geográfico no formato TerraLib)
  
 +// 
 +//
 ** 1. Conectando banco, visualizando no R e criando vistasno TV ** ** 1. Conectando banco, visualizando no R e criando vistasno TV **
  
Linha 79: Linha 80:
   - reconectar o banco o TV e ver que "apareceu" a vista agora chamada "parana" !!!   - reconectar o banco o TV e ver que "apareceu" a vista agora chamada "parana" !!!
  
 +//
 +//
 ** 2. Operações da TL via aRT ** ** 2. Operações da TL via aRT **
  
Linha 92: Linha 95:
 </code> </code>
  
 +//
 +//
 ** 3. Trazendo e manipulando dados e fazendo análises no R ** ** 3. Trazendo e manipulando dados e fazendo análises no R **
  
Linha 124: Linha 129:
 Fazendo agora uma interpolação dos dados de chuva via geoestatística/krigagem Fazendo agora uma interpolação dos dados de chuva via geoestatística/krigagem
  
-  Carregando o pacote geoR+  Carregando o pacote geoR
 <code R> <code R>
 require(geoR) require(geoR)
 </code> </code>
-  Convertendo dados os dados para o formato "geodata" da geoR (conveniente) +  Convertendo dados os dados para o formato "geodata" da geoR (conveniente) 
 <code R> <code R>
 geo <- as.geodata(dados, data.col=2) geo <- as.geodata(dados, data.col=2)
Linha 134: Linha 139:
 plot(geo) plot(geo)
 </code> </code>
-  Estimando parâmetros +  Estimando parâmetros 
 <code R> <code R>
 ml <- likfit(geo, trend="1st", ini=c(1000, 100)) um pouco demorado ml <- likfit(geo, trend="1st", ini=c(1000, 100)) um pouco demorado
Linha 142: Linha 147:
 points(loc0, pch=".", col=2) veja gráfico! points(loc0, pch=".", col=2) veja gráfico!
 </code> </code>
-  e fazendo krigagem+  e fazendo krigagem
 <code R> <code R>
 kc <- krige.conv(geo, loc=loc0, krige=krige.control(obj=ml), bor=geo$borders) kc <- krige.conv(geo, loc=loc0, krige=krige.control(obj=ml), bor=geo$borders)
 </code> </code>
-  visualizando no R+  visualizando no R
 <code R> <code R>
 image(kc, col=terrain.colors(15), coords=parana$coords) ver gráfico image(kc, col=terrain.colors(15), coords=parana$coords) ver gráfico
Linha 183: Linha 188:
  
  
-** Adicionando uma nova coluna (com predições) à tabela de dados **\\+// 
 +// 
 +** 4. Adicionando uma nova coluna (com predições) à tabela de dados **\\
 Para exemplificar isto vamos fazer a krigagem nos pontos onde há dados Para exemplificar isto vamos fazer a krigagem nos pontos onde há dados
 <code R> <code R>
Linha 198: Linha 205:
  
  
-**// Fazer no TV://**+**//Fazer no TV://**
   - reconectar o banco e ver que há agora uma nova tabela no layer de dados   - reconectar o banco e ver que há agora uma nova tabela no layer de dados
  
-** 4. Some plots **+// 
 +// 
 +** 5. Some plots **
  
 Plotting directly from the data-base (not using R objects) Plotting directly from the data-base (not using R objects)
Linha 214: Linha 223:
 <code R> <code R>
 deleteLayer(con, "pr.l_pred") deleteLayer(con, "pr.l_pred")
-<code>+</code>
  
  

QR Code
QR Code software:art:curso:chuvapr (generated for current page)