Sentinel2_preprocessing.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# <div style=\"text-align: right;font-family:Times New Roman\">Приложение Б-1</div>\n",
    "# <div style=\"text-align: center;font-family:Times New Roman\">Предобработка снимков Sentinel-2</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Preprocessing of Sentinel-2 imagery\n",
    "# Предобработка снимков Sentinel-2\n",
    "---\n",
    "In this notebook Sentinel-2 imagery is being preprocessed. Preprocessing includes atmospheric correction using Sen2Cor algorythm, applying superresolution algorythm to 20- and 60-meters bands, cloud masking, reprojecting and mosaic creation.\n",
    "\n",
    "В этом блокноте выполняется предобработка снимков Sentinel-2. Предобработка включает в себя атмосферную коррекцию с применением алгоритма Sen2Cor, применение алгоритма суперразрешения к каналам с разрешением 20 и 60 м, маскировка облаков, перепроецирование и создание мозаики."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Importing dependencies\n",
    "### Импорт необходимых библиотек"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import subprocess\n",
    "import shutil\n",
    "import random\n",
    "import re\n",
    "from typing import List, Tuple\n",
    "from collections import defaultdict\n",
    "from pathlib import Path\n",
    "import numpy as np\n",
    "from glob import glob\n",
    "import geopandas as gpd\n",
    "from osgeo import gdal\n",
    "import rasterio as rio\n",
    "from rasterio.warp import calculate_default_transform, reproject, Resampling\n",
    "import rasterio.merge\n",
    "import rasterio.fill\n",
    "from rasterio.io import MemoryFile\n",
    "from rasterio.enums import Resampling\n",
    "from rasterio.windows import Window\n",
    "from rasterio.transform import Affine\n",
    "import fiona\n",
    "import shapely\n",
    "from shapely.geometry import box\n",
    "import zipfile"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Unzipping all archives with images\n",
    "### Распаковка всех архивов со снимками"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [],
   "source": [
    "#path to archives / путь к архивам\n",
    "zippath = 'E:\\\\Work\\\\Sentinel\\\\'\n",
    "archives = glob(zippath+'*')\n",
    "#where to unpack / конечный путь распаковки\n",
    "wpath = 'F:\\\\Work\\\\CorrProizv\\\\LVL1\\\\'\n",
    "for archive in archives:\n",
    "    with zipfile.ZipFile(archive, 'r') as zip_ref:\n",
    "        zip_ref.extractall(wpath)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Applying [Sen2Cor](https://step.esa.int/main/snap-supported-plugins/sen2cor/) atmospheric correction algorythm\n",
    "### Применение алгоритма атмосферной коррекции [Sen2Cor](https://step.esa.int/main/snap-supported-plugins/sen2cor/)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true,
    "jupyter": {
     "outputs_hidden": true
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VMN_20210915T184413.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VNQ_20210916T052032.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VPQ_20210915T081105.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20180630T072621_N9999_R049_T41VML_20210916T045925.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20180723T055641_N9999_R091_T44VNP_20210915T192826.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20180826T071621_N9999_R006_T41VMK_20210915T093843.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMN_20210915T070556.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMR_20210915T195148.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190715T054641_N9999_R048_T44VPN_20210915T042924.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T43VFJ_20210915T012154.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLN_20210915T182127.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLP_20210915T005950.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLQ_20210915T073007.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VNQ_20210915T032204.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T45VUK_20210915T024042.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40VFR_20210914T225005.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFS_20210915T003945.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFT_20210915T020246.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VLL_20210914T231151.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VMJ_20210916T043910.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VML_20210915T014030.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMM_20210915T001800.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMN_20210915T021927.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VDG_20210914T212730.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VFH_20210914T220628.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLN_20210914T222817.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLP_20210914T204329.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VMQ_20210914T233339.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190814T054641_N9999_R048_T44VPP_20210915T074756.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190814T054641_N9999_R048_T45VUJ_20210914T214401.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190818T070621_N9999_R106_T41VPL_20210914T174935.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190820T074611_N9999_R135_T41WMN_20210914T132020.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNJ_20210914T153753.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNL_20210914T125849.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VPK_20210914T143947.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WNM_20210914T141333.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WPM_20210914T191419.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T42VUR_20210914T123742.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWM_20210915T140446.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWP_20210914T115544.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXN_20210915T142700.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXP_20210915T175658.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXQ_20210915T153835.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXR_20210915T160235.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFH_20210914T102732.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFJ_20210914T111038.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T44VLN_20210914T113238.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200605T060641_N9999_R134_T44VMQ_20210914T100459.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200609T054641_N9999_R048_T44VPQ_20210914T082206.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200609T054641_N9999_R048_T45VVJ_20210914T084258.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200610T065631_N9999_R063_T42VUP_20210914T073838.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VEJ_20210915T164842.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VFJ_20210914T133734.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200628T061641_N9999_R034_T44VMQ_20210914T162133.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200703T070631_N9999_R106_T41VNK_20210915T053945.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200703T070631_N9999_R106_T42VUP_20210914T071651.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200708T061631_N9999_R034_T43VFJ_20210914T065417.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200709T072621_N9999_R049_T40VFQ_20210914T063312.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200709T072621_N9999_R049_T41VLJ_20210914T050403.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNF_20210914T034243.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNG_20210914T032628.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VUM_20210914T012244.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVM_20210915T131934.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVN_20210914T040112.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVP_20210914T061220.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVQ_20210914T024337.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVR_20210914T030501.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWN_20210914T044450.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWP_20210914T042245.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWQ_20210914T054847.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWR_20210914T052551.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VXQ_20210914T020454.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WVS_20210914T014405.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WWS_20210914T010059.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VNH_20210914T022220.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPG_20210913T235742.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPH_20210913T233407.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VWN_20210913T215201.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXM_20210914T004022.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXN_20210913T221031.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXP_20210913T212857.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXQ_20210913T191848.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCG_20210913T210750.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCH_20210913T225443.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCJ_20210913T204528.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDG_20210913T231520.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDJ_20210914T001753.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDK_20210913T223144.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200716T071621_N9999_R006_T41VNG_20210913T171218.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVL_20210915T125706.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVM_20210913T160650.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VXR_20210913T145941.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCK_20210913T152148.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCL_20210913T154448.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VDL_20210913T162838.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200719T072621_N9999_R049_T41VMJ_20210913T143930.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200802T070631_N9999_R106_T41VPL_20210913T141855.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200803T063631_N9999_R120_T43VCF_20210913T125627.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200804T060641_N9999_R134_T44VMP_20210913T123421.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200805T053651_N9999_R005_T45VVH_20210913T121305.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VUK_20210913T135813.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVL_20210913T133800.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVM_20210913T115207.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VWM_20210913T131744.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200807T061641_N9999_R034_T43VEG_20210913T103018.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200808T054641_N9999_R048_T44VPP_20210913T100855.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VDH_20210913T084024.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEG_20210913T113034.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEH_20210913T110922.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEJ_20210913T090157.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEK_20210915T062100.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VFK_20210913T092345.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2A_MSIL2A_20200829T065631_N9999_R063_T42VVP_20210913T081938.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180705T072729_N9999_R049_T41WMM_20210915T091557.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180714T061629_N9999_R034_T43VFK_20210915T171207.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180715T072619_N9999_R049_T41VML_20210915T100116.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNH_20210915T115154.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNL_20210915T104413.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41WNM_20210915T110616.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180724T061629_N9999_R034_T44VMP_20210915T190617.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFS_20210915T085334.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFT_20210915T064455.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T41VLL_20210915T083139.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T41WMN_20210916T054759.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20180828T070609_N9999_R106_T41VPJ_20210915T102232.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190702T062639_N9999_R077_T43VDG_20210915T145002.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190703T055649_N9999_R091_T44VNP_20210915T045336.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPF_20210915T121406.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPJ_20210915T040715.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPK_20210915T030003.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUN_20210915T051743.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUP_20210915T034512.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190715T063639_N9999_R120_T43VDJ_20210915T162517.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190725T063639_N9999_R120_T43VDL_20210915T151438.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190802T055649_N9999_R091_T44VNN_20210914T235559.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VCG_20210915T173546.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VEK_20210914T210540.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190815T074619_N9999_R135_T40WFS_20210914T202143.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T40VFQ_20210914T160028.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T40WFS_20210914T181115.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLJ_20210914T150753.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLL_20210914T183122.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMJ_20210914T185302.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMK_20210914T195925.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VNK_20210914T164557.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMM_20210914T193713.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMN_20210914T170532.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190820T065629_N9999_R063_T41VPH_20210914T172740.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190820T065629_N9999_R063_T42VUM_20210915T123502.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20190829T072619_N9999_R049_T40VFR_20210914T104948.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUP_20210915T112829.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUQ_20210915T060139.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200607T055639_N9999_R091_T45VUH_20210914T092151.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200607T073619_N9999_R092_T41VLK_20210914T094423.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200608T070629_N9999_R106_T41VMH_20210914T090442.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200610T060639_N9999_R134_T44VNQ_20210914T075955.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200617T055639_N9999_R091_T44VMN_20210914T121732.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPF_20210913T181456.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPG_20210913T194132.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VUL_20210913T202449.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVP_20210913T183538.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVQ_20210913T175303.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VWQ_20210913T200226.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VXR_20210913T185752.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VCF_20210915T134227.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VDF_20210913T165044.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VEF_20210913T173207.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200808T063629_N9999_R120_T42VWN_20210913T105049.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200808T063629_N9999_R120_T43VDK_20210913T094611.SAFE',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\S2B_MSIL2A_20200814T065629_N9999_R063_T42VVQ_20210913T075830.SAFE']"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#importing data / импорт данных\n",
    "wpath = 'F:\\\\Work\\\\CorrProizv\\\\Imagery\\\\'\n",
    "paths = glob(wpath+'*')\n",
    "paths"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#creating env for running Sen2Cor / создание среды для запуска Sen2Cor\n",
    "my_env = os.environ.copy()\n",
    "my_env['PATH'] = r'C:\\ProgramData\\Sen2Cor-02.09.00-win64\\bin;' + r'C:\\ProgramData\\Sen2Cor-02.09.00-win64;' + my_env['PATH']\n",
    "my_env['SEN2COR_HOME'] = r'C:\\ProgramData\\Sen2Cor-02.09.00-win64'\n",
    "my_env['SEN2COR_BIN'] = r'C:\\ProgramData\\Sen2Cor-02.09.00-win64\\lib\\python2.7\\site-packages\\sen2cor'\n",
    "my_env['LC_NUMERIC'] = 'C'\n",
    "my_env['GDAL_DATA'] = r'C:\\ProgramData\\Sen2Cor-02.09.00-win64\\share\\gdal'\n",
    "my_env['GDAL_DRIVER_PATH']= 'disable'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# running Sen2Cor for every image / применение Sen2Cor к каждому снимку\n",
    "for path in paths:\n",
    "    cmd = r'L2A_Process.bat \"'+ path + '\" '\n",
    "    try:\n",
    "        result = subprocess.check_output(cmd, env=my_env, shell=True)\n",
    "        shutil.rmtree(path)\n",
    "        print(path)\n",
    "    except subprocess.CalledProcessError as e:\n",
    "        print(e)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Applying superresolution algorythm\n",
    "### Применение алгоритма суперразрешения"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "#importing python library for superresolution algorythm / импорт библиотеки, реализующей алгоритм суперразрешения\n",
    "#https://github.com/simonreise/s2-superres\n",
    "\n",
    "import sys\n",
    "sys.path.insert(1, r'C:\\Users\\mosko\\Downloads\\s2-superres\\src\\s2_superres')\n",
    "import s2_superres"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "out_dir = r'F:\\Work\\CorrProizv\\superres\\\\'\n",
    "\n",
    "#applying superresolution to every image / применяем алгоритм суперразрешения к каждому снимку\n",
    "for path in paths:\n",
    "    pref = Path(path).stem\n",
    "    if (Path(out_dir + pref + '_superresolution.tif').is_file()):\n",
    "        print('Pass: ' + path)\n",
    "    else:\n",
    "        s2_superres.Superresolution(input_dir = path, output_dir = out_dir, copy_original_bands = True).start()\n",
    "        print('DONE: ' + path)\n",
    "        print('!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Cloud masking, errors fixing and reprojection\n",
    "### Маскировка облаков, исправление технических ошибок и перепроецирование"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true,
    "jupyter": {
     "outputs_hidden": true
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VMN_20210915T184413.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VNQ_20210916T052032.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VPQ_20210915T081105.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20180630T072621_N9999_R049_T41VML_20210916T045925.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20180723T055641_N9999_R091_T44VNP_20210915T192826.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20180826T071621_N9999_R006_T41VMK_20210915T093843.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMN_20210915T070556.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMR_20210915T195148.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190715T054641_N9999_R048_T44VPN_20210915T042924.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190721T060641_N9999_R134_T43VFJ_20210915T012154.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLN_20210915T182127.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLP_20210915T005950.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLQ_20210915T073007.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VNQ_20210915T032204.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190721T060641_N9999_R134_T45VUK_20210915T024042.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T40VFR_20210914T225005.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFS_20210915T003945.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFT_20210915T020246.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VLL_20210914T231151.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VMJ_20210916T043910.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VML_20210915T014030.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMM_20210915T001800.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMN_20210915T021927.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VDG_20210914T212730.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VFH_20210914T220628.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLN_20210914T222817.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLP_20210914T204329.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VMQ_20210914T233339.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190814T054641_N9999_R048_T44VPP_20210915T074756.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190814T054641_N9999_R048_T45VUJ_20210914T214401.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190818T070621_N9999_R106_T41VPL_20210914T174935.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190820T074611_N9999_R135_T41WMN_20210914T132020.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNJ_20210914T153753.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNL_20210914T125849.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VPK_20210914T143947.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WNM_20210914T141333.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WPM_20210914T191419.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190821T071621_N9999_R006_T42VUR_20210914T123742.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWM_20210915T140446.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWP_20210914T115544.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXN_20210915T142700.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXP_20210915T175658.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXQ_20210915T153835.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXR_20210915T160235.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFH_20210914T102732.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFJ_20210914T111038.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20190823T061631_N9999_R034_T44VLN_20210914T113238.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200605T060641_N9999_R134_T44VMQ_20210914T100459.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200609T054641_N9999_R048_T44VPQ_20210914T082206.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200609T054641_N9999_R048_T45VVJ_20210914T084258.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200610T065631_N9999_R063_T42VUP_20210914T073838.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VEJ_20210915T164842.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VFJ_20210914T133734.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200628T061641_N9999_R034_T44VMQ_20210914T162133.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200703T070631_N9999_R106_T41VNK_20210915T053945.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200703T070631_N9999_R106_T42VUP_20210914T071651.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200708T061631_N9999_R034_T43VFJ_20210914T065417.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200709T072621_N9999_R049_T40VFQ_20210914T063312.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200709T072621_N9999_R049_T41VLJ_20210914T050403.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNF_20210914T034243.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNG_20210914T032628.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VUM_20210914T012244.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVM_20210915T131934.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVN_20210914T040112.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVP_20210914T061220.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVQ_20210914T024337.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVR_20210914T030501.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWN_20210914T044450.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWP_20210914T042245.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWQ_20210914T054847.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWR_20210914T052551.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VXQ_20210914T020454.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WVS_20210914T014405.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WWS_20210914T010059.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VNH_20210914T022220.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPG_20210913T235742.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPH_20210913T233407.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VWN_20210913T215201.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXM_20210914T004022.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXN_20210913T221031.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXP_20210913T212857.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXQ_20210913T191848.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCG_20210913T210750.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCH_20210913T225443.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCJ_20210913T204528.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDG_20210913T231520.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDJ_20210914T001753.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDK_20210913T223144.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200716T071621_N9999_R006_T41VNG_20210913T171218.SAFE\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVL_20210915T125706.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVM_20210913T160650.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VXR_20210913T145941.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCK_20210913T152148.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCL_20210913T154448.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VDL_20210913T162838.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200719T072621_N9999_R049_T41VMJ_20210913T143930.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200802T070631_N9999_R106_T41VPL_20210913T141855.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200803T063631_N9999_R120_T43VCF_20210913T125627.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200804T060641_N9999_R134_T44VMP_20210913T123421.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200805T053651_N9999_R005_T45VVH_20210913T121305.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VUK_20210913T135813.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVL_20210913T133800.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVM_20210913T115207.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VWM_20210913T131744.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200807T061641_N9999_R034_T43VEG_20210913T103018.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200808T054641_N9999_R048_T44VPP_20210913T100855.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VDH_20210913T084024.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEG_20210913T113034.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEH_20210913T110922.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEJ_20210913T090157.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEK_20210915T062100.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VFK_20210913T092345.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2A_MSIL2A_20200829T065631_N9999_R063_T42VVP_20210913T081938.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180705T072729_N9999_R049_T41WMM_20210915T091557.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180714T061629_N9999_R034_T43VFK_20210915T171207.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180715T072619_N9999_R049_T41VML_20210915T100116.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNH_20210915T115154.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNL_20210915T104413.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180722T071619_N9999_R006_T41WNM_20210915T110616.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180724T061629_N9999_R034_T44VMP_20210915T190617.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFS_20210915T085334.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFT_20210915T064455.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180827T073719_N9999_R092_T41VLL_20210915T083139.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180827T073719_N9999_R092_T41WMN_20210916T054759.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20180828T070609_N9999_R106_T41VPJ_20210915T102232.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190702T062639_N9999_R077_T43VDG_20210915T145002.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190703T055649_N9999_R091_T44VNP_20210915T045336.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPF_20210915T121406.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPJ_20210915T040715.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPK_20210915T030003.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUN_20210915T051743.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUP_20210915T034512.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190715T063639_N9999_R120_T43VDJ_20210915T162517.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190725T063639_N9999_R120_T43VDL_20210915T151438.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190802T055649_N9999_R091_T44VNN_20210914T235559.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VCG_20210915T173546.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VEK_20210914T210540.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190815T074619_N9999_R135_T40WFS_20210914T202143.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T40VFQ_20210914T160028.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T40WFS_20210914T181115.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLJ_20210914T150753.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLL_20210914T183122.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMJ_20210914T185302.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMK_20210914T195925.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VNK_20210914T164557.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMM_20210914T193713.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMN_20210914T170532.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190820T065629_N9999_R063_T41VPH_20210914T172740.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190820T065629_N9999_R063_T42VUM_20210915T123502.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20190829T072619_N9999_R049_T40VFR_20210914T104948.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUP_20210915T112829.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUQ_20210915T060139.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200607T055639_N9999_R091_T45VUH_20210914T092151.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200607T073619_N9999_R092_T41VLK_20210914T094423.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200608T070629_N9999_R106_T41VMH_20210914T090442.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200610T060639_N9999_R134_T44VNQ_20210914T075955.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200617T055639_N9999_R091_T44VMN_20210914T121732.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPF_20210913T181456.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPG_20210913T194132.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VUL_20210913T202449.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVP_20210913T183538.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVQ_20210913T175303.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VWQ_20210913T200226.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VXR_20210913T185752.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VCF_20210915T134227.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VDF_20210913T165044.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VEF_20210913T173207.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200808T063629_N9999_R120_T42VWN_20210913T105049.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200808T063629_N9999_R120_T43VDK_20210913T094611.SAFE\n",
      "F:\\Work\\CorrProizv\\Imagery\\S2B_MSIL2A_20200814T065629_N9999_R063_T42VVQ_20210913T075830.SAFE\n"
     ]
    }
   ],
   "source": [
    "for path in paths:\n",
    "    #loading bands / загрузка всех каналов\n",
    "    if os.path.isfile('F:\\\\Work\\\\CorrProizv\\\\projected\\\\' + path.split('\\\\')[4][:-5] + '\\\\B1.tif'):\n",
    "        print('Skip', path)\n",
    "    else:\n",
    "        img = 'F:\\\\Work\\\\CorrProizv\\\\superres\\\\' + path.split('\\\\')[4][:-5] + '_superresolution.tif'\n",
    "        with rio.open(img) as bands:\n",
    "            img = bands.read()\n",
    "            meta = bands.profile\n",
    "            bounds = bands.bounds\n",
    "            descs = bands.descriptions\n",
    "            transform = bands.transform\n",
    "        #masking zeros / маскировка нулей\n",
    "        img = np.where(img == 0, 1, img)\n",
    "        #masking 65555 pixels / маскировка значений 65555\n",
    "        mask = np.where(img > 15000, 0, 1)\n",
    "        for i in range(img.shape[0]):\n",
    "            img[i] = rio.fill.fillnodata(img[i], mask[i], max_search_distance = 500)\n",
    "        #masking clouds by Sentinel L1 mask / маскировка облаков при помощи всроенной маски Sentinel L1\n",
    "        try:\n",
    "            shape = glob(path + '\\\\GRANULE\\\\**\\\\QI_DATA\\\\MSK_CLOUDS_B00.gml')[0]\n",
    "            with fiona.open(shape, \"r\") as shapefile:\n",
    "                cloudmask = [feature[\"geometry\"] for feature in shapefile]\n",
    "            with MemoryFile() as memfile:\n",
    "                with memfile.open(\n",
    "                    driver='GTiff',\n",
    "                    height=img.shape[1],\n",
    "                    width=img.shape[2],\n",
    "                    count=img.shape[0],\n",
    "                    dtype=img.dtype,\n",
    "                    compress = 'lzw',\n",
    "                    crs=meta['crs'],\n",
    "                    transform=transform\n",
    "                ) as temp:\n",
    "                    temp.write(img)\n",
    "                    img, transform = rio.mask.mask(temp, cloudmask, invert = True, filled = True)\n",
    "        except:\n",
    "            pass\n",
    "        #masking clouds by algorythm based on Scene Classification of L2 product and bands 1 and 9\n",
    "        #маскировка облаков при помощи алгоритма, основанного на классификации снимка уровня 2 и каналов 1 и 9\n",
    "        l1path = glob(path + r'\\GRANULE\\*\\IMG_DATA\\R20m\\*SCL_20m.jp2')\n",
    "        with rio.open(l1path[0]) as classification:\n",
    "            cls = classification.read()\n",
    "            clsmeta = classification.profile\n",
    "        mask = np.empty_like(img[0])\n",
    "        reproject(\n",
    "            cls, mask,\n",
    "            src_transform = clsmeta['transform'],\n",
    "            dst_transform = meta['transform'],\n",
    "            src_crs = clsmeta['crs'],\n",
    "            dst_crs = meta['crs'],\n",
    "            resampling = Resampling.nearest)\n",
    "        for i in range(len(descs)):\n",
    "            if 'B1 ' in descs[i]:\n",
    "                band1 = i\n",
    "        for i in range(len(descs)):\n",
    "            if 'B9 ' in descs[i]:\n",
    "                band9 = i\n",
    "        clmask = np.where(((img[band1] >= 3200) & (img[band9] >= 6000) & (mask >= 7)), 1, 0)\n",
    "        #delete 1px gaps in mask / удаление пробелов в 1 пиксел\n",
    "        clmask = np.where((rio.fill.fillnodata(np.where(clmask == 1, 0, 1), np.where(clmask == 1, 0, 1), max_search_distance = 2)) == 1, 0, 1)\n",
    "        #emerge mask / увеличение маски\n",
    "        clmask = rio.fill.fillnodata(clmask, clmask, max_search_distance = 22)\n",
    "        img = np.where(clmask > 0, 0, img)\n",
    "        #masking nodata / маскировка\n",
    "        img = np.where(mask == 0, 0, img)\n",
    "        #reprojecting / перепроецирование\n",
    "        if meta['crs'] != 'EPSG:32642':\n",
    "            transform, width, height = calculate_default_transform(\n",
    "                meta['crs'], 'EPSG:32642', meta['width'], meta['height'], *bounds)\n",
    "            proj = np.zeros((img.shape[0], width, height), img.dtype)\n",
    "            img, transform = reproject(\n",
    "                source=img,\n",
    "                destination=proj,\n",
    "                src_transform=meta['transform'],\n",
    "                src_crs=meta['crs'],\n",
    "                dst_transform=transform,\n",
    "                dst_crs='EPSG:32642',\n",
    "                resampling=Resampling.nearest)\n",
    "            img = proj\n",
    "        #save / сохранение\n",
    "        Path('F:\\\\Work\\\\CorrProizv\\\\projected\\\\' + path.split('\\\\')[4][:-5]).mkdir(parents=True, exist_ok=True)\n",
    "        for i in range(img.shape[0]):\n",
    "            pathres = 'F:\\\\Work\\\\CorrProizv\\\\projected\\\\' + path.split('\\\\')[4][:-5] + '\\\\' + descs[i].split(' ')[1] + '.tif'\n",
    "            with rio.open(\n",
    "                pathres,\n",
    "                'w',\n",
    "                driver='GTiff',\n",
    "                height=img.shape[1],\n",
    "                width=img.shape[2],\n",
    "                count=1,\n",
    "                dtype=img.dtype,\n",
    "                compress = 'deflate',\n",
    "                PREDICTOR = 2,\n",
    "                ZLEVEL=9,\n",
    "                crs='EPSG:32642',\n",
    "                transform=transform,\n",
    "                nodata = 0\n",
    "            ) as outfile:\n",
    "                outfile.write(img[i],1)\n",
    "        print(path)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Merging all Sentinel-2 images into one mosaic for whole area\n",
    "### Объединение всех снимков Sentinel-2 в одну мозаику на всю территорию"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true,
    "jupyter": {
     "outputs_hidden": true
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VMN_20210915T184413',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VNQ_20210916T052032',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VPQ_20210915T081105',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20180630T072621_N9999_R049_T41VML_20210916T045925',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20180723T055641_N9999_R091_T44VNP_20210915T192826',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20180826T071621_N9999_R006_T41VMK_20210915T093843',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMN_20210915T070556',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMR_20210915T195148',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190715T054641_N9999_R048_T44VPN_20210915T042924',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T43VFJ_20210915T012154',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLN_20210915T182127',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLP_20210915T005950',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLQ_20210915T073007',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VNQ_20210915T032204',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T45VUK_20210915T024042',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40VFR_20210914T225005',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFS_20210915T003945',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFT_20210915T020246',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VLL_20210914T231151',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VMJ_20210916T043910',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VML_20210915T014030',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMM_20210915T001800',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMN_20210915T021927',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VDG_20210914T212730',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VFH_20210914T220628',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLN_20210914T222817',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLP_20210914T204329',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VMQ_20210914T233339',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190814T054641_N9999_R048_T44VPP_20210915T074756',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190814T054641_N9999_R048_T45VUJ_20210914T214401',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190818T070621_N9999_R106_T41VPL_20210914T174935',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190820T074611_N9999_R135_T41WMN_20210914T132020',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNJ_20210914T153753',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNL_20210914T125849',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VPK_20210914T143947',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WNM_20210914T141333',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WPM_20210914T191419',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T42VUR_20210914T123742',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWM_20210915T140446',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWP_20210914T115544',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXN_20210915T142700',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXP_20210915T175658',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXQ_20210915T153835',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXR_20210915T160235',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFH_20210914T102732',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFJ_20210914T111038',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T44VLN_20210914T113238',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200605T060641_N9999_R134_T44VMQ_20210914T100459',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200609T054641_N9999_R048_T44VPQ_20210914T082206',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200609T054641_N9999_R048_T45VVJ_20210914T084258',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200610T065631_N9999_R063_T42VUP_20210914T073838',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VEJ_20210915T164842',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VFJ_20210914T133734',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200628T061641_N9999_R034_T44VMQ_20210914T162133',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200703T070631_N9999_R106_T41VNK_20210915T053945',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200703T070631_N9999_R106_T42VUP_20210914T071651',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200708T061631_N9999_R034_T43VFJ_20210914T065417',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200709T072621_N9999_R049_T40VFQ_20210914T063312',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200709T072621_N9999_R049_T41VLJ_20210914T050403',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNF_20210914T034243',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNG_20210914T032628',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VUM_20210914T012244',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVM_20210915T131934',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVN_20210914T040112',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVP_20210914T061220',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVQ_20210914T024337',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVR_20210914T030501',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWN_20210914T044450',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWP_20210914T042245',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWQ_20210914T054847',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWR_20210914T052551',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VXQ_20210914T020454',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WVS_20210914T014405',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WWS_20210914T010059',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VNH_20210914T022220',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPG_20210913T235742',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPH_20210913T233407',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VWN_20210913T215201',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXM_20210914T004022',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXN_20210913T221031',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXP_20210913T212857',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXQ_20210913T191848',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCG_20210913T210750',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCH_20210913T225443',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCJ_20210913T204528',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDG_20210913T231520',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDJ_20210914T001753',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDK_20210913T223144',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200716T071621_N9999_R006_T41VNG_20210913T171218',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVL_20210915T125706',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVM_20210913T160650',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VXR_20210913T145941',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCK_20210913T152148',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCL_20210913T154448',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VDL_20210913T162838',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200719T072621_N9999_R049_T41VMJ_20210913T143930',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200802T070631_N9999_R106_T41VPL_20210913T141855',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200803T063631_N9999_R120_T43VCF_20210913T125627',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200804T060641_N9999_R134_T44VMP_20210913T123421',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200805T053651_N9999_R005_T45VVH_20210913T121305',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VUK_20210913T135813',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVL_20210913T133800',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVM_20210913T115207',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VWM_20210913T131744',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200807T061641_N9999_R034_T43VEG_20210913T103018',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200808T054641_N9999_R048_T44VPP_20210913T100855',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VDH_20210913T084024',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEG_20210913T113034',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEH_20210913T110922',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEJ_20210913T090157',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEK_20210915T062100',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VFK_20210913T092345',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200829T065631_N9999_R063_T42VVP_20210913T081938',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180705T072729_N9999_R049_T41WMM_20210915T091557',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180714T061629_N9999_R034_T43VFK_20210915T171207',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180715T072619_N9999_R049_T41VML_20210915T100116',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNH_20210915T115154',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNL_20210915T104413',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41WNM_20210915T110616',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180724T061629_N9999_R034_T44VMP_20210915T190617',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFS_20210915T085334',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFT_20210915T064455',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T41VLL_20210915T083139',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T41WMN_20210916T054759',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180828T070609_N9999_R106_T41VPJ_20210915T102232',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190702T062639_N9999_R077_T43VDG_20210915T145002',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190703T055649_N9999_R091_T44VNP_20210915T045336',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPF_20210915T121406',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPJ_20210915T040715',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPK_20210915T030003',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUN_20210915T051743',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUP_20210915T034512',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190715T063639_N9999_R120_T43VDJ_20210915T162517',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190725T063639_N9999_R120_T43VDL_20210915T151438',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190802T055649_N9999_R091_T44VNN_20210914T235559',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VCG_20210915T173546',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VEK_20210914T210540',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190815T074619_N9999_R135_T40WFS_20210914T202143',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T40VFQ_20210914T160028',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T40WFS_20210914T181115',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLJ_20210914T150753',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLL_20210914T183122',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMJ_20210914T185302',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMK_20210914T195925',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VNK_20210914T164557',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMM_20210914T193713',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMN_20210914T170532',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190820T065629_N9999_R063_T41VPH_20210914T172740',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190820T065629_N9999_R063_T42VUM_20210915T123502',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190829T072619_N9999_R049_T40VFR_20210914T104948',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUP_20210915T112829',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUQ_20210915T060139',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200607T055639_N9999_R091_T45VUH_20210914T092151',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200607T073619_N9999_R092_T41VLK_20210914T094423',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200608T070629_N9999_R106_T41VMH_20210914T090442',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200610T060639_N9999_R134_T44VNQ_20210914T075955',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200617T055639_N9999_R091_T44VMN_20210914T121732',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPF_20210913T181456',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPG_20210913T194132',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VUL_20210913T202449',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVP_20210913T183538',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVQ_20210913T175303',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VWQ_20210913T200226',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VXR_20210913T185752',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VCF_20210915T134227',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VDF_20210913T165044',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VEF_20210913T173207',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200808T063629_N9999_R120_T42VWN_20210913T105049',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200808T063629_N9999_R120_T43VDK_20210913T094611',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200814T065629_N9999_R063_T42VVQ_20210913T075830']"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#getting all images / список всех снимков\n",
    "spaths = glob('F:\\\\Work\\\\CorrProizv\\\\Projected\\\\'+'*')\n",
    "spaths"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "###### Sorting images by cloud cover - less cloudy on top\n",
    "###### Сортировка снимков по облачности - наименее облачные сверху"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "#getting number of nodata pixels for each image / количество значений \"нет данных\" для каждого снимка\n",
    "zeros = []\n",
    "for path in spaths:\n",
    "    path = path + '\\\\B1.tif'\n",
    "    with rio.open(path) as bnd:\n",
    "        img = bnd.read(1)\n",
    "    try:\n",
    "        zeros.append(np.count_nonzero(img==0))\n",
    "    except:\n",
    "        zeros.append(0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true,
    "jupyter": {
     "outputs_hidden": true
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLQ_20210915T073007',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T45VUK_20210915T024042',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFT_20210915T020246',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNG_20210914T032628',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200608T070629_N9999_R106_T41VMH_20210914T090442',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200617T055639_N9999_R091_T44VMN_20210914T121732',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190820T074611_N9999_R135_T41WMN_20210914T132020',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T40WFS_20210914T181115',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VDG_20210914T212730',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VPQ_20210915T081105',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLP_20210915T005950',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T43VFJ_20210915T012154',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40WFS_20210915T003945',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200609T054641_N9999_R048_T44VPQ_20210914T082206',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190815T074619_N9999_R135_T40WFS_20210914T202143',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200607T055639_N9999_R091_T45VUH_20210914T092151',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VXQ_20210914T020454',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VNK_20210914T164557',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200805T053651_N9999_R005_T45VVH_20210913T121305',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200609T054641_N9999_R048_T45VVJ_20210914T084258',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190814T054641_N9999_R048_T45VUJ_20210914T214401',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200808T054641_N9999_R048_T44VPP_20210913T100855',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200628T061641_N9999_R034_T44VMQ_20210914T162133',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T40VFR_20210914T225005',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190829T072619_N9999_R049_T40VFR_20210914T104948',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T41VNF_20210914T034243',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T44VLN_20210914T113238',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200709T072621_N9999_R049_T40VFQ_20210914T063312',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VNQ_20210915T032204',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T40VFQ_20210914T160028',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMR_20210915T195148',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200610T060639_N9999_R134_T44VNQ_20210914T075955',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180724T061629_N9999_R034_T44VMP_20210915T190617',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDG_20210913T231520',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190701T060641_N9999_R134_T44VMN_20210915T070556',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200605T060641_N9999_R134_T44VMQ_20210914T100459',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190721T060641_N9999_R134_T44VLN_20210915T182127',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFT_20210915T064455',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190703T055649_N9999_R091_T44VNP_20210915T045336',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190814T054641_N9999_R048_T44VPP_20210915T074756',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200703T070631_N9999_R106_T41VNK_20210915T053945',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T40WFS_20210915T085334',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20180723T055641_N9999_R091_T44VNP_20210915T192826',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VEF_20210913T173207',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VNQ_20210916T052032',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VMQ_20210914T233339',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190715T054641_N9999_R048_T44VPN_20210915T042924',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200804T060641_N9999_R134_T44VMP_20210913T123421',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLP_20210914T204329',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190802T055649_N9999_R091_T44VNN_20210914T235559',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20170731T060631_N9999_R134_T44VMN_20210915T184413',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T44VLN_20210914T222817',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VFJ_20210914T133734',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEK_20210915T062100',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200716T071621_N9999_R006_T41VNG_20210913T171218',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200709T072621_N9999_R049_T41VLJ_20210914T050403',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLL_20210914T183122',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200808T063629_N9999_R120_T42VWN_20210913T105049',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VWN_20210913T215201',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUQ_20210915T060139',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VLL_20210914T231151',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWN_20210914T044450',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFJ_20210914T111038',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMN_20210914T170532',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPF_20210915T121406',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200719T072621_N9999_R049_T41VMJ_20210913T143930',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VPK_20210914T143947',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMN_20210915T021927',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VEK_20210914T210540',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190823T061631_N9999_R034_T43VFH_20210914T102732',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200708T061631_N9999_R034_T43VFJ_20210914T065417',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VDL_20210913T162838',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200807T061641_N9999_R034_T43VEG_20210913T103018',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200607T073619_N9999_R092_T41VLK_20210914T094423',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCL_20210913T154448',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPH_20210913T233407',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPK_20210915T030003',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200803T063631_N9999_R120_T43VCF_20210913T125627',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDK_20210913T223144',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T43VCK_20210913T152148',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200802T070631_N9999_R106_T41VPL_20210913T141855',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180715T072619_N9999_R049_T41VML_20210915T100116',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VXR_20210913T185752',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200808T063629_N9999_R120_T43VDK_20210913T094611',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41WNM_20210915T110616',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T41WMN_20210916T054759',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VML_20210915T014030',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41WMM_20210915T001800',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180705T072729_N9999_R049_T41WMM_20210915T091557',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VFK_20210913T092345',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41WMM_20210914T193713',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEG_20210913T113034',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WNM_20210914T141333',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T41VPJ_20210915T040715',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VLJ_20210914T150753',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMK_20210914T195925',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180827T073719_N9999_R092_T41VLL_20210915T083139',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180828T070609_N9999_R106_T41VPJ_20210915T102232',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VDJ_20210914T001753',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20180826T071621_N9999_R006_T41VMK_20210915T093843',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VNH_20210914T022220',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190819T072619_N9999_R049_T41VMJ_20210914T185302',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPG_20210913T194132',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190725T072621_N9999_R049_T41VMJ_20210916T043910',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41WPM_20210914T191419',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNL_20210914T125849',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20180630T072621_N9999_R049_T41VML_20210916T045925',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190725T063639_N9999_R120_T43VDL_20210915T151438',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNL_20210915T104413',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190818T070621_N9999_R106_T41VPL_20210914T174935',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200829T065631_N9999_R063_T42VVP_20210913T081938',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEJ_20210913T090157',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180714T061629_N9999_R034_T43VFK_20210915T171207',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCJ_20210913T204528',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200611T062631_N9999_R077_T43VEJ_20210915T164842',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCG_20210913T210750',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T41VNJ_20210914T153753',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190715T063639_N9999_R120_T43VDJ_20210915T162517',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200713T070631_N9999_R106_T41VPG_20210913T235742',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190820T065629_N9999_R063_T41VPH_20210914T172740',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190814T063629_N9999_R120_T43VCG_20210915T173546',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T43VCH_20210913T225443',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VDH_20210913T084024',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190813T061631_N9999_R034_T43VFH_20210914T220628',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20180722T071619_N9999_R006_T41VNH_20210915T115154',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200810T062631_N9999_R077_T43VEH_20210913T110922',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190702T062639_N9999_R077_T43VDG_20210915T145002',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VCF_20210915T134227',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T41VPF_20210913T181456',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200716T062629_N9999_R077_T43VDF_20210913T165044',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVM_20210913T115207',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VUM_20210914T012244',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WVS_20210914T014405',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200814T065629_N9999_R063_T42VVQ_20210913T075830',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVQ_20210913T175303',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VXR_20210913T145941',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VUK_20210913T135813',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VWM_20210913T131744',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200806T064631_N9999_R020_T42VVL_20210913T133800',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200703T070631_N9999_R106_T42VUP_20210914T071651',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXQ_20210913T191848',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200610T065631_N9999_R063_T42VUP_20210914T073838',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXN_20210915T142700',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VVP_20210913T183538',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190820T065629_N9999_R063_T42VUM_20210915T123502',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200605T065629_N9999_R063_T42VUP_20210915T112829',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWP_20210914T042245',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWP_20210914T115544',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVM_20210913T160650',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VWQ_20210913T200226',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUP_20210915T034512',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VWM_20210915T140446',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXR_20210915T160235',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXM_20210914T004022',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXP_20210915T175658',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190822T064631_N9999_R020_T42VXQ_20210915T153835',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXP_20210913T212857',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200714T063631_N9999_R120_T42VXN_20210913T221031',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42WWS_20210914T010059',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWQ_20210914T054847',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVN_20210914T040112',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVR_20210914T030501',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVP_20210914T061220',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVM_20210915T131934',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VVQ_20210914T024337',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200710T065631_N9999_R063_T42VWR_20210914T052551',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20190714T070629_N9999_R106_T42VUN_20210915T051743',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20200717T064631_N9999_R020_T42VVL_20210915T125706',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2A_MSIL2A_20190821T071621_N9999_R006_T42VUR_20210914T123742',\n",
       " 'F:\\\\Work\\\\CorrProizv\\\\Projected\\\\S2B_MSIL2A_20200715T065629_N9999_R063_T42VUL_20210913T202449']"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#sorting images by number of nodata pixels / сортировка снимков по количеству пикселов \"нет даннных\"\n",
    "zerosdict = dict(zip(spaths, zeros))\n",
    "sortedzeros = dict(sorted(zerosdict.items(), key=lambda item: item[1], reverse = True))\n",
    "order = list(sortedzeros.keys())\n",
    "order"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "merge done\n"
     ]
    }
   ],
   "source": [
    "#merging images into one mosaic / объединение снимков\n",
    "bands = ['B1', 'B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'B9', 'B11', 'B12']\n",
    "for band in bands:\n",
    "    #reading n-th band from every image / чтение n-ного канала каждого из снимков\n",
    "    ordfiles = []\n",
    "    for img in order:\n",
    "        img = img + '\\\\' + band + '.tif'\n",
    "        ordfiles.append(img)\n",
    "    files = []\n",
    "    for path in ordfiles:\n",
    "        pathfile = rio.open(path)\n",
    "        files.append(pathfile)\n",
    "    #merging / объединение\n",
    "    final, final_trans = rio.merge.merge(files,method = 'last')\n",
    "    print('merge done')\n",
    "    #filling the gaps (if exist) / заполнение пустот (если они есть)\n",
    "    final = rio.fill.fillnodata(final, final, max_search_distance=250)\n",
    "    print('fill done')\n",
    "    #clipping the mosaic to the extent of KhMAO / обрезка мозаики по границе ХМАО\n",
    "    files = None\n",
    "    memfile = None\n",
    "    temp = None\n",
    "    with fiona.open(\"E:\\\\Work\\\\CorrProizv\\\\clipper.gpkg\", \"r\") as shapefile:\n",
    "        shape = [feature[\"geometry\"] for feature in shapefile]\n",
    "    with MemoryFile() as memfile:\n",
    "        with memfile.open(\n",
    "            driver='GTiff',\n",
    "            height=final.shape[1],\n",
    "            width=final.shape[2],\n",
    "            count=final.shape[0],\n",
    "            dtype=final.dtype,\n",
    "            compress = 'lzw',\n",
    "            crs='EPSG:32642',\n",
    "            transform=final_trans,\n",
    "            BIGTIFF='YES',\n",
    "            nodata = 0\n",
    "        ) as temp:\n",
    "            temp.write(final)\n",
    "            final, final_trans = rio.mask.mask(temp, shape,crop=True, filled=True)\n",
    "    print('clipping done')\n",
    "    #saving / сохранение\n",
    "    with rio.open(\n",
    "        'F:\\\\Work\\\\CorrProizv\\\\final_noadj\\\\'+ band + '.tif',\n",
    "        'w',\n",
    "        driver='GTiff',\n",
    "        height=final.shape[1],\n",
    "        width=final.shape[2],\n",
    "        count=1,\n",
    "        dtype=final.dtype,\n",
    "        compress = 'lzw',\n",
    "        crs='EPSG:32642',\n",
    "        transform=final_trans,\n",
    "        BIGTIFF='YES',\n",
    "        nodata = 0\n",
    "    ) as outfile:\n",
    "        outfile.write(final)\n",
    "    print(band)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Sentinel-2 processing tools from this notebook are now avaliable in a Python library named [Remote Sensing Processor](https://github.com/simonreise/remote-sensing-processor)\n",
    "#### Инструменты предобработки Sentinel-2, использованные в этом блокноте доступны в библиотеке Python [Remote Sensing Processor](https://github.com/simonreise/remote-sensing-processor)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}