1. Density of Starbucks

1.1. Task

You will write a program that will calculate the density of Starbucks locations within a specified area based on real data of Starbucks locations (as of 2018). You will have to work with some existing code that is provided to you, and you will need to write your own functions to ultimately complete the assignment. You will

  • Use a provided function to load data into your program

  • Write a function to convert latitude/longitude units from degrees, minutes, and seconds to decimal

  • Write a function to calculate the surface area of a sphere as defined by latitudes and longitudes

  • Write a function to count the number of starbucks contained within an area defined by latitudes and longitudes

  • Write a function to calculate the density of Starbucks within an area


You should not expect to be able to sit down and just start coding a solution. Programming does not work like that. Expect this assignment to take hours. Expect to get things wrong, and then, expect to get them wrong again — this is normal.

1.2. Provided Files

You are provided with

  • A notebook file called asn1.ipynp containing the starting point of the assignment

    • This file is to be uploaded to Google Colab

    • Alternatively, if you prefer to complete the assignment with an IDE on your own computer, you may download and use the asn1.py file

  • The data file called starbucks2018.csv containing latitude and longitude pairs of Starbucks locations


Do not alter the function details in the provided .ipynb/.py files

  • Do not change the name of the functions

  • Do not remove the function description

  • Do not remove or add to the parameter list

1.3. Part 0 — Read the Assignment

Read the assignment description in its entirety before starting.

1.4. Part 1 — Uploading Files to Colab

After downloading the notebook file above, you will need to upload it to Colab to get started. See the below image to help find how to do this. I recommend saving a copy of this notebook file to your Google drive and then work with that one. You don’t have to, but you will have to re-upload the project every time you want to work on it.


Additionally, you will need to upload the Starbucks location data file to Colab. The way you upload a data file like this is different from uploading the notebook file. See the below image for an example of how to upload this file. Unfortunately, you must re-upload this data file every time you open your Colab project. Not a big deal, but something you will need to remember.


1.5. Part 2 — Read Over Loading Function

The function that loads the data, load_starbucks_data, is already provided for you. It contains ideas we have yet to discuss in class, but it shouldn’t be too difficult to get an intuition about what exactly it is doing if you look over it and play with it a little.

Get used to looking at code that is not yours, using unfamiliar ideas, and trying to figure out what existing code does. This is not a trivial thing, but as a programmer, it is something you will end up spending a lot of time doing.


For the loading function to work, the .csv file must be in the same directory as your Python script. This means in your Colab project (see Part 1). If it is not, this function will not work and you will see an error message like No such file or directory: 'starbucks2018.csv'.

1.6. Part 3 — Degrees to Decimal

Humans like to express latitudes and longitudes in degrees, minutes, and seconds. You could work with those units in Python, but your computations (and, hence, code) will be much cleaner and easier to follow if you convert the data into the single unit “degrees”, using decimals of a degree to represent “arcminutes” (1/60th of a degree) and “arcseconds” (1/60th of a arcminute).

Complete the function convert_degrees_to_decimal such that it converts the provided degrees, arcminutes, and arcseconds to decimal. The parameters are a latitude/longitude in degrees, arcminutes, and arcseonds. The function should return the same latitude/longitude as a single value in decimal degrees (a single value of type float).

If you are not familiar with the conversion, check Wikipedia.

If you had a look at the data file, you will have noticed that the data is already stored as a decimal. This means you do not actually need to use this function to convert the data from the Starbucks location file.

1.7. Part 4 — Subtended Area

Since we want to calculate the density of starbucks, we need to consider the units used. A reasonable measure would be Starbucks per square kilometers. In order to compute this, we must first calculate the area of the “rectangle” defined by two latitude and two longitudes. You are to complete the function subtended_area that takes four parameters. Two latitudes defining the top and bottom of the “rectangle” and two longitudes defining the sides of the “rectangle”. The function will return the area of the “rectangle” in kilometers squared (\(km^{2}\)).

However, in reality, since it is a sphere that the surface area is being calculated on, the latitude and longitudes do not actually define a rectangle we are familiar with. This means that the simple \(length * height\) will not work. Instead we need to calculate it with the following equation.

\(\frac{\pi}{180} \cdot R^{2} \cdot \lvert sin(lat_{1}) - sin(lat_{2}) \rvert \cdot \lvert lon_{1} - lon_{2} \rvert\)

In our case we will use \(R = 6371\) for Earth, which is stored in the provided file as the constant EARTH_RADIUS.


Does Python’s trig functions (eg., math.sin) expect parameters in degrees or radians? Read the relevant documentation to find out.

1.8. Part 5 — Counting Starbucks

In order to calculate the density of Starbucks, the number of Starbucks within the specified area needs to be known. With the data available, the way to do this is to check each Starbucks’ latitude & longitude and check if it falls within the specified “rectangle”. In other words, check if the Starbucks’ latitude falls between the “rectangle’s” latitudes and if the longitude falls between the “rectangle’s” longitudes. See the below image for an example.


The function number_starbucks_within_area() takes the list of Starbucks locations and the latitude and longitudes specifying the “rectangle” as parameters. For simplicity, assume latitude_line_1 < latitude_line_2 and longitude_line_1 < longitude_line_2. The function is already set up to loop over every Starbucks location in a list. Have a look a the loop — even if we haven’t formally discussed this in class, there is a good chance you can make sense of what it’s doing.

You are to complete the body of the loop. Each time through the loop, we’ll be considering a new Starbucks location. The existing code already stores the current Starbucks location’s latitude and longitude values in their respective variables. You are to figure out if this specific location falls within the area defined by the latitudes and longitudes passed to the function as parameters. If the location is within the area, we count it, otherwise, we do not. Keep track of the running total of Starbucks within the area and, when the loop is finished checking each Starbucks location, the function will return the final count.

1.9. Part 6 — Calculate Starbucks Density

Complete the starbucks_per_square_kilometer function that, given a file name and the latitude and longitudes to define a “rectangle”, calculates and returns the density of Starbucks within that “rectangle”.

The function will need to

  • Load the data

  • Calculate the area of the “rectangle”

  • Count the number of Starbucks within the “rectangle”

  • Calculate the density of Starbucks within the “rectangle” — divide the number of Starbucks by the area

  • Return the density

1.10. Part 7 — Using Your Function

Play around with the starbucks_per_square_kilometer function. Try some small “rectangles” and big ones. What area has the highest Starbucks density you can find? The lowest?

After playing with the function a little, record within a text file the smallest and largest densities you found and what the parameters were that you used to get the densities. You do not need to find the largest or smallest possible densities — simply try a few parameters and see what you get.

1.11. Part 8 — Testing

To help ensure that your program is correct, run the provided assertion tests. Each function is followed by a series of commented out assertion tests that will help you test your code. When you are ready to test your functions, simply make them not comments (remove the #) to include them in your running program. There is no guarntee that if your code passes all the tests that you will be correct, but it certainly helps provide peace of mind that things are working as they should.

Realistically you should have been running tests after you complete each of the above parts, but this part is here to remind you. Remember, we are lucky that we get to test our solutions for correctness ourselves; you don’t need to wait for the marker to return your assignment before you have an idea of if it works correctly.

1.12. Some Hints

  • Work on one function at a time

  • Get each function working perfectly before you go on to the next one

  • Test each function as you write it
    • This is a really nice thing about programming; you can call your functions and see what result gets returned

    • Mentally test before you even write — what does this function do? What problem is it solving?

  • If you need help, ask
    • Drop by office hours

1.13. Some Marking Details


Just because your program produces the correct output, that does not necessarily mean that you will get perfect, or even that your program is correct.

Below is a list of both quantitative and qualitative things we will look for:

  • Correctness?

  • Did you follow instructions?

  • Comments?

  • Variable Names?

  • Style?

  • Did you do just weird things that make no sense?

1.14. What to Submit to Moodle

  • Make sure your NAME and STUDENT NUMBER appear in a comment at the top of the program

  • Submit your version of asn1.py to Moodle
    • Do not submit the .ipynb file

    • To get the asn1.py file from Colab, see the image below

  • Also submit your text file describing the areas you found with the highest, and lowest, Starbucks densities and a short description of how you found them
    • Don’t worry about finding the highest or lowest density values, just try a few and pick your highest and lowest


Verify that your submission to Moodle worked. If you submit incorrectly, you will get a 0.


1.15. General FAQ:

  • Does my text file have enough details?
    • Probably. The shorter the better. I really just want to see that you played around a little.

  • Is my area a high/low enough density?
    • I really don’t care how high/low it is. Just try a few things and see what you get.

  • I don’t know how to do X.
  • It’s not working, therefore Python is broken!
    • Probably not; you’re very likely doing something wrong

  • My thing keeps telling me No such file or directory: 'starbucks2018.csv'
    • Then the starbucks file probably isn’t where python is looking.

  • But density will grow larger the smaller I make the area.
    • Correct, you understand basic arithmetic.

  • But I never used the one function!!!!1!
    • Fine, but write the code anyways and make sure it works.

  • But the degrees values don’t specify a cardinal direction!
    • Make use of changing +/- if you need to change hemispheres.

  • Wtf do the functions do that you gave me?
    • Read the descriptions. Try figuring it out. This is actually part of the assignment learning objectives.

  • Some of the code in the functions you gave us look like magic.
    • At this point it may seem that way, but by the end of the semester, they will lose their magic.

  • Do I have enough comments?
    • I don’t know, maybe? If the code is too complex to understand simply by looking at it, add a comment.

  • I know you told me to do it this way, but I did it another way, and I think my way is better.
    • Your way may be better, but I don’t care. Do it the way I told you.

  • Can I work with my friend?
    • No.

  • I know I cheated, I know I know I was cheating, but I’m reeeeaaaaaaaaallllllly sorry [that I got caught]. Can we just ignore it this time?
    • Lol, no.

  • If I submit it at 11:56pm, you’ll still mark it, right? I mean, commmmon!
    • No. 11:55pm and earlier is on time. Anything after 11:55pm is late. Anything late is not marked. It’s rather simple really.

  • Moodle was totally broken, it’s not my fault it’s late.
    • Nice try.

  • I accidentally submitted the wrong code. Here is the right code, but it’s late. But you can see that I submitted the wrong code on time! You’ll still accept it, right?
    • Do you think I was born yesterday? No.

  • Will I really get 0 if I do the submission wrong? Like, what if I submit the .ipynb instead of the .py?
    • Yes, you’ll really get a ZERO.