We use cookies and other tracking technologies to improve your browsing experience on our website, to show you personalized content and targeted ads, to analyze our website traffic, and to understand where our visitors are coming from.
Please recreate the figure, above. You are welcome to reuse the code found on the website. Using this figure as starting point, do the following:
When plotting the earthquakes, include only the earthquake data for earthquakes of magnitude greater than the 75th percentile. Add a point for five of your favourite South Pacific island nations. Ensure the point is correctly associated with the island name and that the map is correctly labelled, has a title, and it is as close to publication quality as you can make it. Your script needs to show all the steps (thoroughly annotated) leading to the final figure.
Marks will also be assigned for the overall aethetic appearance of the map. Feel free to be creative, but ensure the final product remains publication quality. (\30)
Question 2
Successfully completing on of the options available in this task will earn you a bonus of up to 8 or 10% onto your CA mark.
You have until 31 March 2025 to complete it.
A map that is worthy of display will become a large format poster to display in the BCB Department. Your name displayed next to it will immortalise you for continued fame and glory amongst future BCB students.
The winner of each category of map (hypometric and non-hypsometric) will also get a box of Lindt chocolate.
Option 1 [up to 10% bonus]: Create a hypsometric map based on these examples
The maps show the locations of linefish catches along the SA coast as per a DFFE dataset. I do not expect that you add these data points as you don’t have access to this dataset. However, the location of the 58 coastal sections indicated by circles can be plotted using the data provided here. You are also welcome to create a map of any topographically-interesting region on Earth, but be sure to include a few data points of some kind to draw our attention to some interesting features or statistics. Be creative!
Since I think a few of you might actually accomplish this, best add a few improvements to it to make your map even better than mine and stand out from that of your peers. There can be only one winner in each category, and the best one wins (although everyone can benefit from the bonus marks).
Warning: You’ll need a fairly beefy computer to accomplish this task.
Option 2 [up to 8% bonus]: Create an artistic map of your choice
Alternatively, if you cannot access a powerful computer, for a bonus of up to 8% onto your CAM, create any (non-hypsometric) map of your choice of any region on Earth. Make something that you would be proud to display as a large format poster. The map may draw attention to an interesting regional geophysical, ecological, or socio-ecological (etc.) phenomena, or it may simply showcase your unique (but tasteful!) artistic ability. Show me some examples of what you wish to create before you start to avoid wasting your time on something too simple or entirely tasteless. There are many examples of beautiful maps on the internet that you may use as source of inspiration.
Whichever option you choose, please also submit your code together with the final product in a well-described Quarto .html document. Explain each step of the way and describe the rationale for the approach you take.
---title: "BCB744 Bonus Task"format: html: fig-format: svg fig_retina: 2 fig-dpi: 400params: hide_answers: true---# [[Assessment Sheet](BCB744_Task_Bonus_Surname.xlsx)]{.my-highlight} {#sec-assessment}# 11. The Fiji Earthquake Data```{r, echo=FALSE,warning=FALSE,message=FALSE,fig.width=6,fig.asp=0.65,out.width="80%",fig.align='center'}library(tidyverse)library(sf)library(rnaturalearth)library(rnaturalearthdata)world <-ne_countries(returnclass ='sf',scale =10, type ="countries") |>select(continent, sovereignt, iso_a3) NE_proj <-"+proj=natearth +lon_0=170 "world_1 <-ne_countries(returnclass ='sf',scale =10, type ="countries") |>select(continent, sovereignt, iso_a3) |>st_break_antimeridian(lon_0 =170) |>st_transform(NE_proj)library(countrycode)world_1$region <-countrycode(world_1$iso_a3, origin ="iso3c",destination ="un.regionsub.name")sw_pacific <- world_1 |>filter(region %in%c("Australia and New Zealand", "Melanesia", "Micronesia","Indonesia", "Polynesia", "South-eastern Asia")) |>group_by(continent) |>summarise()quakes <-as_tibble(datasets::quakes)margin <-15.0xmin <-min(quakes$long) - margin; xmax <-max(quakes$long) + marginymin <-min(quakes$lat) - margin; ymax <-max(quakes$lat) + marginWGS84_proj <-"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"bbox <-st_sfc(st_point(c(xmin, ymin)), st_point(c(xmax, ymax)),crs = WGS84_proj)bbox_trans <-st_transform(bbox, NE_proj)sw_pacific_cropped <- sw_pacific |>st_crop(bbox_trans)quakes_sf <- quakes |>st_as_sf(coords =c("long", "lat"),crs = WGS84_proj)quakes_sf_trans <-st_transform(quakes_sf, NE_proj)quakes_sf_trans_sub <- quakes_sf_trans |>filter(mag >quantile(mag, 0.75))ggplot() +geom_sf(data = sw_pacific_cropped, colour ="black", fill ="grey70") +geom_sf(data = quakes_sf_trans, aes(colour = mag, size = mag),stat ="sf_coordinates",shape ="*", alpha =0.4) +scale_colour_continuous(type ="viridis") +guides(size ="none") +coord_sf(expand =FALSE) +labs(x =NULL, y =NULL,title ="The Fiji Earthquake Data",subtitle ="All data")```The questions concern the above figure.## Question 1Please recreate the figure, above. You are welcome to reuse the code found on the website. Using this figure as starting point, do the following:When plotting the earthquakes, include **only** the earthquake data for earthquakes of magnitude greater than the 75th percentile. Add a point for five of your favourite South Pacific island nations. Ensure the point is correctly associated with the island name and that the map is correctly labelled, has a title, and it is as close to publication quality as you can make it. Your script needs to show all the steps (thoroughly annotated) leading to the final figure.Marks will also be assigned for the overall aethetic appearance of the map. Feel free to be creative, but ensure the final product remains publication quality. **(\30)**`r if (params$hide_answers) "::: {.content-hidden}"`**Answer**This is the 'base' map, which you will modify:```{r, echo=TRUE,warning=FALSE,message=FALSE,fig.width=6,fig.asp=0.65,out.width="80%",fig.align='center'}library(tidyverse)library(sf)library(rnaturalearth)library(rnaturalearthdata)world <-ne_countries(returnclass ='sf', scale =10, type ="countries") |>select(continent, sovereignt, iso_a3) NE_proj <-"+proj=natearth +lon_0=170 "world_1 <-ne_countries(returnclass ='sf',scale =10, type ="countries") |>select(continent, sovereignt, iso_a3) |>st_break_antimeridian(lon_0 =170) |>st_transform(NE_proj)library(countrycode)world_1$region <-countrycode(world_1$iso_a3, origin ="iso3c",destination ="un.regionsub.name")sw_pacific <- world_1 |>filter(region %in%c("Australia and New Zealand", "Melanesia", "Micronesia","Indonesia", "Polynesia", "South-eastern Asia")) |>group_by(continent) |>summarise()quakes <-as_tibble(datasets::quakes)margin <-15.0xmin <-min(quakes$long) - margin; xmax <-max(quakes$long) + marginymin <-min(quakes$lat) - margin; ymax <-max(quakes$lat) + marginWGS84_proj <-"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"bbox <-st_sfc(st_point(c(xmin, ymin)), st_point(c(xmax, ymax)),crs = WGS84_proj)bbox_trans <-st_transform(bbox, NE_proj)sw_pacific_cropped <- sw_pacific |>st_crop(bbox_trans)quakes_sf <- quakes |>st_as_sf(coords =c("long", "lat"),crs = WGS84_proj)quakes_sf_trans <-st_transform(quakes_sf, NE_proj)# Save to reuse laterquake1 <-ggplot() +geom_sf(data = sw_pacific_cropped, colour ="black", fill ="grey70") +scale_colour_continuous(type ="viridis") +guides(size ="none") +coord_sf(expand =FALSE) +labs(x =NULL, y =NULL,title ="The Fiji Earthquake Data",subtitle ="All data")# Add the earthquakesquake1 +geom_sf(data = quakes_sf_trans, aes(colour = mag, size = mag),stat ="sf_coordinates",shape ="*", alpha =0.4)```This is the final map you will have produced: ```{r, echo=TRUE,warning=FALSE,message=FALSE,fig.width=6,fig.asp=0.65,out.width="80%",fig.align='center'}# Filter the dataquakes_sf_trans_sub <- quakes_sf_trans |># ✓filter(mag >quantile(mag, 0.75))# Add the pointsislands <-read.csv("../data/pacific_nations.csv") # ✓# Convert to sfislands_sf <- islands |># ✓st_as_sf(coords =c("Longitude", "Latitude"),crs = WGS84_proj)# Transform to the new projectionislands_sf_trans <-st_transform(islands_sf, NE_proj) # ✓quake1 +geom_sf(data = quakes_sf_trans_sub, aes(colour = mag, size = mag),stat ="sf_coordinates",shape ="*", alpha =0.4) +geom_sf(data = islands_sf_trans, size =1.2, colour ="red") +# ✓ geom_sf_text(data = islands_sf_trans, aes(label = Country), nudge_x =500000, nudge_y =200000, size =3) +# ✓scale_colour_viridis_c(name ="Magni-\ntude") +# ✓labs(title ="The Fiji Earthquake Data", # ✓subtitle ="Quakes > 75 percentile",x ="Lon (°E)", y ="Lat (°S)")````r if (params$hide_answers) ":::"`## Question 2Successfully completing on of the options available in this task will earn you a *bonus of up* to **8 or 10% onto your CA mark**.You have until 31 March 2025 to complete it.- A map that is worthy of display will become a large format poster to display in the BCB Department. Your name displayed next to it will immortalise you for continued fame and glory amongst future BCB students.- The winner of each category of map (hypometric and non-hypsometric) will also get a box of Lindt chocolate.**Option 1 [up to 10% bonus]: Create a hypsometric map based on these examples**The maps show the locations of linefish catches along the SA coast as per a DFFE dataset. I do not expect that you add these data points as you don't have access to this dataset. However, the location of the 58 coastal sections indicated by circles can be plotted using the [data provided here](../../data/interp_LowRes_coords_58_sites.txt). You are also welcome to create a map of any topographically-interesting region on Earth, but be sure to include a few data points of some kind to draw our attention to some interesting features or statistics. Be creative!Since I think a few of you might actually accomplish this, best add a few improvements to it to make your map even better than mine and stand out from that of your peers. There can be only one winner in each category, and the best one wins (although everyone can benefit from the bonus marks).Warning: You'll need a fairly beefy computer to accomplish this task.**Option 2 [up to 8% bonus]: Create an artistic map of your choice**Alternatively, if you cannot access a powerful computer, for a **bonus of up to 8% onto your CAM**, create any (non-hypsometric) map of your choice of any region on Earth. Make something that you would be proud to display as a large format poster. The map may draw attention to an interesting regional geophysical, ecological, or socio-ecological (etc.) phenomena, or it may simply showcase your unique (but tasteful!) artistic ability. Show me some examples of what you wish to create before you start to avoid wasting your time on something too simple or entirely tasteless. There are many examples of beautiful maps on the internet that you may use as source of inspiration.Whichever option you choose, please also submit your code together with the final product in a well-described Quarto .html document. Explain each step of the way and describe the rationale for the approach you take.Good luck!
