def test_calc_pixel_neighbors_square_diagonal(): x, y = np.meshgrid(np.arange(20), np.arange(20)) cam = CameraGeometry( cam_id='test', pix_id=np.arange(400), pix_type='rectangular', pix_x=u.Quantity(x.ravel(), u.cm), pix_y=u.Quantity(y.ravel(), u.cm), pix_area=u.Quantity(np.ones(400), u.cm**2), ) cam._neighbors = cam.calc_pixel_neighbors(diagonal=True) assert set(cam.neighbors[21]) == {0, 1, 2, 20, 22, 40, 41, 42}
def test_calc_pixel_neighbors_square_diagonal(): x, y = np.meshgrid(np.arange(20), np.arange(20)) cam = CameraGeometry( cam_id='test', pix_id=np.arange(400), pix_type='rectangular', pix_x=u.Quantity(x.ravel(), u.cm), pix_y=u.Quantity(y.ravel(), u.cm), pix_area=u.Quantity(np.ones(400), u.cm**2), ) cam._neighbors = cam.calc_pixel_neighbors(diagonal=True) assert set(cam.neighbors[21]) == {0, 1, 2, 20, 22, 40, 41, 42}
def test_calc_pixel_neighbors_square_diagonal(): """ check that neighbors for square-pixel cameras are what we expect, namely that the diagonals are included if requested. """ x, y = np.meshgrid(np.arange(20), np.arange(20)) cam = CameraGeometry( camera_name="test", pix_id=np.arange(400), pix_type="rectangular", pix_x=u.Quantity(x.ravel(), u.cm), pix_y=u.Quantity(y.ravel(), u.cm), pix_area=u.Quantity(np.ones(400), u.cm**2), ) cam._neighbors = cam.calc_pixel_neighbors(diagonal=True) assert set(cam.neighbors[21]) == {0, 1, 2, 20, 22, 40, 41, 42}