On Qua, 2010-06-23 at 12:34 -0700, Pedro Venâncio wrote: > > > O v.surf.rst permite criar validação cruzada, com o parâmetro cvdev > (ainda não testei), este script permite determinar iterativamente os > valores óptimos de tensão e smooth? Que vantagens tem em relação à > validação cruzada do próprio v.surf.rst? > > Em relação aos MDT hidrológicamente consistentes, qual é a metodologia > que usa? A rede hidrográfica tem informação altimétrica associada? > Algumas aproximações que já vi quando não se dispõe de rede > hidrográfica 3D, é criar o MDT, extrair a informação altimétrica do > MDT para a rede hidrográfica (eventualmente através do v.drape), > subtrair um valor constante a essa altimetria (por exemplo 0,2 - 0,3m) > e depois refazer o MDT, incluindo a rede hidrográfica como se de > curvas de nível se tratassem. > > Como funciona o v.breach? E, já agora, como o coloco a funcionar? > > Hoje fiz um teste para criar um MDT com as minhas curvas de nível, com > o v.surf.rst. Demorou cerca de 5h, mas o resultado, numa análise > visual, e com recurso a alguns perfis topográficos a atravessar alguns > vales e cumeadas, parece bastante bom. O tempo é que é demasiado... > Vou tentar generalizar a informação (v.generalize) e refazer o MDT > pelo v.surf.rst, e comparar o tempo e os resultados. > Estou neste momento também a fazer uma tentativa com o > r.surf.contour.. Está a demorar mais ainda que o v.surf.rst, e fiquei > com uma dúvida, converti as curvas de nível para raster com um pixel > de 2m, para ficarem semelhantes às originais em vector, e estou a usar > o r.surf.contour com uma resolução de 25m. Este procedimento estará > correcto? Deveria usar a mesma resolução na rasterização das curvas de > nível e na interpolação por r.surf.contour? >
O script serve para correres o v.surf.rst em modo de validação cruzada, variando um determinado intervalo de valores nos parâmetros e escrevendo as estatísticas dos resíduos num ficheiro de texto. Executa-se assim uma análise de sensibilidade à variação dos parâmetros, para um determinado conjunto de dados. Em relação ao v.breach, é usado para retirar os artefactos introduzidos em linhas de água com vales apertados. Tem o inconveniente de necessitar de polilinhas que representem as linhas de água. O produto é um conjunto de pontos, com cotas, que se adicionam aos dados vectoriais iniciais. Em relação ao tempo necessário para o v.surf.rst, experimenta criar vários Mapsets, copiar os dados vectoriais, e executar vários v.surf.rst (vários threads) com regiões de análise (g.region) diferentes. No final juntas os vários resultados com uma combinação de r.mapcalc e r.patch. Se os teus dados têm uma distribuição homogénea podes sempre baixar os valores de segmax (nº de pontos num determinado segmento) e npmin ( nº de pontos usados para interpolação). A computação será mais rápida. Mas se os segmentos forem visíveis terás de aumentar o npmin ou até usar o r.random.cells para densificar com dependência espacial as amostras nas áreas em défice a partir de um mdt provisório. Em relação à resolução, salvo raras excepções, nas computações será sempre usado a resolução definida para a região de análise, independentemente da resolução dos dados usados. Em caso de dados com resolução diferente da região será usado um resampling "on the fly". Para nosso controlo sobre os processos recomendo preparar os dados para a resolução da região de análise. Luís _______________________________________________ Portugal mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/portugal
