Le portail francophone de la géomatique

Pardon pour cette imprécision,

en fait je voudrais la distance du point au segment du polygone qui lui est le plus proche
....par contre le réponse ne vient toujours pas.....

bonjour,

finalement j'ai changé de méthode, je suis parti sur une grille, ce qui me permet de faire des distance de centroides de polygones à centreoides de polygone.

le probleme c'est que je n'arrive pas à automatiser ma démarche.

j'ai deux tables: sel et selvide, qui sont issues d'une troisième qui s'appelle grille
sel contient les case "pleines" et selvide les cases "vides"

je voudrais donc calculer la distance en tre les cases de sel et celle de selvide. voici ma routine:

'------------------------------------------------
dim dist as float
dim nbsel, nbselvide, i,j as integer

nbsel=tableinfo(sel, tab_info_nrows)
nbselvide=tableinfo(selvide, tab_info_nrows)

'sel.col3 et sel.col4 sont les X et Y des centroides des objets de la table "sel"
'selvide.col3 et selvide.col4 sont les X et Y des centroides des objets de la table "selvide"

for i=1 to nbselvide
fetch rec i from selvide
                     for j=1 to nbsel
    fetch rec j from sel

'dist=((x2-x1)²+(y2-y1)²)^(1/2)
dist=(((sel.col3-selvide.col3)*(sel.col3-selvide.col3))+((sel.col4-selvide.col4)*(sel.col4-selvide.col4)))^(1/2)

    update selvide set distance=abs(dist)
    next
next
'----------------------------------------------------------------------------------------

et la c'est le drame, je me retrouve avec une valeur de "0" dans ma colonne distance


d'ou cela peut-il venir?

merci pour votre aide

merci pour tes conseils Pierre

c'est vrai, autant pour moi

je n'ai pas encore trouvé mais je vais chercher

j'ai effectué la formule suivante : 1.25+(sqr(sel.col3-selvide.col3)+sqr(sel.col4-selvide.col4))^(1/2)
mon champ se met a jour avec la valeur 1.25

le probleme vient donc de mes colonnes. mais je ne comprend pas pourquoi cela ne fonctionne pas alors qu'en remplacant par des valeurs (au lieu de mes colonnes) cela fonctionne
un autre petit souci, j'ai la même valeurs partout, du coup un problème dans l'indentation.......?

merci

bon alors en fait, j'ai appelé mes colonnes dans un "note" et je m'aperçois que sel.col3 et selvide.col3 ont la même valeur, idem pour sel.col4 et selvide.col4. du coup c'est normal que je trouve 0 au final!!!!!
je vais essayer de trouver le problème

..............comme quoi, juste en ecrivant on trouve des solutions........

Apparement non, car la table "sel" possede 830 enregistrement et la table "selvide" 5090 pour un total de 5920 pour ma table "grille".
sur ce point la c'est tout bon

ce que je voudrai, c'est la distance d'un polygone de la table "selvide" au plus proche polygone de la table "sel" (par l'intermédiaire des centroïdes), et cela pour chaque polygone de la table "selvide".

je m'y suis peut etre mal pris dans ma démarche, mais j'essaye dans un premier temps de calculer une distance entre deux objets de 2 tables différentes, et c'est la que je coince....

en espérant avoir été plus clair...

Merci

Bonjour Damien,

en effet, cela fait quelque jours que j'essaye une autre méthode, qui justement est celle que tu m'indiques
voici mon code:

Dmin=1000000                    'bornage du dmin
i=0
fetch first from sel                'positionnement au premier enregistrement de la table sel
i=i+1
     do while not EOT(sel)
           di=sqr((selvide.X-sel.X)*(selvide.X-sel.X)+(selvide.Y-sel.Y)*(selvide.Y-sel.Y))

                             if di<=dmin then
                                dmin=di
                             end if

select*from selvide where rowid=i
update selection set distance=dmin      'l'enregistrement doit etre selectionné pour pouvoir etre mis a jour
fetch next from sel
     loop
note "fini"


avec ce code, j'arrive à mettre a jour mon premier enregistrement de la table selvide avec la bonne valeur de distance la plus courte entre le premier objet de la table sel et le premier objet de la table selvide.

maintenant je dois trouver comment faire la manip pour tout les objets de la table selvide


j'avance petit à petit

merci

je me suis encore une fois, mal exprimé:

j'arrive à mettre a jour mon premier enregistrement de la table selvide avec la bonne valeur de distance la plus courte entre le premier objet de la table selvide et TOUT les objets de la table sel.


j'avais fait des tests avec des "print" pour vérifier mes valeurs et cela avait bien fonctionné

milles excuses

Dernière modification par Stephane_84 (Tue 04 August 2009 13:51)

Bonjour,

alors voilà où j'en suis:


browse * from sel        'affichage des tables
browse * from selvide

dmin=1000000

for i=1 to nbselvide
    for j=1 to nbsel
        di=sqr((selvide.X-sel.X)*(selvide.X-sel.X)+(selvide.Y-sel.Y)*(selvide.Y-sel.Y))
            if di<=dmin then
               dmin=di
            end if
    next
        select*from selvide where rowid=i update selection set distance=dmin
next

en théorie ca devrait fonctionner, mais en pratique non!
j'ai des valeurs fausses dans le champ distance

mais pourquoi??

Que dire ...

je ne voit pas le fetch sur les tables ... dans ce cas pas de changement normalement on reste sur le premier enregistrement ou le dernier sélectionner ...Sans oublier l'initialisation de dmin ...

Ensuite "select*from selvide where rowid=i update selection set distance=dmin" peut être remplacé par
update selvide set distance=dmin where rowid = id ça doit marcher sur 1 seul enregistrement ...

Bonjour,
Bravo.
Je suppose que "dim=1000000" est une faute de frappe au lieu de "dmin=1000000"

Bonne continuation.