tile_dem_file.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas as gpd\n",
    "from subprocess import Popen\n",
    "import os\n",
    "from oggm import utils"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "ifile = '/home/mowglie/disk/OGGM_INPUT/download_cache/data.pgc.umn.edu/elev/dem/setsm/REMA/mosaic/v1.1/100m/REMA_100m_dem.tif'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf = gpd.read_file('REMA_Tile_Index_Rel1.1.shp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWUAAAD8CAYAAACvm7WEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAEl9JREFUeJzt3X+s3XV9x/Hne61UjdMWuMGm7V1rbLaVxUw4gS5bFiMOCroVF3TlD6lI7DYw02TJbEcim/oHaDI2Jv5oBrEYYmWoaRMhXUGI+6fArSJQtPaCEtpUmJQfS0xkxff+OJ8L33u9vfeec8+553PvfT6Sk/s97+/n+/18Tnv64svn++NGZiJJqsNvDXoAkqTXGMqSVBFDWZIqYihLUkUMZUmqiKEsSRUxlCWpIoayJFXEUJakiiwd9ABqceaZZ+batWsHPQxJFTt48OAvMnOon30YysXatWsZGRkZ9DAkVSwinup3H05fSFJFDGVJqoihLEkVMZQlqSKGsiRVxFCWpIoYypJUEUNZkirizSNa9NZu/860bX52/XvnYCSSR8qSVBVDWZIqYihLUkUMZUmqiCf6tKDM5KRdTfudjicYFx+PlCWpIoayJFXEUJakijinrHlrUPO8c2myz+g888I26yPliFgTEfdFxOMRcSgiPl7qp0fE/og4Un6uKPWIiJsiYjQiHomIcxr72lraH4mIrY36uRHxaNnmpoiIbvuQpJr1YvriJPD3mbkB2AhcExEbgO3AvZm5Hri3vAe4GFhfXtuAL0E7YIHrgPOB84DrxkK2tPloY7tNpd5RH5JUu1mHcmYez8zvl+X/BX4ErAI2A7tKs13ApWV5M3Bbth0AlkfESuAiYH9mnsjM54H9wKay7s2ZeSAzE7htwr466UOSqtbTOeWIWAu8E3gAOCszj5dVPwfOKsurgKcbmx0ttanqRyep00Ufxxs1ImIb7SNphoeHZ/YhNTCLYQ55Jvr15+BcdR16dvVFRLwJ+Cbwicx8qbmuHOFmr/qaTDd9ZObOzGxlZmtoaKhPI5OkmetJKEfE62gH8u2Z+a1SfmZsyqD8fLbUjwFrGpuvLrWp6qsnqXfThyRVrRdXXwRwC/CjzPyXxqq9wNgVFFuBPY36FeUKiY3Ai2UKYh9wYUSsKCf4LgT2lXUvRcTG0tcVE/bVSR+SVLVezCn/MfAh4NGIeLjU/hG4HrgjIq4CngI+WNbdBVwCjAK/BK4EyMwTEfEZ4KHS7tOZeaIsXw18FXgDcHd50WkfklS7aE/FqtVq5cjIyKCHocKTenPPE33Ti4iDmdnqZx/eZi1JFTGUJakihrIkVcQHEvWRD5PRfDJX8/j+G5iaR8qSVBFDWZIqYihLUkUMZUmqiCf6ZsEbHKTOzeTfzWI+GeiRsiRVxFCWpIoYypJUEeeU59jE+bTFPHfW5Py81OaRsiRVxFCWpIoYypJUEUNZkipiKEtSRQxlSaqIoSxJFTGUJaki3jwyYN5MIqnJI2VJqoihLEkVMZQlqSLOKc/CxPnfXjxUx9+ALS1uHilLUkUMZUmqiKEsSRUxlCWpIp7omwf87b9abBbzTVUeKUtSRQxlSaqIoSxJFXFOuYf6cTPJTE3XV21zcv72amlyPTlSjohbI+LZiHisUTs9IvZHxJHyc0WpR0TcFBGjEfFIRJzT2GZraX8kIrY26udGxKNlm5siIrrtQ5Jq1qvpi68CmybUtgP3ZuZ64N7yHuBiYH15bQO+BO2ABa4DzgfOA64bC9nS5qON7TZ104ck1a4noZyZ3wNOTChvBnaV5V3ApY36bdl2AFgeESuBi4D9mXkiM58H9gObyro3Z+aBzEzgtgn76qQPSapaP+eUz8rM42X558BZZXkV8HSj3dFSm6p+dJJ6N30cb9SIiG20j6QZHh7u4KPNzGTzuIOaS+3Xtc7ODUu9NSdXX5Qj3Kytj8zcmZmtzGwNDQ31aWSSNHP9DOVnxqYMys9nS/0YsKbRbnWpTVVfPUm9mz4kqWr9DOW9wNgVFFuBPY36FeUKiY3Ai2UKYh9wYUSsKCf4LgT2lXUvRcTGctXFFRP21UkfklS1nswpR8TXgXcBZ0bEUdpXUVwP3BERVwFPAR8sze8CLgFGgV8CVwJk5omI+AzwUGn36cwcO3l4Ne0rPN4A3F1edNqHJNUu2lOxarVaOTIy0vd+PDEm9cYgboiKiIOZ2epnH95mLUkVMZQlqSKGsiRVxAcSzbFBPrRImq9qe6BWP3mkLEkVMZQlqSKGsiRVxFCWpIoYypJUEUNZkipiKEtSRQxlSaqIDyQq5uqBRNPxZhJpZnwgkSSp7wxlSaqIoSxJFXFOuahlTnkmnHeWftNczDE7pyxJi4yhLEkVMZQlqSKGsiRVxN88ImlBmHgCfL7+thKPlCWpIoayJFXEUJakinjzSDGfbh6ZjDeUSNOb7TyzN49I0iJjKEtSRQxlSaqIoSxJFTGUJakihrIkVcRQlqSKGMqSVBEfSDQPeaNIb83khgL/zDVXFvSRckRsiojDETEaEdsHPR5Jms6CDeWIWALcDFwMbAAuj4gNgx2VJE1twYYycB4wmplPZubLwG5g84DHJElTWshzyquApxvvjwLnNxtExDZgG8Dw8PDcjUxS3/mQ+3koM3dmZiszW0NDQ4MejiQt6FA+BqxpvF9dapJUrYUcyg8B6yNiXUScBmwB9g54TJI0pQU7p5yZJyPiY8A+YAlwa2YeGvCwJGlKCzaUATLzLuCuQY9DdfPGENVkIU9fSNK8YyhLUkUMZUmqyIKeU5a0eM3kXEGNN5h4pCxJFTGUJakihrIkVcRQlqSKGMqSVBFDWZIqYihLUkUMZUmqiDePSFo0arxZZCKPlCWpIoayJFXEUJakijinPA91My/mg9yl+cEjZUmqiKEsSRUxlCWpIoayJFXEUJakihjKklQRQ1mSKmIoS1JFvHlE0oI0Hx4+NBmPlCWpIoayJFXEUJakijinvEhMNr/mQ4q0kMzXOeSJPFKWpIoYypJUEUNZkipiKEtSRQxlSarIrEI5Ij4QEYci4tcR0ZqwbkdEjEbE4Yi4qFHfVGqjEbG9UV8XEQ+U+jci4rRSX1bej5b1a7vtQ5JqN9sj5ceAvwS+1yxGxAZgC3A2sAn4YkQsiYglwM3AxcAG4PLSFuAG4MbMfDvwPHBVqV8FPF/qN5Z23fYhSVWbVShn5o8y8/AkqzYDuzPzV5n5U2AUOK+8RjPzycx8GdgNbI6IAN4N3Fm23wVc2tjXrrJ8J3BBad9RH7P5nJI0V/p188gq4EDj/dFSA3h6Qv184Azghcw8OUn7VWPbZObJiHixtO+0j98QEduAbQDDw8Mz/GgLx8SL7b2ZRBq8aUM5Iu4B3jrJqmszc0/vhzR3MnMnsBOg1WrlgIcjSdOHcma+p4v9HgPWNN6vLjVOUX8OWB4RS8vRcrP92L6ORsRS4C2lfad9SFL1+nVJ3F5gS7lyYh2wHngQeAhYX660OI32ibq9mZnAfcBlZfutwJ7GvraW5cuA75b2HfXRp88pST01qznliHg/8O/AEPCdiHg4My/KzEMRcQfwOHASuCYzXynbfAzYBywBbs3MQ2V3nwR2R8RngR8At5T6LcDXImIUOEE7ZOmyD0mqWrQPOtVqtXJkZGTQwxgoT/RpPhnEU+Ei4mBmtqZv2T3v6JOkihjKklQRQ1mSKuJvHtGrupmjcx5ac2Gh/FaRmfBIWZIqYihLUkUMZUmqiHPKmpWZzPV1M+88k9++3at5RufFVROPlCWpIoayJFXEUJakihjKklQRH0hU+EAi9ZonELtX680iPpBIkhYZQ1mSKmIoS1JFvHlE6gHnj9UrHilLUkUMZUmqiKEsSRUxlCWpIp7ok3pgJk+10+RqvVFkUDxSlqSKGMqSVBFDWZIq4pyy1APOH6tXPFKWpIoYypJUEUNZkirinLKkOeV1yVPzSFmSKmIoS1JFDGVJqoihLEkVmdWJvoj4PPDnwMvAE8CVmflCWbcDuAp4Bfi7zNxX6puAfwOWAP+RmdeX+jpgN3AGcBD4UGa+HBHLgNuAc4HngL/KzJ9104ek/vIk3uzN9kh5P/AHmfkO4CfADoCI2ABsAc4GNgFfjIglEbEEuBm4GNgAXF7aAtwA3JiZbweepx22lJ/Pl/qNpV23fUhS1WYVypn5X5l5srw9AKwuy5uB3Zn5q8z8KTAKnFdeo5n5ZGa+TPvIeHNEBPBu4M6y/S7g0sa+dpXlO4ELSvuO+pjN55SkudLLOeWPAHeX5VXA0411R0vtVPUzgBcaAT9WH7evsv7F0r7TPiSpetPOKUfEPcBbJ1l1bWbuKW2uBU4Ct/d2eP0VEduAbQDDw8MDHo00/ziH3HvThnJmvmeq9RHxYeB9wAWZmaV8DFjTaLa61DhF/TlgeUQsLUfDzfZj+zoaEUuBt5T2nfYx2WfbCewEaLVaOVkbSZpLs5q+KFc5/APwF5n5y8aqvcCWiFhWrqpYDzwIPASsj4h1EXEa7RN1e0uY3wdcVrbfCuxp7GtrWb4M+G5p31Efs/mckjRXZvvsiy8Ay4D97XNvHMjMv8nMQxFxB/A47WmNazLzFYCI+Biwj/blardm5qGyr08CuyPis8APgFtK/RbgaxExCpygHbJ02YckVS1em3FY3FqtVo6MjAx6GJqnFutD7hfbnHJEHMzMVj/78ClxUg/0Kpz6Fe4TxzeTfhZb4NbC26wlqSKGsiRVxFCWpIp4oq/wRJ+k6czFiT6PlCWpIoayJFXEUJakihjKklQRQ1mSKmIoS1JFDGVJqoihLEkVMZQlqSLe0VdExP8ATw16HDN0JvCLQQ+iA463vxxvfzXH+zuZOdTPzgzleSgiRvp9q2cvOd7+crz9NdfjdfpCkipiKEtSRQzl+WnnoAfQIcfbX463v+Z0vM4pS1JFPFKWpIoYynMoIj4fET+OiEci4tsRsbyxbkdEjEbE4Yi4qFHfVGqjEbG9UV8XEQ+U+jci4rRSX1bej5b1a2fRxwci4lBE/Doixp19rnG8vdDPfZf93xoRz0bEY43a6RGxPyKOlJ8rSj0i4qYylkci4pzGNltL+yMRsbVRPzciHi3b3BQR0W0fZf2aiLgvIh4v34WP1zzmiHh9RDwYET8s4/3nUu/796+bPiaVmb7m6AVcCCwtyzcAN5TlDcAPgWXAOuAJYEl5PQG8DTittNlQtrkD2FKWvwz8bVm+GvhyWd4CfGMWffw+8LvA/UCr8TmqHG8P/n76tu9GH38KnAM81qh9Dthelrc3vheXAHcDAWwEHij104Eny88VZXlFWfdgaRtl24u76aMxtpXAOWX5t4GflL+bKsdcam8qy68DHiht5uLfS0d9nPI7MuigWqwv4P3A7WV5B7CjsW4f8Eflta9R31FeQfti9rGAf7Xd2LZleWlpF532MWGs9zM+lKse7yz+Tvq27wn9rGV8KB8GVpbllcDhsvwV4PKJ7YDLga806l8ptZXAjxv1V9t12scUY98D/Nl8GDPwRuD7wPn9/v7RxXf8VH/GTl8Mzkdo/9ceYBXwdGPd0VI7Vf0M4IXMPDmhPm5fZf2LpX2nfUxlvo13pvq576mclZnHy/LPgbOmGc9U9aOT1Lvp4zeU/+1+J+2jz2rHHBFLIuJh4FlgP+0j235//7r5jk9q6alWqDsRcQ/w1klWXZuZe0qba4GTwO1zObZT+BSwOiI2lvdvAd4YEfvHxqu5k5kZEX29JKqbPiLiTcA3gU9k5ktl2rfr/XWqkz4y8xXgD6N9zubbwO/1c2y9Zij3WGa+Z6r1EfFh4H3ABVn+fwY4BqxpNFtdapyi/hywPCKWlv/yNtuP7etoRCylHbLPTdHHXwP/lJkXlfHtKJ9jqkAe5HhP1UcvTNVnPz0TESsz83hErKR9hDfVeI4B75pQv7/UV0/Svps+XhURr6MdyLdn5rfmw5gBMvOFiLiP9lRCv79/3XzHJ9fr+TJfU84lbgIeB4Ym1M9m/EmFJ2mfUFhaltfx2kmFs8s2/8n4kwpXl+VrGH9S4Y5u+2iM737GzylXPd5Z/P30bd8T+lnL+DnlzzP+hNbnyvJ7GX9C68FSPx34Ke0TZivK8ull3cSTZpd000djbAHcBvzrhHqVYwaGgOVl+Q3Af9M+COr796/TPk75/Rh0UC2mFzBKe27p4fL6cmPdtbTnvg5Tzj6X+iW0z3g/QXsKZKz+tvJlHi1fhmWl/vryfrSsf9ss+ng/7bmxXwHPMP4ER3Xj7dHfUd/2Xfb/deA48H/lz/Yq2vOL9wJHgHt4LawCuLmM5VHG/4fxI+XPbBS4slFvAY+Vbb7AazeIddxHWf8nQAKPNL63l9Q6ZuAdwA/KeB8DPjWH/1467mOyl3f0SVJFvPpCkipiKEtSRQxlSaqIoSxJFTGUJakihrIkVcRQlqSKGMqSVJH/B+VDSG30oaBOAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gdf.plot();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1516, 1516)"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(gdf), len(gdf.tile.unique())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(len(gdf)):\n",
    "    sel = gdf.iloc[[i]]\n",
    "    \n",
    "    name = '{}_100m_v1.1'.format(sel.iloc[0].tile)\n",
    "    odir = os.path.join('/home/mowglie/disk/TMP_Data/REMA_100m_v1.1', name)\n",
    "    utils.mkdir(odir, reset=True)\n",
    "    ofile = os.path.join(odir, name + '_reg_dem.tif')\n",
    "    index_dir = os.path.join(odir, 'index')\n",
    "    shp_file = os.path.join(index_dir, name + '_index.shp')\n",
    "    utils.mkdir(index_dir)\n",
    "    sel.to_file(shp_file)\n",
    "    command = 'gdalwarp -cutline {} -crop_to_cutline -of GTiff {} {}'.format(shp_file, ifile, ofile)\n",
    "    Popen(command, shell=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "hide_input": false,
  "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.5.2"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": false,
   "sideBar": true,
   "skip_h1_title": true,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}