Le portail francophone de la géomatique

Bonjour a tous

Suite a quelques problemes de connexion pour recuperer le document je
vous invite a essayer l'adresse suivante :

ftp://lct.teledetection.fr/pub/sigmap/dgal/

Attention : Jacques Paris vient de m'envoyer une version amelioree du
document, je prend le temps de le lire avant de le mettre a disposition.

Cordialement

De: Eric Mauviere
J'ai lu avec un tres grand interet le document
elabore par Sebastien Roddier et Jacques Paris.
J'ai apprecie leur approche progressive et l'exhaustivite du document,
nourri de nombreux exemples.

Je voudrai apporter ici 3 cas concrets supplementaires,
qui correspondent a des problemes frequemment rencontres dans ma vie de statisticien.

1) lissage geographique d'un taux
Prenons l'exemple du taux des jeunes de moins de 20 ans dans la population totale,
etudie par commune.
Le lissage evoque ici, bien connu desormais, consiste a calculer le taux sur une zone
un peu plus large que la commune a laquelle il sera affecte, en prenant en compte
la premiere couronne des communes adjacentes. Il reprend le principe de la moyenne mobile,
utilise pour lisser des courbes.

Si le calcul est possible en SQL simple, il suppose de creer une copie du fond de carte
de travail, avec les deux variables constitutives du taux, puis de faire une jointure cartesienne
du fonds source et de sa copie, avec une condition intersects .

Mapinfo offre une alternative plus elegante encore,
qui permet de faire le calcul en deux lignes, sans avoir a creer de copie,
grace a la fonction Mettre a jour colonne :

add column data (_varc) from data set to sum((10^9)*var1+var2) where intersects
add column data (_varliss) from data set to 100*round(_varc/10^9,1)/(_varc -10^9*round(_varc/10^9,1))

Commentaire : data est le fonds source (ex : region PACA par commune), var1 le premier element du taux (ex : nb des moins de 20 ans),
var2 le 2eme element du taux (ex population totale de la commune).
add column permet de faire une jointure d'une table sur elle-meme, ce que le SQL de Mapinfo n'autorise pas.
on a dans la premiere ligne groupe le calcul de la somme mobile des elements du taux, en faisant l'hypothese
que ces elements etaient inferieurs a 10^9.

2) calcul d'une somme cumulee

Soit data la table de travail, var le champ numerique dont on veut produire la somme cumulee.
On suppose que data est deja triee dans l'ordre desire pour le cumul.

Alter Table data ( add _rowid Integer, _test smallint )
Update data Set _rowid = rowid, _test=1
Commit Table data
Commit Table data As c:tempdatac.tab
Open Table c:tempdatac.tab
Select data._rowid,data.var,sum(datac.var) csom_var from data, datac where data._test=datac._test and data._rowid>=datac._rowid group by 1 into csum noselect

Commentaire : on numerote les lignes de la table source et l'on cree un champ test=1
avant de faire une copie de la table.
La jointure cartesienne de la table source et de sa copie est toleree par Mapinfo des lors qu'une condition de jointure geographique
ou qu'une egalite entre deux champs lui est fournie (role de la variable test ).

3) determination automatique de la commune la plus peuplee d'un canton donne
l'exemple commune/canton pris ici s'etend a tout couple de zonages emboites.

Select cant,max(val(str(psdc99) & com)) test from FOND group by 1 into temp1 noselect
Select cant,right(str(test),5) compole from temp1 into tab_poles

Commentaire : cant est le code canton, com le code commune, psdc99 la variable population du fonds de carte FOND .
La table resultante tab_poles indique pour chaque canton present dans FOND le code de sa commune la plus peuplee.

Eric Mauviere
http://www.geoclip.net