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)
_______________________________________________
Spanish mailing list
http://lists.osgeo.org/mailman/listinfo/spanish
http://es.osgeo.org
http://twitter.com/osgeoes
