No conozco ArcPy, pero trabajo con archivos de forma y ráster en Python durante años
- Para procesar archivos shape en Python, hay muchos módulos como osgeo / ogr , Fiona , Pysal o Pyshp ( shpUtils es uno de ellos y no el más utilizado), y otros, vea Pypi: GIS y ejemplos en gis.stackexchange y muchos ejemplos en la web (no solo en inglés). La mayoría de ellos son mucho más antiguos que ArcPy (o arcgisscripting) ...
- para procesar ráster puede usar osgeo / gdal , el estándar
- Para procesar geometrías geoespaciales, hay bien proporcionado
- Para trazar las geometrías puede usar matplotlib y posiblemente descartes , "extensión" de matplotlib para áreas, pero también muchos, muchos otros módulos, vea Pypi: Trazado y módulos como mayavi para representación 3D (también matplotlib)
- También hay módulos como mapnik que le brindan directamente las posibilidades de 1) leer un archivo shape y 4) trazar con el módulo Pycairo .
Después de eso, es como un SIG:
- usa los módulos 1) para abrir, guardar los archivos de forma y llevar a cabo los tratamientos con otros módulos como numpy o scipy, si lo desea.
- Puede utilizarlo para manipular y analizar los objetos geométricos (búfer, etc.).
- puede usar matplotlib para trazar las geometrías, pero matplotlib no sabe lo que quiere trazar. Es su trabajo con los módulos 1) o 3) especificar qué trazar (atributos, etc.) y cómo.
Si quiero visualizar una determinada columna de mi archivo de forma, ¿cómo puedo implementar esto en el código?
Por lo tanto, debe aprender matplotib y los otros módulos. Tienes que aprender ArcPy, es lo mismo ... (hay muchos tutoriales excelentes en la web, especialmente para matplolib, y es más fácil que ArcPy porque es Python puro).
Algunos ejemplos solo con Python
Mapa geológico (archivo de forma poligonal) con colores basados en un atributo
Puntos 3D (archivo de forma PointZ) con color basado en un atributo
Puntos 3D (archivo de forma de punto con z como atributo) y línea 3D (archivo de forma PolyLineZ) en un DEM y en un ráster desplegado sobre la superficie del DEM.
Perfil topográfico con valores z y colores basados en atributos (formaciones geológicas = sección transversal) del archivo de forma original (archivo de forma de polilínea)
DEM (GeoTIFF) con el módulo Mayavi2
DEM (ESRI ascii grid, .asc) y Point shapefiles (con z como atributo) con el módulo visvis
Perforaciones (búfer 3D de una polylineZ con colores basados en un atributo (formaciones geológicas), con una superficie de cuadrícula calculada con los módulos numpy y matplotlib a partir de un archivo de forma de puntos (con z como atributo), visualizado con el visvis del módulo