posté par farid

JE VOUDRAI AVOIR UN DOCUMENT EXPLICATIF SUR LES SYSTEMES DE PROJECTION (UTM ET LAMBER).

posté par slim

La Transverse universelle de Mercator (en anglais Universal Transverse Mercator ou UTM) est un type de projection conforme de la surface de la Terre. L’Allemagne l’utilise sous le nom de Projection de Gauss-Kruger. Cette projection est une projection cylindrique où l’axe du cylindre croise perpendiculairement l’axe des pôles de l’ellipsoïde terrestre au centre de l’ellipsoïde.

L'UTM est également un système de référence géospatiale permettant d'identifier n'importe quel point sur notre planète .

En pratique, pour couvrir la surface de la Terre, les géomaticiens ont la découpé en 60 fuseaux de 6 degrés en séparant l’hémisphère Nord et l’hémisphère Sud. Soit au total 120 zones (60 pour le Nord et 60 pour le Sud). Ils ont développé alors un cylindre tangent à l’ellipsoïde le long d’un méridien pour obtenir une représentation plane.

Ce méridien est une ligne imaginaire tracée sur la Terre (ou la sphère céleste), joignant les pôles, et correspondant à

l'ensemble des points partageant une même longitude.

Si les latitudes peuvent être mesurées à partir de l'équateur, il n'existe pas de référence naturelle équivalente pour fixer l'origine des longitudes. Il est donc nécessaire de définir un méridien d'origine, où les points ont par définition une longitude nulle.

Actuellement, le méridien d'origine pour la plupart des systèmes géodésiques est voisin du méridien de Greenwich, qui passe par l'observatoire de Greenwich, en Angleterre. Jusqu'au début du XXe siècle, différents pays utilisèrent d'autres méridiens d'origine comme le méridien de Paris en France.

Ce fameux méridien de Paris passe par le centre de l'Observatoire de Paris. Il est situé à 2° 20' 14.025" à l'est de celui de Greenwich. Il est aussi connu sous le nom de Méridienne de France qui est sa matérialisation sur le territoire Français.

De ce fait, il est primordial d’ajouter aux anciennes coordonnées géographiques 2° 20' 14.025" afin de les convertir dans un système de projection (UTM).

posté par ibrahim74

salut ,

un document sur la topo en général dont les differents systeme de coordonnées UTM LAMBERT .....ETC

http://www.zshare.net/download/243960749d3bef/

posté par MedVal

j'essaye d'apporter un petit appui Incha Allah:

La projection Lambert est :la projection de l'ellipsoide (modele mathematique de la terre) sur un cone dont le sommet est sur l'axe des poles et ce cone est tangent au parallele origine de latitude phi0 ou parallel central.

Connaissant les coordonnees geographiques (lambda,phi) d'un lieu,on se propose de determiner les coordonnees X,Y Lambert.

Sachant la valeur phi0 (latitude du parrallel origine )

1. Calcul de R0= N0 cotg (phi0)

R0 reduit = N0 cotg(phi0) k ou k est le coefficient de reduction d'echelle.

N0= a/ (1-e*2 sin*2(phi0))**0.5 avec a et e parametres de l'ellipsoide.

2.Calcul de l0,l sachant phi (latitude)

li= ln|tg(50 grades+ phi/2)|-e/2 ln |(1+e sin(phi))/(1-e sin(phi)|

3.Ri=Ro reduit . exp [lo-lo]sin (phi0)

4.yo(i)=R0+ constante en Y-R(i)

5. dl =|l-3 grade| , gama = dl sin (phi0)

6.xi=R(i) sin (gama)

7. a= x(i) tang (gama/2)

pour V.L.U +135 m en X

8.Xi Lambert=constante en X+xi pour V.L.U +90 m en Y

Yi Lambert=Yo(i) +a

ou Yi Lambert=constante en Y +R0-R(i) cos (gama)

Pour le cas de l'Algerie vous avez les constantes suivantes:

cte Lambert Nord Algerie Lambert Sud Algerie

lo 3 grade Est Grrenwich 3 grade Est Greenwich

phi0 40 grade Nord 37 grades Nord

cteX 500 Km 300 Km

cteY 300 Km 300 Km

K 0.999625544 0.999625769

Limites d'emploi:

Lambert Nord Algerie: phi [37.5 grade,42 grade]

Lambert Sud Algerie : phi [ 34.5 grade ,39,5 grade]

Exemple ou Jeux d'essai:

Considerons l'ellipsoide de Clark 1880.A

ou e= 0.0824834 excentricite

a= 6378249.145

phi0= 40 grade Lambert Nord Algerie

k=0.999625544 coefficient de reduction d'echelle.

PHI A =40,981701 Nord, latitude du point

l A= 3,4845337 grade Est longitude du point.

N = 6385758,602 la grande normale.

R0=8785951,503

RI=8687944.738

Y0(I)= 398006.7600

Xlambert=538866.7433 m

YLambert=398093.6984 m

a l'epoque ou j'etais etudiant j'ai fait ce petit programme sur calculatrice HP41 CV que vous pouvez telecharger(calculatrice HP41CV) ,J'espere qu'il vous sera un resume explicatif:

f lbl Tran ( debut du programme de nom : Tran)

GRAD

0.824834

sto 00 ( valeur de e)

6378249.145

sto 02 (valeur de a)

PHI0=?

prompt

sto 03 ( introduire la valeur de phi0)

sto 08

XEQ01

rcl 09

sto 01 ( valeur de longitude origine)

K=?

prompt

sto 04

PHI a =? (latitude du point)

prompt

sto 05

LAMDA a =? (longitude du point)

prompt

sto 06

rcl 03

sin

x*2

rcl 00

x*2

*

CHS

1

+

v- (racine carree)

rcl 02

x>

posté par MedVal

je vois que j'ai depasse la limite de texte permise donc la suite du programe a ete troncuee ,desole. Mais on va remedier a ca ,voila la suite du programme:

v- ( racine carree)

rcl 02

x>< y ( permutation des registres x et y) [ figurent dans le precedent message]

/

rcl 03

tan

1/x

*

rcl 04

*

sto 07 ( R0 reduit)

R0 r+

ARCL.X

AVIEW

R/S

rcl 05

sto 08

XEQ 01

rcl 09

sto 11

rcl o1

x>

posté par MedVal

desole j'ai redige la totalite du programme mais ce que je vois des toncatures je le complete sur la meme ligne alors

sto11 rcl01 x><y (du message precedent) - rcl03 sin * e*(expon-) rcl07 * sto12 RI= ARCL.X AVIEW R/S CHS 300000 (300095 pou v.l.u) + rcl07 + sto13 (Y0(I)) YI= ARCL.X AVIEW R/S rcl 06 3 - ABS (valeur absolue) rcl 03 sin * sto14(gama) sin rcl12 * sto15 500000 ( 500135 pour v.l.u) + X= ARCL.X AVIEW R/S rcl14 2 / tan rcl 15 * rcl 13 + Y= ARCL.X AVIEW RTN ( ligne 109 du prgm-) f lbl 01 (sous prgm- 01) rcl08 2 / 50 + tan ABS LN sto99 rcl08 sin rcl00 * 1 + enter rcl08 sin rcl00 * 1 + enter rcl08 sin rcl00 * CHS 1 + / ABS LN rcl00 2 / * chs enter rcl99 + sto09 RTN ( fin du prgm a la ligne 146).