diff options
| author | Jovina Dsouza | 2014-07-25 00:00:07 +0530 |
|---|---|---|
| committer | Jovina Dsouza | 2014-07-25 00:00:07 +0530 |
| commit | fd5a671b82455b88fd313d8d0bee2793ab27739a (patch) | |
| tree | e35148d9f2c5d1ee88f62f903b2ca46292b6f568 /Surveying_Volume_3 | |
| parent | c8733e4b6b4bffcddf7eb45ff1c72ccc837aa3af (diff) | |
| download | Python-Textbook-Companions-fd5a671b82455b88fd313d8d0bee2793ab27739a.tar.gz Python-Textbook-Companions-fd5a671b82455b88fd313d8d0bee2793ab27739a.tar.bz2 Python-Textbook-Companions-fd5a671b82455b88fd313d8d0bee2793ab27739a.zip | |
adding book
Diffstat (limited to 'Surveying_Volume_3')
| -rwxr-xr-x | Surveying_Volume_3/Chapter1.ipynb | 83 | ||||
| -rwxr-xr-x | Surveying_Volume_3/Chapter1_1.ipynb | 83 | ||||
| -rwxr-xr-x | Surveying_Volume_3/Chapter3.ipynb | 104 | ||||
| -rwxr-xr-x | Surveying_Volume_3/Chapter4.ipynb | 970 | ||||
| -rwxr-xr-x | Surveying_Volume_3/Chapter5.ipynb | 398 | ||||
| -rwxr-xr-x | Surveying_Volume_3/Chapter6.ipynb | 629 | ||||
| -rwxr-xr-x | Surveying_Volume_3/README.txt | 10 | ||||
| -rwxr-xr-x | Surveying_Volume_3/screenshots/chapter1.png | bin | 0 -> 112592 bytes | |||
| -rwxr-xr-x | Surveying_Volume_3/screenshots/chapter4.png | bin | 0 -> 122796 bytes | |||
| -rwxr-xr-x | Surveying_Volume_3/screenshots/chapter5.png | bin | 0 -> 213035 bytes |
10 files changed, 2277 insertions, 0 deletions
diff --git a/Surveying_Volume_3/Chapter1.ipynb b/Surveying_Volume_3/Chapter1.ipynb new file mode 100755 index 00000000..015050ae --- /dev/null +++ b/Surveying_Volume_3/Chapter1.ipynb @@ -0,0 +1,83 @@ +{ + "metadata": { + "name": "S3-C1" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Electronic Distance Measurement" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.1,Page 22" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nT=273+25.0#temperature \np=752.0;#pressure mm Hg\nNo=294.0e-6;\nns=1.000284;\nD1=1438.254;#recorded distance in m\nh=263.42-243.25;#height difference in m\nR=6370e3;#radius of earth in m\n\n#calculation\nn=1+No*(273/T)*(p/760);\nD=D1*ns/n;\ncg=-h**2/2/D;\nHm=263.42/2+243.25/2;\nD=D+cg;\nch=-D*Hm/R\nEL=D+ch;\n\n#result\nprint \"equivalent length in m\",round(EL,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "equivalent length in m 1438.081\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.2,Page 22" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nT=273+18.8;#temperature\np=713;\ne=3;\nc=299792.5e3;#speed of light in km/s\nf=11e6;#frequency in Hz\n\n#calculation\nn=1+(103.49/T*(p-e)+86.26/T*(1+5748.0/T)*e)/1e6;\nV=c/n;\nl=V/f;\n\n#result\nprint \"wavelenght of light in m\",round(l,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "27.247 wavelenght of light in m\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.3,Page 23" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nT=273.0+30.0;#temperature\nl=0.85;#wavelength in m\np=752.4;\nc=299792.5e3;#speed of light in km/s\nf=24e6;#frequency in Hz\n\n#calculation\nno=1+(287.604+4.8864/l**2+0.068/l**4)/1e6;\nns=1+(no-1)*273/T*p/760;\nV=c/ns;\nl=V/f;\n\n#result\nprint \"wavelength of light in m\",round(l,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "wavelenght of light in m 12.488\n" + } + ], + "prompt_number": 1 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Surveying_Volume_3/Chapter1_1.ipynb b/Surveying_Volume_3/Chapter1_1.ipynb new file mode 100755 index 00000000..015050ae --- /dev/null +++ b/Surveying_Volume_3/Chapter1_1.ipynb @@ -0,0 +1,83 @@ +{ + "metadata": { + "name": "S3-C1" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Electronic Distance Measurement" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.1,Page 22" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nT=273+25.0#temperature \np=752.0;#pressure mm Hg\nNo=294.0e-6;\nns=1.000284;\nD1=1438.254;#recorded distance in m\nh=263.42-243.25;#height difference in m\nR=6370e3;#radius of earth in m\n\n#calculation\nn=1+No*(273/T)*(p/760);\nD=D1*ns/n;\ncg=-h**2/2/D;\nHm=263.42/2+243.25/2;\nD=D+cg;\nch=-D*Hm/R\nEL=D+ch;\n\n#result\nprint \"equivalent length in m\",round(EL,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "equivalent length in m 1438.081\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.2,Page 22" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nT=273+18.8;#temperature\np=713;\ne=3;\nc=299792.5e3;#speed of light in km/s\nf=11e6;#frequency in Hz\n\n#calculation\nn=1+(103.49/T*(p-e)+86.26/T*(1+5748.0/T)*e)/1e6;\nV=c/n;\nl=V/f;\n\n#result\nprint \"wavelenght of light in m\",round(l,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "27.247 wavelenght of light in m\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.3,Page 23" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nT=273.0+30.0;#temperature\nl=0.85;#wavelength in m\np=752.4;\nc=299792.5e3;#speed of light in km/s\nf=24e6;#frequency in Hz\n\n#calculation\nno=1+(287.604+4.8864/l**2+0.068/l**4)/1e6;\nns=1+(no-1)*273/T*p/760;\nV=c/ns;\nl=V/f;\n\n#result\nprint \"wavelength of light in m\",round(l,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "wavelenght of light in m 12.488\n" + } + ], + "prompt_number": 1 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Surveying_Volume_3/Chapter3.ipynb b/Surveying_Volume_3/Chapter3.ipynb new file mode 100755 index 00000000..4c4c37ac --- /dev/null +++ b/Surveying_Volume_3/Chapter3.ipynb @@ -0,0 +1,104 @@ +{ + "metadata": { + "name": "S3-C3" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Trilateration" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.1,Page 47" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos\nAB=25145.32;\nR=6370.0e3;#radius of earth\nha=325.14;\n\n#calculation\ntheta=AB*cos(3+9.0/60+40.0/3600)/R;\nAB_dash=AB/sin(pi/2+theta/2)*sin(pi/2-theta/2-(3+9.0/60+40.0/3600)*pi/180);\nCD=AB_dash-AB_dash*ha/R;\nS=CD+CD**3/24.0/R**2;\n\n#result\nprint \"sea level length in m\",round(S,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "sea level length in m 25108.53\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": "Example 3.2,Page 50" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,asin\nAB=31325.14\nR=6370.0e3;#radius of earth\nha=1582.15;\nh=2669.17\n\n#calculation\nAB_dash=AB-h**2/2/AB;\ntheta=2*asin(AB_dash/2/R);\nAB_dash=AB-(h*sin(theta/2)+h**2/AB/2)\nCD=AB_dash-AB_dash*ha/R;\nS=CD+CD**3/24.0/R**2;\n\n#result\nprint \"sea level length in m\",round(S,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "sea level length in m 31197.163\n" + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.3,Page 53" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin, acos\nAB=1525.456;\nBC=2176.945;\nCD=1697.435;\nAD=2401.435;\nAC=3073.845;\nBD=2483.115;\n\n#calculation\n#alpha=A and beta=B and those are angles\nA1=acos((CD**2+AC**2-AD**2)/(2*CD*AC));\nA1=A1*180/pi;\nA2=acos((AD**2+BD**2-AB**2)/(2*AD*BD));\nA2=A2*180/pi;\nA3=acos((AB**2+AC**2-BC**2)/(2*AB*AC));\nA3=A3*180/pi;\nA4=acos((BC**2+BD**2-CD**2)/(2*BC*BD));\nA4=A4*180/pi;\nB1=acos((CD**2+BD**2-BC**2)/(2*CD*BD));\nB1=B1*180/pi;\nB2=acos((AD**2+AC**2-CD**2)/(2*AD*AC));\nB2=B2*180/pi;\nB3=acos((AB**2+BD**2-AD**2)/(2*AB*BD));\nB3=B3*180/pi;\nB4=acos((AC**2+BC**2-AB**2)/(2*AC*BC));\nB4=B4*180/pi;\ne1=360-A1-A2-A3-A4-B1-B2-B3-B4;#error\ne2=A1+B1-A3-B3;#error\ne3=A2+B2-A4-B4;#error\n#angle update\nA1=A1+e1/8-e2/4;\nA3=A3+e1/8+e2/4;\nB1=B1+e1/8-e2/4;\nB3=B3+e1/8+e2/4;\nA2=A2+e1/8-e3/4;\nB2=B2+e1/8-e3/4;\nA4=A4+e1/8+e3/4;\nB4=B4+e1/8+e3/4;\n#updating sides\nAD=1525.456*sin(B3*pi/180)/sin(A2*pi/180);\nBD=1525.456*sin(A3*pi/180+B3*pi/180)/sin(A2*pi/180);\nAC=1525.456*sin(A4*pi/180+B3*pi/180)/sin(B4*pi/180);\nBC=1525.456*sin(A3*pi/180)/sin(B4*pi/180);\nCD=BC*sin(A4*pi/180)/sin(B1*pi/180)\n\n#result\n\nprint \"equation for B2 is wrong\"\nprint \"corrected length of AD in m \",round(AD,3)\nprint \"corrected length of AC in m \",round(AC,3)\nprint \"corrected length of BD in m \",round(BD,3)\nprint \"corrected length of BC in m \",round(BC,3)\nprint \"corrected length of CD in m \",round(CD,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "27.5965629411\nequation for B2 is wrong\ncorrected length of AD in m 2401.251\ncorrected length of AC in m 3073.446\ncorrected length of BD in m 2413.581\ncorrected length of BC in m 2176.758\ncorrected length of CD in m 1697.416\n" + } + ], + "prompt_number": 25 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.4,Page 59" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,log,sin\n#angles found by cosine law\n#alpha=A and beta=B and those are angles\nA1=45.801596;\nA2=40.605250;\nA3=50.143258;\nA4=43.077646;\nB1=48.779868;\nB2=44.141587;\nB3=49.733152;\nB4=37.737035;\n\n#calculation\ne1=360-A1-A2-A3-A4-B1-B2-B3-B4;#error\n#angle update\nA1=A1+e1/8;\nA3=A3+e1/8;\nB1=B1+e1/8;\nB3=B3+e1/8;\nA2=A2+e1/8;\nB2=B2+e1/8;\nA4=A4+e1/8;\nB4=B4+e1/8;\nE2=log(sin(A1*pi/180),10)*log(sin(A2*pi/180),10)*log(sin(A3*pi/180),10)*log(sin(A4*pi/180),10)-log(sin(B1*pi/180),10)*log(sin(B2*pi/180),10)*log(sin(B3*pi/180),10)*log(sin(B4*pi/180),10);\nc3=E2/17.1;\nc4=E2/17.0;\nA1=A1-c3;\nB1=B1+c3;\n\n#result\nprint \"corrected angle A1 in degrees\",A1,\"corrected angle B1 in degrees\",B1", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "corrected angle A1 in degrees 45.7991704812 corrected angle B1 in degrees 48.7774455188\n" + } + ], + "prompt_number": 2 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Surveying_Volume_3/Chapter4.ipynb b/Surveying_Volume_3/Chapter4.ipynb new file mode 100755 index 00000000..f2b5c7ec --- /dev/null +++ b/Surveying_Volume_3/Chapter4.ipynb @@ -0,0 +1,970 @@ +{ + "metadata": { + "name": "S3-C4" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Principles of Field Astronomy" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.1,Page 81" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ntheta=25+14.0/60;\nR=6370.0;\n\n#calculation\nAB=cos(theta*pi/180)*(16+10.0/60)\ndis=2*pi*R*AB/360.0;\n\n#result\nprint \"distance of AB in m \",round(dis,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "distance of AB in m 1625.863\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.2,Page 81" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nimport numpy as np\nfrom scipy import linalg\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nLongA=52+24.0/60;\nLongB=55+30.0/60;\nlatA=18+10.0/60;\nlatB=15.0;\nR=6370.0;#radius of earth\n\n#calculation\nP=LongB-LongA;\nPA=90-latA;\nPB=90-latB;\nAB=acos(cos(PB*pi/180)*cos(PA*pi/180)+sin(PB*pi/180)*sin(PA*pi/180)*cos(P*pi/180));\ndis=AB*R;\n#solving for A\na = np.array([[0.5,0.5], [0.5,-0.5]])\nb = np.array([atan(cos((PB/2-PA/2)*pi/180)/tan(P*pi/180)/cos((PB/2+PA/2)*pi/180))*180/pi,atan(sin((PB/2-PA/2)*pi/180)/tan(P*pi/180)/sin((PB/2+PA/2)*pi/180))*180/pi])\nx = np.linalg.solve(a, b)\nA=x[0];\nA=deg_to_dms(A);\n\n#result\nprint \"distance of AB in km \",round(dis,2)\nprint \"direction of B to A in deg min sec\",A", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "distance of AB in km 482.72\ndirection of B to A in deg min sec [117, 8, 30.72]\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.3,Page 83" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nb=40.0;\np=6.0;\n\n#calculation\na=pi/2-asin(cos(b*pi/180)*cos(p*pi/180));\nBC=a*180/pi-b;\nBC=BC*1.853*60;\nB=asin(sin(b*pi/180)/sin(a))\nB=deg_to_dms(B*180/pi);\n\n#result\nprint \"distance BC in km\",round(BC,2)\nprint \"angle of B deg min sec\",B", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "distance BC in km 41.43\nangle of B deg min sec [82, 53, 56.83]\n" + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.4,Page 105" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ncoal=90.0-28.0-24.0/60;#coaltitude\ncola=90.0-48.0-30.0/60;#colatitude\n\n\n#calculation\ndelta=pi/2-acos((cos(coal*pi/180)*cos(cola*pi/180)+sin(coal*pi/180)*sin(cola*pi/180)*cos(50*pi/180)));\nH=acos(cos(coal*pi/180)/(sin(cola*pi/180)*cos(delta))-tan(delta)/tan(cola*pi/180));\nH=deg_to_dms(H*180/pi); \ndelta=deg_to_dms(delta*180/pi); \n\n#result\nprint \"H in deg min sec\",H\nprint \"declination in deg mi sec\",delta", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "H in deg min sec [99, 7, 55.94]\ndelta in deg mi sec [46, 57, 38.26]\n" + } + ], + "prompt_number": 25 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.5,Page 107" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ndelta=22+45.0/60;\ntheta=55.0;\nH=45+15.0/60;\n\n#calculation\nalpha=asin((cos(H*pi/180)+tan(theta*pi/180)*tan(delta*pi/180))*cos(theta*pi/180)*cos(delta*pi/180));\nalpha=alpha*180/pi;\nA=acos(sin(delta*pi/180)/(cos(theta*pi/180)*cos(alpha*pi/180))-tan(alpha*pi/180)*tan(theta*pi/180));\nA=deg_to_dms(360-A*180/pi);\n\n#result\nprint \"azimuth in deg mi sec\",A", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "azimuth in deg mi sec [244, 40, 8.11]\n" + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.6,Page 109" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ncoal=42.0;\ncode=100+40.0/60;#co-declination\n\n#calculation\nalpha=pi/2-acos((cos(coal*pi/180)*cos(code*pi/180)+sin(coal*pi/180)*sin(code*pi/180)*cos(35*pi/180)));\nA=acos((cos(code*pi/180)-cos(coal*pi/180)*cos(pi/2-alpha))/(sin(coal*pi/180)*sin(pi/2-alpha)));\nA=deg_to_dms(A*180/pi);\nalpha=deg_to_dms(alpha*180/pi);\n\n#result\nprint \"azimuth in deg min sec\",A\nprint \"alpha in deg min sec\",alpha", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "azimuth in deg mi sec [142, 1, 26.46]\nalpha in deg min sec [23, 38, 48.21]\n" + } + ], + "prompt_number": 33 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.7,Page 113" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ntheta=33+42.0/60+34.0/3600;\ndelta=81+55.0/60+14.0/3600;\n\n#calculation\nH1=acos(tan(theta*pi/180)/tan(delta*pi/180));\nH1=deg_to_dms(360-H1*180/pi);\nalpha=asin(sin(theta*pi/180)/sin(delta*pi/180));\nalpha=deg_to_dms(alpha*180/pi);\nA=asin(cos(delta*pi/180)/cos(theta*pi/180));\nA=deg_to_dms(A*180/pi);\n\n\n#result\nprint \"azimuth in deg min sec\",A\nprint \"alpha in deg min sec\",alpha\nprint \"hour angle in deg min sec\",H1", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "azimuth in deg min sec [9, 43, 37.11]\nalpha in deg min sec [34, 5, 36.7]\nhour angle in deg min sec [275, 26, 3.85]\n" + } + ], + "prompt_number": 41 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.8,Page 114" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nalpha=0.0;\ntheta=25+45.0/60;\n\n#calculation\ndelta=asin(sin(alpha)/sin(theta*pi/180));\n\n#result\nprint \"declination in degrees\",delta", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "declination in degrees 0.0\n" + } + ], + "prompt_number": 42 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.9,Page 116" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ncola=90.0-49.0;#colatitude\ncode=90.0+19.0;#codeclination\n\n#calculation\nA1=acos((cos(code*pi/180)-cos(cola*pi/180)*cos(pi/2))/sin(cola*pi/180)/sin(pi/2));\nA1=deg_to_dms(360-A1*180/pi);\nH=acos((0-cos(cola*pi/180)*cos(code*pi/180))/sin(cola*pi/180)*sin(code*pi/180));\nH=deg_to_dms(H*180/pi)\n\n#result\nprint \"azimuth in deg min sec\",A1\nprint \"hour angle in deg min sec\",H", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "azimuth in deg min sec [240, 14, 52.45]\nhour angle in deg min sec [69, 15, 38.0]\n" + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.10,Page 120" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\n#part1\ndelta=38+15.0/60;\ntheta=25+10.0/60;\n\n#calculation\nz=delta-theta;\nalpha=90-z;\nz=deg_to_dms(z);\nalpha=deg_to_dms(alpha);\n\n#result\nprint \"zenith distance in deg min sec\",z\nprint \"altitude in deg min sec\",alpha\n\n#part2\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\ndelta=22+40.0/60;\ntheta=25+10.0/60;\n\n#calculation\nz=-delta+theta;\nalpha=90-z;\nz=deg_to_dms(z);\nalpha=deg_to_dms(alpha);\n\n#result\nprint \"zenith distance in deg min sec\",z\nprint \"altitude in deg min sec\",alpha\n\n#part3\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md+1)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ndelta=70+20.0/60;\ntheta=25+10.0/60;\n\n#calculation\nz=delta-theta;\nz=deg_to_dms(z);\n\n\n#result\nprint \"zenith distance in deg min sec\",z\n\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "zenith distance in deg min sec [13, 4, 60.0]\naltitude in deg min sec [76, 55, 0.0]\nzenith distance in deg min sec [2, 30, 0.0]\naltitude in deg min sec [87, 30, 0.0]\nzenith distance in deg min sec [45, 10, -0.0]\n" + } + ], + "prompt_number": 20 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.11,Page 122" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\ntheta=42+50.0/60;\ndelta=83+40.0/60;\n\n#calculation\nz=180-delta-theta;\nalpha=90-z;\nz=deg_to_dms(z);\nalpha=deg_to_dms(alpha);\n\n#result\nprint \"zenith distance in deg min sec\",z\nprint \"altitude in deg min sec\",alpha", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "zenith distance in deg min sec [53, 29, 60.0]\naltitude in deg min sec [36, 30, 0.0]\n" + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.12,Page 133" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#part1\ndol=110-82-30.0/60; #differnce in longitude\nST=18+35.0/60+10.0/3600;#standard time\n\n#calculation\nLMT=ST+dol/15.0;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec\",LMT\n\n#part2\ndol=82+30.0/60+30; #differnce in longitude\nST=18+35.0/60+10.0/3600;#standard time\n\n#calculation\nLMT=ST-dol/15.0;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec\",LMT\n\n#part1\ndol=82+30.0/60-30; #differnce in longitude\nST=18+35.0/60+10.0/3600;#standard time\n\n#calculation\nLMT=ST-dol/15.0;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec\",LMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec [20, 25, 10.0]\nLMT in hr min sec [11, 5, 10.0]\nLMT in hr min sec [15, 5, 10.0]\n" + } + ], + "prompt_number": 66 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.13,Page 134" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#part1\nLMT=8+30.0/60+15.0/3600;\nLong=45+30.0/60;\n\n#calculation\nGMT=LMT+Long/15.0;\nGMT=deg_to_dms(GMT)\n\n#result\nprint \"GMT in hr min sec (AM)\",GMT\n\n\n#part2\nLMT=6+40.0/60+10.0/3600;\nLong=55+30.0/60;\n\n#calculation\nGMT=LMT-Long/15.0;\nGMT=deg_to_dms(GMT)\n\n#result\nprint \"GMT in hr min sec (PM)\",GMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "GMT in hr min sec (AM) [11, 32, 15.0]\nGMT in hr min sec (PM) [2, 58, 10.0]\n" + } + ], + "prompt_number": 71 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.14,Page 134" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#part1\nGMT=20+30.0/60+15.0/3600;\nLong=82+30.0/60;\n\n#calculation\nLMT=GMT+Long/15.0-24;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec (next day)\",LMT\n\n#part1\nGMT=20+30.0/60+15.0/3600;\nLong=120.0;\n\n#calculation\nLMT=GMT-Long/15.0;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec (same day)\",LMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec (next day) [2, 0, 15.0]\nLMT in hr min sec (same day) [12, 30, 15.0]\n" + } + ], + "prompt_number": 74 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.15,Page 135" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#part1\nRA=6+15.0/60+20.0/3600;\nHA=8+10.0/60+30.0/3600;\n\n#calculation\nLST=RA+HA;\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST in hr min sec\",LST\n\n#part2\nRA=8+40.0/60+15.0/3600;\nHA=3+50.0/60+20.0/3600;\n\n#calculation\nLST=RA-HA;\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST in hr min sec\",LST\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LST in hr min sec [14, 25, 50.0]\nLST in hr min sec [4, 49, 55.0]\n" + } + ], + "prompt_number": 76 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.16,Page 135" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#part1\nLHA=6+30.0/60+10.0/3600;\n\n#calculation\nLAT=LHA+12;\nLAT=deg_to_dms(LAT);\n\n#result\nprint \"LAT in hr min sec\",LAT;\n\n#part2\nLHA=18+40.0/60+20.0/3600;\n\n#calculation\nLAT=LHA+12-24;\nLAT=deg_to_dms(LAT);\n\n#result\nprint \"LAT in hr min sec\",LAT;", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LAT in hr min sec [18, 30, 10.0]\nLAT in hr min sec [6, 40, 20.0]\n" + } + ], + "prompt_number": 78 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.17,Page 136" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=60.0;\nLHA=5+30.0/60+20.0/3600;\n\n#calculation\nGMT=LHA+12-Long/15.0;\nGMT=deg_to_dms(GMT);\n\n#result\nprint \"GMT in hr min sec\",GMT\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "GMT in hr min sec [13, 30, 20.0]\n" + } + ], + "prompt_number": 79 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.18,Page 138" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#parta\nGMT=0;\nET=10.0/60+1.8/3600;\n\n#calculaion\nGAT=GMT+ET;\nGAT=deg_to_dms(GAT);\n\n#result\nprint \"GAT in hr min sec\",GAT\n\n#partb\nGMT=0;\nET=-13.0/60-28.5/3600;\n\n#calculaion\nGAT=GMT+ET+24;\nGAT=deg_to_dms(GAT);\n\n#result\nprint \"GAT in hr min sec\",GAT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "GAT in hr min sec [0, 10, 1.8]\nGAT in hr min sec [23, 46, 31.5]\n" + } + ], + "prompt_number": 81 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.19,Page 138" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nET1=-3.0/60-51.4/3600; #ET at april 2\nET2=-3.0/60-33.5/3600; #ET at april 3\n\n#calculation\ndET=(ET2-ET1)*18.0/24; #change in ET\nET=ET1+dET;\nET=deg_to_dms(ET);\n\n#result\nprint \"ET (-ve) in hr min s\",ET", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "ET in hr min s [0, 3, 37.97]\n" + } + ], + "prompt_number": 83 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.20,Page 138" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLAT=15+12.0/60+40.0/3600;\nLong=20+3.0/60;\nGMN=5.0/60+10.65/3600;\n\n#calculation\nGAT=LAT+Long/15.0;\ne1=(GAT-12)*0.22/3600+GMN;\nLAT=GAT+e1-Long/15.0;\nLAT=deg_to_dms(LAT);\n\n#result\nprint \"LAT in hr min sec\",LAT\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LAT in hr min sec [15, 17, 51.65]\n" + } + ], + "prompt_number": 84 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.21,Page 139" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=56+35.0/60;\nLMT=7+15.0/60+25.0/3600;\nGMN=3.0/60+54.0/3600;\n\n#calculation\nGMT=LMT-Long/15.0;\ne1=GMN+(12-GMN)*0.25/3600;\nLAT=GMT+Long/15.0+e1;\nLAT=deg_to_dms(LAT);\n\n#result\nprint \"LAT in hr min sec\",LAT\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LAT in hr min sec [7, 19, 21.98]\n" + } + ], + "prompt_number": 86 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.22,Page 142" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nST=7+15.0/60+30.0/3600;\n\n#calculation\ntr=9.8296/3600*ST;\nST=ST-tr;\nST=deg_to_dms(ST);\n\n#result\nprint \"solar mean time in hr min sec\",ST", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "solar mean time in hr min sec [7, 14, 18.65]\n" + } + ], + "prompt_number": 87 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.23,Page 142" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nST=7+45.0/60+50.0/3600;\n\n#calculation\ntr=9.8565/3600*ST;\nST=ST+tr;\nST=deg_to_dms(ST);\n\n#result\nprint \"solar mean time in hr min sec\",ST", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "solar mean time in hr min sec [7, 47, 6.52]\n" + } + ], + "prompt_number": 90 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.24,Page 143" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding LST\n\n#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=140+35.0/60+20.0/3600;\nGST=13+15.0/60+30.0/3600;\n\n#calculation\nLST=GST+9.8565/3600*Long/15.0;\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST in deg min sec\",LST\n\n#part2\nLong=160+45.0/60+30.0/3600;\nGST=13+15.0/60+30.0/3600;\n\n#calculation\nLST=GST-9.8565/3600*Long/15.0;\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST in deg min sec\",LST", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LST in deg min sec [13, 17, 2.38]\nLST in deg min sec [13, 13, 44.37]\n" + } + ], + "prompt_number": 93 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.25,Page 145" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=75.0;\nGST=15+55.0/60+13.0/3600;\nLMT=11.0;\n\n#calculation\nLST=GST+Long/15.0*9.8565/3600+LMT+9.8565/3600*LMT;\nLST=LST-24;\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST of next day in deg min sec\",LST", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LST of next day in deg min sec [2, 57, 50.7]\n" + } + ], + "prompt_number": 95 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.26,Page 145" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=75.0;\nGST=5+25.0/60+15.0/3600;\nLMT=5.0+40.0/60;\n\n#calculation\nLST=GST-Long/15.0*9.8565/3600+LMT+9.8565/3600*(LMT);\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST of next day in deg min sec\",LST", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LST of next day in deg min sec [11, 5, 21.57]\n" + } + ], + "prompt_number": 96 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.27,Page 146" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=75.0;\nGST=11+45.0/60+10.0/3600;\nLST=26+35.0/60+42.0/3600;\n\n#calculation\nLMM=GST+Long/15*9.8565/3600;\nLMT=LST-LMM-(LST-LMM)*9.8296/3600;\nLMT=deg_to_dms(LMT-12)\n\n#result\nprint \"LMT in hr min sec (PM)\",LMT \n\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec (PM) [2, 47, 16.96]\n" + } + ], + "prompt_number": 97 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.28,Page 147" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=90.0;\nGST=9+15.0/60+14.0/3600;\nLST=31+35.0/60+12.0/3600;\n\n#calculation\nLMM=GST-Long/15*9.8565/3600;\nLMT=LST-LMM-(LST-LMM)*9.8296/3600;\nLMT=deg_to_dms(LMT-12)\n\n#result\nprint \"LMT in hr min sec (PM)\",LMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec (PM) [10, 17, 17.46]\n" + } + ], + "prompt_number": 100 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.29,Page 149" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=150.0;\nGST=12+25.0/60+15.0/3600;\nLST=30+15.0/60+10.0/3600;\nLMN=18+15.0/60+10.0/3600;\n\n#calculation\nLMM=GST+Long/15*9.8565/3600;\nLMT=LMN-LMM-(LMN-LMM)*9.8296/3600;\nLMT=deg_to_dms(LMT)\n\n#result\nprint \"LMT in hr min sec (PM)\",LMT\n\n#part2\nLMM=GST+Long/15*9.8565/3600;\nLMT=LST-LMM-(LST-LMM)*9.8296/3600;\nLMT=deg_to_dms(LMT-12)\n\n#result\nprint \"LMT in hr min sec (AM)\",LMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec (PM) [5, 47, 19.38]\nLMT in hr min sec (AM) [5, 45, 21.42]\n" + } + ], + "prompt_number": 102 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.30,Page 150 " + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\ntheta=35+15.0/60+20.0/3600;\ndelta=88+15.0/60+45.0/3600;\nRA=1+45.0/60+15.0/3600;\n\n#calculation\nH=acos(tan(theta*pi/180)/tan(delta*pi/180));\nH=H*180/pi/15;\nLST=RA+H;\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST in hr min sec\",LST", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LST in hr min sec [7, 40, 20.12]\n" + } + ], + "prompt_number": 103 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.31,Page 151" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nRA=22+25.0/60+10.0/3600;\nST=14+45.0/60;#sidereal time\nLong=90.0;\nGMT=27+15.0/60;\nLMN=15+21.0/60+15.0/3600;\n\n#calculation\nMT=12+24-GMT+Long/15;#mean time interval\nacc=9.8565/3600*MT;\nHA=LMN-acc-MT;\nLST=HA+24-RA;\nLST=deg_to_dms(LST);\n\n#result\nprint \"LST in hr min sec\",LST\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LST in hr min sec [2, 8, 39.62]\n" + } + ], + "prompt_number": 105 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.32,Page 152" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=60;\nGMT=11+15.0/60+20.0/3600;\n\n#calculation\nLMT=GMT-Long/15*9.8296/3600;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec \",LMT\n\n#part2\nLong=45;\nGMT=11+15.0/60+20.0/3600;\n\n#calculation\nLMT=GMT+Long/15*9.8296/3600;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec \",LMT\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec [11, 14, 40.68]\nLMT in hr min sec [11, 15, 49.49]\n" + } + ], + "prompt_number": 107 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": "Example 4.33,Page 153" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding LMT\n\n#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=75;\nGMT=6+18.0/60+20.0/3600;\nLST=10+25.0/60+15.0/3600;\n\n#calculation\nLMT=GMT+Long/15*9.8296/3600;\nLST=LST-9.8296/3600*LST;\nLMT=LMT+LST;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec \",LMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec [16, 42, 41.72]\n" + } + ], + "prompt_number": 111 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.34,Page 154" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nGST=7+35.0/60+40.0/3600;\n\n#calculation\nGMT=24-GST-(24-GST)*9.8296/3600;\nGMT=deg_to_dms(GMT);\n\n#result\nprint \"GMT in hr min sec\",GMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "GMT in hr min sec [16, 21, 38.74]\n" + } + ], + "prompt_number": 112 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": "Example 4.34,Page 154" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding LMT\n\n#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=120.0;\nGMT=12+3.0/60+24.6/3600;\n\n#calculation\nLMT=GMT-17.8/24*Long/15.0/3600;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT of LAN in hr min sec\",LMT\n\n#part2\nLong=45;\nGMT=12+3.0/60+24.6/3600;\n\n#calculation\nLMT=GMT+17.8/24*Long/15.0/3600;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT of LAN in hr min sec\",LMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT of LAN in hr min sec [12, 3, 18.67]\nLMT of LAN in hr min sec [12, 3, 26.83]\n" + } + ], + "prompt_number": 114 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.35,Page 156" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#one of the 2 solution\nf0=5+1.9/60;\nn=0.25;\ndel0=0;\ndel1=-0.1;\nd2=23.0;#del1/2\n\n#calculation\nfn=f0+n*d2/60+n*(n-1)/2*(del1+del0)/60;\nfn=deg_to_dms(fn)\n\n#result\nprint \"sun declination in deg min sec\",fn", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "sun declination in deg min sec [5, 7, 39.56]\n" + } + ], + "prompt_number": 117 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.36,Page 157" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nalt=23+40.0/60;\nazi=145.0;\nlat=50.0;\n\n#calculation\ndelta=pi/2-acos(cos(pi/2-lat*pi/180)*cos(pi/2-alt*pi/180)+sin(pi/2-lat*pi/180)*sin(pi/2-alt*pi/180)*cos(azi*pi/180));\nH=acos((cos(pi/2-alt*pi/180)-cos(pi/2-lat*pi/180)*cos(pi/2-delta))/sin(pi/2-lat*pi/180)*sin(pi/2-delta))\nH=deg_to_dms(360-H*180/pi);\ndelta=deg_to_dms(delta*180/pi);\n\n#result\nprint \"there is a calculation mistake in calculating H in the book\"\nprint \"declination in deg min sec\",delta\nprint \"hour angle in deg min sec\",H", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "declination in deg min sec [-10, 3, 51.85]\nhour angle in deg min sec [325, 4, 44.11]\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.37,Page 158" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nalt=25+30.0/60;\nazi=45.0;\nlat=42.0;\n\n#calculation\ndelta=pi/2-acos(cos(pi/2-lat*pi/180)*cos(pi/2-alt*pi/180)+sin(pi/2-lat*pi/180)*sin(pi/2-alt*pi/180)*cos(azi*pi/180));\nH=acos((cos(pi/2-alt*pi/180)-cos(pi/2-lat*pi/180)*cos(pi/2-delta))/sin(pi/2-lat*pi/180)*sin(pi/2-delta))\nH=deg_to_dms(H*180/pi);\ndelta=deg_to_dms(delta*180/pi);\n\n#result\nprint \"declination in deg min sec\",delta\nprint \"hour angle in deg min sec\",H", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "declination in deg min sec [49, 40, 22.06]\nhour angle in deg min sec [93, 58, 30.79]\n" + } + ], + "prompt_number": 122 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.38,Page 158" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding hour angle,azimuth\n\n#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\ndelta=21.0*pi/180+25.0/60*pi/180;\nlat=25+40.0/60;\n\n#calculation\nH=acos((0-cos(pi/2-lat*pi/180)*cos(pi/2-delta))/sin(pi/2-lat*pi/180)*sin(pi/2-delta));\nA=acos(cos(pi/2-delta)/sin(pi/2-lat*pi/180));\nH=deg_to_dms(360-H*180/pi);\nA=deg_to_dms(A*180/pi);\n\n#result\nprint \"hour angle in deg min sec\",H\nprint \"azimuth in deg min sec\",A", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "hour angle in deg min sec [260, 35, 53.78]\nazimuth in deg min sec [66, 6, 4.59]\n" + } + ], + "prompt_number": 130 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.39,Page 161" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding latitude\n\n#initialisation of variable\nimport numpy as np\nfrom scipy import linalg\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\n#calculation\nA=np.array([[1,1],[1,-1]]);\nb=np.array([[90],[0]]);\nx=np.linalg.solve(A,b);\n\n#result\nprint \"latitude in degrees\",x[0]", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "latitude in degrees [ 45.]" + }, + { + "output_type": "stream", + "stream": "stdout", + "text": "\n" + } + ], + "prompt_number": 131 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.40,Page 161" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding altitude\n\n#initialisation of variable\nimport numpy as np\nfrom scipy import linalg\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\ntheta=53+20.0/60;\ndelta=53+20.0/60;\n\n#calculation\nalpha=theta+delta-90;\nalpha=deg_to_dms(alpha);\n\n#result\nprint \"altitude in deg min sec\",alpha", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "altitude in deg min sec [16, 40, 0.0]\n" + } + ], + "prompt_number": 132 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.41,Page 162" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding GAT\n\n#initialisation of variable\nimport numpy as np\nfrom scipy import linalg\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nGMT=18+30.0/60;\nET=1.0/60+25.4/3600-0.67*6.5/3600;\n\n#calculation\nGAT=GMT+ET;\nGAT=deg_to_dms(GAT)\n\n#result\nprint \"GAT in hr min sec\",GAT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "GAT in hr min sec [18, 31, 21.05]\n" + } + ], + "prompt_number": 133 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.42,Page 163" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding LMT\n\n#initialisation of variable\nimport numpy as np\nfrom scipy import linalg\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=30.0;\nGAT=13+15.0/60+10.0/3600;\nET=6.0/60+15.35/3600+0.3/3600*1.25278;\n\n#calculation\nLMT=GAT+ET-Long/15.0;\nLMT=deg_to_dms(LMT);\n\n#result\nprint \"LMT in hr min sec\",LMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LMT in hr min sec [11, 21, 25.73]\n" + } + ], + "prompt_number": 134 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.43,Page 163" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding LHA\n\n#initialisation of variable\nimport numpy as np\nfrom scipy import linalg\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nLong=45.0;\nE=11+55.0/60+5.0/3600-1.5/6*17.0/3/3600;\n\n#calculation\nGMT=14+40.0/60+Long/15.0;\nGHA=GMT+E;\nLHA=GHA-24-Long/15;\nLHA=deg_to_dms(LHA);\n\n#result\nprint \"LHA of the sun in hr min sec\",LHA\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "LHA of the sun in hr min sec [2, 35, 3.58]\n" + } + ], + "prompt_number": 135 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 4.44,Page 164" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding GMT\n\n#initialisation of variable\nimport numpy as np\nfrom scipy import linalg\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n if sd==60:\n m=m+1;\n sd=0;\n \n sd=round(sd,2)\n return [d, m, sd]\nE=11+55.0/60+24.0/3600+0.5/3600*(2+40.0/60+21.2/3600);\nGHA=8+35.0/60+45.2/3600;\n\n#calculation\nGMT=GHA+24-E;\nGMT=deg_to_dms(GMT);\n\n#result\nprint \"GMT in hr min sec\",GMT", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "GMT in hr min sec [20, 40, 19.86]\n" + } + ], + "prompt_number": 136 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Surveying_Volume_3/Chapter5.ipynb b/Surveying_Volume_3/Chapter5.ipynb new file mode 100755 index 00000000..4d4ee9e6 --- /dev/null +++ b/Surveying_Volume_3/Chapter5.ipynb @@ -0,0 +1,398 @@ +{ + "metadata": { + "name": "S3-C5" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Uses Of Field Astronomy in surveying" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.1,Page 174" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\n#printing result in degree minute and seconds respectively\nl1=11.5;\nl2=13.5;\nr1=8.5;\nr2=6.5;\nalpha=3+15.0/60+28.0/3600;\nOB=121+45.0/60+18.0/3600;\nOA=43+25.0/60+20.51/3600;\n\n#calculation\ngamma=(l1+l2)/4-(r1+r2)/4;\ne=gamma*tan(alpha*pi/180)/3600; #correction\nCH=OB-OA-e;\nCH=deg_to_dms(CH);\n\n#result\nprint \"corrected horizontal angle in deg,min,sec respectively\",CH", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "corrected horizontal angle in deg,min,sec respectively [78, 19, 57.35]\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.2,Page 184" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nalpha=30+32.0/60+18.0/3600;#latitude\nd=16.0/60+2.85/3600; #diameter of sun\n\n#calculation\nC1=-58.0/3600/tan(alpha*pi/180); \nC2=8.8/3600*cos(alpha*pi/180); \nC3=d; \nCL=alpha+C1+C2+C3;\nCL=deg_to_dms(CL);\n\n#result\nprint \"corrected latutude in deg,min,sec respectively\",CL", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "corrected latutude in deg,min,sec respectively [30, 46, 50.12]\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.3,Page 184" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nalpha=40+52.0/60+10.0/3600;#latitude\n\n#calculation\nC1=-58.0/3600/tan(alpha*pi/180);\nCL=alpha+C1;\nCL=deg_to_dms(CL);\n\n#result\nprint \"corrected latutude in deg,min,sec respectively\",CL", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "corrected latutude in deg,min,sec respectively [40, 51, 2.97]\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.4,Page 197" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nLMT=21+23.0/60+05.0/3600; #local chronometer time\nLong=65.0+19.0/60; #longitude\nGST=13+15.0/60+20.0/3600;\nRA=9+32.0/60+15.0/3600;\nLong2=82.0+30.0/60; #longitude of India\n\n#calculation\ne1=Long/15*9.8565/3600; #error\nSIT=RA+24-GST+e1; #sidereal time interval after LMM\ne2=SIT*9.8296/3600; #error\nMI=SIT-e2; #mean interval after LMM\nLMT=LMT-(Long2-Long)/15.0;\nCE=MI-LMT;\nCE=deg_to_dms(CE);\n\n#result\nprint \"chronometer error in hours,min,sec respectively\",CE", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "chronometer error in hours,min,sec respectively [0, 0, 2.56]\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.5,Page 198" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nMST=12+32.0/60+15.0/3600; #mean sidereal time\nRA=15+45.0/60+10.0/3600;\ntheta=55+14.0/60+20.0/3600;#latitude\ndelta=15+24.0/60+30.0/3600;#declination\nalpha=35+44.0/60+10.0/3600;#zenith deistance\n\n#calculation\nc=90-theta;\np=90-delta;\nz=90-alpha;\nH=acos(cos(z*pi/180)/sin(c*pi/180)/sin(p*pi/180)-1/(tan(p*pi/180)*tan(c*pi/180)))\nH=H/15*180/pi;\nLST=RA-H;\nCE=MST-LST;\nCE=deg_to_dms(CE);\n\n#result\nprint \"chronometer error in hours,min,sec respectively\",CE", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "chronometer error in hours,min,sec respectively [0, 0, 12.94]\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.6,Page 199" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nLMTe=6+34.0/60+18.0/3600;# LMT east\nLMTw=8+58.0/60+2.0/3600; # LMT west\nRA=16+11.0/60+25.0/3600;\nLong=125+33.0/60;\nGST=8+25.0/60+14.0/3600;\n\n#calculation\ne1=Long/15*9.8565/3600; #error\nSIT=RA-GST+e1; #sidereal time interval after LMM\ne2=SIT*9.8296/3600;\nMI=SIT-e2; #mean time interval after LMM\nLMTav=(LMTe+LMTw)/2; #mean LMT\nCE=LMTav-MI;\nCE=deg_to_dms(CE);\n\n#result\nprint \"chronometer error in slower side in hours,min,sec respectively\",CE", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "chronometer error in slower side in hours,min,sec respectively [0, 0, 6.9]\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.7,Page 204" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nLMM=15+9.0/60+5.21/3600;# mean LMT\nGMT=10+9.0/60+3.76/3600;\nLong=75.0;#longitude\nalpha=42+30.0/60+42.0/3600;\ntheta=34+48.0/60+12.0/3600;\ndelta=15+36.0/60+48.0/3600;\n\n#calculation\nH=acos(sin(alpha*pi/180)/cos(theta*pi/180)/cos(delta*pi/180)-(tan(delta*pi/180)*tan(theta*pi/180)))\nH=H/15*180/pi;\nGAT=12+H-Long/15;\nLMT=GAT+Long/15-5.0/60-40.0/3600;\nCE=LMM-LMT;\nCE=deg_to_dms(CE);\n\n#result\nprint \"chronometer error in slower side in hours,min,sec respectively\",CE", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "[0, 0, 1.45] chronometer error in slower side in hours,min,sec respectively\n" + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.8,Page 219" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nd=5+1.0/60+50.0/3600;\ndel1=75+14.0/60+20.0/3600;\ndel2=70+12.0/60+30.0/3600;\n\n#calculation\nk=cos(del1*pi/180)/cos(del2*pi/180);\nA2=pi/2-atan((cos(d*pi/180)-k)/sin(d*pi/180));\nA2=A2*180/pi;\nA2=120+15.0/60+10.0/3600-A2;\nCR=360-A2;\nA2=deg_to_dms(A2);\nCR=deg_to_dms(CR);\n\n#result\nprint \"azimuth of angle R in degree,minites,seconds respectively\",A2\nprint \"true bearing of CR in degree,minites,seconds respectively\",CR", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "azimuth of angle R in degree,minites,seconds respectively [100, 27, 40.0]\ntrue bearing of CR in degree,minites,seconds respectively [259, 32, 20.0]\n" + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.9,Page 223" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\na=26.0/60+51.0/3600;\np=56.0/60+5.1/3600;#polar distance\n\n#calculation\nH=acos(a/p);\nA=p*sin(H)/cos(30.75694*pi/180);\nCR=25+35.0/60+40.0/3600-A;\nCR=deg_to_dms(CR);\n\n#result\nprint \"azimuth of angle CR in degree,minites,seconds respectively\",CR", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "azimuth of angle CR in degree,minites,seconds respectively [24, 38, 22.01]\n" + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.10,Page 227 " + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nLong=75.0;#longitude\nGST=11+40.0/60+32.4/3600;\nRA=12+25.0/60+18.35/3600;\nGMT=15+45.0/60+25.3/3600;\ndelta=22+6.0/60+32.5/3600;\n\n#calculation\ne1=Long/15*9.8565/3600;\nLSTofLMM=GST-e1;\nLMT=GMT+Long/15;\nSIT=LMT+LMT*9.8565/3600;#sidereal time interval\nLHA=SIT+LSTofLMM;\nH=RA+24-LHA;\nH=H*15;\nB=atan(tan(delta*pi/180)*tan(H*pi/180));\nB=B*180/pi;\nA=atan(tan(H*pi/180)*cos(B*pi/180)/sin((B-32-15.0/60)*pi/180))\nA=A*180/pi;\nTB=360+A-135-15.0/60-20.0/3600;\nTB=deg_to_dms(TB);\n\n#result\nprint \"true bearing TB in degree,minites,seconds respectively\",TB\nprint \"there is slight difference in the answers due to rounding off error in the book\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "true bearing TB in degree,minites,seconds respectively [313, 17, 36.07]\nthere is slight difference in the answers due to rounding off error in the book\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.11,Page 237" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nz=51+47.0/60+18.0/3600;#zenith distance\np=88+57.0/60+57.0/3600;#polar distance\nc=61+27.0/60+55.0/3600;#co-latitude\n\n#calculation\ns=(z+p+c)/2;\nA=2*atan(sqrt(sin((s-z)*pi/180)/sin(s*pi/180)*sin((s-c)*pi/180)/sin((s-p)*pi/180)));\nA=A*180/pi;\nTB=360-A-165-18.0/60-20.0/3600;\nTB=deg_to_dms(TB);\n\n#result\nprint \"true bearing TB in degree,minites,seconds respectively\",TB", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "true bearing TB in degree,minites,seconds respectively [80, 59, 47.52]\n" + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.12,Page 241" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#finding latitude\n\n#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nz2=90-40-13.0/60-15.0/3600;\ndel2=12+15.0/60+30.0/3600;#declination of star\n\n#calculation\ntheta=z2+del2;\ntheta=deg_to_dms(theta);\n\n#result\nprint \"altitude in degree,minites,seconds respectively\",theta", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "altitude in degree,minites,seconds respectively [62, 2, 15.0]\n" + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.13,Page 244" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nalpha1=30+45.0/60+25.0/3600;\nalpha2=40+48.0/60+30.0/3600;\n\n#calculation\ne1=-58/3600/tan(alpha1*pi/180); #error 1\ne2=-58/3600/tan(alpha2*pi/180); #error 2\ntheta=(alpha1+alpha2+e1+e2)/2;\ntheta=deg_to_dms(theta)\n\n#result\nprint \"latitude in degree,minites,seconds respectively\",theta", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "latitude in degree,minites,seconds respectively [34, 21, 48.1]\n" + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.14,Page 258" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nZP=37+29.0/60+40.0/3600;#colatitde\nZM=56+24.0/60+50.0/3600;#coaltitude\nPM=67+54.0/60+24.0/3600;#codeclination\n\n#calculation\nA1=acos((cos(PM*pi/180)-cos(ZP*pi/180)*cos(ZM*pi/180))/(sin(ZP*pi/180)*sin(ZM*pi/180)));\nA1=A1*180/pi;\nA=360-A1;\nA=deg_to_dms(A);\n\n#result\nprint \"azimuth of sun in degree,minites,seconds respectively\",A", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "azimuth of sun in degree,minites,seconds respectively [262, 53, 12.16]\n" + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.15,Page 259" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\ntheta=54+30.0/60;#latitude\ndelta=62+12.0/60+21.0/3600;#declination\n\n#calculation\nalpha=asin(sin(theta*pi/180)/sin(delta*pi/180));\nA1=acos(tan(theta*pi/180)/tan(alpha));\nA1=A1*180/pi;\nTB=360-A1-65-18.0/60-42.0/3600;\nTB=deg_to_dms(TB);\nalpha=deg_to_dms(alpha*180/pi);\nH=atan(tan(theta*pi/180)/tan(delta*pi/180));\nH=deg_to_dms(H*180/pi);\n\n#result\nprint \"true bearing in degree,minites,seconds respectively\",TB\nprint \"altitude in degree,minites,seconds respectively\",alpha\nprint \"hour angle in degree,minites,seconds respectively\",H", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "true bearing in degree,minites,seconds respectively [241, 16, 19.55]\naltitude in degree,minites,seconds respectively [66, 58, 7.13]\nhour angle in degree,minites,seconds respectively [36, 27, 49.32]\n" + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.16,Page 261" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nalpha=44+12.0/60+30.0/3600;\nd=15.0/60+45.86/3600;#diameter correction\nLong=7+20.0/60+15.0/3600;#longitude\n\n#calculation\nalpha=alpha+d-58/3600/tan(alpha)+8.8/3600*cos(alpha);\nGAT=Long/15;\ne2=6.82/3600*GAT;\ndelta=22+18.0/60+12.8/3600+e2;\ntheta=delta+90-alpha;\ntheta=deg_to_dms(theta);\n\n#result\nprint \"altitude in degree,minites,seconds respectively\",theta", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "[67, 49, 51.7] altitude in degree,minites,seconds respectively\n" + } + ], + "prompt_number": 47 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.17,Page 262" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nGMT=16+22.0/60+55.0/3600;\nET=3.0/60+43.0/3600;\nc=90-42-20.0/60;\np=90-18-45.0/60-50.0/60;\nz=90-43-38.0/60;\n\n#calculation\nH=acos(cos(z*pi/180)/sin(c*pi/180)/sin(p*pi/180)-1/tan(c*pi/180)*1/tan(p*pi/180));\nH=H*180/pi;\nLAT=12-H/15;\nLMT=LAT-ET;\nLong=GMT-LMT;\nLong=Long*15;\nLong=deg_to_dms(Long);\n\n#result\nprint \"Longitude in degree,minites,seconds respectively\",Long", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Longitude in degree,minites,seconds respectively [114, 50, 53.21]\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.18,Page 263" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,tan,sqrt,sin,cos,acos,atan,asin\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nalpha=21+35.0/60+30.0/3600;#mean observed altitude\nC=(4.5+5.5-3.5-2.5)/4*15.0/3600;\nc=44+30.0/60;#colatitude\nz=68+26.0/60+34.0/3600;#coaltitude\np=94+4.0/60+15.0/3600;#codeclination\ns=(c+p+z)/2;\n\n#calculation\ncr=-58/3600/tan(alpha);#correction refraction\ncp=8.8/3600*cos(alpha);#correction parallax\nalpha=alpha+C+cr+cp; #corrected altitude\nA=2*atan(sqrt(sin((s-z)*pi/180)/sin(s*pi/180)*sin((s-c)*pi/180)/sin((s-p)*pi/180)));\nA=A*180/pi;\nMh=(121+45.0/60+20.0/3600+122+47.0/60)/2;#mean horizontal angle\nAZ=360-Mh-A;\nAZ=deg_to_dms(AZ);\n\n#result\nprint \"Azimuth from north(clockwise) in degree,minites,seconds respectively\",AZ", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "[117, 0, 19.45] Azimuth from north(clockwise) in degree,minites,seconds respectively\n" + } + ], + "prompt_number": 55 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Surveying_Volume_3/Chapter6.ipynb b/Surveying_Volume_3/Chapter6.ipynb new file mode 100755 index 00000000..3fb19c9c --- /dev/null +++ b/Surveying_Volume_3/Chapter6.ipynb @@ -0,0 +1,629 @@ +{ + "metadata": { + "name": "S3-C6" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Photogrammetry" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.1,Page 281" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nDa=184.32;\nDb=95.84;\nAx=-115.0\nBy=-115.0;\n\n#calculation\nphi=atan(Ax/By);\nAB=sqrt(Ax**2+By**2);\ntheta=acos((Da**2+AB**2-Db**2)/2/Da/AB);\nalpha=phi-theta;\nxc=Da*cos(alpha)-115.0;\nyc=-Da*sin(alpha);\n\n#result\nprint \"the coordiantes in mm x is\",round(xc),\"y is\",round(yc,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "the coordiantes in mm x is 64.0 y is -43.72\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.2,Page 290" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nf=0.152;\nH=1800;#elevation of topmost point\nh=300;#elevation of ground\n\n#calculation\nS=f/(H-h);\n\n#result\nprint \"scale of photograph in 1 in\",round(1/S)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "scale of photograph in 1 in 9868.0\n" + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.3,Page 290" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nf=0.150 #focal length\nh1=1500.0;#elevation A\nh2=1200.0;#elevation B\nh3=1000.0;#elevation C\nH=3000.0;#height\n\n#calculation\nhav=1.0/3*(h1+h2+h3);\nS1=f/(H-h1);\nS2=f/(H-h2);\nS3=f/(H-h3);\nSav=f/(H-hav);\n\n#result\nprint \"scale of point 1 in 1 in\",round(1/S1);\nprint \"scale of point 2 in 1 in\",round(1/S2);\nprint \"scale of point 3 in 1 in\",round(1/S3);\nprint \"average scale in 1 in\",round(1/Sav);", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "scale of point 1 in 1 in 10000.0\nscale of point 2 in 1 in 12000.0\nscale of point 3 in 1 in 13333.0\naverage scale in 1 in 11778.0\n" + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.4,Page 292" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nab=188.0;\nAB=120;\nSm=1.0/20000;\n\n#calculation\nS=ab/AB*Sm;\n\n#result\nprint \"scale of photograph in 1 in\",round(1/S)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "scale of photograph in 1 in 12766.0\n" + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.5,Page 295" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nxa=45.35;\nxb=-40.16;\nf=152.4;#focal length\nH=1500.0;#actual height\nha=200.0;#height A\nhb=150.0;#height B\nya=38.41;\nyb=-45.65;\n\n#calculation\nXa=xa*(H-ha)/f;\nYa=ya*(H-hb)/f;\nXb=xb*(H-ha)/f;\nYb=yb*(H-hb)/f;\nAB=sqrt((Xb-Xa)**2+(Yb-Ya)**2);\n\n#result\nprint \"distance of AB in m\",round(AB,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "distance of AB in m 1042.361\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.6,Page298" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nd=62.4;#displacement\nH=250.0;#height of datum\nr=115.4;#image distance\n\n#calculation\nh=H*d/r;\n\n#result\nprint \"height of chimney in m\",round(h,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "height of chimney in m 135.182\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.7,Page 299" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nra=88.25;#image distance of A\nrb=81.23;#image distance of B\nrc=68.14;#image distance of C\nH=2000.0;\nha=255;\nhb=200;\nhc=145;\nf=0.1524;\n\n#calculation\naa=ra*ha/H;\nbb=rb*hb/H;\ncc=rc*hc/H;\nS=f/H;\n\n#result\nprint \"relief distance of A in mm\",round(aa,2)\nprint \"relief distance of B in mm\",round(bb,2)\nprint \"relief distance of C in mm\",round(cc,2)\nprint \"scale of photograph in 1 in\",round(1/S)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "relief distance of A in mm 11.25\nrelief distance of B in mm 8.12\nrelief distance of C in mm 4.94\nscale of photograph in 1 in 13123.0\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.8,Page 300" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nf=0.1524;#focal length\nS=0.08251/1000;#scale\n\n#calculation\nH=f/S;\n\n#result\nprint \"flying height in m\",round(H,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "flying height in m 1847.049\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.9,Page 300" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nimport numpy as np\n\n#solving the quadratic polynomial in H\n#0=0.4064-365.929H-289685.07\ncoeff=[0.4064, -365.929,-289685.926];\nH=np.roots(coeff);\n\n#result\nprint \"height required in m\",round(H[0],2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "height required in m 1407.02\n" + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.10,Page 309" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nS=1.0/10000;#scale\nA=500.0;#area\npw=0.3;\nl=0.23;\nw=0.23;\n\n#calculation\na=(1-0.6)*(1-pw)*l*w/S**2/1000/1000;\nN=A/a;\n\n#result\nprint \"no. of photographs taken\",round(N)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "no. of photographs taken 338.0\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.11,Page 310" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nL=25.0e3;\nk=0.23e4;#l/s=w/s;\npl=0.6;\npw=0.3;\nW=20.0e3;\n\n#calculation\nN=((L/((1-pl)*k)+1))*((W/((1-pw)*k)+1)+1);\n\n#result\nprint \"no. of photographs taken\",round(N)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "no. of photographs taken 406.0\n" + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.12,Page 310" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nf=0.1524;#focal length\nS=1.0/10000;#scale\npw=0.3;#side lap\nw=0.23;#format width\npl=0.6;\nl=0.23;\n\n#calculation\nW=(1-pw)/S*w;\nH=f/S+300;\nN2=30/W+1;\nN2=round(N2)\nL=(1-pl)*1/S*l/1000;\nT=3600*L/240.0;\nAd=T*240e3/60.0/60.0;#adjusted ground distance\nN1=40.0e3/Ad+1;\nN1=round(N1)\nN=N1*N2;\n\n#result\nprint \"height over datum in m\",H\nprint \"no. of flight strips\",round(N2-1)\nprint \"length of each photograph cover in km\",round(L,3)\nprint \"exposure time in s\",round(T)\nprint \"no. of photographs taken\",round(N)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "height over datum in m 1824.0\nno. of flight strips 0.0\nlength of each photograph cover in km 0.92\nexposure time in s 14.0\nno. of photographs taken 44.0\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.12b,Page 317" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nt=3.0/180*pi;\nya=82.25;\nxa=-62.45;\ns=220;\nf=152.4;#focal length\nH=2500.0e3;\nh=500.0e3;\n\n#calculation\ntheta=s-180;\nya_dash=xa*sin(theta*pi/180)+ya*cos(theta*pi/180)+f*tan(t)\nS=(f/cos(t)-ya_dash*sin(t))/(H-h);\n\n#result\nprint \"scale of photograph in 1 in\",round(1/S)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "scale of photograph in 1 in 13246.0\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.13,Page 319" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nt=3*pi/180;\nxa=-62.45;\nxb=78.25;\nf=152.4;#focal length\nH=2500.0;#actual height\nhb=800#height B\nha=500.0;#height A\nya=82.25;\nyb=-41.15;\ns=220.0;\n\n#calculation\ntheta=s-180;\nya1=xa*sin(theta*pi/180)+ya*cos(theta*pi/180)+f*tan(t); # ya'\nxa1=xa*cos(theta*pi/180)-ya*sin(theta*pi/180); #xa'\nxb1=xb*cos(theta*pi/180)-yb*sin(theta*pi/180); #xb'\nyb1=xb*sin(theta*pi/180)+yb*cos(theta*pi/180)+f*tan(t); # yb'\nXa=xa1*(H-ha)/(f/cos(t)-ya1*sin(t));\nXb=xb1*(H-hb)/(f/cos(t)-yb1*sin(t));\nYa=ya1*cos(t)*(H-ha)/(f/cos(t)-ya1*sin(t));\nYb=yb1*cos(t)*(H-hb)/(f/cos(t)-yb1*sin(t));\nAB=sqrt((Xb-Xa)**2+(Yb-Ya)**2);\n\n#result\nprint \"distance of AB in m\",round(AB,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "distance of AB in m 2307.753\n -100.708756369 78.25 82.25 -41.15 0.642787609687\n" + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.14,Page 324" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nri=95.0;\nf=152.4;\nt=3*pi/180;\nl=50*pi/180;\n\n#calculation\ndt=ri**2*sin(t)*cos(l)**2/(f-ri*sin(t)*cos(l));\n \n#result\nprint \"tilt displacement of the image in mm\",round(dt,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "tilt displacement of the image in mm 1.31\n" + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.15,Page 351" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nd=230.0;#square side\nf=152.4#focal length\npl=0.6;\n\n#calculation\nk=(1-pl)*d/f;\nV=k/0.15;\n\n#result\nprint \"vertical exaggeration is\",round(V,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "vertical exaggeration is 4.02\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.16,Page 360" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nxa=51.23;\nxb=91.48;\nya=48.33;\nyb=-51.63;\nf=152.4;#focal length\nB=425.0;#actual height\nhb=842.86#height B\nha=820.97;#height A\nr1=10.42;\nr2=9.67;\nb1=89.12;#b'\nb=89.43;\nra=11.62;\nrb=14.53;\n\n#calculation\nC=0.5*((b1-r1)+(b-r2))\npa=C+ra;\npb=C+rb;\nXa=B*xa/pa;\nXb=xb*B/pb;\nYa=ya*B/pa;\nYb=yb*B/pb;\nAB=sqrt((Xb-Xa)**2+(Yb-Ya)**2);\n\n#result\nprint \"distance of AB in m\",round(AB,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "distance of AB in m 492.28\n" + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.17,Page 363" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nC=79.0;\nra=11.42;#elevarion in image\nrb=15.65;#elevarion in image\nhb=651;#height of B\nH=1500;#height\n\n#calculation\ndelp=ra-rb;#pa=ra+c and pb=rb+c so ra-rb=pa-pb\npa=ra+C;\nha=hb+delp/pa*(H-hb);\n\n#result\nprint \"height of A in m\",round(ha,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "height of A in m 611.282\n" + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.18,Page 364" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nB=741.0;\nf=152.4;\npa=94.32;\nha=325;\n\n#calculation\nH=ha+B*f/pa;\n\n#result\nprint \"height in m\",round(H,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "height in m 1522.29\n" + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.19,Page 364" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nH=1632.0;\nf=152.4;#focal length\npa=82.75;\nha=283;\n\n#calculation\nB=pa/f*(H-ha);\n\n#result\nprint \"width of air base in m\",round(B,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "width of air base in m 732.479\n" + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.20,Page 377" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nalpha=50.0;\nbeta=46.0;\nf=300.0;\nxa=24.0;\nxb=30.0;\nbeta=46.0;\n\n#calculation\ndela=xa/f;\ndelb=xb/f;\nA=alpha+dela*180/pi; #angle A\nB=beta-delb*180/pi; #angle B\nD=180-A-B;\nAD=1300.0*sin(B*pi/180)/sin(D*pi/180);\nY=6/(sqrt(xa**2+f**2))*AD;\nRD=60.12+Y;\n\n#result\nprint \"distance of AD in m\",round(AD,2)\nprint \"RL of D in m\",round(RD,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "distance of AD in m 843.34\nRL of D in m 76.93\n" + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.21,Page 378" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nf=152.4;\nL=120;#length\nx1=40.0;\nx2=-90.0;\n\n#calculation\nX=f*L/(x1-x2);\nY=L*x1/(x1-x2);\nh=X*(30-20)/f;\n\n#result\nprint \"the coordinates of D in m is X=\",round(X,2),\"Y =\",round(Y,2);\nprint \"elevation of D in m\",round(h,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "the coordinates of D in m is X= 140.68 Y = 36.92\nelevation of D in m 9.23\n" + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.22,Page 380" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nf=150.4;\nxc=-32.43;\nxd=9.52;\n\n#calculation\nthc=atan(xc/f);\nthd=atan(xd/f);\nth=thd-thc;\nth=th*180/pi;\nAz=325+15.0/60+th;\nAz=deg_to_dms(Az);\n\n#result\nprint \"Azimuth of D in deg,min,sec respectively\",Az", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Azimuth of D in deg,min,sec respectively [341, 2, 23.9]\n" + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.23,Page 382" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nBC=66.0;\nAC=81.6;\nxb=3.0;\nya=1.25;\nxa=3.3;\ntheta=23+43.0/60;\n\n#calcualtions\nf=(xa+xb)/2/tan(theta*pi/180)+sqrt((xa+xb)**2/4/(tan(theta*pi/180))**2+xa*xb);\naa=atan(ya/sqrt(xa**2+f**2))\nVa=AC*tan(aa);\nab=atan(-1.87/sqrt(xa**2+f**2));\nVb=-BC*tan(ab);\n\n#result\nprint \"focal length in cm\",round(f,2)\nprint \"horizontal distance in m\",round(Vb+Va,2)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "focal length in cm 15.0\nhorizontal distance in m 14.68\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.24,Page 383" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nEab=300.0;#average elevation\nf=152.4;\nxa=28.4;\nxb=-22.5;\nya=24.5;\nyb=38.4;\nHa=2322.0;\nha=400.0;\nhb=200.0;\nab=61.05;\nAB=810;\n\n#calculation\nHa=300+AB/ab*f;\nXa=round((Ha-ha)/f*xa,2);\nXb=round((Ha-ha)/f*xb,2);\nYa=round((Ha-hb)/f*ya,2);\nYb=round((Ha-hb)/f*yb,2);\nAB=sqrt((Xa-Xb)**2+(Ya-Yb)**2);\nH=300+810/AB*(Ha-Eab);\nXa=(H-ha)/f*xa;\nXb=(H-ha)/f*xb;\nYa=(H-hb)/f*ya;\nYb=(H-hb)/f*yb;\nAB1=sqrt((Xa-Xb)**2+(Ya-Yb)**2);\n\n#result\nprint \"Xb is calculated wrong in the book that resulted in the error\"\nprint \"length AB in m\", round(AB,2)\nprint \"corrected length AB in m\", round(AB1,2)\nprint \"flying height in m\",round(H,3)\n\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "length AB in m 670.47\ncorrected length AB in m 816.12\nflying height in m 2742.807\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.25,Page 386" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nAB=300.0;\nab=102.4;\nf=152.4;#focal length\nhab=320.0;\nd=7.8;\nr=75.4;\n\nH=hab+AB/ab*f;\nh=d*H/r;\n\n#result\nprint \"height difference in m\",round(h,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "height difference in m 79.29\n" + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.26,Page 386" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nf=152.4;\nb=74.25;\nht=100.0;\nH=700.0;#flying height\n\n#calculation\nB=b*H/f;\npb=f*B/H;\npt=f*B/(H-ht);\ndelp=pt-pb;\nht=delp/pt*(H);\n\n#result\nprint \"error due to parallax in mm\",round(delp,2)\nprint \"height of chimney in m\",round(ht,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "error due to parallax in mm 12.38\nheight of chimney in m 100.0\n" + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.27,Page 387" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\ndef deg_to_dms(deg):\n d = int(deg)\n md = abs(deg - d) * 60\n m = int(md)\n sd = (md - m) * 60\n sd=round(sd,2)\n return [d, m, sd]\nB=180.0;\nf=120.0;\npa=54.32\npb=46.35;\n\n#calculation\ndelH=B*f/pa/pb*(pa-pb);\n\n#result\nprint \"height difference in m\",round(delH,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "height difference in m 68.38\n" + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.28,Page 387" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nfrom math import pi,sqrt,acos,asin,atan,cos,sin,tan\nL1=30000.0;\npl=0.6;\nk=12000.0*0.2;#=l/S and w/S\npw=0.3;\nW1=24000;\n\n#calculation\nN=round(((L1/((1-pl)*k)+1)+1))*round(((W1/((1-pw)*k)+1)+1));\nNf=N/33-1;#flight strips\ngd=(1-pl)*k;#grounf distance\nI=gd/(200e3)*60.0*60.0;#exposure interval\nad=round(I)/60.0/60*200e3;#actual distance\n\n#result\nprint \"no. of photographs taken\",round(N)\nprint \"no. of flight strips\",Nf\nprint \"ground distance in m\",round(gd,2)\nprint \"exposure interval in s\", round(I)\nprint \"actual distance in m\",round(ad,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "no. of photographs taken 528.0\nno. of flight strips 15.0\nground distance in m 960.0\nexposure interval in s 17.0\nactual distance in m 944.44\n" + } + ], + "prompt_number": 35 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Surveying_Volume_3/README.txt b/Surveying_Volume_3/README.txt new file mode 100755 index 00000000..2c27e6c8 --- /dev/null +++ b/Surveying_Volume_3/README.txt @@ -0,0 +1,10 @@ +Contributed By: Tarun Kumar Das +Course: btech +College/Institute/Organization: College of Engineering +Department/Designation: Industrial Engineering +Book Title: Surveying Volume 3 +Author: A K Arora +Publisher: Lakshmi Publications, Chennai +Year of publication: 2011 +Isbn: 9788189401276 +Edition: 1
\ No newline at end of file diff --git a/Surveying_Volume_3/screenshots/chapter1.png b/Surveying_Volume_3/screenshots/chapter1.png Binary files differnew file mode 100755 index 00000000..14135a76 --- /dev/null +++ b/Surveying_Volume_3/screenshots/chapter1.png diff --git a/Surveying_Volume_3/screenshots/chapter4.png b/Surveying_Volume_3/screenshots/chapter4.png Binary files differnew file mode 100755 index 00000000..ad0a4bc9 --- /dev/null +++ b/Surveying_Volume_3/screenshots/chapter4.png diff --git a/Surveying_Volume_3/screenshots/chapter5.png b/Surveying_Volume_3/screenshots/chapter5.png Binary files differnew file mode 100755 index 00000000..bd0c856e --- /dev/null +++ b/Surveying_Volume_3/screenshots/chapter5.png |
