geom_corr2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from osgeo import gdal, osr, gdalconst\n",
    "from osgeo import ogr\n",
    "import geopandas as gpd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import re\n",
    "import copy\n",
    "import os\n",
    "import cv2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2011181.0325.061.2017324101641.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0200.061.2017334143411.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2012182.0335.061.2017334143816.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2013181.0305.061.2017299063350.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2014181.0235.061.2017312140235.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0205.061.2017321044912.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2015181.0345.061.2017321045213.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2016182.0220.061.2017326170448.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2017181.0325.061.2017315065028.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2018181.0255.061.2018181133528.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n",
      "The georeference of ./hdfs/MOD021KM.A2019181.0225.061.2019181131027.hdf is complete.\n"
     ]
    }
   ],
   "source": [
    "import os\n",
    "from osgeo import gdal\n",
    "\n",
    "\n",
    "def swath_georeference(in_file, geo_file, target_dataset_name, out_file, out_geo_range, x_res_out, y_res_out):\n",
    "    dataset = gdal.Open(in_file)\n",
    "    subdatasets = dataset.GetSubDatasets()\n",
    "    dataset_pos = 0\n",
    "    target_pos = 0#-1\n",
    "    for subdataset in subdatasets:\n",
    "        subdataset_name = subdataset[0]\n",
    "        if subdataset_name.endswith(target_dataset_name):\n",
    "            target_pos = dataset_pos\n",
    "        dataset_pos += 1\n",
    "        \n",
    "        if target_pos != -1:\n",
    "            target_dataset = dataset.GetSubDatasets()[target_pos][0]\n",
    "            vrt_file = os.path.splitext(in_file)[0] + '.vrt'\n",
    "            vrt_out = gdal.Translate(vrt_file, target_dataset,\n",
    "                                     format='vrt', unscale='true')\n",
    "            georeferenced_data = gdal.Warp(out_file+\"_\"+str(target_pos+1)+\".tif\", vrt_out, multithread=True, outputBounds=out_geo_range,\n",
    "                                           format='GTiff', geoloc=True, dstSRS=\"EPSG:4326\",\n",
    "                                           xRes=x_res_out, yRes=y_res_out, dstNodata=0.0, outputType=gdal.GDT_Float32)\n",
    "            target_pos += 1\n",
    "            print('The georeference of {} is complete.'.format(in_file))\n",
    "            os.remove(vrt_file)\n",
    "        else:\n",
    "            print('{} has no the target dataset.'.format(in_file))\n",
    "\n",
    "\n",
    "def main():\n",
    "    file_prefix = '.hdf'\n",
    "    input_directory = './hdfs/'\n",
    "    output_directory = './hdf-output/'\n",
    "    if not os.path.exists(output_directory):\n",
    "        os.mkdir(output_directory)\n",
    "    file_list = os.listdir(input_directory)\n",
    "    for file_i in file_list:\n",
    "        if file_i.endswith(file_prefix):\n",
    "            file_name = input_directory + file_i\n",
    "            out_file = output_directory + os.path.splitext(file_i)[0]# + '_geo.tiff'\n",
    "            dataset_name = 'Image_Optical_Depth_Land_And_Ocean'\n",
    "            geo_range = None  # [90.0, 55.0, 120.0, 25.0]\n",
    "            x_res = 0.01\n",
    "            y_res = 0.01\n",
    "            swath_georeference(file_name, file_name, dataset_name, out_file, geo_range, x_res, y_res)\n",
    "\n",
    "\n",
    "if __name__ == '__main__':\n",
    "    main()\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "MOD021KM.A2011181.0325.061.2017324101641_5.tif\n",
      "1352.6150066102987 -3.090690017977832e-06 290.61072989740563\n",
      "1353 0 290.61047957125885\n",
      "2230.32484992828 -1.1335507037525168e-06 373.30844880284707\n",
      "2230 0 373.3084940191718\n",
      "MOD021KM.A2011181.0325.061.2017324101641_8.tif\n",
      "1519.7290279418173 -3.1262600652439585e-06 328.3977316575475\n",
      "1520 0 328.3979096452941\n",
      "2760.9612789306243 -1.7878082111906037e-06 301.85815196723627\n",
      "2761 0 301.85805650708585\n",
      "2643.8940579156847 -1.633568624626066e-06 304.9740656387131\n",
      "2644 0 304.9739884043758\n",
      "2385.6331045348275 -1.330045569161508e-06 198.88942906284754\n",
      "2386 0 198.88934570015746\n",
      "1943.8356449396815 -8.830357103446287e-07 114.77932946516616\n",
      "1944 0 114.7792787784193\n",
      "MOD021KM.A2012182.0200.061.2017334143411_5.tif\n",
      "3691.882807551365 -3.1559880276182373e-06 127.08707693834334\n",
      "3692 0 127.0870440569455\n",
      "2229.8390512830792 -1.150619908984426e-06 139.94600886285636\n",
      "2230 0 139.94620252824242\n",
      "MOD021KM.A2012182.0200.061.2017334143411_8.tif\n",
      "3758.685771563983 -3.2712335071216977e-06 179.00998344655483\n",
      "3759 0 179.00974877091858\n",
      "2898.291676378588 -1.9450183493967655e-06 141.43935539994462\n",
      "2898 0 141.4394847446934\n",
      "2691.2201899284987 -1.676535395626647e-06 111.42371820734252\n",
      "2691 0 111.42341518845318\n",
      "2419.054652211602 -1.3543978588427082e-06 70.81818975844243\n",
      "2419 0 70.81865622821215\n",
      "1980.1115797619848 -9.003841542058954e-07 43.94728396133161\n",
      "1980 0 43.947181850280586\n",
      "MOD021KM.A2012182.0335.061.2017334143816_5.tif\n",
      "3637.708245938404 -3.148083193559912e-06 249.88991146835514\n",
      "3638 0 249.89002686180487\n",
      "2214.492439049174 -1.147737889530109e-06 338.64057784056456\n",
      "2214 0 338.6408286000668\n",
      "MOD021KM.A2012182.0335.061.2017334143816_8.tif\n",
      "3705.728954541795 -3.2630705238782872e-06 271.0904272571181\n",
      "3706 0 271.0907406357393\n",
      "2867.2600968053653 -1.9401761638699877e-06 255.9230685390307\n",
      "2867 0 255.92335974215666\n",
      "2652.9403909705425 -1.6723363423210483e-06 290.5563420657706\n",
      "2653 0 290.5562697564381\n",
      "2384.4810792259036 -1.351003305757396e-06 206.3032123571453\n",
      "2384 0 206.3032008456072\n",
      "1944.1722691256334 -8.981278526082626e-07 130.63168574813596\n",
      "1944 0 130.63172070345846\n",
      "MOD021KM.A2013181.0305.061.2017299063350_5.tif\n",
      "3677.5973384114013 -3.202957518052357e-06 232.04694053290063\n",
      "3678 0 232.04693381550848\n",
      "2260.546387088984 -1.16656351297853e-06 294.3365471401645\n",
      "2261 0 294.33641206513335\n",
      "MOD021KM.A2013181.0305.061.2017299063350_8.tif\n",
      "3793.3946499366198 -3.4077718728587367e-06 267.4618712196547\n",
      "3793 0 267.4619309219455\n",
      "2937.4131790975316 -2.040639343923505e-06 243.39331482373882\n",
      "2937 0 243.39336740311686\n",
      "2683.758397753801 -1.7031455627911338e-06 262.1917796477237\n",
      "2684 0 262.1917267145001\n",
      "2404.6880461685464 -1.3672812719117272e-06 188.88391626858342\n",
      "2405 0 188.88391435053595\n",
      "1971.1441082906238 -9.181712681644569e-07 115.9633042442873\n",
      "1971 0 115.96319317477632\n",
      "MOD021KM.A2014181.0235.061.2017312140235_5.tif\n",
      "1338.581311811958 -3.251257468856124e-06 230.85711218002834\n",
      "1339 0 230.8572697382467\n",
      "2266.5765391070536 -1.1821679712204956e-06 277.6600603822605\n",
      "2267 0 277.66004181351764\n",
      "MOD021KM.A2014181.0235.061.2017312140235_8.tif\n",
      "1342.6834931190263 -3.221042650960449e-06 272.792764593021\n",
      "1343 0 272.79272566325176\n",
      "1309.4332237800872 -2.061597989589939e-06 241.51373041079847\n",
      "1309 0 241.51371601625243\n",
      "2711.9720196970993 -1.7363527473854699e-06 231.3458738534429\n",
      "2712 0 231.3459117678939\n",
      "2424.339693156819 -1.3874208979703215e-06 153.02402301534983\n",
      "2424 0 153.02417957805463\n",
      "1991.5362187373107 -9.362627231995815e-07 91.13193763572204\n",
      "1992 0 91.13216881756102\n",
      "MOD021KM.A2015181.0205.061.2017321044912_5.tif\n",
      "1292.5296560731229 -3.337586214945851e-06 253.8240659132979\n",
      "1293 0 253.82382159355032\n",
      "2314.264342410562 -1.2125879336641829e-06 291.11508377261345\n",
      "2314 0 291.1150005565553\n",
      "MOD021KM.A2015181.0205.061.2017321044912_8.tif\n",
      "1347.772567934467 -3.3766867124243956e-06 303.77220769923946\n",
      "1348 0 303.77223256193\n",
      "1295.2518872985336 -2.2581994004288614e-06 265.8083950319641\n",
      "1295 0 265.8084949552238\n",
      "1053.5394917502776 -1.7622962251240107e-06 226.02183497204823\n",
      "1054 0 226.02170913492756\n",
      "760.2896513322395 -1.4024686446486196e-06 135.38817750521775\n",
      "760 0 135.38785791805915\n",
      "573.7322059538108 -9.603503087329123e-07 81.95336280046014\n",
      "574 0 81.95342741988917\n",
      "MOD021KM.A2015181.0345.061.2017321045213_5.tif\n",
      "1338.9207406632206 -3.343362606982042e-06 271.3077603859171\n",
      "1339 0 271.30778283028803\n",
      "2304.721449978678 -1.2146880288745687e-06 359.07337498626094\n",
      "2305 0 359.07347904276645\n",
      "MOD021KM.A2015181.0345.061.2017321045213_8.tif\n",
      "1369.0712530219025 -3.3825281950105677e-06 307.2139436851064\n",
      "1369 0 307.2142434951515\n",
      "1331.7026912836582 -2.262095977715213e-06 283.6006870547553\n",
      "1332 0 283.60068229567315\n",
      "2741.0724987987946 -1.7653444538070437e-06 309.0540248461357\n",
      "2741 0 309.0540852957476\n",
      "895.9557383482554 -1.4048929370478692e-06 208.07404092553966\n",
      "896 0 208.07377991871869\n",
      "690.1898193086505 -9.620132862450659e-07 120.06398849989628\n",
      "690 0 120.06372521650661\n",
      "MOD021KM.A2016182.0220.061.2017326170448_5.tif\n",
      "1236.6266681969569 -3.428039394932924e-06 150.29047478771196\n",
      "1237 0 150.29070791471315\n",
      "2320.6605608466834 -1.253583102612631e-06 162.10155076848557\n",
      "2321 0 162.10144021485777\n",
      "MOD021KM.A2016182.0220.061.2017326170448_8.tif\n",
      "1290.459993987626 -3.513793442072028e-06 201.20268358172692\n",
      "1290 0 201.2028795262837\n",
      "1237.2391658848458 -2.417604242023576e-06 162.7842960542045\n",
      "1237 0 162.7842247073423\n",
      "836.2088225055048 -1.929116176605095e-06 136.33599149487824\n",
      "836 0 136.33585384786065\n",
      "788.4358514215585 -1.4836097057933832e-06 87.30803532583175\n",
      "788 0 87.30790654612237\n",
      "640.4932081787491 -1.0478079373380835e-06 53.42241934205344\n",
      "640 0 53.42224903206316\n",
      "MOD021KM.A2017181.0325.061.2017315065028_5.tif\n",
      "1300.3833768287002 -3.463300939146519e-06 294.3867160162516\n",
      "1300 0 294.3867815047614\n",
      "2336.2449423163675 -1.2525152070495153e-06 374.1859160424517\n",
      "2336 0 374.1861317497879\n",
      "MOD021KM.A2017181.0325.061.2017315065028_8.tif\n",
      "1476.5401712418684 -3.711281741743264e-06 333.5886967198285\n",
      "1477 0 333.5887622511127\n",
      "1318.4609321082523 -2.5896540652773614e-06 304.6247177793456\n",
      "1318 0 304.6250513995553\n",
      "997.9607271870586 -2.08971328072388e-06 303.26847746170756\n",
      "998 0 303.2684274736339\n",
      "735.0626350087597 -1.525471624291823e-06 173.25319166776816\n",
      "735 0 173.2531499704309\n",
      "563.6735011499628 -1.0742870508952711e-06 99.3801636064436\n",
      "564 0 99.38008898093827\n",
      "MOD021KM.A2018181.0255.061.2018181133528_5.tif\n",
      "1266.8613505277392 -3.5523168436500767e-06 255.38605495688856\n",
      "1267 0 255.3859371985095\n",
      "2370.110247525251 -1.2805405973845016e-06 317.00141733470645\n",
      "2370 0 317.00136319601114\n",
      "MOD021KM.A2018181.0255.061.2018181133528_8.tif\n",
      "1314.5498517101876 -3.844276747603632e-06 294.79195526294575\n",
      "1315 0 294.791834110379\n",
      "1264.6648923037314 -2.790017969936016e-06 265.88319638682555\n",
      "1265 0 265.8833613639505\n",
      "3040.8146222072824 -2.178247480720522e-06 269.26987456449183\n",
      "3041 0 269.26998074466854\n",
      "2557.2436722598814 -1.5405349433518345e-06 170.54795803192908\n",
      "2557 0 170.54806588741775\n",
      "2183.328745333007 -1.1229634858572917e-06 101.85303022936162\n",
      "2183 0 101.85264096484967\n",
      "MOD021KM.A2019181.0225.061.2019181131027_5.tif\n",
      "1220.5388665350372 -3.620474617300573e-06 228.49134864724567\n",
      "1221 0 228.49113859783455\n",
      "2391.8741003462324 -1.3144685850921571e-06 253.69330089603284\n",
      "2392 0 253.6933709439287\n",
      "MOD021KM.A2019181.0225.061.2019181131027_8.tif\n",
      "1267.7310246833492 -4.043364226152953e-06 273.43146055280715\n",
      "1268 0 273.43142635805424\n",
      "1219.9762766818158 -3.0053316181612945e-06 238.9936652570275\n",
      "1220 0 238.99392669491115\n",
      "1245.1183462824913 -2.3093714275506617e-06 212.99675972246663\n",
      "1245 0 212.99659873356183\n",
      "1070.578435778038 -1.5976527205679667e-06 121.42505233877523\n",
      "1071 0 121.4249508009289\n",
      "939.6078787656569 -1.1878126156737459e-06 70.88575200007351\n",
      "940 0 70.88560501832322\n"
     ]
    }
   ],
   "source": [
    "folder_path = \"./hdf-output/\"\n",
    "for filename in os.listdir(folder_path):\n",
    "    #print(filename[-6:-4])\n",
    "    if filename[-6:-4] == \"_5\" or filename[-6:-4] == \"_8\": \n",
    "        # MODDIS前7个波段的数据\n",
    "        print(filename)\n",
    "        hdf_name = \"./hdfs/\" + filename[:-6] + \".hdf\"\n",
    "        hdf_data = gdal.Open(hdf_name)\n",
    "        solar_zenith_data = gdal.Open(hdf_data.GetSubDatasets()[14][0])\n",
    "        # 太阳天顶角 它的角度范围是0-18000, 所以要除200 并化到0-pi\n",
    "        zenith = np.deg2rad(\n",
    "            solar_zenith_data.GetRasterBand(1).ReadAsArray() / 200)\n",
    "\n",
    "        data = gdal.Open(folder_path + filename)\n",
    "        rad_off = data.GetMetadataItem(\"radiance_offsets\")\n",
    "        rad_off = np.array(\n",
    "                rad_off.replace(' ', '').split(\",\"), dtype=np.float64)\n",
    "        rad_scales = data.GetMetadataItem(\"radiance_scales\")\n",
    "        rad_scales = np.array(\n",
    "                rad_scales.replace(' ', '').split(\",\"), dtype=np.float64)\n",
    "        ref_off = data.GetMetadataItem(\"reflectance_offsets\")\n",
    "        ref_off = np.array(\n",
    "            ref_off.replace(' ', '').split(\",\"), dtype=np.float64)\n",
    "        ref_scale = data.GetMetadataItem(\"reflectance_scales\")\n",
    "        ref_scale = np.array(\n",
    "            ref_scale.replace(' ', '').split(\",\"), dtype=np.float64)\n",
    "        band_size = data.RasterCount\n",
    "        driver = gdal.GetDriverByName(\"GTiff\")\n",
    "        out_ds = driver.Create(\"./rad-corr/\" + filename, \n",
    "                                data.RasterXSize, data.RasterYSize,\n",
    "                                band_size*2, gdal.GDT_Int64)\n",
    "        for i in range(1,band_size+1):\n",
    "            band = data.GetRasterBand(i).ReadAsArray()\n",
    "            corr_band = (band - rad_off[i-1]) * rad_scales[i-1]\n",
    "            # 辐射\n",
    "            out_band = out_ds.GetRasterBand(i)\n",
    "            # 反射\n",
    "            ref_band = out_ds.GetRasterBand(i + band_size)\n",
    "            # 重采样太阳天顶角的数据\n",
    "            z = cv2.resize(zenith, dsize=corr_band.shape,\n",
    "                            interpolation=cv2.INTER_LINEAR)\n",
    "            # 计算反射率\n",
    "            ref = ref_scale[i-1] * (band - ref_scale[i-1]) / np.cos(z).T\n",
    "            # 将反射率映射到扩大1000背, 因为类型是Int32\n",
    "            ref = ref * 1000\n",
    "            out_band.WriteArray(corr_band)\n",
    "            ref_band.WriteArray(ref)\n",
    "            print(ref.max(), ref.min(), ref.mean())\n",
    "            ref2 = out_ds.GetRasterBand(i + band_size).ReadAsArray()\n",
    "            print(ref2.max(), ref2.min(), ref2.mean())\n",
    "            pass\n",
    "        out_ds.SetProjection(data.GetProjection())\n",
    "        out_ds.SetGeoTransform(data.GetGeoTransform())\n",
    "        out_ds.FlushCache()\n",
    "        out_ds = None\n",
    "        pass"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}