Processing math: 100%

CS440/ECE448 Fall 2018

Assignment 2: Configuration Space Planning

Due date: Monday September 30, 11:59pm

Credits (Fall 2018): Micky Abir, Chris Benson, Krishna Kothapalli, and JD Lee
Credits (Fall 2019): Zih-Siou Hung, Devansh Shah

In this assignment you will write code that transforms a 2D planning problem for a robotic arm into a configuration space, and then searches for a path in that space. See Section 25.4 of the textbook and Lecture 7 for background information.

Contents

General guidelines

Basic instructions are the same as in MP 1. Your code must be in Python 3.6 or 3.7 and will be submitted to gradescope. Your code may import extra modules, but only ones that are part of the standard python library . Extra credit is capped at 10% (so the score for the MP will not exceed 110%).

Two changes from MP 1:

  • In addition to the standard python libraries and pygame, you may import numpy.
  • The extra credit portion will be submitted separately on gradescope.

For general instructions, see the main MP page and the course policies.

You will need to re-use your bfs code from MP 1.

Problem Statement

You are given a two-link arm in 2D space as explained in class. This arm has two links of length L1 and L2 respectively. The arm is free to rotate about its base that is pivoted on the ground and link-2 can rotate about the joint where it connects with link-1. Let's use α for the angle between the link-1 and the ground (equivalent to θ1). Let's use β for the angle between link-2 and link-1 (equivalent to θ2). The robot is shown at the top of the page with visualization of θ1 and θ2. Note that the angles are measured counter-clockwise.

For each planning problem, you will be given:

  • The starting angles for the robotic arm.
  • The location (x,y) and radius r of a circular goal.
  • The locations (x,y) and radius r of various circular obstacles.

More details can be found in Map Configuration.

You need to find a shortest path for the robotic arm from its starting position to the goal so that the tip (open end of the arm) touches or is inside the goal.

The tricky bit is that the robotic arm may not pass through any of the circular objects including a goal and an obstacle. Also, any part of the arm should not get out of the given window. So configurations like the following are NOT allowed:

Note that the blue and the red circle denote a goal and an obstacle, respectively.

You will do your path planning in two steps:

  • Compute a configuration space map (Maze) that shows which joint angles are blocked by obstacles or by the window.
  • Use your code from MP 1 (search.py) to compute a shortest path in this configuration space map. Use the bfs search algorithm for search. Make sure to return the empty path if there is no path from start to goal in the maze

For compatibility with MP 1, you will digitize the angles α and β when creating your configuration space map so that the map will only be an approximation to real life. Also, the arm is allowed to only change one of α or β in one step. Each step will change one of the angles by a fixed number called the angular resolution that is measured in degrees. The angular resolution will be an input parameter to your method, so that it can be adjusted for your tests and for our evaluation of your code. There is a tradeoff here: finer resolution will allow the arm to get through smaller gaps between obstacles but it will also slow down your search.

Part 0: Map configuration

First of all, you need to understand the 2D space that describes a robotic arm, a goal, and obstacles. The following example map, named 'Test1', would help you to understand map configuration specifications.
[Test1]
Window    : (300, 200)      # (Width, Height)   
ArmBase   : (150, 190)      # (x-coordinate, y-coordinate)
ArmLinks  : [
                (100, 40, (20, 160)),    # (length, initial angle, (min angle, max angle)
                (80, 0, (-150, 150)),
            ]
Obstacles : [
                (70, 50, 15),     # (x-coordinate, y-coordinate, radius)
                (140, 30, 17),
                (115, 75, 17)
            ]
Goals     : [
                (110, 40, 10)      # (x-coordinate, y-coordinate, radius)
            ]

The window size for the given example map is 300x200 pixel. The arm base is placed at (150, 200). Note that (0, 0) is the left-top corner and (300, 200) is the right-bottom corner. There are two arm links. The length of the first arm link is 100 and its initial angle (α) is 40. The minimum and maximum angle for this arm link is 20 and 160. The next arm link can be explained in the same manner.

Look carefully at the conventions for α and β in the diagram at the top of this page. Both angles are measured counter-clockwise. α is relative to the ground and β is relative to the direction of link-1.

There are three obstacles, each of them is described in tuple (x-coordinate, y-coordinate, radius). Goals can be described in the same convention. Technically, there can be more than one goal. However, for simplicity, you can assume that there is only one goal for this MP. The generated 2D map for above configuration is shown in following:

You can play with the above map with the following command if your geometry functions are correctly implemented:

python3 mp2.py --human --map Test1 --config test_config.txt

Once the window pops up, you can rotate the arm using the following keys:

  • z / x: increase / decrease α
  • a / s: increase / decrease β
  • q / w: increase / decrease γ (if you have 3 arm links)

Part 1: Geometry

The first part of this MP is to work out the geometrical details.

The first problem in geometry is to compute (start, end) position of each arm link for the given α and β. The base position, initial angle, and length of link-1 is given in configuration, so the end position of link-1 can be computable. The end position becomes the base position of link-2. Then, the end position of link-2 can be also computable for the given length and β. We kindly provide two helper classes: Arm and ArmLink. One arm can have at least 1 arm link and at most 3 arms links. You do not have to modify these two classes. To solve MP2, the following Arm methods would be useful:

arm.py

  • getBase(): Return (x, y) of the arm base
  • getArmAngle(): Return relative angles of all arm links. If there are two arm links, the return value would be (alpha, beta)
  • getArmLimit(): Return (min angle, max angle) of all arm links
  • setArmAngle(angles): Set angles(alpha, beta, gamma) for all arm links
  • getEnd(): Returns (x, y) of the arm tip
  • getArmPos(): Returns (start, end) of all arm links

The relationship among arm links are already defined in Arm. What you need to define is how to compute the end position of the arm link for the given length and angle. To this end, you may need to implement the following function:

geometry.py

  • computeCoordinate(start, length, angle): Return the end position for the given start position and length. Please round down the trigonometry function to integer.

Once you correctly define the above function, above red colored getEnd() and getArmPos() functions would perform as expected.

The next geometry problem is then to check whether

  • the arm tip (the open end of the arm) touches or is inside the goal.
  • the arm runs into obstacles. Note that the goal should be also considered as an obstacle unless the robotic arm tip touches it.
  • all parts of the arm remain in the given window.

To do this, you need to implement the following functions:

geometry.py

  • doesArmTouchGoals(armEnd, goals): Return True if the given arm tip position is anywhere in the given goals, including the boundary. Or return False.
  • doesArmTouchObstacles(armPos, obstacles): Return True if the given arm touches any obstacle, including the boundary. Or return False.
  • isArmWithinWindow(armPos, window): Return True if all parts of the arm are within the window. Or return False.

Once all geometry related functions are implemented, you can check whether they function correctly through manual robotic arm control. You will get a 'SUCCESS' message when the arm tip touches anywhere of the goal.

Part 2: Transformation to Maze

The second part of this MP is to transform the generated 2D map to the maze in MP1. Rows and columns in the maze would match with alpha and beta, joint angles of the robotic arm. The number of rows and columns of the maze can be decided by 1) degree granularity (angular resolution) and 2) minimum / maximum value of the angle. The number of positions in the row or column is maxangleminanglegranularity+1. For example, if range of alpha is (45, 135) and degree granularity is 2, then the number of rows becomes 46=135452+1. When the output of the division is not an integer, round it down.

Once the maze size is determined, then you need to fill the starting point, walls, goals in the maze. Initial alpha and beta would be the starting point. To contruct walls and goals, you need to use geometry functions defined in Part 1.

The naive approach is to do bruteforce checking for all pairs of alpha and beta. Or you can apply some smart short-cut techniques to bypass unnecessary checking. For example, if the first link already hits an obstacle, you may not need to check beta for that alpha. Although you use your own shortcut algorithm, the output maze should look same as the result of the bruteforce algorithm.

Your main effort for Part 2 would be to implement the following function.

transform.py

  • transformToMaze(arm, goals, obstacles, window, granularity): This function will return the maze instance created based on the given input arguments.
  • The input arguments are:
    • arm (Arm): arm instance
    • goals (list): [(x, y, r)] of goals, where x and y are coordinates, and r is radius
    • obstacles (list): [(x, y, r)] of obstacles, where x and y are coordinates, and r is radius
    • window (tuple): (width, height) of the window
    • granularity (int): unit of increasing/decreasing degree for angles
  • Hint: use utils.py.

The following images are examples of output maze for Test1 above. The left image has the finer granularity and the right image has the coarser granularity. As you can see, the maze resolution is determined depending on the input granularity. As can be seen, you should be able to handle any level of granularity. Granularity is a natural number for simplicity.

mp2.py does not generate the maze map as above. If you would like to visualize the map, you should run mp2.py to export the maze to a text file, which is an input file of mp1.py. Then, you can run your mp1.py to create the maze map. That is,

python3 mp2.py --map Test1 --config test_config.txt --save-maze [MAZEFILE_PATH]
python3 mp1.py [MAZEFILE_PATH]

Note that maze.py in MP1 and MP2 is slightly different, so please do not use MP1 maze.py for MP2 and vice versa.

Part 3: Searching the path in Maze

Now we are back to MP1. You can exploit your previous search functions to find the path to exit the maze with minor modifications. In MP2, you do not have to worry about multiple goals (dots). Your search function should be terminated after your agent reaches any of goals, Otherwise, it would take ages for agents to visit all the goals in the maze. The final state of the robotic arm for Test1 is shown in the left image. You can leave the footprint of the arm rotation in every n move with "--trajectory n" option. Accordingly, the explored shortest path for the maze is described in the right image.

Extra Credit

Generalize your computation to robots with one and three links. This means that your configuration space map will need to be three-dimensional. We have included a sample problem involving a 3-link (Test1, Test2) robot. Can your algorithm solve this problem?

In the geometry section above, rough approximation to decide a contact between the robotic arm and objects will make the maze inaccurate. Can you do a more exact analysis of when the link can contact a circular object?

Provided Code Skeleton

The code you need to get started is in tar file and zip file. You will only have to modify following files:

  • geometry.py
  • transform.py
  • search.py
  • maze.py modify and submit this only for the extra credit part

Do not modify other provided code. It will make your code not runnable.

To understand how to run the MP, read the provided README.md or run python3 mp2.py -h into your terminal. The following command will display a 2D map and let you control the robotic arm manually (q and w for link-1, a and w for link-2, z and x for link-3).

python3 mp2.py --human --map [MAPNAME] --config [CONFIGFILENAME]

You can also save your output picture using --save-image option. If you want to export your maze to the text file, use --save-maze option.

The following command will run your bfs search method on the maze.

python3 mp1.py --method bfs --map [MAZEFILE_PATH]

Tips

  • Think about how to accurately compute the distance between a line segment and a point.
  • When constructing the maze, start with the bruteforce algorithm and then think about the shortcuts.
  • When loading the maze in mp1.py, if the maze map size is too big, control the size using --scale option to make a proper size.

Deliverables

This MP will be submitted via gradescope.

Please upload only search.py, transform.py, geometry.py one by one to gradescope for MP2.

For extra credit, upload search.py, transform.py, geometry.py, and maze.py to MP2_extra.

After submission, gradescope will run preliminary tests to determine whether or not your submission appears valid. These preliminary tests are worth 0 points, and are for your information only. Passing all of these preliminary tests does not guarantee that your implementations are correct.

Here's what the autograder outputs if you simply submit the provided skeleton code:

If you don't see a similar series of test results in a reasonable amount of time, the autograder has probably crashed (e.g. because you imported a module that wasn't allowed for this MP) or has gone into an infinite loop. Fix your code and resubmit. If the autograder is still running, a resubmission will cause it to stop work and restart on your new submission.