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 Ambos lados da revisão seguinte | ||
|
software:art:curso:chuvapr [2007/12/12 00:43] paulojus |
software:art:curso:chuvapr [2007/12/12 00:44] paulojus |
||
|---|---|---|---|
| Linha 124: | Linha 124: | ||
| 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 134: | ||
| 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 142: | ||
| 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 198: | Linha 198: | ||
| - | **// 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 | ||
