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