Sistemas de Informação Geográfica

12
Gdal Dataset.ReadAsArray () trava o Python

Estou usando o Python 2.6.5 (32 bits) com o Numpy 1.3 e o Gdal 1.9.1 instalados no Windows 7 de 64 bits. Estou tentando ler um conjunto de dados de varredura Imagine de 800 MB Imagine (.img) em uma matriz Numpy para fazer uma álgebra de varredura, mas assim que executo o código a seguir, o...

12
Buffer em metros / km usando camadas WGS84 usando QGIS

Estou usando o QGIS 1.7.4-Wroclaw e trabalhando com uma camada SHP no WGS84. Quero fazer um buffer para um certo número de quilômetros dessa camada usando a ferramenta de buffer ftools. Entendo que a ferramenta de buffer sempre usa as unidades de camada, que para WGS84 são graus decimais. Todos...

12
Lendo geodatabase de arquivo usando R?

O geodatabase de arquivo (fgdb) inclui várias tabelas de geodatabase de arquivo. Tanto quanto eu sei, eles existem como dbfarquivos, mas estão dentro de um arquivo Database.gdb. No ArcCatalog, o caminho do arquivo é semelhante C:\Users\...\Database.gdb\Stats_AA. Como ler todos esses dbfarquivos...

12
Como obter interativamente as coordenadas de um ponto no QGIS

Eu tenho um shapefile com uma única camada do tipo point. Está no CRS do Google (EPSG: 900913). Quando uso a ferramenta 'Identify Feature', ele fornece muitas informações sobre o ponto clicado, mas não consigo ver as coordenadas. Existe uma maneira de descobrir a latência de um ponto...

12
Removendo linhas no shapefile em R

Importei um shapefile para o R e juntei-o a uma tabela. Meu shapefile contém todos os IDs do censo, enquanto minha tabela contém apenas os IDs selecionados. Agora estou tentando excluir todas as linhas que não obtiveram correspondência. É assim que meu conjunto de dados se parece (estou tentando...

12
Como exibir linhas sobrepostas?

Eu tenho várias camadas de linhas que representam rotas de ônibus e muitas delas se sobrepõem (ou seja, perto de terminais, locais de transferência etc.). Existe uma função no Qgis que exibe os dois recursos próximos / paralelos? No link de exemplo, as linhas Vermelha e Roxa são duas rotas...

12
Por que o script do ArcPy é lento?

Eu tenho um script arcpy simples para atualizar um campo em um shapefile de ponto com informações do recurso de polígono em que ele está. Demora 9 minutos para fazer 100 pontos no arcpy, mas uma junção espacial no arcmap é instantânea. Tenho certeza de que existe uma maneira rápida e estabelecida...

12
Como converter um shapefile para WKT?

Estou tentando converter um shapefile no formato simples Texto conhecido (WKT) no exemplo abaixo (em http://en.wikipedia.org/wiki/Well-known_text ). MULTIPOLYGON (((40 40, 20 45, 45 30, 40 40)), ((20 35, 45 20, 30 5, 10 10, 10 30, 20 35), (30 20, 20 25, 20 15, 30 20))) Sei que uso o QGIS para...