diff --git a/README.html b/README.html index e13a2cd..9266ca7 100644 --- a/README.html +++ b/README.html @@ -69,7 +69,7 @@ - + @@ -553,7 +553,7 @@
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"
\\n\"+\n \"
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
The Kernel crashed while executing code in the the current cell or a previous cell. Please review the code in the cell(s) to identify a possible cause of the failure. Click <a href='https://aka.ms/vscodeJupyterKernelCrash'>here</a> for more info. View Jupyter <a href='command:jupyter.viewOutput'>log</a> for further details.
+
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
Earthquake Intensity Data was obtained from USGS and was overlayed with shapefiles of Morocco admin regions taken from UNOCHA and from GADM[https://gadm.org/download_country.html]. This can be used as a baseline map upon which other indicators can be visualized.
+Earthquake Intensity Data was obtained from USGS and was overlayed with shapefiles of Morocco admin regions taken from UNOCHA and from GADM[https://gadm.org/download_country.html]. This can be used as a baseline map upon which other indicators can be visualized. The MMI is a range that allows us to estimate the amount of damage done to buildings.
The impact of the earthquake in this situation is measured using the Modified Mercelli Intensity Scale (mmi).
+The impact of the earthquake in this situation is measured using the Modified Mercalli Intensity Scale (mmi).
The effect of an earthquake on the Earth’s surface is called the intensity. The intensity scale consists of a series of certain key responses such as people awakening, movement of furniture, damage to chimneys, and finally - total destruction. Although numerousintensity scales have been developed over the last several hundred years to evaluate the effects of earthquakes, the one currently used in the United States is the Modified Mercalli (MM) Intensity Scale. It was developed in 1931 by the American seismologists Harry Wood and Frank Neumann. This scale, composed of increasing levels of intensity that range from imperceptible shaking to catastrophic destruction, is designated by Roman numerals. It does not have a mathematical basis; instead it is an arbitrary ranking based on observed effects.
The Modified Mercalli Intensity value assigned to a specific site after an earthquake has a more meaningful measure of severity to the nonscientist than the magnitude because intensity refers to the effects actually experienced at that place.
The lower numbers of the intensity scale generally deal with the manner in which the earthquake is felt by people. The higher numbers of the scale are based on observed structural damage. Structural engineers usually contribute information for assigning intensity values of VIII or above.
-Source: USGS Modified Mercelli Intensity Scale
+Source: USGS Modified Mercalli Intensity Scale
Alternatively, another variable that can be used is also the distance from the epicenter.
+
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
- Last updated on Oct 20, 2023.
+ Last updated on Oct 24, 2023.
The building footprints are obtained from satellite imagery which is obstructed by cloud cover. Given Morocco is a dessert area, the results are somewhere near 80% accurate. Data from OpenStreetMaps is crowdsourced.
- -Text(-0.35, -0.05, 'Source: USGS')
+
percentage_tourist_damages
+
{'Al Haouz': {'No': 71.30801687763713, 'Yes': 28.69198312236287},
+ 'Chichaoua': {'No': 80.0, 'Yes': 20.0},
+ 'Taroudannt': {'No': 96.42857142857143, 'Yes': 3.5714285714285716},
+ 'Marrakech': {'No': 99.88331388564761, 'Yes': 0.11668611435239207}}
+
percentage_population=round(100*financial_services['population_without_children'].sum()/total_population_affected_region['population_wihtou_children'].sum())
+print(f'{percentage_population}% of people in the affected regions lost access to financial services after the earthquake. These include {financial_service_types}')
+
40% of people in the affected regions lost access to financial services after the earthquake. These include ['bank' 'post_office' 'bureau_de_change' 'atm' 'money_transfer']
+
fig, ax = plt.subplots()
+affected_regions.plot(cmap='Reds', column='PARAMVALUE', ax=ax)
+financial_services.to_crs(epsg=4326).plot(column='population_without_children', cmap='Blues', ax=ax, legend=True)
-eq_adm4.plot(column='PARAMVALUE', ax=ax, cmap = 'Reds', alpha = 0.5)
-affected_settlement.plot(ax=ax)
-financial_services.plot(color='purple', ax=ax, label = 'financial services')
-tourist_attractions.plot(color = 'green', ax=ax)
+plt.title('Total population within 5km radius of financial services', fontsize=12, ha='center')
-legend_labels = ['financial services', 'tourist attractions']
-handles = [Line2D([0], [0], marker='o', color='w', label=label, markerfacecolor=color_dict[label]) for label in legend_labels]
-ax.legend(handles=handles, loc='upper left', frameon=False)
+ax.spines['top'].set_visible(False)
+ax.spines['bottom'].set_visible(False)
+ax.spines['right'].set_visible(False)
+ax.spines['left'].set_visible(False)
-plt.axis('off')
+ax.set_xticks([])
+ax.set_yticks([])
-plt.title('Settlements in worst affected Admin 1 regions overlayed with Points of Inetrest', fontsize=20, ha='center')
-subtitle_text = "Source: Settlement layer from Atlas AI, earthquake intensity from USGS and Points of Interest from OpenStreetMaps"
-ax.text(-0.35, -0.05, subtitle_text, transform=ax.transAxes, fontsize=15, ha='left')
+subtitle_text = "Source: Population estimates from Meta's HRSL and financial points from OSM"
+ax.text(-0.35, -0.05, subtitle_text, transform=ax.transAxes, fontsize=10, ha='left')
Text(-0.35, -0.05, 'Source: Settlement layer from Atlas AI, earthquake intensity from USGS and Points of Interest from OpenStreetMaps')
+Text(-0.35, -0.05, "Source: Population estimates from Meta's HRSL and financial points from OSM")
-
+