def test_polar_point3d_3():
ppp = Point("1", 10, 20, 30)
ooo = PolarObservation("1", "station", Angle("0", "DMS"), None, None, 1.0)
sss = Station(ppp, ooo)
oo2 = PolarObservation("2", None, Angle("90", "DMS"), Angle("90", "DMS"), Distance(20, "SD"), 3.5)
pp2 = Calculation.polarpoint(sss, oo2)[0]
assert compare(pp2, Point('2', 30.0000, 20.0000, 27.5))
def test_observation():
# pnt = [Point('1', 1000, 2000, 50), Point('2', 1500, 2000, 60)]
obs = [
PolarObservation('1', 'station', None, None, None, 1.54),
PolarObservation('2', None, Angle(60.9345, 'GON'),
Angle(89.855615, 'DEG'), Distance(501.105, 'SD'),
1.80)
]
assert compare(obs[1].horiz_dist(), 501.103)
def test_polar_points():
# polar points
p101pol = Point("101", 13456.25, 12569.75)
p201pol = Point("201", 13102.13, 11990.13)
p202pol = Point("202", 13569.11, 12788.66)
p203pol = Point("203", 13861.23, 12001.54)
o101pol = PolarObservation('101', "station")
s101pol = Station(p101pol, o101pol)
o201pol = PolarObservation("201", None, Angle("112-15-15", "DMS"))
o202pol = PolarObservation("202", None, Angle("288-06-30", "DMS"))
o203pol = PolarObservation("203", None, Angle("45-21-12", "DMS"))
o9pol = PolarObservation("9", None, Angle("145-10-16", "DMS"), None, Distance(206.17, "HD"))
o10pol = PolarObservation("10", None, Angle("201-30-47", "DMS"), None, Distance(219.38, "HD"))
z101pol = Calculation.orientation(s101pol, [[p201pol, o201pol], [p202pol, o202pol], [p203pol, o203pol]])[0]
assert compare(z101pol.get_angle('DMS'), '99-10-05')
s101pol.o.hz = z101pol
p9pol = Calculation.polarpoint(s101pol, o9pol)[0]
p10pol = Calculation.polarpoint(s101pol, o10pol)[0]
assert compare(p9pol, Point("9", 13270.4141, 12480.4691))
assert compare(p10pol, Point("10", 13267.5785, 12681.6903))
pA1pol = Point("A1", 153.867, 456.430)
pT1pol = Point("T1", -237.865, -297.772)
pT2pol = Point("T2", -1549.927, 669.6126)
pT3pol = Point("T3", 1203.064, -220.0314)
oA1pol = PolarObservation('A1', "station")
sA1pol = Station(pA1pol, oA1pol)
oT1pol = PolarObservation("T1", None, Angle("73-02-35", "DMS"))
oT2pol = PolarObservation("T2", None, Angle("142-43-39", "DMS"))
oT3pol = PolarObservation("T3", None, Angle("348-24-26", "DMS"))
oP1pol = PolarObservation("P1", None, Angle("112-43-47", "DMS"), None, Distance(673.699, "HD"))
oP2pol = PolarObservation("P2", None, Angle("84-0-44", "DMS"), None, Distance(788.105, "HD"))
zA1pol = Calculation.orientation(sA1pol, [[pT1pol, oT1pol], [pT2pol, oT2pol], [pT3pol, oT3pol]])[0]
assert compare(zA1pol.get_angle('DMS'), '134-24-16')
sA1pol.o.hz = zA1pol
pP1pol = Calculation.polarpoint(sA1pol, oP1pol)[0]
pP2pol = Calculation.polarpoint(sA1pol, oP2pol)[0]
assert compare(pP1pol, Point("P1", -466.8905, 194.6468))
assert compare(pP2pol, Point("P2", -335.8415, -161.0610))
def test_traverse2():
# closed at both ends and known bearings at both ends
pKtra = Point("K", 599767.21, 148946.70)
pVtra = Point("V", 599733.75, 149831.76)
p1015tra = Point("1015", 598642.17, 148436.26)
p1016tra = Point("1016", 600136.60, 148588.85)
p1017tra = Point("1017", 600264.30, 149325.79)
p1018tra = Point("1018", 598258.90, 149496.78)
p1019tra = Point("1019", 600092.33, 150676.80)
oKtra = PolarObservation('K', "station")
sKtra = Station(pKtra, oKtra)
oK_1017tra = PolarObservation("1017", None, Angle("61-28-18", "DMS"))
oK_1016tra = PolarObservation("1016", None, Angle("142-53-28", "DMS"))
oK_1015tra = PolarObservation("1015", None, Angle("254-23-32", "DMS"))
oK_1tra = PolarObservation("1", None, Angle("17-14-18", "DMS"), None, Distance(139.82, "HD"))
sKtra.o.hz = Calculation.orientation(sKtra,
[[p1017tra, oK_1017tra], [p1016tra, oK_1016tra], [p1015tra, oK_1015tra]])[0]
assert compare(sKtra.o.hz.get_angle('DMS'), '351-12-05')
s1tra = Station(None, PolarObservation('1', "station"))
o1_Ktra = PolarObservation("K", None, Angle("79-28-20", "DMS"), None, Distance(139.85, "HD"))
o1_2tra = PolarObservation("2", None, Angle("236-13-46", "DMS"), None, Distance(269.32, "HD"))
# o1_501tra = PolarObservation("501", None, Angle("204-58-10", "DMS"), None, Distance(59.12, "HD"))
s2tra = Station(None, PolarObservation('2', "station"))
o2_3tra = PolarObservation("3", None, Angle("82-18-45", "DMS"), None, Distance(169.40, "HD"))
o2_1tra = PolarObservation("1", None, Angle("217-58-34", "DMS"), None, Distance(269.36, "HD"))
s3tra = Station(None, PolarObservation('3', "station"))
o3_2tra = PolarObservation("2", None, Angle("262-18-44", "DMS"), None, Distance(169.45, "HD"))
o3_Vtra = PolarObservation("V", None, Angle("41-18-10", "DMS"), None, Distance(345.90, "HD"))
# o3_502tra = PolarObservation("502", None, Angle("344-28-25", "DMS"), None, Distance(55.46, "HD"))
oVtra = PolarObservation('V', "station")
sVtra = Station(pVtra, oVtra)
oV_3tra = PolarObservation("3", None, Angle("257-44-08", "DMS"), None, Distance(345.94, "HD"))
oV_1018tra = PolarObservation("1018", None, Angle("346-24-11", "DMS"))
oV_1019tra = PolarObservation("1019", None, Angle("112-12-06", "DMS"))
oV_1017tra = PolarObservation("1017", None, Angle("222-50-58", "DMS"))
sVtra.o.hz = Calculation.orientation(sVtra,
[[p1018tra, oV_1018tra], [p1019tra, oV_1019tra], [p1017tra, oV_1017tra]])[0]
assert compare(sVtra.o.hz.get_angle('DMS'), '270-47-46')
plist = Calculation.traverse([[sKtra, None, oK_1tra], [s1tra, o1_Ktra, o1_2tra], [s2tra, o2_1tra, o2_3tra],
[s3tra, o3_2tra, o3_Vtra], [sVtra, oV_3tra, None]])
assert compare(plist[0], Point('1', 599787.74865, 149085.0079))
assert compare(plist[1], Point('2', 599718.9691, 149345.3886))
assert compare(plist[2], Point('3', 599802.5094, 149492.7786))
def test_traverse1():
# closed at one end and known bearings at one end (free traverse)
p5241otra = Point("5241", 646414.44, 211712.77)
p5245otra = Point("5245", 646938.71, 212635.92)
p5246otra = Point("5246", 646380.61, 212793.97)
p5247otra = Point("5247", 646381.14, 212476.49)
o5247otra = PolarObservation('5247', "station")
s5247otra = Station(p5247otra, o5247otra)
o5241otra = PolarObservation("5241", None, Angle("245-23-41", "DMS"))
o5245otra = PolarObservation("5245", None, Angle("141-56-11", "DMS"))
o5246otra = PolarObservation("5246", None, Angle("67-47-14", "DMS"))
s5247otra.o.hz = Calculation.orientation(s5247otra,
[[p5241otra, o5241otra],
[p5245otra, o5245otra], [p5246otra, o5246otra]])[0]
assert compare(s5247otra.o.hz.get_angle('DMS'), '292-06-34')
o5247_111otra = PolarObservation("111", None, Angle("241-26-57", "DMS"), None, Distance(123.42, "HD"))
s111otra = Station(None, PolarObservation('111', "station"))
o111_5247otra = PolarObservation("5247", None, Angle("225-39-00", "DMS"))
o111_112otra = PolarObservation("112", None, Angle("92-38-43", "DMS"), None, Distance(142.81, "HD"))
s112otra = Station(None, PolarObservation('112', "station"))
o112_111otra = PolarObservation("111", None, Angle("227-16-34", "DMS"))
o112_113otra = PolarObservation("113", None, Angle("69-16-28", "DMS"), None, Distance(253.25, "HD"))
s113otra = Station(None, PolarObservation('113', "station"))
o113_112otra = PolarObservation("112", None, Angle("102-56-44", "DMS"))
o113_114otra = PolarObservation("114", None, Angle("205-46-21", "DMS"), None, Distance(214.53, "HD"))
s114otra = Station(None, PolarObservation('114', "station"))
o114_113otra = PolarObservation("113", None, Angle("104-23-11", "DMS"))
o114_115otra = PolarObservation("115", None, Angle("305-54-29", "DMS"), None, Distance(234.23, "HD"))
s115otra = Station(None, PolarObservation("115", 'station'))
plist = Calculation.traverse([[s5247otra, None, o5247_111otra], [s111otra, o111_5247otra, o111_112otra],
[s112otra, o112_111otra, o112_113otra], [s113otra, o113_112otra, o113_114otra],
[s114otra, o114_113otra, o114_115otra], [s115otra, None, None]])
assert compare(plist[0], Point('111', 646394.9860, 212353.8491))
assert compare(plist[1], Point('112', 646302.1362, 212245.3429))
assert compare(plist[2], Point('113', 646077.3941, 212128.60997))
assert compare(plist[3], Point('114', 646131.5459, 211921.02697))
assert compare(plist[4], Point('115', 646103.4028, 211688.4938))
def test_traverse4():
# loop traverse
pA5tra = Point("A5", 646333.5695, 276616.4171)
pA4tra = Point("A4", 646284.6886, 276659.2165)
oA5tra = PolarObservation('A5', "station")
sA5tra = Station(pA5tra, oA5tra)
oA5_A4tra = PolarObservation("A4", None, Angle("311-12-21", "DMS"))
oA5_A6tra = PolarObservation("A6", None, Angle("136-36-49", "DMS"), Angle("100-13-22", "DMS"),
Distance(58.405, "SD"))
# oA5_A11tra = PolarObservation("A11", None, Angle("146-26-58", "DMS"), Angle("90-00-00", "DMS"),
# Distance(306.234, "SD"))
sA5tra.o.hz = Calculation.orientation(sA5tra, [[pA4tra, oA5_A4tra]])[0]
assert compare(sA5tra.o.hz.get_angle('DMS'), '359-59-57')
sA6tra = Station(None, PolarObservation('A6', "station"))
oA6_A5tra = PolarObservation("A5", None, Angle("316-36-49", "DMS"), Angle("79-56-53", "DMS"),
Distance(58.378, "SD"))
oA6_A7tra = PolarObservation("A7", None, Angle("136-43-52", "DMS"), Angle("101-55-08", "DMS"),
Distance(43.917, "SD"))
sA7tra = Station(None, PolarObservation('A7', "station"))
oA7_A6tra = PolarObservation("A6", None, Angle("316-43-52", "DMS"), Angle("78-21-53", "DMS"),
Distance(43.878, "SD"))
oA7_A8tra = PolarObservation("A8", None, Angle("141-48-30", "DMS"), Angle("102-29-18", "DMS"),
Distance(48.459, "SD"))
sA8tra = Station(None, PolarObservation('A8', "station"))
oA8_A7tra = PolarObservation("A7", None, Angle("321-48-30", "DMS"), Angle("77-49-44", "DMS"),
Distance(48.401, "SD"))
oA8_A9tra = PolarObservation("A9", None, Angle("153-25-18", "DMS"), Angle("100-01-16", "DMS"),
Distance(47.098, "SD"))
sA9tra = Station(None, PolarObservation('A9', "station"))
oA9_A8tra = PolarObservation("A8", None, Angle("333-25-18", "DMS"), Angle("80-17-11", "DMS"),
Distance(47.040, "SD"))
oA9_A10tra = PolarObservation("A10", None, Angle("153-59-32", "DMS"), Angle("97-46-19", "DMS"),
Distance(58.077, "SD"))
sA10tra = Station(None, PolarObservation('A10', "station"))
oA10_A9tra = PolarObservation("A9", None, Angle("333-59-32", "DMS"), Angle("82-27-53", "DMS"),
Distance(58.045, "SD"))
oA10_A11tra = PolarObservation("A11", None, Angle("154-05-41", "DMS"), Angle("97-06-32", "DMS"),
Distance(58.188, "SD"))
sA11tra = Station(None, PolarObservation('A11', "station"))
oA11_A10tra = PolarObservation("A10", None, Angle("334-05-41", "DMS"), Angle("83-09-29", "DMS"),
Distance(58.151, "SD"))
# oA11_A5tra = PolarObservation("A5", None, Angle("326-26-13", "DMS"), Angle("90-00-00", "DMS"),
# Distance(306.233, "SD"))
plist = Calculation.traverse(
[[sA5tra, None, oA5_A6tra], [sA6tra, oA6_A5tra, oA6_A7tra], [sA7tra, oA7_A6tra, oA7_A8tra],
[sA8tra, oA8_A7tra, oA8_A9tra], [sA9tra, oA9_A8tra, oA9_A10tra], [sA10tra, oA10_A9tra, oA10_A11tra],
[sA11tra, oA11_A10tra, None]])
assert compare(plist[0], Point('A6', 646373.0539, 276574.6449))
assert compare(plist[1], Point('A7', 646402.5093, 276543.3546))
assert compare(plist[2], Point('A8', 646431.7631, 276506.1698))
assert compare(plist[3], Point('A9', 646452.5118, 276464.6981))
assert compare(plist[4], Point('A10', 646477.7451, 276412.9820))
assert compare(plist[5], Point('A11', 646502.9712, 276361.0454))
def test_traverse3():
# closed at both ends and known bearing at one end
pA5tra = Point("A5", 646333.5695, 276616.4171)
pA11tra = Point("A11", 646502.8710, 276361.2386)
pA4tra = Point("A4", 646284.6886, 276659.2165)
oA5tra = PolarObservation('A5', "station")
sA5tra = Station(pA5tra, oA5tra)
oA5_A4tra = PolarObservation("A4", None, Angle("311-12-21", "DMS"))
oA5_A6tra = PolarObservation("A6", None, Angle("136-36-49", "DMS"), Angle("100-13-22", "DMS"),
Distance(58.405, "SD"))
sA5tra.o.hz = Calculation.orientation(sA5tra, [[pA4tra, oA5_A4tra]])[0]
assert compare(sA5tra.o.hz.get_angle('DMS'), '359-59-57')
sA6tra = Station(None, PolarObservation('A6', "station"))
oA6_A5tra = PolarObservation("A5", None, Angle("316-36-49", "DMS"), Angle("79-56-53", "DMS"),
Distance(58.378, "SD"))
oA6_A7tra = PolarObservation("A7", None, Angle("136-43-52", "DMS"), Angle("101-55-08", "DMS"),
Distance(43.917, "SD"))
sA7tra = Station(None, PolarObservation('A7', "station"))
oA7_A6tra = PolarObservation("A6", None, Angle("316-43-52", "DMS"), Angle("78-21-53", "DMS"),
Distance(43.878, "SD"))
oA7_A8tra = PolarObservation("A8", None, Angle("141-48-30", "DMS"), Angle("102-29-18", "DMS"),
Distance(48.459, "SD"))
sA8tra = Station(None, PolarObservation('A8', "station"))
oA8_A7tra = PolarObservation("A7", None, Angle("321-48-30", "DMS"), Angle("77-49-44", "DMS"),
Distance(48.401, "SD"))
oA8_A9tra = PolarObservation("A9", None, Angle("153-25-18", "DMS"), Angle("100-01-16", "DMS"),
Distance(47.098, "SD"))
sA9tra = Station(None, PolarObservation('A9', "station"))
oA9_A8tra = PolarObservation("A8", None, Angle("333-25-18", "DMS"), Angle("80-17-11", "DMS"),
Distance(47.040, "SD"))
oA9_A10tra = PolarObservation("A10", None, Angle("153-59-32", "DMS"), Angle("97-46-19", "DMS"),
Distance(58.077, "SD"))
sA10tra = Station(None, PolarObservation('A10', "station"))
oA10_A9tra = PolarObservation("A9", None, Angle("333-59-32", "DMS"), Angle("82-27-53", "DMS"),
Distance(58.045, "SD"))
oA10_A11tra = PolarObservation("A11", None, Angle("154-05-41", "DMS"), Angle("97-06-32", "DMS"),
Distance(58.188, "SD"))
sA11tra = Station(pA11tra, PolarObservation('A11', "station"))
oA11_A10tra = PolarObservation("A10", None, Angle("334-05-41", "DMS"), Angle("83-09-29", "DMS"),
Distance(58.151, "SD"))
plist = Calculation.traverse(
[[sA5tra, None, oA5_A6tra], [sA6tra, oA6_A5tra, oA6_A7tra], [sA7tra, oA7_A6tra, oA7_A8tra],
[sA8tra, oA8_A7tra, oA8_A9tra], [sA9tra, oA9_A8tra, oA9_A10tra], [sA10tra, oA10_A9tra, oA10_A11tra],
[sA11tra, oA11_A10tra, None]])
assert compare(plist[0], Point('A6', 646373.0353, 276574.6808))
assert compare(plist[1], Point('A7', 646402.4768, 276543.4174))
assert compare(plist[2], Point('A8', 646431.7153, 276506.2621))
assert compare(plist[3], Point('A9', 646452.4490, 276464.8193))
assert compare(plist[4], Point('A10', 646477.6637, 276413.1392))