Le portail francophone de la géomatique

include "C:\Program Files\MapInfo\MapBasic\MAPBASIC.DEF"

declare sub main     

sub main

set coordsys Earth Projection 8, 9999, 4, 97.295, -263.237, 310.866, 1.599931, -0.838952, -3.14118, 13.325864, 0, "m", 165, 0, 0.9996, 500000, 10000000 Bounds (-7746230.6469, 1712.61611073) (8746230.6469, 19998287.3839) 

'les variables
dim xdol as float    'coordonnées pt A
dim ydol as float
dim xpt as float            'coordonnées pt B
dim ypt as float
dim xclst as float        'coordonnées pt + près
dim yclst as float
dim distmin as float
dim distcalc as float

dim nb_ptdoline as integer
dim nb_ptligne as integer
dim n as integer
dim j as integer

dim dol as object     'objet point de la table A
dim opt as object     'objet point de la table B

'ouverture des tables A, B et création de la table de stockage des lignes droites
open table "C:\Etude\A09-52 - Doline\SELECTION DE DOLINE.tab" as ptdoline
open table "C:\Etude\A09-52 - Doline\LIMITE_PLATEAU_ptp.tab" as ptligne
Create Table "CHEMIN_COURT" (OBJ_ID char(10), LONGUEUR float) file "C:\Etude\CHEMIN_COURT.TAB" TYPE NATIVE Charset "WindowsLatin1"

'fenêtre carte pour les chemin les plus courts
Create Map For CHEMIN_COURT CoordSys Earth Projection 8, 9999, 4, 97.295, -263.237, 310.866, 1.599931, -0.838952, -3.14118, 13.325864, 0, "m", 165, 0, 0.9996, 500000, 10000000 Bounds (-7746230.6469, 1712.61611073) (8746230.6469, 19998287.3839)
Map from CHEMIN_COURT
Set Map Layer 1 Editable On

nb_ptdoline = tableinfo(ptdoline, TAB_INFO_NROWS)
nb_ptligne = tableinfo(ptligne, TAB_INFO_NROWS)

print "nombre de dolines à traiter: " & nb_ptdoline
print "nombre de points a comparer: " & nb_ptligne

distmin = 1000000
n = 1
j = 1

for n = 1 to nb_ptdoline
  print "doline traitée : " & n & " sur " & nb_ptdoline
    select * from ptdoline where rowid = n
    dol = selection.obj
    xdol = objectgeography(dol, OBJ_GEO_POINTX)
    print xdol
  ydol = objectgeography(dol, OBJ_GEO_POINTY)
  print ydol
  
  for j = 1 to nb_ptligne
      select * from ptligne where rowid = j
      opt = selection.obj
      xpt = objectgeography(opt, OBJ_GEO_POINTX)        
      print xpt
      ypt = objectgeography(opt, OBJ_GEO_POINTY)
    print ypt
    
    'calcul de la distance 
    distcalc = sqr(abs(xdol - xpt) * abs(xdol - xpt) + abs(ydol - ypt) * abs(ydol - ypt))
    print distcalc
    
    if distcalc < distmin 
        then
        distmin = distcalc
        xclst = xpt
      yclst = ypt
    end if 
      
      j = j + 1
  next
  
  note distmin
  create line into window frontwindow() (xdol, ydol) (xclst, yclst)
  commit table CHEMIN_COURT 
  
    n = n + 1  
next

end sub