Sentinel2b_test.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys\n",
    "import ee\n",
    "%run Sentinel2b.ipynb\n",
    "#ee.Authenticate()\n",
    "ee.Initialize()\n",
    "\n",
    "modulename = 'ee'\n",
    "if modulename not in sys.modules:\n",
    "    # import GEE and Authenticate, token or log in will be asked from web browser\n",
    "    import ee\n",
    "    # TO DO : UNCOMMENT THIS IF IT IS THE FIRST TIME AUTHINTICATE GEE ON THIS MACHINE\n",
    "    # ee.Authenticate()\n",
    "    ee.Initialize()\n",
    "# else:\n",
    "    # google earth engine already imported and authenticated\n",
    "\n",
    "\n",
    "modulename = 'ipynb_MapVisualisation'\n",
    "if modulename not in sys.modules:\n",
    "    %run MapVisualisation.ipynb\n",
    "    # adding an identifier to sys.modules to avoiding loading the same file multiple times\n",
    "    sys.modules['ipynb_MapVisualisation'] = None\n",
    "# else\n",
    "    # Utils modules has already been loaded somewhere else\n",
    "\n",
    "\n",
    "modulename = 'ipynb_Sentinel2b'\n",
    "if modulename not in sys.modules:\n",
    "    %run Sentinel2b.ipynb\n",
    "    # adding an identifier to sys.modules to avoiding loading the same file multiple times\n",
    "    sys.modules['ipynb_Sentinel2b'] = None\n",
    "# else\n",
    "    # Utils modules has already been loaded somewhere else\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "%run MapVisualisation.ipynb\n",
    "# These are the coordinates of Spain\n",
    "my_map = folium.Map(location=[42, -3], zoom_start=11)\n",
    "# Add custom basemaps\n",
    "basemaps['Google Maps'].add_to(my_map)\n",
    "basemaps['Google Satellite Hybrid'].add_to(my_map)\n",
    "\n",
    "\n",
    "def update_cld_shdw_mask(img):\n",
    "    print('test')\n",
    "\n",
    "\n",
    "countries = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017')\n",
    "AOI = countries.filter(ee.Filter.eq('country_na', 'Spain'))\n",
    "\n",
    "AOI = ee.Geometry.Polygon(\n",
    "        [[[-3.2242859153727754, 40.999456684741496],\n",
    "          [-3.2242859153727754, 40.928422009540796],\n",
    "          [-3.121975734708713, 40.928422009540796],\n",
    "          [-3.121975734708713, 40.999456684741496]]])\n",
    "\n",
    "\n",
    "AOI = ee.Geometry.Polygon(\n",
    "        [[[-2.97404756802236, 41.99472813320794],\n",
    "          [-2.97404756802236, 41.8107480842418],\n",
    "          [-2.66643038052236, 41.8107480842418],\n",
    "          [-2.66643038052236, 41.99472813320794]]])\n",
    "\n",
    "# cloudmask01 = getCloudMask(median)\n",
    "\n",
    "START_DATE = '2019-01-01'\n",
    "END_DATE = '2019-12-31'\n",
    "\n",
    "\n",
    "\n",
    "masks = { 'gsw': 1}#,'forestMask': {'buffer':10, 'startDate':'2000-01-01', 'endDate':'2022-12-31',}} # # , 'lmask': 30,\n",
    "\n",
    "\n",
    "s2 = Sentinel2(AOI, START_DATE, END_DATE, 50, masks)\n",
    "# s2.byMonth(2020)\n",
    "\n",
    "medianS2 = s2.getMedian()\n",
    "bandNames = medianS2.bandNames()\n",
    "print('Band names medianS2: ', bandNames.getInfo())\n",
    "med = s2.getCollection().first()\n",
    "\n",
    "# Subset layers and prepare them for display.\n",
    "clouds = med.select('clouds').selfMask()\n",
    "shadows = med.select('shadows').selfMask()\n",
    "\n",
    "\n",
    "dark_pixels = med.select('dark_pixels')  # .selfMask()\n",
    "probability = med.select('probability')\n",
    "cloudmask = med.select('cloudmask')\n",
    "cloudmaskMasked = cloudmask  # .selfMask()\n",
    "cloud_transform = med.select('cloud_transform')\n",
    "\n",
    "medianS2 = medianS2.select(['B4', 'B3', 'B2'])\n",
    "\n",
    "task = ee.batch.Export.image.toDrive(**{\n",
    "    'image': medianS2,\n",
    "    'description': 'medRGB1_GB1_2019',\n",
    "    'folder': 'earth_engine_demos',\n",
    "    'scale': 10,\n",
    "    'region': AOI.getInfo()['coordinates'],\n",
    "    'maxPixels': 1549491660\n",
    "})\n",
    "task.start()\n",
    "\n",
    "\n",
    "# Add Land Mask to the map\n",
    "my_map.add_ee_layer(medianS2, {\n",
    "                    'bands': ['B4', 'B3', 'B2'], 'min': 0, 'max': 2500, 'gamma': 1.1}, 'median')\n",
    "my_map.add_ee_layer(AOI, {}, 'AOI')\n",
    "\n",
    "\n",
    "\n",
    "# Add a layer control panel to the map.\n",
    "my_map.add_child(folium.LayerControl())\n",
    "plugins.Fullscreen().add_to(my_map)\n",
    "\n",
    "# Display the map.\n",
    "display(my_map)\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.11.1"
  },
  "orig_nbformat": 4,
  "vscode": {
   "interpreter": {
    "hash": "cab1d7ef7b90e69a2393a883ac82077044fd5f1d4df4dae9ffefa7c49ee44033"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}