Positionnement GPS photo

y = R*LOG[TAN(Pi/4+Lat/2)]

' Positionnement sur un planisphère (en projection de Mercator)
DIM p$,wpp,hpp,wp,hp,rp,pi,fact,eq,gn,x,y,lon$,lat$,xgps%,ygps%
p$ = "C:\GRAPH\JPG\Planisphere.jpg"

wpp = 5152: hpp = 2621: ' dimensions totales de la carte, en pixels
pi = 4*ATN(1)

rp = wpp/(2*pi): ' rayon du globe, en pixels
fact = 3: ' facteur de réduction de l'image dans le picture
wp = wpp/fact: hp = hpp/fact

eq = hpp/2: ' équateur
gn = wpp/2-130: ' méridien de Greenwich

WIDTH 0,wp+16: HEIGHT 0,hp+40
PICTURE 1: WIDTH 1,wp: HEIGHT 1,hp: STRETCH_ON 1
FILE_LOAD 1,p$
2D_TARGET_IS 1: 2D_PEN_WIDTH 3
2D_LINE 0,eq,wpp,eq: ' équateur
2D_LINE gn,0,gn,hpp: ' greenwich
' Coordonnées GPS données dans les Exifs du fichier photo JPG (appareil Panasonic)
lat$ = "Lat N 31d 37m 20.45s": ' par exemple (environs Marrakech)
lon$ = "Long W  7d 59m 31.03s"

PositGPS(rp,lat$,lon$)
2D_PEN_COLOR 255,0,0: 2D_PEN_WIDTH 5
2D_LINE gn+xgps%,eq+ygps%-100,gn+xgps%,eq+ygps%+100: ' ligne V
2D_LINE gn+xgps%-100,eq+ygps%,gn+xgps%+100,eq+ygps%: ' ligne

END
' ==============================================================================
SUB PositGPS(R,lat$,lon$)
' Positionnement d'un point sur un planisphère total ou partiel en projection de Mercator
' Paramètres: R = rayon de l'image du globe terrestre, en pixels.
'            lon$, lat$: longitude et latitude donnés, de la forme par exemple
'            Lat N 31d 37m 20.45s et Long W  7d 59m 31.03s
'            (forme donnée dans les données Exifs des photos numériques).
' Retour:    coordonnées en pixels du point donné, par rapport au point zéro
'            (intersection Équateur méridien de Greenwich), dans xgps% et ygps%
DIM_LOCAL pi,k%,l$,a$,d,m,s,v
pi = 4*ATN(1)
k% = INSTR(lat$,"d"): a$ = MID$(lat$,k%-4,1): ' Latitude N ou S
l$ = MID$(lat$,k%-2,14): ' extraction des degrés, minutes, secondes
k% = INSTR(l$,"d"): d = VAL(MID$(l$,k%-2,2))
k% = INSTR(l$,"m"): m = VAL(MID$(l$,k%-2,2))
k% = INSTR(l$,"s"): s = VAL(MID$(l$,k%-5,5))
v = d + m/60 + s/3600: v = v*pi/180: ' en degrés décimaux, puis radians
ygps% = R*LOG(TAN(pi/4+v/2))
IF a$ = "N" THEN ygps% = -1*ygps%: ' -> latitude, en pixels
k% = INSTR(lon$,"d"): a$ = MID$(lon$,k%-4,1): ' Longitude E ou W
l$ = MID$(lon$,k%-2,14): ' extraction des degrés, minutes, secondes
k% = INSTR(l$,"d"): d = VAL(MID$(l$,k%-2,2))
k% = INSTR(l$,"m"): m = VAL(MID$(l$,k%-2,2))
k% = INSTR(l$,"s"): s = VAL(MID$(l$,k%-5,5))
v = d + m/60 + s/3600: v = v*pi/180: ' en degrés décimaux, puis radians
xgps% = R*SIN(v)
IF a$ = "W" THEN xgps% = -1*xgps%: ' -> longitude, en pixels
END_SUB
' ==============================================================================