[Gfoss] spatialite

pigreco pigrecoinfinito a gmail.com
Dom 27 Maggio 2018 18:35:47 CEST


Maurizio Trevisani wrote
> Non l'ho testato e potrebbe essermi sfuggito qualche errore:
> 
> prova un po'
> 
> -------------------------
> 
> select 'Creazione indice spaziale su ', 'nodes_all','geom',
> coalesce(checkspatialindex('nodes_all','geom'),CreateSpatialIndex('nodes_all','geom'));
> 
> select dropgeotable('_templines');
> create table _templines(pk_uid integer primary key autoincrement, id
> integer);
> select addgeometrycolumn('_templines','geom',(SELECT cast(srid as
> integer) FROM geometry_columns WHERE lower(f_table_name) =
> lower('strade') AND lower(f_geometry_column) = lower('geom')),
> 'multilinestring', 'xy');
> 
> insert into _templines (id, geom)
> SELECT a.pk AS id,
> casttomulti(ST_Split( a.geom, (select
> casttomultipoint(st_collect(b.geom)) from "points_snapped" as b where
> b.rowid in (select rowid from spatialindex where f_table_name =
> 'points_snapped' and f_geometry_column = 'geom' and search_frame =
> a.geom)))) as geom
> from "strade" a;
> 
> select dropgeotable('lines_split');
> select ElementaryGeometries( '_templines' , 'geom' , 'lines_split' ,
> 'out_pk' , 'out_multi_id' ,  1, ) as num, ' lines splitted' as label;
> select createspatialindex('lines_split','geom');
> select dropgeotable('_templines');
> 
> select UpdateLayerStatistics('lines_split');
> 
> -------------------------
> Ciao,
> Maurizio

Ciao e grazie.

Ho appena provato lo script [0] e ottengo la geotabella 'lines_split' ma
lanciando il check geometry:
count(*)|Geometry Type("geom")|Srid("geom")|CoordDimensio("geom")
849       |Null                             | Null            |Null

quindi come posso proseguire?

[0]  https://gist.github.com/pigreco/a00164ded1a67e4b067b3a8c167a80fa


-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://gfoss-geographic-free-and-open-source-software-italian-mailing.3056002.n2.nabble.com/


Maggiori informazioni sulla lista Gfoss