Le portail francophone de la géomatique

Salut et bienvenue,

Ce type de traitement peut être fait avec la FMB (Fenêtre MapBasic) Options --> Voir Fenêtre MapBasic
Cette fenêtre couplée avec un petit Log appelé Scriptgen peut t'aider.
N'hésite pas à faire des recherches sur le forum avec comme mot clef "Scritpgen" "générateur de scripts".
A+
Joël

Merci Joël pour ta promptitude.
Comme je l'écrivais dans mon post, je ne suis pas un chantre de la
programmation. toutefois je vais aller fouiller le "Scriptgen" sur le forum.
JMR

Salut à vous deux,

Je ne comprend pas tout
MapBasic c'est de la programmation pure et dure (production d'outil compilé en .mbx)
La FMB c'est utilisé quelques requêtes MapBasic au travers d'une fenêtre de commandes
On va dire --> la programmation en dose homéopathique
Scriptgen est juste là pour dupliquer les commandes pour arriver à une automatisation.
Je te mets en pièce jointe ce qui est faisable avec la FMB
Voir "Calculer un angle par rapport à une droite Horizontale"
PS:Pour en revenir à scriptgen, ce petit logiciel permets de dupliquer des commandes et donc être utile pour tout logiciel ayant une fenêtre de commande --> Autocad, Arcgis, Mapinfo, Fenêtre de commande Dos etc...........
et aussi un gros avantage, il est facilement portabilisable donc sur une clef USB.
A+
Joël

Dernière modification par Spacejo (Tue 05 February 2013 14:55)

Merci Joël de rebondir sur les tenants et aboutissants de MapBasic et
Scriptgen, mais vous êtes tous les deux beaucoup plus calés que moi sur le
sujet. Alors je lis consciencieusement ... mais je ne sais toujours pas
comment résoudre mon "problème" de lignes et rosaces.
Concrètement, il me faudrait extraire l'azimut de chaque lignes... avec la
contrainte que ce ne sont pas toujours des droites (puisque ce sont des
failles / linéaments et que la Nature n'est pas que pure Géométrie !).

Je sais que ArcGis a une routine qui fait ce genre d'opération... mais
peut-être pas MapInfo / Discover.
Discover fait des rosaces... mais il lui faut une table avec des azimuts.

Merci pour la réponse.
JM Regnoult

Salut,

Pour ce qui concerne tes failles (pas de problème car il te faudra extraire les points de début et de fin).
Il existe des outils pour ça ---> En voilà 2 --> http://www.directionsmag.com/files/view … nts/139018
---> http://www.routeware.dk/toolbox.php ---> Voir "Add Nodes"
Maintenant si tu suis le pdf (l'angle est mesuré par rapport à une horizontale) --> en changeant la formule tu auras l'angle par rapport à une verticale.
1-Extraire les points de début et de fin de tes failles
2-Crée le code pour une routine (calcul de l'angle avec la verticale grâce aux points extraits--> pour une faille)
3-Automatisation de cette routine grâce à une liste de tes points (début et fin) pour chaque faille (avec Scriptgen)
4-Traçage de tes rosaces avec les éléments contenus dans ta table (champs longueur et azimuth)
Voilà le plan d'attaque
A+
Joël

Dernière modification par Spacejo (Wed 06 February 2013 22:06)

Salut,

Le pdf est celui que j'ai mis en pièce jointe ci dessus.
A+
Joël

Bonjour Joël,

Ton aide a été efficace car j'ai réussi à extraire mes longueurs (facile), puis les coordonnées des extrémités de mes lignes (plus compliqué car il s'agissait de polylignes que j'ai dues "éclater" et j'ai ratatouillé pour le Code). Enfin, je me sens un peu plus familier avec la "FMB" (qui jusqu'à lors restait pour moi une fenêtre "magique")

Et maintenant je cale sur les angles...

Cordialement.
JMR

Merci pour le suivi, Joël,

Grâce à tes premières explications, j'ai réussi à extraire mes points de début et fin de segments.

Maintenant je bloque sur l'extraction automatique de l'azimuth des segments.
Explication:
- voici le code que j'ai utilisé dans la FMB pour calculer l'azimuth :

Select * from MaTable where (Y2-y1) > 0 and (X2-x1) < 0  Into azimut
'==> NW => angles de 90 à 180°
Update MaTable Set azimut = 180 - atn(abs(Y2-y1)/abs(X2-x1))*57.2958
Select * from MaTable where (Y2-y1) < 0 and (X2-x1) < 0  Into azimut2
'==> SW => angles de 180 à 270°
Update MaTable Set azimut2 = 180 + atn(abs(Y2-y1)/abs(X2-x1))* 57.2958
Select * from MaTable where (Y2-y1) < 0 and (X2-x1) > 0  Into azimut3
'==> SE => angles de 270 à 360°
Update MaTable Set Azimut3 = 360 - atn(abs(Y2-y1)/abs(X2-x1))* 57.2958

Explication sur les paramètres:
- (x1,y1) coord. début segment et  (X2,Y2) coord. fin segment;
- azimut : écriture anglaise sans "h" = direction par rapport au Nord géographique du segment.  "Azimut" est le nom d'une colonne dans MaTable. J'ai créé une colonne pour chaque quadrangle (azimut, azimut2, azimut3).

Préalables:
- MI 10.5.2 sous Win7.
- des essais ont été faits avec des cartes en Long-lat DD WGS84 et en UTM-WGS84.
- tous les segments sont dans le même sens.

Remarque :
- pourquoi MapInfo emplit-il toutes les lignes des colonnes "azimut", alors que je lui impose des conditions (delta y) - (delta x) etc . ?  Mon code est-il mal rédigé !?

Résultats:
- je trouve des azimuths corrects, sauf pour certains segments aberrants (environ 10% ?!).  Par contre lorsque j'utilise l'outil "Calculer la direction" de MapCad dans MI 10.5.2, je trouve que tous les azimuths sont corrects en mesurant  les segments un par un (tu parles d'un boulot!).

Interprétation possible (selon moi!): l'origine des erreurs provient de ce que je ne maîtrise pas le point d'origine des segments (puisque le calcul a été fait automatiquement par MI) et que Mapinfo "cafouille" à ce moment.

Solution : As-tu une suggestion ? (moi je suis sec).

Au plaisir de te relire.
JM Regnoult

Bonsoir Joël,

Pas de repos pour les braves :-(

Ton Code ne change pas tellement les choses par rapport à ce que j'ai pu déjà faire grâce à tes conseils. Pour être plus clair, tu trouveras ci-joint une copie d'écran et un fichier texte du tableau de mes segments.
La colonne la plus à droite correspond à des mesures faites "à la main" (outil MapCad) ... donc justes, par rapport au Nord géographique !
Le "malheur" est qu'il n'y a pas de règle absolue qui se dégage du résultat dans le tableau (ci-joint).
Comme je le disais dans mon précédent post, ça doit venir du fait que les points de début des segments ne sont pas tous orientés de la même façon (bien que j'ai pris soin de vérifier que tous les segments sont dans le même sens).

Merci pour ton aide active ... ça fait du bien d'être  épaulé.
JMR

PS : désolé , mais je n'ai pas réussi à attacher mes fichiers en ligne ... comme recommandé sur le site ! Je fais glisser mon icône sur le message ou son bord, mais rien ne se passe ...  Comment fait-on ? Merci


Merci pour m'avoir expliqué la "jointure". J'espère que les 2 fichiers (1 petite image et 1 Excel) sont là !

Dernière modification par data_public (Mon 18 February 2013 10:19)

Salut

J'avais une question --> pourquoi 2 nouveaux champs azimut?
La différence entre mon code et le tien c'est que mes "update" ce font uniquement sur mes sélection succéssives.
Pour ton cas à chaque fois tu "update" toute ta table donc tes conditions et donc tes requêtes ne servent pas.

Select * from MaTable where (Y2-y1) > 0 and (X2-x1) < 0  Into Selection1
Update Selection1 Set azimut = 180 - atn(abs(Y2-y1)/abs(X2-x1))*57.2958
Select * from MaTable where (Y2-y1) < 0 and (X2-x1) < 0  Into Selection2
Update Selection2 Set azimut2 = 180 + atn(abs(Y2-y1)/abs(X2-x1))* 57.2958
Select * from MaTable where (Y2-y1) < 0 and (X2-x1) > 0  Into Selection3
Update Selection3 Set azimut3 = 360 - atn(abs(Y2-y1)/abs(X2-x1))* 57.2958

A essayer avec tes 3 champs "azimut"
PS:Je regarde tes fichiers pour mieux comprendre ta problèmatique
Peux tu m'envoyer le fichier .xlsx en .dbf --> plus facile pour moi je suis resté en excel 2003
A+
Joël

Dernière modification par Spacejo (Mon 18 February 2013 10:44)

Salut,

Après une petite analyse, je me suis aperçu que tu as des entités qui ont été digitalisées du nord vers le sud (donc à l'envers).
J'avais une question --> pourquoi avoir besoin de 4 colonnes "azimut".
A+
Joël