[Spanish] obtener la extension de los fenomenos(features)

[Pedro-Juan Ferrer Matoses]

Respuesta de sabado por la tarde

Tengo un script en bash que utiliza ogr2ogr que CASI lo hace... pero
extrañamente me deja algunas geometrías en blanco...

lo pongo aquí a ver si a alguien se le ocurre que le pasa

Nota: el método es gocho, pero válido
********************************************************************
#!/bin/bash

fidlist=`ogrinfo -ro -al -geom=SUMMARY
Poligonos_Peninsula_IB_CE_ML_ED50_LL.shp|grep -E ":[0-9]+$"|grep -o -E
"[0-9]+$"`

echo '<OGRVRTDataSource>
    <OGRVRTLayer name="salida">
        <SrcDataSource>salida.csv</SrcDataSource>
        <GeometryType>wkbMultiPolygon</GeometryType>
        <LayerSRS>WGS84</LayerSRS>
        <GeometryField encoding="WKT" field="GEOM"/>
    </OGRVRTLayer>
</OGRVRTDataSource>' > salida.vrt

for i in $fidlist
do
geom=`ogrinfo -fid $i -q Poligonos_Peninsula_IB_CE_ML_ED50_LL.shp
Poligonos_Peninsula_IB_CE_ML_ED50_LL|grep MULTIPOLYGON`
echo "ID,GEOM" > salida.csv
echo $i,'"'$geom'"' >> salida.csv
extent=`ogrinfo -ro -so -al -geom=SUMMARY -fid $i salida.vrt|grep
Extent|awk '{print $2 $3 $5 $ 6}'`
echo $i $extent >> extents.txt
done

********************************************************************

Un saludo

-- 
Pedro-Juan Ferrer Matoses
Valencia (España)