Mais

Usando a ferramenta Vector Ruggedness Measure (VRM) da ArcScripts na versão mais recente do ArcGIS for Desktop?


Preciso usar a ferramenta Vector ruggedness measure (VRM) disponível em http://arcscripts.esri.com/details.asp?dbid=15423. No entanto, estou usando ArcGIS 10.2 for Desktop e essas ferramentas são feitas para 9.0, 9.1 e 9.3 apenas. Tentei executar esta ferramenta mesmo assim e obtive o erro:

Traceback (última chamada mais recente): Arquivo "C: Users jc221340 Documents project Processed data My_tools Downloaded ruggedness.py", linha 20, em gp.AddToolbox ("C: / Program Files / ArcGIS / ArcToolbox / Toolboxes / Spatial Analyst Tools.tbx ") RuntimeError: Objeto: Toolbox C: / Arquivos de programas / ArcGIS / ArcToolbox / Toolboxes / Spatial Analyst Tools.tbx não existe

Falha ao executar (Robustez (VRM)).

Como posso usar esta ferramenta?


Como você está usando a v10.2, ele está procurando o caminho da caixa de ferramentas do Spatial Analyst, que mudou ligeiramente os locais do diretório desde a v9:

Para atualizar, abra o robustez.py arquivo usando o Bloco de notas ou algum outro editor de texto e altere o caminho da linha 20 para algo assim:

C:  Arquivos de programas (x86)  ArcGIS  Desktop10.2  ArcToolbox  Toolboxes  Spatial Analyst Tools.tbx

Eu também encontrei problemas ao executar VRM no ArcGIS 10.1. A resposta acima foi útil, mas não o suficiente. A ferramenta VRM gerou arquivos de aspecto, inclinação, x, y, z ... etc, e então mostrou a conclusão do processo. Os métodos a seguir, como geração de camada de robustez e exclusão de arquivos intermediários, foram ignorados.

O problema foi resolvido usando "Benthic Terrain Modeler" em vez de calcular o VRM.


Eu tive o mesmo problema. Mandei um email para o criador da ferramenta, Mark Sappington, ele atualizou esta ferramenta para versões mais recentes do ArcGIS e o novo script Python e caixa de ferramentas ArcGIS podem ser baixados aqui: http://www.arcgis.com/home/item.html?id= 9e4210b3ee7b413bbb1f98fb9c5b22d4

Isso funcionou com a minha versão 10.2.

Boa sorte!


Você pode replicar facilmente essa métrica usando álgebra raster. Este exemplo foi escrito em pseudocódigo de álgebra raster e deve ser fácil de se adaptar a qualquer software GIS que suporte análise raster.

Calcule a inclinação e o aspecto usando uma constante para converter em radianos

inclinação = inclinação (dem, unidades = graus) * 0,0175 aspecto = aspecto (dem) * 0,0175

Aplique as funções cosseno e seno. Para aspecto, normalmente o software GIS atribui -1 às áreas planas. Por causa disso, você pode precisar usar uma instrução con para lidar com a conversão de valores negativos para zero.

sin.slp = cos (inclinação) cos.slp = sin (inclinação) sin.asp = con (aspecto == -1, 0, sin (aspecto) * sin.asp) cos.asp <- con (aspecto == - 1, 0, cos (aspecto) * sin.asp)

Aplique a função de soma focal para a escala desejada (vizinhança focal)

escala = 5 x.sum = focal (sin.asp, janela = escala, função = soma) y.sum = focal (cos.asp, janela = escala, função = soma) z.sum = focal (cos.slp, janela = escala, função = soma)

O Vector Ruggedness Measure (VRM) é 1 - a raiz quadrada das somas quadradas dos rasters focais dividida pela escala quadrada.

vrm = 1 - (sqrt (sqr (x.sum) + sqr (y.sum) + sqr (z.sum)) / sqr (escala))