Mapping DHS with leaflet.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "#  <span style=\"color:#3386cd; font-size:28px\">Mapping DHS with leaflet </span>\n",
    "\n",
    "*Version 1.0 August 24th, 2016*\n",
    "#### &#9733;<i>Thomas Roca, PhD - AFD Research Direction &#9733;</i>\n",
    "---\n",
    "\n",
    "The Demographic Health Surveys (DHS) are household surveys covering most of developping countries (and more) in the areas of population, health, and nutrition. These surveys are supported by <a href=\"https://www.usaid.gov/\" target=\"_blank\">USAID</a> and the <a href=\"https://www.census.gov/\" target=\"_blank\"> U.S. Census Bureau</a> (Kudos!)\n",
    "More information about DHS: http://dhsprogram.com/\n",
    "\n",
    "This script uses geolocated dataset to map DHS Data at the subnational level using <a href=\"http://leafletjs.com/\" target=\"_blank\">leaflet</a>.\n",
    "More information about geotagged DHS data: http://spatialdata.dhsprogram.com/home/\n",
    "\n",
    "\n",
    "\n",
    "![](http://spatialdata.dhsprogram.com/global/images/logo.png)\n",
    "\n",
    "---\n",
    "\n",
    "Packages used: fiona, json"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Meta data as presented on : http://spatialdata.dhsprogram.com/data/#/common/download\n",
    "\n",
    "Countries: AL, AO, AM, AZ, BD, BJ, BO, BT, BR, BF, BU, KH, CM, CF, TD, CO, KM, CG, CD, CI, DR, EC, EG, ES, ER, ET, GA, GH, GU, GN, GY, HT, HN, IA, ID, JO, KK, KE, KY, LS, LB, MD, MW, MV, ML, MR, MX, MB, MA, MZ, NM, NP, NC, NI, NG, PK, PY, PE, PH, RW, ST, SN, SL, ZA, LK, SD, SZ, TJ, TZ, TH, GM, TL, TG, TT, TN, TR, TM, UG, UA, UZ, VN, YE, ZM, ZW\n",
    "\n",
    "### Indicators:\n",
    "- Households with electricity (i7C3CB1)\n",
    "- Women with secondary or higher education (iA999FA)\n",
    "- Women who are literate (i797B5DA)\n",
    "- Total fertility rate (i133C8F8)\n",
    "- Married women currently using any method of contraception (i209E190)\n",
    "- Married women currently using any modern method of contraception (i209E191)\n",
    "- Median age at first marriage: 25-49 (i349C430)\n",
    "- Median age at first sexual intercourse: 25-49 (i36848B0)\n",
    "- Unmet need total (iCE37FE2)\n",
    "- Infant mortality rate (i42FFDB2)\n",
    "- Under-five mortality rate (i42FFDB4)\n",
    "- Place of delivery: Health facility (i49B3AD0)\n",
    "- Received all vaccinations (i4E75F0A)\n",
    "- Treatment of diarrhea: Either ORS or RHS (i533BDCB)\n",
    "- Median duration of exclusive breastfeeding (i5716309)\n",
    "- Children stunted (iCC96C1A)\n",
    "- Children underweight (iCC97002)\n",
    "- Children wasted (iCC973EA)\n",
    "- Men receiving an HIV test and receiving test results in the last 12 months (i42759491)\n",
    "- Women receiving an HIV test and receiving test results in the last 12 months (i42759492)\n",
    "- HIV prevalence among men (i41075E91)\n",
    "- HIV prevalence among women (i41075E92)\n",
    "- HIV prevalence among general population (i41075E93)\n",
    "- Children under 5 who slept under an insecticide treated net (ITN) (iC6D4E13)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Exemple of map produced: [link](http://afd.countrydashboards.com/dataviz/local/Map_i42FFDB4_BFA.html)\n",
    "\n",
    "[![image](http://afd.countrydashboards.com/dataviz/local/BFA.png)](http://afd.countrydashboards.com/dataviz/local/Map_i42FFDB4_BFA.html)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1. Get global shapefile:\n",
    "\n",
    "Download the [25 most commons indicators](http://spatialdata.dhsprogram.com/data/#/common/download) for all countries such as:\n",
    "- Geography Level: subnational\n",
    "- Data format: Shapefile\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2. Convert the shapefile into geojson"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "folder='C:/Users/rocat/OneDrive/Bureau/DHS/Global/shps/'\n",
    "file_name='sdr_subnational_data_quickstats'\n",
    "\n",
    "#We use fiona package to convert the shapefile into json\n",
    "with fiona.open(folder+file_name+'.shp') as source:\n",
    "    records = list(source)\n",
    "    geojson = {\"type\": \"FeatureCollection\",\"features\": records}\n",
    "\n",
    "with open(file_name+\".json\", \"w\") as f:\n",
    "    f.write(json.dumps(geojson))\n",
    "    f = open(file_name+\".json\",'r')\n",
    "    filedata = f.read()\n",
    "    f.close()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Parse the global file DHS by country"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Load global json"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import json\n",
    "filename='C:/Users/rocat/OneDrive/Bureau/DHS/Global/shps/sdr_subnational_data_quickstats.json'\n",
    "json_data=open(filename).read()\n",
    "data = json.loads(json_data)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create national json and js files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of region available:  750\n"
     ]
    }
   ],
   "source": [
    "from string import Template\n",
    "from IPython.display import HTML\n",
    "\n",
    "# ISO 3166-1 alpha-2 list of countries\n",
    "isolist=['AL','AO','AM','AZ','BD','BJ','BO','BT','BR','BF','BU','KH','CM','CF','TD','CO','KM','CG','CD','CI','DR','EC','EG','ES','ER','ET','GA','GH','GU','GN','GY','HT','HN','IA','ID','JO','KK','KE','KY','LS','LB','MD','MW','MV','ML','MR','MX','MB','MA','MZ','NM','NP','NC','NI','NG','PK','PY','PE','PH','RW','ST','SN','SL','ZA','LK','SD','SZ','TJ','TZ','TH','GM','TL','TG','TT','TN','TR','TM','UG','UA','UZ','VN','YE','ZM','ZW']\n",
    "\n",
    "# len(data[\"features\"]) will count the overall number of indicators (one by DHS region for each country)\n",
    "print(\"Number of region available: \", len(data[\"features\"]))\n",
    "\n",
    "# loop over isolist\n",
    "for iso3 in isolist:\n",
    "    \n",
    "    #create en empty list\n",
    "    vars()['json_'+str(iso3)]=[]\n",
    "    \n",
    "    #loop over features lenght\n",
    "    for features in range (0,len(data[\"features\"])-1):  \n",
    "        iso=data[\"features\"][features][\"properties\"][\"ISO\"]\n",
    "        content=data[\"features\"][features]    \n",
    "        if iso3==iso: \n",
    "            vars()['json_'+str(iso3)].append(content)\n",
    "        \n",
    "        #write extra info used in leaflet\n",
    "        string='{\"type\": \"FeatureCollection\", \"features\":'+ str(vars()['json_'+str(iso3)]) + '}'\n",
    "        #we replace None by string None\n",
    "        string=string.replace('None', '\"None\"')\n",
    "        #some information are coded using single quote, we replace the single quote by double\n",
    "        string=string.replace(\"'\",'\"')\n",
    "        #Save as json to be able to play with it in Python\n",
    "        with open('C:/Users/rocat/OneDrive/Bureau/DHS/local/'+iso3+\".json\", \"w\") as f: f.write(string)\n",
    "        \n",
    "        #Save as js to map them with leaflet\n",
    "        string='var statesData ='+string\n",
    "        with open('C:/Users/rocat/OneDrive/Bureau/DHS/local/'+iso3+\".js\", \"w\") as f: f.write(string)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## List missing values\n",
    "\n",
    "In this section we read all json files and trace indicators for which values are missing (nb. missing values are coded 9999)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import json\n",
    "folder='C:/Users/rocat/OneDrive/Bureau/DHS/local/'\n",
    "isolist=['AL','AO','AM','AZ','BD','BJ','BO','BT','BR','BF','BU','KH','CM','CF','TD','CO','KM','CG','CD','CI','DR','EC','EG','ES','ER','ET','GA','GH','GU','GN','GY','HT','HN','IA','ID','JO','KK','KE','KY','LS','LB','MD','MW','MV','ML','MR','MX','MB','MA','MZ','NM','NP','NC','NI','NG','PK','PY','PE','PH','RW','ST','SN','SL','ZA','LK','SD','SZ','TJ','TZ','TH','GM','TL','TG','TT','TN','TR','TM','UG','UA','UZ','VN','YE','ZM','ZW']\n",
    "\n",
    "# We use the indicators dictionnary that contains: the indicator code (e.g. i7C3CB1) the definition and \n",
    "# arbitrary min and max we set to define the range\n",
    " \n",
    "ind_dict={\"i7C3CB1\":\"Households with electricity,0,100\",\"iA999FA\":\"Women with secondary or higher education,0,100\",\\\n",
    "          \"i797B5DA\":\"Women who are literate,0,100\",\"i133C8F8\":\"Total fertility rate,0,10\",\\\n",
    "          \"i209E190\":\"Married women currently using any method of contraception,0,100\",\\\n",
    "          \"i209E191\":\"Married women currently using any modern method of contraception,0,100\",\\\n",
    "          \"i349C430\":\"Median age at first marriage:25-49,12,35\",\\\n",
    "          \"i36848B0\":\"Median age at first sexual intercourse:25-49,12,35\",\\\n",
    "          \"iCE37FE2\":\"Unmet need total,0,100\",\\\n",
    "          \"i42FFDB2\":\"Infant mortality rate,0,300\",\\\n",
    "          \"i42FFDB4\":\"Under-five mortality rate,0,150\",\\\n",
    "          \"i49B3AD0\":\"Place of delivery: Health facility,0,300\",\\\n",
    "          \"i4E75F0A\":\"Received all vaccinations,0,100\",\\\n",
    "          \"i533BDCB\":\"Treatment of diarrhea: Either ORS or RHS,0,100\",\\\n",
    "          \"i5716309\":\"Median duration of exclusive breastfeeding,0,100\",\\\n",
    "          \"iCC96C1A\":\"Children stunted,0,100\",\\\n",
    "          \"iCC97002\":\"Children underweight,0,100\",\\\n",
    "          \"iCC973EA\":\"Children wasted,0,100\",\\\n",
    "          \"i42759491\":\"Men receiving an HIV test and receiving test results in the last 12 months,0,100\",\\\n",
    "          \"i42759492\":\"Women receiving an HIV test and receiving test results in the last 12 months,0,100\",\\\n",
    "          \"i41075E91\":\"HIV prevalence among men,0,100\",\\\n",
    "          \"i41075E92\":\"HIV prevalence among women,0,100\",\\\n",
    "          \"i41075E93\":\"HIV prevalence among general population,0,100\",\\\n",
    "          \"iC6D4E13\":\"Children under 5 who slept under an insecticide treated net,0,100\"}\n",
    "\n",
    "#create en empty list to host the missing values list\n",
    "string=[]\n",
    "\n",
    "# Loop over countries\n",
    "for iso in isolist:\n",
    "    #open json files\n",
    "    json_data=open(folder+iso+'.json').read()\n",
    "    data0 = json.loads(json_data)\n",
    "    \n",
    "    #loop over dictionnary values\n",
    "    for code,info_ind in ind_dict.items():\n",
    "        item=info_ind.split(',')\n",
    "        definition=item[0]\n",
    "        \n",
    "        #Read the features\n",
    "        for features in range (0,len(data0[\"features\"])-1):\n",
    "            Region=data0[\"features\"][features][\"properties\"][\"DHSREGEN\"]\n",
    "            if data0[\"features\"][features][\"properties\"][code]==9999 :\n",
    "                #Append the list of indicators\n",
    "                string.append(iso+','+Region+','+code+','+definition+',no data\\n')\n",
    "                \n",
    "# Finally Write the missing values list into a csv file \n",
    "with open('C:/Users/rocat/OneDrive/Bureau/DHS/local/missing.csv', \"w\") as f: f.write(''.join(map(str,string)))\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create Maps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import json\n",
    "from string import Template\n",
    "from IPython.display import HTML\n",
    "\n",
    "isolist=['AL','AO','AM','AZ','BD','BJ','BO','BT','BR','BF','BU','KH','CM','CF','TD','CO','KM','CG','CD','CI','DR','EC','EG','ES','ER','ET','GA','GH','GU','GN','GY','HT','HN','IA','ID','JO','KK','KE','KY','LS','LB','MD','MW','MV','ML','MR','MX','MB','MA','MZ','NM','NP','NC','NI','NG','PK','PY','PE','PH','RW','ST','SN','SL','ZA','LK','SD','SZ','TJ','TZ','TH','GM','TL','TG','TT','TN','TR','TM','UG','UA','UZ','VN','YE','ZM','ZW']\n",
    "\n",
    "Indicator_list=[\"i7C3CB1\",\"i797B5DA\",\"i209E190\",\"i209E191\",\"i349C430\",\"i36848B0\",\"iCE37FE2\",\"i42FFDB2\",\"i42FFDB4\",\\\n",
    "                \"i49B3AD0\",\"i4E75F0A\",\"i533BDCB\",\"i5716309\",\"iCC96C1A\",\"iCC97002\",\"iCC973EA\",\"i42759491\",\"i42759492\",\\\n",
    "                \"i41075E91\",\"i41075E92\",\"i41075E93\",\"iC6D4E13\"]\n",
    "\n",
    "Indicator_list=[\"i42FFDB4\"]\n",
    "\n",
    "\n",
    "ind_dict={\"i7C3CB1\":\"Households with electricity,0,100\",\"iA999FA\":\"Women with secondary or higher education,0,100\",\\\n",
    "          \"i797B5DA\":\"Women who are literate,0,100\",\"i133C8F8\":\"Total fertility rate,0,10\",\\\n",
    "          \"i209E190\":\"Married women currently using any method of contraception,0,100\",\\\n",
    "          \"i209E191\":\"Married women currently using any modern method of contraception,0,100\",\\\n",
    "          \"i349C430\":\"Median age at first marriage:25-49,12,35\",\\\n",
    "          \"i36848B0\":\"Median age at first sexual intercourse:25-49,12,35\",\\\n",
    "          \"iCE37FE2\":\"Unmet need total,0,100\",\\\n",
    "          \"i42FFDB2\":\"Infant mortality rate,0,300\",\\\n",
    "          \"i42FFDB4\":\"Under-five mortality rate,0,300\",\\\n",
    "          \"i49B3AD0\":\"Place of delivery: Health facility,0,100\",\\\n",
    "          \"i4E75F0A\":\"Received all vaccinations,0,100\",\\\n",
    "          \"i533BDCB\":\"Treatment of diarrhea: Either ORS or RHS,0,100\",\\\n",
    "          \"i5716309\":\"Median duration of exclusive breastfeeding,0,100\",\\\n",
    "          \"iCC96C1A\":\"Children stunted,0,100\",\\\n",
    "          \"iCC97002\":\"Children underweight,0,100\",\\\n",
    "          \"iCC973EA\":\"Children wasted,0,100\",\\\n",
    "          \"i42759491\":\"Men receiving an HIV test and receiving test results in the last 12 months,0,100\",\\\n",
    "          \"i42759492\":\"Women receiving an HIV test and receiving test results in the last 12 months,0,100\",\\\n",
    "          \"i41075E91\":\"HIV prevalence among men,0,100\",\\\n",
    "          \"i41075E92\":\"HIV prevalence among women,0,100\",\\\n",
    "          \"i41075E93\":\"HIV prevalence among general population,0,100\",\\\n",
    "          \"iC6D4E13\":\"Children under 5 who slept under an insecticide treated net,0,100\"}\n",
    "\n",
    "# We created a GPS dictionnary to store: the gps location of the capital city to center the map so as the zoom level\n",
    "gps_dict={\"AL\":\"ALB,41.1533,20.1683,7\", \"AO\":\"AGO,-11.2027,17.8739,5\", \"AM\":\"ARM,40.0691,45.0382,7\", \"AZ\":\"AZE,40.1431,47.5769,7\",\\\n",
    "          \"BD\":\"BGD,23.685,90.3563,7\", \"BJ\":\"BEN,9.30769,2.31583,6\", \"BO\":\"BOL,-16.2902,-63.5887,5\", \"BT\":\"BTN,27.5142,90.4336,8\",\\\n",
    "          \"BR\":\"BRA,-14.235,-51.9253,4\", \"BF\":\"BFA,12.2383,-1.56159,6\", \"BU\":\"null,-32.87,26.12,6\", \"KH\":\"KHM,12.5657,104.991, 7\",\\\n",
    "          \"CM\":\"CMR,7.36972,12.3547, 5\", \"CF\":\"CAF,6.61111,20.9394,6\", \"TD\":\"TCD,15.4542,18.7322,5\", \"CO\":\"COL,4.57087,-74.2973,5\",\\\n",
    "          \"KM\":\"COM,-11.875,43.8722,8\", \"CG\":\"COG,-0.228021,15.8277,6\", \"CD\":\"COD,-4.03833,21.7587,5\", \"CI\":\"CIV,7.53999,-5.54708,6\",\\\n",
    "          \"DR\":\"null,-32.87,26.12,7\", \"EC\":\"ECU,-1.83124,-78.1834,6\", \"EG\":\"EGY,26.8206,30.8025,5\", \"ES\":\"ESP,40.4637,-3.74922,6\",\\\n",
    "          \"ER\":\"ERI,15.1794,39.7823,7\", \"ET\":\"ETH,9.145,40.4897,5\", \"GA\":\"GAB,-0.803689,11.6094,6\", \"GH\":\"GHA,7.94653,-1.02319,6\",\\\n",
    "          \"GU\":\"GUM,13.4443,144.794,9\", \"GN\":\"GIN,9.94559,-9.69664,6\", \"GY\":\"GUY,4.86042,-58.9302,6\", \"HT\":\"HTI,18.9712,-72.2852,8\",\\\n",
    "          \"HN\":\"HND,15.2,-86.2419,7\", \"IA\":\"null,-32.87,26.12,7\", \"ID\":\"IDN,-0.789275,113.921,4\", \"JO\":\"JOR,30.5852,36.2384,7\",\\\n",
    "          \"KK\":\"null,-32.87,26.12,7\", \"KE\":\"KEN,-0.023559,37.9062,6\", \"KY\":\"CYM,19.5135,-80.567,9\", \"LS\":\"LSO,-29.61,28.2336,8\",\\\n",
    "          \"LB\":\"LBN,33.8547,35.8623,8\", \"MD\":\"MDA,47.4116,28.3699,7\", \"MW\":\"MWI,-13.2543,34.3015,6\", \"MV\":\"MDV,3.20278,73.2207,7\",\\\n",
    "          \"ML\":\"MLI,17.5707,-3.99617,5\", \"MR\":\"MRT,21.0079,-10.9408,5\", \"MX\":\"MEX,23.6345,-102.553,5\", \"MB\":\"null,-32.87,26.12,7\",\\\n",
    "          \"MA\":\"MAR,31.7917,-7.09262,5\", \"MZ\":\"MOZ,-18.6657,35.5296,5\", \"NM\":\"null,-32.87,26.12,7\", \"NP\":\"NPL,28.3949,84.124,7\",\\\n",
    "          \"NC\":\"NCL,-20.9043,165.618,7\", \"NI\":\"NIC,12.8654,-85.2072,7\", \"NG\":\"NGA,9.082,8.67528,6\", \"PK\":\"PAK,30.3753,69.3451,5\",\\\n",
    "          \"PY\":\"PRY,-23.4425,-58.4438,6\", \"PE\":\"PER,-9.18997,-75.0152,5\", \"PH\":\"PHL,12.8797,121.774,5\", \"RW\":\"RWA,-1.94028,29.8739,8\",\\\n",
    "          \"ST\":\"STP,0.18636,6.61308,8\", \"SN\":\"SEN,14.4974,-14.4524,7\", \"SL\":\"SLE,8.46056,-11.7799,7\", \"ZA\":\"ZAF,-30.5595,22.9375,5\",\\\n",
    "          \"LK\":\"LKA,7.87305,80.7718,7\", \"SD\":\"SDN,12.8628,30.2176,5\", \"SZ\":\"SWZ,-26.5225,31.4659,8\", \"TJ\":\"TJK,38.861,71.2761,6\",\\\n",
    "          \"TZ\":\"TZA,-6.36903,34.8888,5\", \"TH\":\"THA,15.87,100.993,5\", \"GM\":\"GMB,13.4432,-15.3101,8\", \"TL\":\"TLS,-8.87422,125.728,8\",\\\n",
    "          \"TG\":\"TGO,8.61954,0.824782,7\", \"TT\":\"TTO,10.6918,-61.2225,8\", \"TN\":\"TUN,33.8869,9.5375,6\", \"TR\":\"TUR,38.9637,35.2433,6\",\\\n",
    "          \"TM\":\"TKM,38.9697,59.5563,7\", \"UG\":\"UGA,1.37333,32.2903,7\", \"UA\":\"UKR,48.3794,31.1656,6\", \"UZ\":\"UZB,41.3775,64.5853,6\",\\\n",
    "          \"VN\":\"VNM,14.0583,108.277,6\", \"YE\":\"YEM,15.5527,48.5164,6\", \"ZM\":\"ZMB,-13.1339,27.8493,6\", \"ZW\":\"ZWE,-19.0154,29.1549,6\"}\n",
    "\n",
    "\n",
    "for indicator in Indicator_list:\n",
    "    # save indicator definition\n",
    "    for code,info_ind in ind_dict.items():\n",
    "        if code==indicator: \n",
    "            item=info_ind.split(',')\n",
    "            definition=item[0]\n",
    "            min_ind=int(item[1])\n",
    "            max_ind=int(item[2])\n",
    "\n",
    "    #Set color range\n",
    "    ind_scale=(max_ind-min_ind)/8\n",
    "    cat0=0\n",
    "    for val in range(1,7+1):\n",
    "        vars()['cat'+str(val)]=round(ind_scale+vars()['cat'+str(val-1)])\n",
    "\n",
    "\n",
    "    # save GPS coordinate to center the map\n",
    "    for iso in isolist:    \n",
    "        for iso2,info in gps_dict.items():\n",
    "            if iso==iso2:\n",
    "                item=info.split(',')\n",
    "                iso3=item[0]\n",
    "                if iso3==\"null\": iso3=iso2\n",
    "                Longitude=item[1]\n",
    "                Latitude=item[2]\n",
    "                zoomlevel=item[3]\n",
    "                #print(iso2,iso3,Latitude,Longitude)\n",
    "\n",
    "        Input = {'iso3':iso3,'iso':iso, 'definition':definition, 'iso2':iso2, 'Latitude':Latitude,'Longitude':Longitude, \\\n",
    "                 'zoomlevel':zoomlevel, 'indicator': indicator,\\\n",
    "                 'cat1':cat1, 'cat2':cat2, 'cat3':cat3,'cat4':cat4,'cat5':cat5,'cat6':cat6, 'cat7':cat7 }   \n",
    "        #html file content\n",
    "        html='''\n",
    "            <!DOCTYPE html>\n",
    "            <html><head>\n",
    "                    <title>Leaflet Layers Control Example</title>\n",
    "                    <meta charset=\"utf-8\" />\n",
    "                    <meta name=\"viewport\" content=\"width=device-width, initial-scale=1.0\">\n",
    "                    <link rel=\"stylesheet\" href=\"https://npmcdn.com/leaflet@1.0.0-rc.3/dist/leaflet.css\" />\n",
    "                    <style>\n",
    "                        #map {\n",
    "                            width: 800px;\n",
    "                            height: 500px;\n",
    "                        }\n",
    "\n",
    "                        .info {\n",
    "                            padding: 6px 8px;\n",
    "                            font: 14px/16px Arial, Helvetica, sans-serif;\n",
    "                            background: white;\n",
    "                            background: rgba(255,255,255,0.8);\n",
    "                            box-shadow: 0 0 15px rgba(0,0,0,0.2);\n",
    "                            border-radius: 5px;\n",
    "                        }\n",
    "                        .info h4 {\n",
    "                            margin: 0 0 5px;\n",
    "                            color: #777;\n",
    "                        }\n",
    "\n",
    "                        .legend {\n",
    "                            text-align: left;\n",
    "                            line-height: 18px;\n",
    "                            color: #555;\n",
    "                        }\n",
    "                        .legend i {\n",
    "                            width: 18px;\n",
    "                            height: 18px;\n",
    "                            float: left;\n",
    "                            margin-right: 8px;\n",
    "                            opacity: 0.7;\n",
    "                        }\n",
    "                    </style>\n",
    "                </head>\n",
    "                <body>\n",
    "                    <div id=\"map\"></div>\n",
    "\n",
    "                    <script src=\"https://npmcdn.com/leaflet@1.0.0-rc.3/dist/leaflet.js\"></script>\n",
    "\n",
    "                    <script type=\"text/javascript\" src=\"$iso.js\"></script>\n",
    "                    <script type=\"text/javascript\">\n",
    "\n",
    "                        var map = L.map('map').setView([$Longitude,$Latitude], $zoomlevel);\n",
    "\n",
    "                        L.tileLayer('https://api.tiles.mapbox.com/v4/{id}/{z}/{x}/{y}.png?access_token=pk.eyJ1IjoibWFwYm94IiwiYSI6ImNpandmbXliNDBjZWd2M2x6bDk3c2ZtOTkifQ._QA7i5Mpkd_m30IGElHziw', {\n",
    "                            maxZoom: 18,\n",
    "                            attribution: 'Map data &copy; <a href=\"http://openstreetmap.org\">OpenStreetMap</a> contributors, ' +\n",
    "                                '<a href=\"http://creativecommons.org/licenses/by-sa/2.0/\">CC-BY-SA</a>, ' +\n",
    "                                'Dataviz from <a href=\"https://twitter.com/Thomas_Roca\">@Thomas_Roca </a> with Python',\n",
    "                            id: 'mapbox.light'\n",
    "                        }).addTo(map);\n",
    "\n",
    "\n",
    "                        // control that shows state info on hover\n",
    "                        var info = L.control();\n",
    "\n",
    "                        info.onAdd = function (map) {\n",
    "                            this._div = L.DomUtil.create('div', 'info');\n",
    "                            this.update();\n",
    "                            return this._div;\n",
    "                        };\n",
    "\n",
    "                        info.update = function (props) {\n",
    "\n",
    "                           this._div.innerHTML = (props ?\n",
    "                           '<h3>'+props.CNTRYNAMEE+'</h3>' +\n",
    "                            '<span style=\"color:#336DA9\"><b>$definition: '+\n",
    "                            '<span style=\"color:#0e7cea\">' + props.$indicator +'</span></b><br><span style=\"font-size:12px; color:black\">'+\n",
    "                            'Survey Region: <b>' + props.DHSREGEN + '</b><br>' +\n",
    "                            'Survey Year: <b> '+props.SVYYEAR + '</b><br>' +\n",
    "                            '</b> <i>Source: <b>DHS program</i></b></a>'\n",
    "                            : '<b>$definition</b>');\n",
    "                        };\n",
    "\n",
    "                        info.addTo(map);\n",
    "\n",
    "\n",
    "                        // get color depending on value\n",
    "                        function getColor(d) {\n",
    "                            return d==9999 ? '#EEEEEE' :\n",
    "                                   d > $cat7 ? '#1a334c' :\n",
    "                                   d > $cat6  ? '#084594' :\n",
    "                                   d > $cat5  ? '#2171b5' :\n",
    "                                   d > $cat4  ? '#4292c6' :\n",
    "                                   d > $cat3   ? '#6baed6' :\n",
    "                                   d > $cat2   ? '#9ecae1' :\n",
    "                                   d > $cat1   ? '#c6dbef' :\n",
    "                                              '#eff3ff';\n",
    "                        }\n",
    "\n",
    "                        function style(feature) {\n",
    "                            return {\n",
    "                                weight: 2,\n",
    "                                opacity: 1,\n",
    "                                color: 'white',\n",
    "                                dashArray: '3',\n",
    "                                fillOpacity: 0.7,\n",
    "                                fillColor: getColor(feature.properties.$indicator)\n",
    "                            };\n",
    "                        }\n",
    "\n",
    "                        function highlightFeature(e) {\n",
    "                            var layer = e.target;\n",
    "\n",
    "                            layer.setStyle({\n",
    "                                weight: 5,\n",
    "                                color: '#666',\n",
    "                                dashArray: '',\n",
    "                                fillOpacity: 0.7\n",
    "                            });\n",
    "\n",
    "                            if (!L.Browser.ie && !L.Browser.opera) {\n",
    "                                layer.bringToFront();\n",
    "                            }\n",
    "\n",
    "                            info.update(layer.feature.properties);\n",
    "                        }\n",
    "\n",
    "                        var geojson;\n",
    "\n",
    "                        function resetHighlight(e) {\n",
    "                            geojson.resetStyle(e.target);\n",
    "                            info.update();\n",
    "                        }\n",
    "\n",
    "                        function zoomToFeature(e) {\n",
    "                            map.fitBounds(e.target.getBounds());\n",
    "                        }\n",
    "\n",
    "                        function onEachFeature(feature, layer) {\n",
    "                            layer.on({\n",
    "                                mouseover: highlightFeature,\n",
    "                                mouseout: resetHighlight,\n",
    "                                click: zoomToFeature\n",
    "                            });\n",
    "                        }\n",
    "\n",
    "                        geojson = L.geoJson(statesData, {\n",
    "                            style: style,\n",
    "                            onEachFeature: onEachFeature\n",
    "                        }).addTo(map);\n",
    "\n",
    "                        map.attributionControl.addAttribution('Data &copy; <a href=\"http://spatialdata.dhsprogram.com/\" target=\"_blank\">DHS Spatial Data Repository</a>');\n",
    "\n",
    "\n",
    "                        var legend = L.control({position: 'bottomright'});\n",
    "\n",
    "                        legend.onAdd = function (map) {\n",
    "\n",
    "                            var div = L.DomUtil.create('div', 'info legend'),\n",
    "                                grades = [$cat1,$cat2,$cat3,$cat4,$cat5,$cat6,$cat7],\n",
    "                                labels = [],\n",
    "                                from, to;\n",
    "\n",
    "                            for (var i = 0; i < grades.length; i++) {\n",
    "                                from = grades[i];\n",
    "                                to = grades[i + 1];\n",
    "\n",
    "                                labels.push(\n",
    "                                    '<i style=\"background:' + getColor(from + 1) + '\"></i> ' +\n",
    "                                    from + (to ? '&ndash;' + to : '+'));\n",
    "                            }\n",
    "\n",
    "                            div.innerHTML = labels.join('<br>');\n",
    "                            return div;\n",
    "                        };\n",
    "\n",
    "                        legend.addTo(map);\n",
    "\n",
    "                    </script>\n",
    "            </body>\n",
    "            </html>\n",
    "            '''\n",
    "\n",
    "\n",
    "        #open an handle\n",
    "        f = open(\"C:/Users/rocat/OneDrive/Bureau/DHS/local/Map_\"+indicator+\"_\"+iso3+\".html\",'w')\n",
    "        #write the content within the handle\n",
    "        file_content = Template(html).safe_substitute(Input)\n",
    "        f.write(file_content)\n",
    "        f.close()\n",
    "\n",
    "#=====================================================================================================================\n",
    "#                                         ! Done\n",
    "#====================================================================================================================="
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.4.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}