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

Sat, 22 Jan 2011 08:55:15 -0800 

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)

_______________________________________________
Spanish mailing list
http://lists.osgeo.org/mailman/listinfo/spanish
http://es.osgeo.org
http://twitter.com/osgeoes

Responder a