def testDPA(self):
with open("CandidateESCLoc_East.txt", "r") as myfile:
data = myfile.readlines()
esc_loc_x = [x / 1000.0 for x in eval(data[0])]
esc_loc_y = [y / 1000.0 for y in eval(data[1])]
centers = zip(esc_loc_x, esc_loc_y)
with open("DPA1.txt", "r") as myfile:
data = myfile.readlines()
dpa_loc_x = [x / 1000.0 for x in eval(data[0])]
dpa_loc_y = [y / 1000.0 for y in eval(data[1])]
dpa_loc = zip(dpa_loc_x, dpa_loc_y)
poly = Polygon(dpa_loc)
coverage_file = "detection-coverage/ITMDetectionCoverage_90deg.txt"
min_ctr_dist = 0
cover = antennacover.min_antenna_area_cover_greedy(
centers, poly, coverage_file, min_center_distance=min_ctr_dist)
testName = "DPA1"
printcover.printAntennaCover(testName, poly, centers, cover,
coverage_file, min_ctr_dist)
def testEstuaryAnneal(self):
"""
Test simulated annealing.
"""
interference_contour = [(20, 55), (35, 65), (40, 60), (45, 65),
(50, 55)]
possible_centers = [(20, 46), (25, 30), (30, 20), (40, 15), (50, 30),
(60, 50)]
antennacover.NDIVISIONS = 200
min_ctr_dist = 0
coverage_file = "detection-coverage/ITMDetectionCoverage_60deg.json"
testName = "Estuary"
poly = excessarea.generate_bounding_polygon(possible_centers,
interference_contour)
cover = antennacover.min_antenna_area_cover_greedy(
possible_centers,
poly,
coverage_file,
min_center_distance=min_ctr_dist)
printcover.printAntennaCover(testName, poly, possible_centers, cover,
coverage_file, min_ctr_dist)
annealr = simannealer.SimAnneal(poly, coverage_file, cover)
annealr.anneal()
testName = "EstuaryAnneal"
improved_cover = annealr.get_result()
printcover.printAntennaCover(testName, poly, possible_centers,
improved_cover, coverage_file,
min_ctr_dist)
def testEastCoastAnneal(self):
with open("InterfContour_EastCoast.txt", "r") as myfile:
data = myfile.readlines()
esc_loc_x = [x / 1000.0 for x in eval(data[0])]
esc_loc_y = [x / 1000.0 for x in eval(data[1])]
ship_loc_x = [x / 1000.0 for x in eval(data[2])]
ship_loc_y = [x / 1000.0 for x in eval(data[3])]
centers = []
for i in range(0, len(esc_loc_x)):
center = (esc_loc_x[i], esc_loc_y[i])
centers.append(center)
interference_contour = []
for i in range(0, len(ship_loc_x)):
p = (ship_loc_x[i], ship_loc_y[i])
interference_contour.append(p)
testName = "EastCoast"
min_ctr_dist = 0
coverage_file = "DetectionCoverage_60deg.txt"
poly = excessarea.generate_bounding_polygon(centers,
interference_contour)
cover = antennacover.min_antenna_area_cover_greedy(
centers, poly, coverage_file, min_center_distance=min_ctr_dist)
printcover.printAntennaCover(testName, poly, centers, cover,
coverage_file, min_ctr_dist)
testName = "EastCoastAnneal"
annealr = simannealer.SimAnneal(poly, coverage_file, cover)
annealr.anneal()
improved_cover = annealr.get_result()
printcover.printAntennaCover(testName, poly, centers, improved_cover,
coverage_file, min_ctr_dist)
def testVBAnneal(self):
# Convert all units to Km
esc_loc_x = [
x / 1000.0 for x in [
1771380, 1769310, 1769790, 1768380, 1767390, 1764690, 1762020,
1759920, 1753110, 1741950, 1752210, 1757010, 1761870, 1768230,
1772820, 1777110, 1781610, 1786920, 1793220
]
]
esc_loc_y = [
y / 1000.0 for y in [
1827030, 1817070, 1806990, 1797090, 1787100, 1776840, 1767270,
1756950, 1746690, 1735050, 1727220, 1717290, 1707360, 1697370,
1687320, 1677450, 1667400, 1657350, 1647360
]
]
ship_loc_x = [
x / 1000.0 for x in [
1847012, 1844913, 1845660, 1834150, 1823280, 1811715, 1807512,
1806671, 1810710, 1807769, 1817910, 1822503, 1827218, 1823623,
1828432, 1842183, 1846928, 1852378, 1858591
]
]
ship_loc_y = [
y / 1000.0 for y in [
1843636, 1833617, 1823583, 1811442, 1799284, 1787072, 1777140,
1767066, 1759078, 1749183, 1741311, 1731358, 1721401, 1709309,
1699318, 1691518, 1681523, 1671542, 1661589
]
]
centers = []
for i in range(0, len(esc_loc_x)):
center = (esc_loc_x[i], esc_loc_y[i])
centers.append(center)
interference_contour = []
for i in range(0, len(ship_loc_x)):
p = (ship_loc_x[i], ship_loc_y[i])
interference_contour.append(p)
testName = "VirginiaBeach"
min_ctr_dist = 60
coverage_file = "detection-coverage/ITMDetectionCoverage_60deg.json"
poly = excessarea.generate_bounding_polygon(centers,
interference_contour)
cover = antennacover.min_antenna_area_cover_greedy(
centers, poly, coverage_file, min_center_distance=min_ctr_dist)
printcover.printAntennaCover(testName, poly, centers, cover,
coverage_file, min_ctr_dist)
testName = "VirginiaBeachAnneal"
annealr = simannealer.SimAnneal(poly, coverage_file, cover)
annealr.anneal()
improved_cover = annealr.get_result()
printcover.printAntennaCover(testName, poly, centers, improved_cover,
coverage_file, min_ctr_dist)
def testSFAnneal(self):
esc_loc_x = m_to_km([
-2300850, -2297160, -2284680, -2283390, -2284800, -2289540,
-2287620, -2287740, -2287620, -2291760, -2289540, -2283720,
-2279730, -2254320, -2252430, -2253120, -2256900, -2273160,
-2273970, -2273910
])
esc_loc_y = m_to_km([
1986840, 1977120, 1966620, 1957680, 1947570, 1937730, 1926720,
1917720, 1907880, 1897830, 1887360, 1876560, 1867620, 1852470,
1843620, 1833720, 1824660, 1817640, 1807710, 1797930
])
ship_loc_x = m_to_km([
-2414875, -2401190, -2405002, -2406089, -2407670, -2402495,
-2402056, -2400759, -2390693, -2394883, -2393517, -2394351,
-2393072, -2396507, -2393492, -2394127, -2399780, -2396297,
-2387368, -2387021
])
ship_loc_y = m_to_km([
2019979, 2007266, 2001267, 1992909, 1982826, 1970152, 1959565,
1950068, 1937330, 1927338, 1917048, 1908072, 1899767, 1892623,
1883306, 1873335, 1864772, 1852304, 1839577, 1829662
])
centers = []
for i in range(0, len(esc_loc_x)):
center = (esc_loc_x[i], esc_loc_y[i])
centers.append(center)
interference_contour = []
for i in range(0, len(ship_loc_x)):
p = (ship_loc_x[i], ship_loc_y[i])
interference_contour.append(p)
testName = "SanFrancisco"
min_ctr_dist = 60
coverage_file = "detection-coverage/ITMDetectionCoverage_60deg.json"
poly = excessarea.generate_bounding_polygon(centers, line_endpoints)
cover = antennacover.min_antenna_area_cover_greedy(
centers, poly, coverage_file, min_center_distance=min_ctr_dist)
printcover.printAntennaCover(testName, poly, centers, cover,
coverage_file, min_ctr_dist)
testName = "SanFranciscoAnneal"
annealr = simannealer.SimAnneal(poly, coverage_file, cover)
annealr.anneal()
improved_cover = annealr.get_result()
printcover.printAntennaCover(testName, poly, centers, improved_cover,
coverage_file, min_ctr_dist)
def testEstuary1(self):
"""
Test for a deep estuary.
"""
interference_contour = [(20, 55), (35, 65), (40, 60), (45, 65),
(50, 55)]
possible_centers = [(20, 46), (25, 30), (30, 20), (40, 15), (50, 30),
(60, 50)]
min_ctr_dist = 0
poly = excessarea.generate_bounding_polygon(possible_centers,
interference_contour)
cover = antennacover.min_antenna_area_cover_greedy(
possible_centers,
poly,
"detection-coverage/ITMDetectionCoverage_90deg.json",
min_center_distance=min_ctr_dist)
testName = "Estuary"
printcover.printAntennaCover(
testName, poly, possible_centers, cover,
"detection-coverage/ITMDetectionCoverage_90deg.json", min_ctr_dist)
with open (protection_region, "r") as myfile:
data=myfile.readlines()
# Units are in Km. This should be converted to json format.
esc_loc_x = [x/1000.0 for x in eval(data[0])]
esc_loc_y = [x/1000.0 for x in eval(data[1])]
ship_loc_x = [x/1000.0 for x in eval(data[2])]
ship_loc_y = [x/1000.0 for x in eval(data[3])]
possible_centers = []
for i in range(0,len(esc_loc_x)):
center = (esc_loc_x[i],esc_loc_y[i])
possible_centers.append(center)
interference_contour = []
for i in range(0,len(ship_loc_x)):
p = (ship_loc_x[i],ship_loc_y[i])
interference_contour.append(p)
bounding_polygon = excess_area.generate_bounding_polygon(possible_centers,interference_contour)
testName = output_file
cover = antennacover.min_antenna_area_cover_greedy(possible_centers, bounding_polygon, coverage_file, min_center_distance=min_ctr_dist,tol=tol)
printcover.printAntennaCover(output_file,bounding_polygon, possible_centers, cover,coverage_file,min_ctr_dist)
if do_anneal:
annealr = simannealer.SimAnneal(bounding_polygon, possible_centers, coverage_file,cover,steps=args.anneal,tol=tol)
annealr.anneal()
cover = annealr.get_result()
printcover.printAntennaCover(output_file + "Anneal", bounding_polygon, possible_centers, cover,coverage_file,min_ctr_dist)
s[1] for s in sorted_detection_coverage_files
]
else:
# We are given a single antenna cover file.
detection_coverage_files = [detection_coverage_dir]
print "Antenna Detection Coverage files ", detection_coverage_files
min_cover_lobes = None
for detection_coverage_file in detection_coverage_files:
# Read the aperture angle from the file.
antenna_cover_patterns = antennacover.read_detection_coverage(
detection_coverage_file, coverage_units="m")
cover = antennacover.min_antenna_area_cover_greedy(
candidate_locs,
dpa_polygon,
detection_coverage_file,
min_center_distance=0,
tol=.005,
coverage_units="m")
cover_lobes = None
for c in cover:
center = c[0]
index = c[1]
angle = c[2]
lobe = antennacover.translate_and_rotate(
antenna_cover_patterns, center, index, angle)
if cover_lobes is None:
cover_lobes = lobe
else:
cover_lobes = cover_lobes.union(lobe)
#print "--- AREA ---", cover_lobes.area