Extrair valores raster

[Lilian]

Oi pessoal,
Queria saber se da para fazer isso no Qgis:
Tenho um conjunto de pontos.
E uma série de camadas raster (.asc)
Queria extrair as informações (uma tablela) de onde seja listado os valores dos rasters para cada ponto.

É possível?

Abraços,
Lilian Costa

[Sidney - Geosaber]

Use o complemento Point Sampling Tool, que coletará as informações de polígonos ou células de um raster e adicionará como atributos ao layer de pontos.

[Lilian]

Que beleza!
Vou tentar ahora!
Obrigada,
Lilian

Em 11 de junho de 2012 18:10, Sidney - Geosaber <sidney...@gmail.com> escreveu:

Use o complemento Point Sampling Tool, que coletará as informações de polígonos ou células de um raster e adicionará como atributos ao layer de pontos.

--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR

[Lilian]

:(

que pena,

não consegui.

Apareceram 2 janelas de erros:

Um erro ocorreu enquanto executava o seguinte código Python:

Traceback (most recent call last):
  File "/home/lili/.qgis//python/plugins/pointsamplingtool/doPointSamplingTool.py", line 98, in updateFieldsList
    self.inData.addItem(str(self.sampItems[i][0].name()) + " : " + str(self.sampItems[i][j][0]) + " (source point)")
UnicodeEncodeError: 'ascii' codec can't encode character u'\xed' in position 5: ordinal not in range(128)

Versão do Python:
2.6.5 (r265:79063, Apr 16 2010, 13:28:26)
[GCC 4.4.3]


Versão do QGIS
1.8.0-Lisboa Lisboa, exported

Caminho para o Python: ['/home/lili/.qgis/python/plugins/elevation', '/usr/share/qgis/python', '/home/lili/.qgis//python', '/home/lili/.qgis//python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.6', '/usr/lib/python2.6/plat-linux2', '/usr/lib/python2.6/lib-tk', '/usr/lib/python2.6/lib-old', '/usr/lib/python2.6/lib-dynload', '/usr/lib/python2.6/dist-packages', '/usr/lib/python2.6/dist-packages/PIL', '/usr/lib/python2.6/dist-packages/gst-0.10', '/usr/lib/pymodules/python2.6', '/usr/lib/python2.6/dist-packages/gtk-2.0', '/usr/lib/pymodules/python2.6/gtk-2.0', '/usr/local/lib/python2.6/dist-packages', '~/.qgis/python', '/home/lili/.qgis/python/plugins/mmqgis/forms', '/usr/share/qgis/python/plugins/fTools/tools']

---

Um erro ocorreu enquanto executava o seguinte código Python:

Traceback (most recent call last):
  File "/home/lili/.qgis//python/plugins/pointsamplingtool/pointSamplingTool.py", line 61, in run
    dialoga = doPointSamplingTool.Dialog(self.iface)
  File "/home/lili/.qgis//python/plugins/pointsamplingtool/doPointSamplingTool.py", line 88, in __init__
    self.updateFieldsList()
  File "/home/lili/.qgis//python/plugins/pointsamplingtool/doPointSamplingTool.py", line 98, in updateFieldsList
    self.inData.addItem(str(self.sampItems[i][0].name()) + " : " + str(self.sampItems[i][j][0]) + " (source point)")
UnicodeEncodeError: 'ascii' codec can't encode character u'\xed' in position 5: ordinal not in range(128)

Versão do Python:
2.6.5 (r265:79063, Apr 16 2010, 13:28:26)
[GCC 4.4.3]


Versão do QGIS
1.8.0-Lisboa Lisboa, exported

Caminho para o Python: ['/home/lili/.qgis/python/plugins/elevation', '/usr/share/qgis/python', '/home/lili/.qgis//python', '/home/lili/.qgis//python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.6', '/usr/lib/python2.6/plat-linux2', '/usr/lib/python2.6/lib-tk', '/usr/lib/python2.6/lib-old', '/usr/lib/python2.6/lib-dynload', '/usr/lib/python2.6/dist-packages', '/usr/lib/python2.6/dist-packages/PIL', '/usr/lib/python2.6/dist-packages/gst-0.10', '/usr/lib/pymodules/python2.6', '/usr/lib/python2.6/dist-packages/gtk-2.0', '/usr/lib/pymodules/python2.6/gtk-2.0', '/usr/local/lib/python2.6/dist-packages', '~/.qgis/python', '/home/lili/.qgis/python/plugins/mmqgis/forms', '/usr/share/qgis/python/plugins/fTools/tools']

[Felipe Dantas]

Sidney, para poder extrair os valores de altitude de uma SRTM, atravez desse plungin que você mencionou (Point Sampling Tool) a raster precisa está em algum formato especifico ou o geotiif basta ?

[Leonardo Lepesqueur]

Olá Sidney!
Tudo bem?
Tenho uma imagem raster com três valores de pixel 0 1 2. Como eu faria para contar quantas células de cada valor eu possuo? Isto é um trabalho de cálculo de área sobre um endvi.
Agradeço o help.

[Sidney - Geosaber]

Léo,

O melhor caminho para você é poligonizar as classes e calcular a área no vetor (supondo que você já reclassificou seu grid para 3 classes de valores, como você mencionou, use a ferramenta do SAGA > Shapes - Grid > Vectorising Grid Classes).

Felipe,

Eu usei um SRTM da Embrapa, convertido o SRC para UTM (ou pode deixar em Geográficas mesmo), mas sempre trabalho com formato GeoTIFF.

Lilian,

Eu peguei um SRTM (Brasil em Relevo - Embrapa), recortei uma pequena área (GDALTools) e gerei pontos aleatórios (fTools).

Rodei o complemento Point Sample Tools e funcionou corretamente, criando um novo shape de pontos com o valor dos pixels na tabela de atributos.

Como eu sempre reclamo, o Python no QGIS, só dá páu!

Você pode usar o SAGA (via SEXTANTE) para fazer o mesmo procedimento, também funcionou perfeitamente aqui no meu Win7x64, no mesmo grupo de ferramentas que eu indiquei ao Leonardo: SAGA > Shapes - Grid > Add Grid Values to Point ou Grid Values to Point (para o SAGA gerar os pontos, sem você ter um shape prévio).

Boa sorte à todos,

Sidney Goveia

http://www.geosaber.com.br

[Felipe Dantas]

Obrigado Sidney pelas informações. São pessoas assim como vocês e todos que fazem parte desse grupo que vão fazer e já fazem a diferença neste mundão de meu Deus.
Pondo em prática aqui os conhecimentos adquiridos !

--

--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 

---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

--

Felipe Ramos Dantas

Graduando em Tecnologia  em Geoprocessamento - IFPI

Contatos: tim (86) 9921-5742 /  claro (86) 9411-4372

Respeite o ciclista.

Uma bike a mais, um carro a menos.
..... __@
...._ \ >_
...(_)/ (_)___________________________

_________
░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

[Leonardo Lepesqueur]

Beleza Sidney. Que bom contar com pessoas que tem esta atitude colaborativa. Tenha certeza que voce esta fazendo do mundo  um lugar melhor.
Vou baixar o SAGA  tentar executar o procedimento. Depois eu conto como foi.
Abraços.

[Sidney - Geosaber]

Use o SAGA diretamente no QGIS via SEXTANTE, é muito mais ágil e você pode criar resultados temporários (camada virtual - útil para não ficar guardando arquivos intermediários) e o final você escolhe para salvar o arquivo e deixar gravado.

Veja mais como fazer a instalação do SAGA e OTB no SEXTANTE em sistemas Windows 7: http://www.geosaber.com.br/2013/04/opus-sig-sextante-no-quantum-gis.html

Sidney Goveia

http://www.geosaber.com.br