diff --git a/Karto/README.md b/Karto/README.md index 8b4629d..5efca94 100644 --- a/Karto/README.md +++ b/Karto/README.md @@ -28,3 +28,6 @@ useful commands: regions: - uni(on) / in(tersection) / su(btract) + +Questions: +- OFFSET in assignment3 gives quite different offsets for Fig. 2 diff --git a/Karto/assignment3/KTZ004_2019_03_ats.py b/Karto/assignment3/KTZ004_2019_03_ats.py index fc67baa..6d7bf3d 100755 --- a/Karto/assignment3/KTZ004_2019_03_ats.py +++ b/Karto/assignment3/KTZ004_2019_03_ats.py @@ -229,25 +229,25 @@ Visas toliau ekspluatuojamas plotas 22960.319 === Iskasto grunto kiekio skaiciavimas === Iskasos Nr.1 plotas (P1) pagal virsutine riba (0.1 m2 tikslumu) -******* +%.1f""" % P1_virsutine + """ Iskasos Nr.1 plotas (P2) pagal slaito apacia (0.1 m2 tikslumu) -******* +%.1f""" % P2_apatine + """ Iskasos Nr.1 turis (V1) paskaiciuotas pagal 2 formule (0.1 m3 tikslumu) -******* +%.1f""" % fig1_v1 + """ Iskasos Nr.1 turis (V2) paskaiciuotas pagal 3 formule (0.1 m3 tikslumu) -******* +%.1f""" % fig1_v2 + """ Galutinis iskasos Nr.1 turis (V) paskaiciuotas pagal 4 formule (1 m3 tikslumu) -******* +%.f""" % fig1_v + """ Iskasos Nr.2 plotas (P3) pagal virsutine riba (0.1 m2 tikslumu) -******* +%.1f""" % P3_virsutine + """ Iskasos Nr.2 plotas (P4) pagal slaito apacia (0.1 m2 tikslumu) -******* +%.1f""" % P4_apatine + """ Iskasos Nr.2 turis (V1) paskaiciuotas pagal 2 formule (0.1 m3 tikslumu) -******* +%.1f""" % fig2_v1 + """ Iskasos Nr.2 turis (V2) paskaiciuotas pagal 3 formule (0.1 m3 tikslumu) -******* +%.1f""" % fig2_v2 + """ Galutinis iskasos Nr.2 turis (V) paskaiciuotas pagal 4 formule (1 m3 tikslumu) -******* +%.0f""" % fig2_v + """ === Apskaiciuotos kainos (0.01Eu tikslumu) === Visa rekultivacijos (R) kaina %.2f""" % (rek_kaina * rek_plotas) + """ diff --git a/Karto/assignment3/KTZ004_2019_03_brez.dwg b/Karto/assignment3/KTZ004_2019_03_brez.dwg index 3438809..62298cf 100644 Binary files a/Karto/assignment3/KTZ004_2019_03_brez.dwg and b/Karto/assignment3/KTZ004_2019_03_brez.dwg differ diff --git a/Karto/assignment3/measure.py b/Karto/assignment3/measure.py index 64471d2..dc53d5c 100755 --- a/Karto/assignment3/measure.py +++ b/Karto/assignment3/measure.py @@ -2,7 +2,7 @@ from decimal import Decimal as Dec from deg import Deg -from math import tan, pi +from math import tan, pi, sqrt def fmt(deg): if deg > 180: @@ -130,6 +130,41 @@ K2K7 = Deg.guess('69.4792556').frac + R K2K8 = Deg.guess('130-50-06.1').frac + L K2K9 = Deg.guess('84-08-54.3').frac + K +H1 = 4.82 +SK1 = 41.08 +A1 = H1 * tan(SK1 * pi / 180) + +H2 = 20.38 +SK2 = 17.69 +A2 = H2 * tan(SK2 * pi / 180) + + +# Area +P1_virsutine = Dec('63425.860') +P2_apatine = Dec('56139.919') + +P3_virsutine = Dec('59232.387') +P4_apatine = Dec('50411.889') + +def f2_v1(h, pv, pa): + return h/3*(float(pv+pa)+sqrt(pv*pa)) + +def f3_v2(h, pv, pa): + return h/2*float(pv+pa) + +def f4_v(v1, v2): + return (v1 + v2) / 2 + +fig1_v1 = f2_v1(H1, P1_virsutine, P2_apatine) +fig1_v2 = f3_v2(H1, P1_virsutine, P2_apatine) +fig1_v = f4_v(fig1_v1, fig1_v2) + +fig2_v1 = f2_v1(H2, P3_virsutine, P4_apatine) +fig2_v2 = f3_v2(H2, P3_virsutine, P4_apatine) +fig2_v = f4_v(fig2_v1, fig2_v2) + +# Kainos + rek_kaina = Dec('871.29') rek_plotas = Dec('97895.626') @@ -145,14 +180,6 @@ visa_kaina = ( ekspl_kaina * ekspl_plotas ) -H1 = 4.82 -SK1 = 41.08 -A1 = H1 * tan(SK1 * pi / 180) - -H2 = 20.38 -SK2 = 17.69 -A2 = H2 * tan(SK2 * pi / 180) - if __name__ == '__main__': print("""Sklypu pradiniu tasku (1 ir 15) koordinates: Atkarpos: