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

[Pedro-Juan Ferrer Matoses]

Uy que lio de saltos de línea, mejor lo mando en HTML

********************************************************************
#!/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
********************************************************************
-- 
Pedro-Juan Ferrer Matoses
Valencia (España)
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: http://lists.osgeo.org/pipermail/spanish/attachments/20110122/a9bca303/attachment-0001.html