Discussion :

[benbannock]

Bonjour à tous,

Voici mon petit problème.

J'utilise une fonction ( .Reduce(x) ) sur un type geometry. Ce sont des multipolygones. La fonction marche et ne renvoie pas d'erreur avec ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 update Zone_VE set geom_simpl = geom.Reduce(2000);

Le problème s'est que je récupère une géométrie de type multipoint, , je perds donc tout mes multipolygones....

J'ai lu dans la doc que tout les types de geometry était pris en compte, mais ca devrait normalement conserver le type de géométrie.

Si quelqu'un à une idée sur le problème ou même (soyons fou) une solution je suis preneur,

Merci d'avance !

[SQLpro]

Jouez sur le niveau de réduction...
Exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
DECLARE @GEO GEOMETRY
SET @GEO = 'POLYGON ((0 0, 1 1, 2 2 , 3 3, 3 0, 0 0))'
SELECT @GEO.Reduce(1).STAsText()

Résultat :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

POLYGON ((0 0, 3 3, 3 0, 0 0))

A +

[benbannock]

Merci tout d'abord de la réponse rapide,

Bon j'ai essayé avec 1 en paramètre :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

update Zone_VE set geom_simpl = geom.Reduce(1);

et du coup j'ai un résutat plutôt bizarre :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
------------------------ | --------------------------------------------...
Geometry original:       | Résultat :
MULTIPOLYGON (((13.55157 | GEOMETRYCOLLECTION (POINT (13.55157470703125...
MULTIPOLYGON (((-74.1210 | GEOMETRYCOLLECTION (POINT (-74.12109375 83.0...
MULTIPOLYGON (((140.9739 | GEOMETRYCOLLECTION (POINT (140.9739685058593...
MULTIPOLYGON (((58.18933 | GEOMETRYCOLLECTION (LINESTRING (59.375610351...
MULTIPOLYGON (((19.82496 | GEOMETRYCOLLECTION (POINT (19.82496643066406...
MULTIPOLYGON (((-1.83583 | GEOMETRYCOLLECTION (POINT (-1.83583450317382...
MULTIPOLYGON (((35.47556 | GEOMETRYCOLLECTION (POLYGON ((34.90316390991...
MULTIPOLYGON (((8.911651 | GEOMETRYCOLLECTION (POINT (8.911651611328125...

Je n'arrive toujours pas à conserver mes multipolygones ....
flûte!


En fait, c'est bon ca marche avec :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

update Zone_VE set geom_simpl = geom.Reduce(0.00001);

Merci encore , j'espère que la généralisation de mes entités sera efficace même avec ce niveau de réduction !

[SQLpro]

Vous n'avez pas des polygones "à plat" ?

A +

[benbannock]

Eeuuuuh des polygones "à plat" ??

J'ai des multipolygones, souvent une entité avec ses extensions regroupé sous un seul et unique identifiant. exemple : La France avec toutes ses îles...

Ca change quelque chose?


Sinon j'ai été un peu vite en besogne. Dans mes multipolygones, certains polygones sont très petits (genre les micro-îles), avec la fonction .Reduce(x) ces derniers suivant leurs taille sont transformés en Linestring ou carrément en point s'ils sont vraiment petits.

Il est donc normal que je récupère un autre type de géométrie en sortie.

J'envisage donc deux solutions hypothétiques :
- L'une serait d'utiliser la fonction .Reduce(x) en l’empêchant de travailler sur des polygones à trois points (ce qui empêcherait d'avoir des linestring ou des point) mais cette solution me parait impossible.
-La seconde serait une fonction qui supprimerait mes points et mes linestring dans mes geometries en sortie qui sont des GEOMETRY COLLECTION. En gros la transformer en Multipolygone.

Si quelqu'un sait comment faire pour la seconde solution....

[SQLpro]

Dans une GEOMETRY COLLECTION, chaque objet géométrique a une position ordinale.

Exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
DECLARE @GEO GEOMETRY
SET @GEO = 'GEOMETRYCOLLECTION(POINT(0 0), POLYGON ((0 0, 1 1, 2 2 , 3 3, 3 0, 0 0)))'
DECLARE @N INT
SET @N = @GEO.STNumGeometries()
WHILE @N >= 1
BEGIN
   SELECT @GEO.STGeometryN(@N).STAsText()
   SET @N -= 1   
END

Peut être devriez vous venir à mon cours : http://www.orsys.fr/pdf-auto/pdfCours/TIA.pdf

A +

[benbannock]

OK, pour la position, du coup j'ai voulu faire ma deuxième solution, j'arrive au curseur suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
declare @idx_tot integer
declare @num_entite integer
declare @champ1 integer
declare	@idx integer
declare	@type varchar(50);
 
declare Cursor_geom cursor   FOR
select id_z_ve, Geom_simpl.STNumGeometries() from Zone_VE 
open Cursor_geom
fetch next FROM Cursor_geom INTO @num_entite, @champ1 
while @@fetch_status = 0
	begin
	set @idx_tot = @champ1
	set	@idx = 1
		while @idx <= 5--@idx_tot
			begin
				set @type = exec('select geom_simpl.STGeometryN('+@idx+').STGeometryType() from Zone_VE')
				if @type <> 'Polygon'
					if @type <> 'MultiPolygon'
					exec('update Zone_VE set geom_simpl =  geom_simpl.STDifference(select geom_simpl.STGeometryN('+@idx+') 
										from Zone_VE )
										 where id_z_ve = '+@num_entite+' ')
				--print @idx
				set @idx = @idx+1
			end
		fetch next FROM Cursor_geom INTO  @num_entite,@champ1 
    end
close Cursor_geom
deallocate Cursor_geom

Évidemment ce dernier ne fonctionne pas puisque la ligne suivante n'est pas correcte :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

set @type = exec('select geom_simpl.STGeometryN('+@idx+').STGeometryType() from Zone_VE')

Je voudrais attribuer le résultat de la requête à la variable @type, comment puis-je faire ??...