Mapas estructurales en Qgis - Python

1. Revisión de la información

1.1 Colección de la información

1.2 Limpieza de datos

2. Manipulación de datos

2.1 Importar librerías

import pandas as pd from shapely.geometry import Point import geopandas as gpd

2.2 Leer la información

#Importar el archivo corregido para ser leido en pandas data = pd.read_csv(r"/work/Mapas_volve/Horizonte/Hugin_Fm_Base_DEPTH_copia.dat") #Vamos a hacer un chequeo de la informacion cargada data.info()

<class 'pandas.core.frame.DataFrame'>
RangeIndex: 186756 entries, 0 to 186755
Data columns (total 5 columns):
 #   Column  Non-Null Count   Dtype  
---  ------  --------------   -----  
 0   at1     186756 non-null  float64
 1   at2     186756 non-null  float64
 2   x       186756 non-null  float64
 3   y       186756 non-null  float64
 4   z       186756 non-null  float64
dtypes: float64(5)
memory usage: 7.1 MB

2.3 Transformar la información

# Combinar las columnas X, Y y Z al objeto geometrico Point() de la libreria shapely data['coordinates'] = data.apply(lambda x: Point((float(x.x), float(x.y), float(x.z))), axis=1) data.info()

<class 'pandas.core.frame.DataFrame'>
RangeIndex: 186756 entries, 0 to 186755
Data columns (total 6 columns):
 #   Column       Non-Null Count   Dtype  
---  ------       --------------   -----  
 0   at1          186756 non-null  float64
 1   at2          186756 non-null  float64
 2   x            186756 non-null  float64
 3   y            186756 non-null  float64
 4   z            186756 non-null  float64
 5   coordinates  186756 non-null  object 
dtypes: float64(5), object(1)
memory usage: 8.5+ MB

# Convertir dataframe a geodataframe Hugin_base = gpd.GeoDataFrame(data, geometry='coordinates') Hugin_base.info()

<class 'geopandas.geodataframe.GeoDataFrame'>
RangeIndex: 186756 entries, 0 to 186755
Data columns (total 6 columns):
 #   Column       Non-Null Count   Dtype   
---  ------       --------------   -----   
 0   at1          186756 non-null  float64 
 1   at2          186756 non-null  float64 
 2   x            186756 non-null  float64 
 3   y            186756 non-null  float64 
 4   z            186756 non-null  float64 
 5   coordinates  186756 non-null  geometry
dtypes: float64(5), geometry(1)
memory usage: 8.5 MB

# Convertir el Geodataframe a .shp Hugin_base.to_file(r'Horizon_shp/Hugin_base.shp', driver='ESRI Shapefile')

/shared-libs/python3.7/py-core/lib/python3.7/site-packages/ipykernel_launcher.py:2: UserWarning: Column names longer than 10 characters will be truncated when saved to ESRI Shapefile.

3. Visualización en QGIS

3.1 Instalar QGIS

3.2 Cargar información a QGIS

3.3 Configurar el CRS

3.4 Crear mapa de contornos

Luego de esto vamos a crear mapa de contornos utilizando el plugin "Countour Plugin". Este Plugin lo vamos a encontrar en el menú principal en Complementos -> Administrar e Instalar Complementos. Luego se abrirá una ventana (figura 6) donde vamos a buscar el plugin Contour plugin y le damos en Instalar Complemento.

Al instalar el complemento lo vamos a ver en la barra de tareas de QGIS (Circulo azul figura 7). Es con este plugin vamos a crear el mapa de contornos.

Luego de esto vamos a configurar el mapa de contornos de acuerdo a la figura 8

Como resultado tendremos un objeto o capa del mapa de contornos y un objeto de líneas de contorno, figura 9.

Ahora que tenemos una representación del subsuelo de la formación Hugin Base, vamos a visualizarlo en un mapa a través de las "Composición de Impresión" de Qgis. Para esto vamos al menú principal en Proyecto -> Nueva Composición de impresión como se muestra en la figura 10. Se abrirá una ventana para modificar y crear el mapa