Executar string com Python

Esses dias me deparei com um problema em Python, queria executar uma string. Eu tinha uma string que era um teste de condicionais e uma lista que seria os nomes das variáveis que estavam sendo usadas nas condicionais. Para executar os testes de condicionais que estavam como string, utilizei o comando eval() e para declarar as novas variáveis que iriam receber um array, usei o exec().

Por exemplo:
>>> import numpy as np
>>> d = "f"     #Criei uma variavel que recebe a string
>>> array = np.array([1, 2, 3, 4])    #criei o array do numpy
>>> condicionais = "np.where( a == 1, 10, np.where (a >= 1, 20, 30))"    #O teste condicionais em string
>>> exec (d+"=array")  #O comando exec vai fazer com que a variável  f receba o array
>>> print f      #Então quando mandar escrever a variável f aparecerá os valores do array
array([1, 2, 3, 4])
>>> eval(condicionais)   #Para que as condicionais seja executada utilizá-se o eval, esse comando avaliará a string
array([10, 20, 20, 20])  #Resultado das condicionais

Bom trabalho!

Trabalhando com Python + R

Trabalhar com gráficos estatísticos no Python ficou mais fácil a partir da biblioteca rpy2. Neste post irei demostrar como instalar o rpy2 no Python 2.7 e plotar alguns gráficos. Faz alguns dias que comecei a migrar para Ubuntu e para minha felicidade foi fácil instalar essa biblioteca, sendo que no Windows estava encontrando problemas para utilizá-la.

1- Vamos instalar o rpy2 no Python, digitando o código abaixo no terminal.

sudo apt-get install python-rpy2

2- Execute o python no terminal

python

3- Os próximos comandos serão executados em ordem no Python pelo terminal

>>> from rpy2 import robjects
>>> r = robjects.r
>>> y = r.rnorm(10) # Gerando 10 números aleatórios
>>> r.boxplot(y) #Plotando o boxplot
>>> r.barplot(y) #Plotando o barplot

Os resultados estão abaixo:

BoxPlot

BarPlot

Lendo dados de um shapefile com OGR + Python

 Nesse post vamos ler um arquivo shapefile e escrever os registros de um determinado campo da tabela. Para isso iremos utilizar a biblioteca OGR que serve para leitura de dados geoespaciais do tipo vetorial, seu código está inserido na biblioteca GDAL e pode ser baixado do site OSGeo. Essa biblioteca suporta vários formatos e entre ele está a extensão .shp.
 O código abaixo está com os comentários na cor verde.

O shapefile é referente aos limites municípais do RS, a leitura foi do campo onde contém os nomes do municípios. O resultado será o contador mais os nomes dos municípios.

 Qualquer dúvida deixem um comentário!!!
Obrigado.

Array de números em Python

Este post tem como objetivo demonstrar como iniciar  e  acessar um Array. Um Array é uma estrutura de dados que contém vários elementos de mesmo tipo de dado, ele pode ser acessado pelo seu índice.

Na primeira linha foi inciada a variável "a" que armazenará 4 linhas e 4 colunas. Na terceira e quarta linha foi criado um loop  que percorerrá esse array. A quinta linha grava a posição da variável "a [ i ] [ j ]" e recebe a soma dessa posição. "i+j". Para finalizar, na sexta linha o programinha escreve qual linha e coluna que está o valor calculado anteriormente.
O resultado fica assim:

Contando o número de vértices

Existem muitos recursos do ArcGis que são executados na tabela pelo Field Calculator, alguns clássicos são o cálculo da área de um polígono, comprimento e centróide. Esse post demonstrará como executar um comando que contará o número de vértices de um polígono ou linha.

Crie uma coluna na tabela Options ->Add Field

Nome: Vertices

Type: Long Integer

Precision: 10 (número de algarismos)

 

Depois que a coluna foi criada pressione o botão direito sobre a coluna e selecione Field Calculator.

Com o Field Calculator aberto selecione a opção Advanced.

Copie e cole os comandos descritos abaixos no Fiel Calculator.

Dim pPointCount As IPointCollection

Set pPointCount = [SHAPE]

pPointCount.PointCount 

Transformando linhas da tabela em shapefile – SAGA-GIS

Este post demonstrará como transformar as linhas de uma tabela em shapefiles ou melhor criar novos shapefiles através das linhas da tabela. Esse exemplo é útil para quando criamos um shapefile para vetorizar vários polígonos e por algum motivo queremos trabalhar apenas com um desses polígonos ou trabalhar com esses polígonos separados.

A ferramenta está em Shapes – Tools -> SplitTable/Shapes by Atribute

Seleciona o shapefile que será dividido e a partir de qual atributo (coluna) da tabela será criado os novos vetores. No meu caso será o vetor do RS e dividirei pelos nomes dos municípios.

Foi criado os novos vetores de acordo com os nomes que estavam na tabela.

Dividindo polígono com ArcGis

Adicione o shapefile no programa.

 

Abra a edição do shapefile com a ferramenta Editor -> Start Editing.

Com edição aberta selecione na janela Task a opção Cut Polygon Features

1 - Selecione o polígono com a ferramenta Tools->Select Features

2 – Posteriormente selecione a ferramenta de edição o lápis (Stech Tool) e crie a linha que servirá para dividir o polígono, para finalizar pressione F2.