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
![ingrese la descripción de la imagen aquí](https://i.stack.imgur.com/hm9DT.jpg)
Mapa geológico (archivo de forma poligonal) con colores basados en un atributo
![ingrese la descripción de la imagen aquí](https://i.stack.imgur.com/Rqiav.jpg)
Puntos 3D (archivo de forma PointZ) con color basado en un atributo
![ingrese la descripción de la imagen aquí](https://i.stack.imgur.com/b6NWR.jpg)
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.
![ingrese la descripción de la imagen aquí](https://i.stack.imgur.com/CoZQn.jpg)
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)
![ingrese la descripción de la imagen aquí](https://i.stack.imgur.com/9olre.png)
DEM (GeoTIFF) con el módulo Mayavi2
![ingrese la descripción de la imagen aquí](https://i.stack.imgur.com/FYEXG.png)
DEM (ESRI ascii grid, .asc) y Point shapefiles (con z como atributo) con el módulo visvis
![ingrese la descripción de la imagen aquí](https://i.stack.imgur.com/nXZJD.jpg)
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