[Spanish] obtener la extension de los fenomenos(features)
Pedro-Juan Ferrer Matoses
pferrer at osgeo.org
Sat Jan 22 11:53:49 EST 2011
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)
More information about the Spanish
mailing list