frame_1frame_2frame_3frame_4frame_5frame_5

Retour sommaire de la programmation

Programme qui résoud des problèmes

d'optiques et de photogrammetrie





 

Ce programme converti les

Latitudes et longitudes

en

coordonées Lambert


 

Ce programme converti les

coordonées Lambert

en

Latitudes et longitudes



 

Ce programme est un

petit bonheur du jour !

10 MODE6:SETF4:CLS:PRINT"(DV5.CAS10/04/99) 1Li//Qua 2Cer 3CNi 4ko 5kL 6m 7ETer 8Exct";:T$=INPUT$(1):T=VAL(T$):CLS
20 ONT GOSUB100,3000,4000,5000,6000,7000,8000,9000:GOTO10
100 S=0:AG=0:AD=0:B=0:H=0:AP=0:PQ=0
140 CLS:PRINT "LimPQ//BaseAB,AngGetD";
150 INPUT " Surface ";S
155 IFS=0THEN150
160 INPUT "AngGauche ";AG
165 IFAG=0THEN160
166 IFAG=100THENAG=100-10^-8
170 INPUT "AngDroite ";AD
175 IFAD=0THEN170
176 IFAD=100THENAD=100-10^-8
180 INPUT "BaseAB ";B
185 IFB=0THEN180
186 ON ERROR GOTO 210
190 PQ=SQR(B^2-2*S*(1/TANAG+1/TANAD)):H=(2*S)/(B+PQ):AP=H/SINAG:BQ=H/SINAD
200 SETF3:PRINT "LimPQ=";PQ:PRINT "HauteurH=";H:PRINT "CoteAP(Gauche)=";AP:PRINT "CoteBQ(Droite)=";BQ
205 SETN:GOTO100
210 PRINT "S>Striangle,DonneesFausses":RESUME205:SETN
215 RETURN
3000 CLS:PRINT"Probl sur Cercle";:PRINT" Quand inconnues Taper EXE":SETF4
3010 RA=0:AO=0:TT=0:TG=0:DV=0:CF=0:FL=0:ST=0:SG=0:AT=0:STT=0
3020 INPUT"Rayon ";RA:INPUT"Angle au centre ";AO:AO=AO/2:IFAO=100THENAO=100-10^-8
3021 IFAO=100-10^-8THENPRINT"Tangentes paralleles Pas de contre fleche"
3022 IFRA*AO<>0THEN3290
3030 INPUT "Angle tgtes(au sommet) ";AT:AT=AT/2
3033 IFAO<>0THENAT=100-AO
3034 IFAT<>0THENAO=100-AT
3035 IFRA*AT<>0THEN3290
3040 INPUT"Biss.(CtreFl) ";CF
3050 IFCF<>0ANDAO<>0THEN RA=(CF*COSAO)/(1-COSAO)
3051 IFCF<>0AND RA<>0THENAO=ACS(RA/(RA+CF))
3060 IFCF<>0AND AO<>0THEN3290
3070 INPUT"Corde ";TT:TT=TT/2
3080 IFRA=0AND AO<>0THENRA=TT/SINAO
3090 IFAO=0AND RA<>0THENAO=ASN(TT/RA)
3100 IFRA*AO<>0THEN3290
3110 INPUT"Fleche ";FL:IFTT+FL=0THEN3190
3115 IFTT<>0ANDFL<>0THENRA=(FL^2+TT^2/4)/2*FL
3120 IFFL<>0ANDAO<>0THENRA=FL/(1-COSAO)
3130 IFFL<>0ANDAO<>0THEN3290
3140 IFTT*FL=0THEN3170
3150 RA=(FL^2+TT^2)/2/FL:AO=ASN(TT/RA):IFRA*AO<>0THEN3290
3160 IFFL*AO<>0THENRA=FL/(1-COSAO)
3170 IFFL*RA<>0THENAO=ACS((RA-FL)/RA)
3180 IFRA*AO<>0THEN3290
3190 INPUT"Tangente ";TG:IFTT*TG<>0THENAO=100-ASN(TT/TG)
3200 IFTG*RA<>0THENAO=ATN(TG/RA)
3210 IFTG*AO<>0THENRA=TG/TANAO
3220 IFRA*AO<>0THEN3290
3230 INPUT"Surf.Secteur ";ST:IFST*AO<>0THENRA=SQR((200*ST)/(PI*AO))
3231 IFST*AO<>0THEN3290
3240 IFST<>0ANDRA<>0THENAO=(400*ST)/(2*PI*RA^2)
3250 IFST<>0ANDRA<>0THEN3290
3260 INPUT"Developpement ";DV:IFST*DV<>0THENRA=2*ST/DV
3261 IFST<>0ANDRA<>0THENAO=(400*ST)/(2*PI*RA^2)
3262 IFST*DV<>0THEN3290
3263 IFRA<>0ANDDV<>0THENAO=100*(DV/RA)/PI:GOTO3290
3264 IFDV*AO<>0THENRA=(DV*100)/(AO*PI):GOTO3290
3270 INPUT"Surf.Segment ";SG:IFSG*AO<>0THENRA= SQR((200*SG)/((PI*AO)-(100*SIN(2*AO))))
3271 IFAO*SG<>0THENGOTO3290
3273 IFST*SG<>0THENSTT=ST-SG:PRINT"Surf Triangle=";STT
3274 IFSTT<>0THENGOTO3000
3275 INPUT "Surf.Triang";STT:IF STT*AO<>0THENRA=SQR(STT/(COSAO*SINAO):GOTO 3290
3280 IFSG*AO=0THENPRINT"Donnees insuffisantes":GOTO 3000:SETN
3290 TG=RA*TANAO:TT=2*RA*SINAO:FL=RA*(1-COSAO): CF=TG*SINAO-FL:DV=PI*RA*AO/100:ST=RA^2*PI*AO/200:AO=2*AO:AT=200-AO: SG=ST-RA^2*SINAO/2:STT=ST-SG
3300 PRINT"Rayon=";RA:PRINT"Angle au centre=";AO: PRINT"Ang.Tgtes(au sommet)=";AT:PRINT"Tangente=";TG:PRINT"Corde=";TT:PRINT"Fleche=";FL
3310 PRINT"Developpement=";DV:PRINT"Biss.(CtreFl)=";CF: PRINT"Surf Secteur=";ST:PRINT"Surf Segment=";SG:PRINT"Surf Triangle=";STT
3320 GOTO 3000:SETN
3330 RETURN
4000 SETF2:CLS:PRINT"Courbe niveau";:INPUT" Z 1er Pt inf(m)";ZI
4010 ONERRORGOTO4100
4040 INPUT"Z 2e Pt sup(m)";ZS
4045 INPUT"Dist entre ces 2 pts(cm)";D
4050 INPUT"Z Courbe(m)";ZC
4070 X=(D/(ZS-ZI))*(ZC-ZI)
4080 PRINT"Dist courbe du pt inf="X;"cm"
4090 GOTO4000
4100 PRINT"Donnees fausses":RESUME4000
4110 RETURN
5000 SETF0:CLS:PRINT"CoefNiv0:ko=-Zmoy(m)*10^3/6370km";:INPUT" Zmoy(m)";ZM
5010 ko=-(ZM/6370)*10^3
5020 PRINT"ko= ";ko;"PPM"
5030 CLS:GOTO5000
5040 RETURN
6000 CLS:SETF1:PRINT"Coef Lambert: kL ";:INPUT"Xmoy(km)";XM:INPUT"Ymoy(km)";YM:INPUT"Xinf(km)";Xinf:INPUT"Yinf(km)";Yinf
6010 INPUT"kXinfYinf(cm)";A
6020 INPUT"kXsupYinf(cm)";B
6030 INPUT"kXinfYsup(cm)";C
6040 INPUT"kXsupYsup(cm)";D
6050 E=XM-Xinf:F=YM-Yinf
6070 G=A+.01*(B-A)*E:H=C+.01*(D-C)*E:I=G+.1*(H-G)*F
6080 PRINT"kXmoyYinf=";G;"cm"
6090 PRINT"kXmoyYsup=";H;"cm"
6100 SETF0:PRINT"kL=";I*10;"PPM
6105 GOTO6000
6110 RETURN
7000 SETF6:CLS:PRINT"Module m=1+k(ppm)/10^6 ";:INPUT"k(PPM)";k
7010 m=1+k*10^-6:PRINT"m=";m
7020 GOTO7000
7030 RETURN
8000 SETF2:CLS:PRINT"Entrees en terre (penteP en ABS)";
8010 INPUT"ZA(proj)";ZA
8020 IFZA=0THEN8010
8030 INPUT"ZB(tn)";ZB
8040 IFZB=0THEN8030
8050 INPUT"P1proj(%)";P1:IFP1=0THENP1=10^-9
8070 INPUT"P2tn(%)";P2:IFP2=0THENP2=10^-10
8080 P=P1/100-P2/100:IFP=0THENP=10^-10
8090 X1=(ZA-ZB)/P:X2=(ZA-ZB)/(P1/100+P2/100)
8100 X1=ABSX1:PRINT"x(cm)P meme sens =";X1;
8105 X2=ABSX2:PRINT" x(cm)P sens dif =";X2
8110 CLS:GOTO8010
9000 SETF4:CLS:PRINT"S=st R=ptpRoche P=Ptloin RP=distC=ASN(SR*SIN(HSP-HSR)/RP)"
9010 INPUT"SR";r:INPUT"HSR";h:INPUT"RP";P:INPUT"HSP";H
9015 C=ASN(r*SIN(H-h)/P):H=H+C: PRINT"HRS=";h+200:PRINT"C(gon)=";C:PRINT"HRP=";H:GOTO9010
65000 CLS:SETN:END:REM FIN DU TRANSFERT

Retour sommaire de la programmation