PERFIL TOPOGRAFICO EN QGIS

Para hacer un perfil topográfico en QGIS se necesita un Modelo Digital de Elevaciones, y un Shape de líneas para trazar tu perfil como en la imagen sig.:

 

Con el plugin qProf se puede graficar sobre la línea las alturas del MDE, seleccionando la capa que contiene la línea, el MDE, el intervalo de muestreo y posteriormente dar clic en Calculate profile

Si todo sale bien aparece la siguiente ventana y dar clic en aceptar

Si se quiere visualizar el perfil antes de guardarlo a un archivo shape o de texto dar clic en la pestaña de Plot profiles.

 

Si se quiere visualizar solo el perfil seleccionar height y si se requiere también las pendientes seleccionar slope

 

Para guardar los datos a un archivo shape seleccionar 2D point shp

para guardar los datos a un archivo de texto seleccionar csv

POLIGONO 3D CON GDAL/OGR

Para hacer un archivo shp de un poligono con ogr2ogr se necsita:

-archivo *.csv en formato WKT como sigue
POLIGONO ((x y z,x y z,xn yn zn,.....))
con las primeras y ultimas coordenadas iguales

-archivo *.vrt en formato XML

<OGRVRTDataSource>
    <OGRVRTLayer name="LT3_geol3d">
        <SrcDataSource relativeToVRT="1" shared="1">C:\PROGSCRIPT\CONVERSION\QGIS\LT3_geol3d.csv</SrcDataSource>
        <GeometryType>wkbPolygon25D</GeometryType>
        <GeometryField>encoding="WKT" field="WKT"</GeometryField>
        <LayerSRS>PROJCS[&quot;NAD_1927_UTM_Zone_15N&quot;,GEOGCS[&quot;GCS_North_American_1927&quot;,DATUM[&quot;North_American_Datum_1927&quot;,SPHEROID[&quot;Clarke_1866&quot;,6378206.4,294.9786982]],PRIMEM[&quot;Greenwich&quot;,0.0],UNIT[&quot;Degree&quot;,0.0174532925199433]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;False_Easting&quot;,500000.0],PARAMETER[&quot;False_Northing&quot;,0.0],PARAMETER[&quot;Central_Meridian&quot;,-93.0],PARAMETER[&quot;Scale_Factor&quot;,0.9996],PARAMETER[&quot;Latitude_Of_Origin&quot;,0.0],UNIT[&quot;Meter&quot;,1.0]]</LayerSRS>
        <Field name="CLAVE" type="String" src="CLAVE" width="10"/>
        <Field name="NOMBRE" type="String" src="NOMBRE" width="10"/>
    </OGRVRTLayer>
</OGRVRTDataSource>


-Comando en la consola de comandos como sigue
ogr2ogr -f "ESRI Shapefile" C:\carpetadestino C:\archivo*.vrt