[Gfoss] spatialite: errore RecoverGeometryColumn

lucadelu · April 14, 2011, 12:17pm

Eccoci di nuovo qui, ho creato una nuova tabella attraverso una select
con ST_Union() e group by, ora vorrei aggiungere la colonna geometrica
alle geometry_columns...

spatialite> select
RecoverGeometryColumn('acque_trento_finali','Geometry',32632,'MULTILINESTRING','XY');
RecoverGeometryColumn(): validation failed
RecoverGeometryColumn('acque_trento_finali','Geometry',32632,'MULTILINESTRING','XY')
------------------------------------------------------------------------------------
0

qualche suggerimento?

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

a.furieri · April 14, 2011, 12:36pm

On Thu, 14 Apr 2011 14:17:00 +0200, Luca Delucchi wrote

Eccoci di nuovo qui, ho creato una nuova tabella attraverso una
select con ST_Union() e group by, ora vorrei aggiungere la colonna geometrica
alle geometry_columns...

spatialite> select
RecoverGeometryColumn('acque_trento_finali','Geometry',32632,
'MULTILINESTRING','XY'); RecoverGeometryColumn(): validation failed
RecoverGeometryColumn('acque_trento_finali','Geometry',32632,
'MULTILINESTRING','XY')

------------------------------------------------------------------------------------

0

qualche suggerimento?

si, un sacco ...
verificare cosa hai veramente nella geometria
generata da ST_Union();

SELECT DISTINCT
  ST_GeometryType(Geometry),
  ST_Srid(Geometry),
  CoordDimension(Geometry)
FROM acque_trento_finali;

vedrai che probabilmente la Union ti ha
generato qualche geometria strana.
se semplicemente vedi che ci sono dei
LINESTRING (senza MULTI) basta che applichi
la "promozione" del tipo geometrico:

UPDATE acque_trento_finali
SET Geometry = CastToMulti(Geometry);

vedrai che poi dovebbe andare meglio

ciao Sandro

lucadelu · April 14, 2011, 1:08pm

Il 14 aprile 2011 14:36, <a.furieri@lqt.it> ha scritto:

vedrai che poi dovebbe andare meglio

chissà come mai ma è andato
Giusto per conoscenza è una cosa che capita anche su postgis? (se
fosse così sono stato sempre fortunato )

ciao Sandro

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

a.furieri · April 14, 2011, 3:37pm

On Thu, 14 Apr 2011 15:08:07 +0200, Luca Delucchi wrote

Il 14 aprile 2011 14:36, <a.furieri@lqt.it> ha scritto:

>
> vedrai che poi dovebbe andare meglio
>

chissÃ come mai ma Ã¨ andato
Giusto per conoscenza Ã¨ una cosa che capita anche su postgis? (se
fosse cosÃ¬ sono stato sempre fortunato )

si, puÃ² capitare anche su postgis (almeno, credo).
le varie ST_Union, ST_Buffer etc etc sono supportate
da GEOS: e sia postgis che splite si limitano a
delegare queste operazioni a GEOS.

Altro esempio "classico": tutti scommetterebbero
che ST_ConvexHull ritorni necessariamente un
POLYGON ... falsissimo ...
puÃ² anche tornare un POINT (se il point-set
comprende un unico punto) ma anche un
LINESTRING (se il point-set comprende
solamente due punti) etc etc
CosÃ¬ come ST_Difference puÃ² tornare un
LINESTRING (se i due poligoni si toccano
lungo un lato) ma anche un POINT (se i due
poligoni si toccano solo su un vertice).

ma sia splite che postgis supportano appunto
una serie di funzioni "casting" che aiutano
a forzare la classe della geometria prodotta
dell'operatore Spatial

ciao Sandro

lucadelu · April 29, 2011, 9:10am

Il 14 aprile 2011 14:36, <a.furieri@lqt.it> ha scritto:

On Thu, 14 Apr 2011 14:17:00 +0200, Luca Delucchi wrote

Eccoci di nuovo qui, ho creato una nuova tabella attraverso una
select con ST_Union() e group by, ora vorrei aggiungere la colonna geometrica
alle geometry_columns...

spatialite> select
RecoverGeometryColumn('acque_trento_finali','Geometry',32632,
'MULTILINESTRING','XY'); RecoverGeometryColumn(): validation failed
RecoverGeometryColumn('acque_trento_finali','Geometry',32632,
'MULTILINESTRING','XY')

------------------------------------------------------------------------------------

0

qualche suggerimento?

si, un sacco ...
verificare cosa hai veramente nella geometria
generata da ST_Union();

SELECT DISTINCT
ST_GeometryType(Geometry),
ST_Srid(Geometry),
CoordDimension(Geometry)
FROM acque_trento_finali;

vedrai che probabilmente la Union ti ha
generato qualche geometria strana.
se semplicemente vedi che ci sono dei
LINESTRING (senza MULTI) basta che applichi
la "promozione" del tipo geometrico:

UPDATE acque_trento_finali
SET Geometry = CastToMulti(Geometry);

vedrai che poi dovebbe andare meglio

l'altra volta aveva funzionato, ma questa volta con un'altra tabella no!
qualche altro suggerimento?

SELECT DISTINCT st_geometrytype(geom), st_srid(geom),
coorddimension(geom) from uso_suolo_trentino2008_new;
st_geometrytype(geom)|st_srid(geom)|coorddimension(geom)
||
MULTIPOLYGON|32632|XY

select count(geom) from uso_suolo_trentino2008_new where not isvalid(geom);
count(geom)
0

select count(geom) from uso_suolo_trentino2008_new where geom is null;
count(geom)
0

select RecoverGeometryColumn('uso_suolo_trentino2008_new','geom',32632,'MULTIPOLYGON','XY');
RecoverGeometryColumn(): validation failed
RecoverGeometryColumn('uso_suolo_trentino2008_new','geom',32632,'MULTIPOLYGON','XY')
0

ciao Sandro

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

a.furieri · April 29, 2011, 9:35am

On Fri, 29 Apr 2011 11:10:43 +0200, Luca Delucchi wrote

l'altra volta aveva funzionato, ma questa volta con un'altra tabella
no! qualche altro suggerimento?

SELECT DISTINCT st_geometrytype(geom), st_srid(geom),
coorddimension(geom)
FROM uso_suolo_trentino2008_new;

MULTIPOLYGON|32632|XY

ok, sono tutti MULTIPOLYGON 2D XY e lo srid Ã¨ 32632

select count(geom)
from uso_suolo_trentino2008_new
where geom is null;

0

e non ci sono NULL

select RecoverGeometryColumn('uso_suolo_trentino2008_new','geom',
32632,'MULTIPOLYGON','XY');

RecoverGeometryColumn(): validation failed

eh no: *deve* funzionare. io non ci vedo nessun errore SQL
non so cosa dirti: bug mi pare quasi impossibile, perchÃ¨
RecoverGeometryColumn per validare esegue esattamente
la query che tu hai giÃ lanciato a mano .. boh !?!?!?!

ciao Sandro

lucadelu · April 29, 2011, 12:06pm

Il 29 aprile 2011 11:35, <a.furieri@lqt.it> ha scritto:

eh no: *deve* funzionare. io non ci vedo nessun errore SQL
non so cosa dirti: bug mi pare quasi impossibile, perchè
RecoverGeometryColumn per validare esegue esattamente
la query che tu hai già lanciato a mano .. boh !?!?!?!

l'unica cosa di strano è quando eseguo

SELECT DISTINCT st_geometrytype(geom), st_srid(geom),
coorddimension(geom) from uso_suolo_trentino2008_new;
st_geometrytype(geom)|st_srid(geom)|coorddimension(geom)
||
MULTIPOLYGON|32632|XY

c'è la prima riga || invece su una tabella corretta non c'è

SELECT DISTINCT st_geometrytype(Geometry), st_srid(Geometry),
coorddimension(Geometry) from uso_suolo_trentino2008;
st_geometrytype(Geometry)|st_srid(Geometry)|coorddimension(Geometry)
POLYGON|32632|XY

penso il problema sia li, ma non so come correggerlo

ciao Sandro

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

a.furieri · April 29, 2011, 12:57pm

On Fri, 29 Apr 2011 14:06:34 +0200, Luca Delucchi wrote

l'unica cosa di strano Ã¨ quando eseguo ...
c'Ã¨ la prima riga || invece su una tabella corretta non c'Ã¨

penso il problema sia li, ma non so come correggerlo

ok, ora Ã¨ chiaro: scusami, ma prima non ci avevo
proprio fatto caso
sono NULL

vedo che tu usi spatialite CLI, che by default
ignora i nulls: se li vuoi visualizzare devi
attivarli esplicitamente:
.nullvalue NULL

concludendo: non hai nessun NULL tra le tue geometrie,
lo avevamo giÃ verificato prima.
perÃ² hai qualche valore assolutamente sballato, cosÃ¬
tanto sballato che non viene neppure riconosciuto come
una Geometry binaria valida: prova ad indentificarli cosÃ¬

SELECT * FROM my_table
WHERE ST_GeometryType(geom) IS NULL;

ciao Sandro

lucadelu · April 29, 2011, 1:08pm

Il 29 aprile 2011 14:57, <a.furieri@lqt.it> ha scritto:

SELECT * FROM my_table
WHERE ST_GeometryType(geom) IS NULL;

ok, trovati... per metterli a posto devo usare sempre la gui?

ciao Sandro

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

a.furieri · April 29, 2011, 1:48pm

On Fri, 29 Apr 2011 15:08:54 +0200, Luca Delucchi wrote

Il 29 aprile 2011 14:57, <a.furieri@lqt.it> ha scritto:

> SELECT * FROM my_table
> WHERE ST_GeometryType(geom) IS NULL;

ok, trovati... per metterli a posto devo usare sempre la gui?

dipende da "cosa" ci hai trovato dentro ....
a naso, l'unica cosa che parrebbe sensata Ã¨
una bella DELETE, ma bisognerebbe verificare

ciao Sandro

lucadelu · April 29, 2011, 2:19pm

Il 29 aprile 2011 15:48, <a.furieri@lqt.it> ha scritto:

dipende da "cosa" ci hai trovato dentro ....
a naso, l'unica cosa che parrebbe sensata è
una bella DELETE, ma bisognerebbe verificare

già provata ma si è portato via un po' troppe features

ciao Sandro

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

a.furieri · April 29, 2011, 2:36pm

On Fri, 29 Apr 2011 16:19:33 +0200, Luca Delucchi wrote

> a naso, l'unica cosa che parrebbe sensata Ã¨
> una bella DELETE, ma bisognerebbe verificare
>

giÃ provata ma si Ã¨ portato via un po' troppe features

riassumendo: quindi il tuo DB ha un sacco di righe
con contengono "roba strana", che spatialite si
rifiuta tassatavimente di riconoscere come
geometrie valide

di grazia, ma cosa Ã¨ successo ?
che procedura e che strumenti hai utilizzato ?

BTW a me a questo punto pare abbastanza chiaro
che cotesto DB Ã¨ fritto ... fai sicuramente
prima a buttarlo e ricaricarlo piuttosto
che tentare di aggiustarlo

ciao Sandro