Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
| Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior | ||
|
software:art:curso:chuvapr [2007/12/12 00:43] paulojus |
software: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> |
