Le portail francophone de la géomatique

IntersectionIl y a au moins une application (INTERSCT.zip) sur DirectionsMag
tools qui fait cela. Mais il faut bien preciser que celle-ci ne conserve
que les zones decoupees qui sont communes aux deux tables.

Jacques Paris

jacques@texte-a-enlever.paris-pc-gis.com

MapInfo Support: http://www.paris-pc-gis.com

Merci pour vos reponses.
En effet, overlay fonctionne tres bien.

Alain Laliberte
Resp. SIG
Falconbridge Ltd

Exemple de code que j'utilise pour faire quelques chose qui ressemble a
votre besoin
(calcul de colonnes au prorata de la somme des surfaces intersectees)
IL suffit de modifier un peu pour recuperer le nom de la commune
(chez-moi ce sont des disques generes pas tampons autour d'un centroide qui
jouent l'equivalent des communes)

La grille est une table ici appelee ~EMPRISE. (elle est generee a partir
d'un sous-ensemble de l'utilitaire gridmaker.mb)
J'utilise une procedure pour pouvoir mettre une jauge d'avancement.

code appellant
****************************************************************************
********************
print Calcul des ratios dans la grille...
'boucle avec barre de progression
iNombreCellule = TableInfo( ~emprise , TAB_INFO_NROWS)
redim aVar(iNombreColonne)           'aVar : tableau declare
comme alias , iNombreColonne est le nombre de colonnes choisies de la table
~EMPRISE
redim fVar(iNombreColonne)            'fVar : tableau de
flottants
redim fTotaux(iNombreColonne)          'fTotaux : tableau de
flottants
for i = 1 to iNombreColonne
aVar(i) = ~disques. + sColonneALisser(i) 'initialisation des alias,
sColonneALisser est un tableau de string -> noms de certaines colones de la
table ~DISQUES
next
O_sel1 = ~emprise + .Obj         'defini comme alias
O_sel2 = ~Disques + .Obj         'defini comme alias
i = 1
Fetch First from ~emprise
i = ~emprise.rowid  'on se place sur le 1ere enregistrement 'non
deleted'
ProgressBar Calcul des ratios dans la grille.
Calling AlimentationGrille Range iNombreCellule
if Not CommandInfo(CMD_INFO_STATUS) then
print Operation interrompue par l'utilisateur
Drop Table ~emprise
Drop Table ~Disques
exit sub
end if

Sub AlimentationGRille
'***************************************************************************
********************

Dim k as integer

Select * from ~Disques where ~disques.obj intersects (select obj from
~emprise where rowid = i) into ~BufferDisques noselect
(On selectionne les disques qui intersectent l'element i de la grille
~EMPRISE, notez la requete imbriquee
ceci accelere grandement les calculs)

OnError goto ERROR3 '(il faut traiter les erreurs car la fonction Overlap
peut echouer dans Map-Info en particulier pour des pb de precision si on
travaille en Lambert)
'(il est conseille d'utiliser l'utilitaire de definition des 'Bounds' dans
la projection de la table si des pb arrivent)
'(Gestionnaire d'outils -> Defini et applique des bornes a une carte (a
partir de la V6.0))

Fetch Rec i from ~emprise       'on se positionne sur
l'enregistrement i
Obj1 = O_sel1                'Obj1, Obj2 defini comme
'object, O_sel1, O_Sel2 comme alias

'remise a 0 des totaux (pas necessaire pour vous)
'-------------------------------
for k = 1 to iNombreColonne
ftotaux(k) = 0
next
'--------------------------------

Fetch First from ~BufferDisques   'on balaye le query resultant de la
requete en debut de procedure
While not EOT( ~BufferDisques )
Obj2 = O_sel2
esurf = Area(overlap(Obj1, Obj2), sq m ) 'surface intersectee avec les
disques

''ici il faut remplacer le traitement pour memoriser la commune ayant la
plus grande surface...

-------
---------------------------------------------------------
if esurf 0 then
DisqueSurf = Area(Obj2, sq m )
for k = 1 to iNombreColonne
fVar(k) = aVar(k)
fTotaux(k) = fTotaux(k) + (esurf / DisqueSurf) * fVar(k)
next
end if

'------
------------------------------------------------------------
APRES_ERREUR:           ' si la fonction overlap echoue on
continue le traitement quand meme apres message d'erreur a l'utilisateur
Fetch next From ~BufferDisques
Wend

' le traitement qui suit consiste a mettre a jour les colonnes de la table
~EMPRISE
-------
--------------------------------------------
'(l'utilisation de run command est rendu necessaire car je ne connais pas a
l'avance le nom des colonnes de la table ~EMPRISE (table de travail)

sCMD = Update  ~emprise  Set  + sColonneALisser(1) +  =  + fTotaux(1)
'sCMD est une string utiliser dans run command
for k = 2 to iNombreColonne
sCMD = sCMD + ,  + sColonneALisser(k) +  =  + fTotaux(k)
next
sCMD = sCMD +  where rowid =  + i
run command sCMD
'------
-------------------------------------------

i = i + 1
onError goto 0
if i > iNombreCellule then
ProgressBar = -1
else
ProgressBar = i
end if

'traitement des erreurs
'------------------------------------
ERROR3:
If Err() then
print **************ERREUR************
print Error()
print ERREUR : ~EMPRISE ligne = + str(~emprise.rowid) + ~DISQUES
ligne = + str(~disques.rowid)
print Les tables ~emprise et ~disques ne seront pas fermees, ->
reperer les lignes en erreur SVP
lerreur = T
Resume APRES_ERREUR
End If

end sub

Cordialement,

FERRATON Alain
alain.ferraton@texte-a-enlever.equipement.gouv.fr