def is_arching_scene(locs, dirs, scene, MAX_ANGLE):
    # for i in range(scene.start, scene.end, 1):
    #     print locs[i], dirs[i]
    if (scene.end - scene.start) < 2:
        return False

    # print scene.end - scene.start

    comb = itertools.combinations(range(scene.start, scene.end + 1, 1), 2)
    for (i, j) in comb:
        if j != i and dirs[i] != dirs[j]:   # non-parallel pairs
            convergent = True
            # only consider pair of locations whose distance is less than 1 meter
            if distance_km(locs[i][0], locs[i][1], locs[j][0], locs[j][1]) < 0.001:
                continue
            l_o = interesection_point(locs[i], locs[j], dirs[i], dirs[j])   # find intersection point
            # print l_o
            if not l_o:
                continue
            for k in range(scene.start, scene.end + 1, 1):
                if k != i and k != j:   # the rest of FOVs
                    if distance_km(l_o[0], l_o[1], locs[k][0], locs[k][1]) > 0.01 * R_100:
                        # print "distance ", distance_km(l_o[0], l_o[1], locs[k][0], locs[k][1])
                        convergent = False
                        break
                    alpha1 = angle_bwn_points(locs[k], l_o)
                    alpha2 = correct_dir(dirs[k])
                    # print alpha1, alpha2, MAX_ANGLE
                    if not within_range(alpha1, alpha2, MAX_ANGLE):
                        convergent = False
                        break

            if convergent:  # exist one convergent center --> return True
                return True
    return False
Esempio n. 2
0
 def cover(self, plat, plng):
     angle = angle_bwn_two_points(self.lat, self.lon, plat, plng)
     distance = distance_km(self.lat, self.lon, plat, plng)
     # print angle, distance
     if abs(angle - self.compass / 2) < 30 and distance < self.R:
         return True
     else:
         return False
Esempio n. 3
0
 def cover(self, plat, plng):
     angle = angle_bwn_two_points(self.lat, self.lon, plat, plng)
     distance = distance_km(self.lat, self.lon, plat, plng)
     # print angle, distance
     if abs(angle - self.compass / 2) < 30 and distance < self.R:
         return True
     else:
         return False
Esempio n. 4
0
 def cover(self, plat, plng):
     angle = angle_bwn_two_points(self.lat, self.lon, plat, plng)
     distance = distance_km(self.lat, self.lon, plat, plng)
     # print angle, distance
     if abs(angle - self.compass/2) < 30 and distance < self.R:
         return True
     else:
         return False
 # print within_circular(34.024734, -118.284988,34.018212,-118.291716,45,1)
Esempio n. 5
0
 def area(self):
     a = distance_km(self.n_box[0, 0], self.n_box[0, 1], self.n_box[0, 0],
                     self.n_box[1, 1]) * distance_km(
                         self.n_box[0, 0], self.n_box[0, 1],
                         self.n_box[1, 0], self.n_box[0, 1])
     return a
Esempio n. 6
0
 def area(self):
     a = distance_km(self.n_box[0,0],self.n_box[0,1],self.n_box[0,0],self.n_box[1,1]) * distance_km(self.n_box[0,0],self.n_box[0,1],self.n_box[1,0],self.n_box[0,1]);
     return a