Não conheço o ArcPy, mas trabalho com shapefiles e raster em Python há anos
- Para processar shapefiles em Python, existem muitos módulos como osgeo / ogr , Fiona , Pysal ou Pyshp (o shpUtils é um deles e não o mais usado) e outros, consulte Pypi: GIS e exemplos em gis.stackexchange e muitos exemplos em a Web (não apenas em inglês). A maioria deles é muito mais antiga que o ArcPy (ou arcgisscripting) ...
- para processar raster você pode usar osgeo / gdal , o padrão
- Para processar geometrias geoespaciais, há formas bem torneadas
- Para plotar as geometrias, você pode usar matplotlib e possivelmente descartes , "extensão" do matplotlib para áreas, mas também muitos outros módulos, consulte Pypi: Plotagem e módulos como mayavi para representação 3D (matplotlib também)
- Existem também módulos como o mapnik, que oferecem diretamente as possibilidades de 1) ler um shapefile e 4) plotar com o módulo Pycairo .
Depois disso, é como um SIG:
- você usa os módulos 1) para abrir, salvar os shapefiles e realizar os tratamentos com outros módulos, como numpy ou scipy, se desejar.
- você pode usar formas para manipulação e análise de objetos geométricos (buffer, etc.).
- você pode usar o matplotlib para plotar as geometrias, mas o matplotlib não sabe o que deseja plotar. É seu trabalho com os módulos 1) ou 3) especificar o que plotar (atributos, etc,) e como.
Se eu quiser visualizar uma determinada coluna do meu shapefile, como posso implementar isso no código?
Portanto, você deve aprender matplotib e os outros módulos. Você precisa aprender o ArcPy, é o mesmo ... (existem muitos tutoriais excelentes na Web, especialmente para o matplolib, e é mais fácil que o ArcPy porque é um Python puro).
Alguns exemplos apenas com Python
Mapa geológico (arquivo de forma de polígono) com cores com base em um atributo
Pontos 3D (shapefile PointZ) com cores baseadas em um atributo
Pontos 3D (ponto shapefile com z como atributo) e linha 3D (PolyLineZ shapefile) em um DEM e em uma varredura colocada na superfície do DEM.
Perfil topográfico com valores z e cores com base em atributos (formações geológicas = seção transversal) do arquivo de forma original (arquivo de forma polilinha)
DEM (GeoTIFF) com o módulo Mayavi2
DEM (ESRI ascii grid, .asc) e Point shapefiles (com z como atributo) com o módulo visvis
Furos (buffer 3D de uma polilinhaZ com cores baseadas em um atributo (formações geológicas), com uma superfície de grade calculada com os módulos numpy e matplotlib a partir de um shapefile de pontos (com z como atributo), visualizado com o módulo visvis