Chapter 36 Distribution of Aus Meteorite Landings with Meteorite Mass

The mass of the Meteorites are indicated by the Radius of the Circles.The mass has been converted to Kgs and then multiplied by 10

top = -45 # north lat
left = 110 # west long
right = 145 # east long
bottom =  -10 # south lat


AusMetLandings = MetLandings %>%
  filter(reclat <=bottom) %>%
  filter ( reclat >= top) %>%
  filter( reclong >= left ) %>%
  filter(reclong <= right)

center_lon = median(AusMetLandings$reclong,na.rm = TRUE)
center_lat = median(AusMetLandings$reclat,na.rm = TRUE)


factpal <- colorFactor(c("red","blue"), 
                       AusMetLandings$fall)

leaflet(AusMetLandings) %>% addProviderTiles("Esri.NatGeoWorldMap") %>%
  addCircles(lng = ~reclong, lat = ~reclat,radius = ~(mass/1e3)*10 ,
             color = ~factpal(fall))  %>%
  # controls
  setView(lng=center_lon, lat=center_lat,zoom = 4) %>%
  
  addLegend("bottomright", pal = factpal, values = ~fall,
            title = "Meteorites landings and fall",
            opacity = 1)

36.1 Peek into the Aus Met Landings Data

AusMetLandings = AusMetLandings %>% arrange(desc(mass))

kable(head(AusMetLandings,6),"html") %>%
  kable_styling(bootstrap_options = c("striped", "hover", "condensed", "responsive")) %>%
  scroll_box(width = "800px")

name	id	nametype	recclass	mass	fall	year	reclat	reclong	GeoLocation
Mundrabilla	16852	Valid	Iron, IAB-ung	24000000	Found	1911	-30.78333	127.5500	(-30.783330, 127.550000)
Youndegin	30374	Valid	Iron, IAB-MG	3800000	Found	1884	-32.10000	117.7167	(-32.100000, 117.716670)
Huckitta	11922	Valid	Pallasite, PMG-an	2300000	Found	1924	-22.36667	135.7667	(-22.366670, 135.766670)
Henbury	11872	Valid	Iron, IIIAB	2000000	Found	1931	-24.56667	133.1667	(-24.566670, 133.166670)
Murnpeowie	16878	Valid	Iron, IC	1143000	Found	1909	-29.58333	139.9000	(-29.583330, 139.900000)
Wolf Creek	24326	Valid	Iron, IIIAB	760000	Found	1947	-19.30000	127.7667	(-19.300000, 127.766670)

36.2 Heatmap of Aus Meteorite Landings

The intensity of the Heatmap is based on the mass of the meteorites.

AusMetLandings %>% leaflet() %>% addProviderTiles("Esri.NatGeoWorldMap") %>% 
  addHeatmap(lng = ~reclong, lat = ~reclat, intensity = ~mass,
             blur = 20, max = 0.05, radius = 15) %>%
  # controls
  setView(lng=center_lon, lat=center_lat,zoom = 4) 

36.3 Distribution of Aus Meteorite Landings with Cluster Markers

The meteorite landings have been clustered and their numbers are being shown in the map.

AusMetLandings %>% leaflet() %>% addProviderTiles("Esri.OceanBasemap") %>% 
  addMarkers(lng = ~reclong, lat = ~reclat,clusterOptions = markerClusterOptions()) %>%
    # controls
  setView(lng=center_lon, lat=center_lat,zoom = 4)