rasterio_vs_gdal.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import rasterio\n",
    "import subprocess\n",
    "import matplotlib.pyplot as plt\n",
    "from rasterio.enums import ColorInterp\n",
    "from rasterio.plot import show\n",
    "from rasterio.windows import Window"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "xoff, yoff, xsize, ysize = 51667529, 27435585, 10, 10\n",
    "source = \"/vsicurl/https://raw.githubusercontent.com/satellogic/telluric/master/tests/data/google.xml\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### gdal_translate"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "args = [\"gdal_translate\"]\n",
    "args.extend((\"-oo\", \"OVERVIEW_LEVEL=1\"))\n",
    "args.extend((\"-srcwin\", \"%s\" % xoff, \"%s\" % yoff, \"%s\" % xsize, \"%s\" % ysize))\n",
    "args.extend((\"-colorinterp\", \"red,green,blue\"))\n",
    "args.extend((source, \"gdal_translate.tif\"))\n",
    "subprocess.check_call(args)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### rasterio"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "with rasterio.open(source, OVERVIEW_LEVEL=1) as src:\n",
    "    out_window = Window(xoff, yoff, xsize, ysize)\n",
    "    out_kwargs = src.profile\n",
    "    out_kwargs.update(\n",
    "        {\n",
    "            \"driver\": \"GTiff\",\n",
    "            \"height\": xsize,\n",
    "            \"width\": ysize,\n",
    "            \"transform\": src.window_transform(out_window)\n",
    "\n",
    "        }\n",
    "    )\n",
    "    with rasterio.open(\"rasterio.tif\", \"w\", **out_kwargs) as out:\n",
    "        out.write(src.read(window=out_window, out_shape=(src.count, ysize, xsize),))\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### compare"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk0AAAE6CAYAAAALNWa8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAHyNJREFUeJzt3XuUZGV97vHn6ct0Mz0jI8IADuhoVC5LEcyEaKJGMShHXSIGIy4lGRPluLI0khU1Go63IB5zWN7WyeWIKKMRvKI5BhUBDUE9ChlhgIFBo4hcZIARkJnumenp7t/5Y781U7R9eXt6v7uqmO9nrb2qq2rX+/6quurXT+3aXdsRIQAAAMytr9MFAAAA9AJCEwAAQAZCEwAAQAZCEwAAQAZCEwAAQAZCEwAAQAZCU4+zHWk5NnP9dWn99xYurVZt93N1p2sBgOnoUfsGQhNmZHttagBXdrqWvdGr4RBAMwr0iI+l5aGaxkMXGuh0AehttgcjYlen6wCw7+imvtOqJSLO7HQtKI8tTV3G9rNt32h71Pa/2P58ejf0UVf+zva9tu+0ffoMt3+r7f9Kt99p+3rbpy6whrWSLkhn/yDNf1u67rZ0/izbN0nakS7/ULpuh+0x2z+0/by2Ma9Mt/uftq9K63zf9uPT9Utsf8L25lT3Hbb/bY4aZ53P9jpJf5pWfU+ad1267qm2v54ew/tsX2z7cQt5fAA0q+2jrzNt/1zSj9PlF9m+K/WMrba/Y/tpbbc70/bP0vVbUh86YjE9Yo5aHvbxnO2DbJ9v+3bbD6UedVL5RwslsaWpi9heIenfJK2QdKWklZL+sG2VtZLepSqofF7Su2cY5gmSbpR0maSDJL1C0mdtr4+I2zJLuVnS5ZJOlHSXpC9Lun/aOu9Ll9/cNu/VkrZIWi3pxZK+ZPuJEbG17XZvk/QFSY+X9HuS3i/pdEl/Iun1km6S9BVJj5X0nDlqnHW+dN+Pl3RUWueHkq6xfYikqyQtk3SJpElJp0o62vaxEbEz58EB0DEfkPQlSdvT+cer6pUPSnqqpOdL+qKko2w/SdJHVPWICyQ9StIzJR2qenrE9Fp2s90n6Wtpvusk/Yek0yR93fZzIuL/1fBYoAM6sqXJ9tm2b7C9wfZlth87y3qPS9dvsn1zW4JfZ/vn6fYbnHaCtn1y27jrbT+7bay/t70xLa9qu3zGseap/xjbP7B9k6utQsOLfUySl6oKTLdKOiEiXqQqALW8Jp1+ICJepyoQTfd2Sf+qKuTcJek+SUOqAkqWiLhG0kXp7E8j4syI+Ltpq30gIk6LiFem86+X9G1Jv5b0X5LGJB0o6WnTbvfxiHiNpPem88el08F0eqOkCyW9TlVonM2s80XERZKuSetdmuq/SFU4e7Skn0q6XXsenyNVNVtgXvSvjnpTRPxpRLwxnf9jVYFnq6Qb0mVHpt9Jq6f8UtUbsbdHxBMlfbemHjG9lnZrVAWmbZKeExGnS/oHVX9z37T3dx+dVnxLU/rIZG1ErG27+NyIeFe6/i9VbTGZ6Yn3GUnnRMTltpdJmmq77m0R8eVp639b0tciImwfo+odx5G2XyLpGZKOVRUgrrT9zYh4aI6xZrs/A5I+K+n0iLje9mMk1fXZ+qp0+uPYcyTlTZKePv36dPqTabUtUdVAnjrD2AfVVGPL99vmfYyqhjXTH4/p816XTh9Mp8vS6WckPU/SyarekYWkK2yfEhGj7QMscL52q9PpUWlp96Q5bod9FP2r67T3nSdLulZ7eki7g9L9e4+kv5T0rXSbH6vacrRxlvFXp9OcHvF9za41zh1t/euWdPr4OW6HLteRLU1tL3ZJGlH1B/JhbB8taSAiLk+32RYRY/OMu60tbLSPe7SkqyJiIj2Bb5A052fLtkdsf8r2Nbavs31yuuqFkm6IiOvTnL+KiMk573C+u9Jp+4vzyBmuPyKdPmXa7Y9WFZgmJP2Wqt9v6+MzL7CW1n2a7TnSvpn6OaoCzGZJh6hq7K1QNH3eiXQ6/Xc+ERGvUrUJ/ShJV6j6eHCmrWk5881U/23p9KsR4daianP9J2e8l8A09K+Oau87L1EVmDao2kJ/cNt1tt2vKrQeqCqo/L2q3vlXaZ3F9oi5Ps5vjXO47aXp51bf/sUct0OX69iO4LbPsX2Hqo+cZto35ymSHrT9lfSiPze9CFrOSZuyP2J7qG3cU2zfIunrkv4sXXy9pJNsL7V9oKrNrIfPM9ZZkr4TEcen9c+1PZLqCtvfsn2t7bcv/tHY7RJVf/yfbPsK25dKOqbt+tZHZn9r+wJJX512+y2q3s0OSPqQqv2SnryXtdyRTn/b9j/ZfsMc696TTg+S9GHt2SdgIV5te5OqfQ/eoj0f6z04w7o587Xqf63tj9l+vqqP/R6UdEr6/X3c9hVp3YMFZKJ/dYVWH3iKqn/1v2za9YdL+qXtL0n6G+0Jmq2eUrJHrFe1r9QySd+1/RlJb1YVhP9pAeOg20REkUXVE2aDqs+G708/b5D0omnrvVPS+2a4/amq9ld5oqoQcLGkP0/XHapqi8KQpE9LevcMt3+upCvazp+V5r9c1QvjzLnGUvWk39hW9+2qtoC8VdLPVe0/s1TSDyS9oMbH7dmq9usZlfQvqoJRqHqX1Kdqx+l7Jd0t6Q3pupB0bLr9X6jaAvPrdJsr0/Wt+7sunX/vPHX0pcfpobT+Jeny29L5501b//3p93yfqtDTWu/l6fpWHWvT+Zen87el88+S9F1Jv5I0rurd2NmSnK5v3c/VmfOtUrX5fGe6/K3p8mNU7Wx/T3qMN0n6qKRlpV4LLL23iP7VVcv013+6rF/S+alH3SnpVW3rHSvpgPRa/2XqKfdK+pykR6fb71WPmKmWmS5XtU/mp1QFrq2q9qF6aacfS5ZFPheLT1Dtp7JujusfJ2njDJc/U9J/tJ0/XdI/zjL+JbOMfaukA2e4/CJJL55rLEk/knTEDOucJunTbeffpWqfgroer/3bfu5LL9hoNVwWFpbmFvoXCwtL+9Kp/55r/8joZO3ZQa7df0paYbu1Y+8JSvvn2D40nVrVFouN6fyT0mWy/QxV775+Zbs/7fCotIPlMUqbcmcbS9WOg29uG++4tsufljaVD0j6A+3Zb6gO59v+nO13q3pXeaSqd0oX1ziHJMn2Sa6+/2n6wneJALOgfwH7rk59T9MHbR+hav+bXyj954ntNZLeGBGvj4hJ22+V9O30wv+RpE+k21+YmpFVbXpu/efKH0n6E9u7VH13xqsiImwPqvpcWao25b42IibmGetsVZtkb3D1nRs/V7Vp9QHbH1bVFEPSNyLi6zU+Nteq+ojtFao2J39B0v+IiJn27VmsZ6r6aGu6ByVdWmA+4JGA/gXso1r7iwAAAGAOHEYFAAAgA6EJAAAgQ5F9mpYMDcd+Iwv9mp6FGejvn3+lGvQ3ME9TH5BOTEzMv9IiTU1Nzb8Sdhvob+Z9y8DA4PwrLcJDv35Q28dGF/oFql1pcMlQDC8dKTpHE31FkgYHyu+22j/QzH2ZnCzfW8bHx4vP0UQfliQ1sOtNU8/jtD9fUQ89eP+WiJj3yBlFXlH7jSzT75340hJD7/boFSuKjt+yfwPzTE0284W8v7p/+jF36ze6bVvxOZrS11f+hbr/suXF55CklSvnOozf4n3ugv9TdPwmDS8d0W8/98Sicyxf3szv/eCDy39n64qGevG2BnrL7bffXnyOLVu2FJ9DkiYb+LvymP33Lz6HJA00EP6/efFns76pnY/nAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMgwUGzmi2NCStHNse9HxW8YGB4vPsaS/3K/hYfM0cF9iaLj4HBMxVXwOSTrkkEOKz7F9ezPP4+3j40XHnyr8em+SJfX3lX0/OTkxUXT8lrHR0eJzNNFXJCmmyr/uR5YuLT5HHHBA8TmaMmA3Ms/y5csbmScHW5oAAAAyEJoAAAAyEJoAAAAyEJoAAAAyEJoAAAAyEJoAAAAyEJoAAAAyEJoAAAAyZIUm2ytsf9n2LbY32X5W6cIAoA70LwB1yf0q6o9JujQiTrW9RFL5r00FgHrQvwDUYt7QZHt/Sc+VtFaSImJcUtljMgBADehfAOqU8/HcEyTdJ+kC29fZPt/2SOG6AKAO9C8AtckJTQOSniHpnyPiOEmjkt4xfSXbZ9heb3v9+M4dNZcJAHtl4f1rfGfTNQLoETmh6U5Jd0bE1en8l1U1oYeJiPMiYk1ErFnSwJHuASDDwvvXkqFGCwTQO+YNTRGxWdIdto9IF71A0s1FqwKAGtC/ANQp97/n3izpwvSfJ7dKel25kgCgVvQvALXICk0RsUHSmsK1AEDt6F8A6sI3ggMAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGTI/UbwBYmpKY1vL3vQ3piYLDp+S19f+Vw5NdzQsfomp5qZpzBPRSPzrFq1qvgcGzduLD6HJN1zzz1Fxx/f+cg5yK37+jQ0VPb4cxHNPIe3b99efI7BwcHic0jSwECRP1cP08R9GW6o309Nle/3jz344OJzSM08j3OxpQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACADoQkAACDDQM5Ktm+TtFXSpKSJiFgz1/oxNaXxHTsXX90cJndNFB1/t6koPsWunePF55Ck/v7+4nMMDGQ9pRZlyZIlxeeQpImJ8s+xoaGh4nNI0s6dZV+PU1H+dbK39qZ/7Srcv6ampoqO3xITk8XnGHAz772Hh4fLTzJZ/vfS19BLJRr423XoYYcVn0OSNm7c2Mg8ORbyF+75EbGlWCUAUA79C8Ci8fEcAABAhtzQFJIus/0j22eULAgAakb/AlCL3I/nnh0Rd9leKely27dExFXtK6RmdIYkDQ3vV3OZALDXFtS/hvdb2okaAfSArC1NEXFXOr1X0lclHT/DOudFxJqIWDM42MyOugAwnwX3ryXN7JwPoPfMG5psj9he3vpZ0gsldc+u7AAwC/oXgDrlfDx3sKSv2m6tf1FEXFq0KgCoB/0LQG3mDU0RcaukpzdQCwDUiv4FoE585QAAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAEAGQhMAAECGgSKDDg5q5cqVJYbeLSKKjt8yPj5efI7R0dHic0iS7eJzjIyMFJ9j2ZJlxeeQpM2bNxefY9WqVcXnkKT+/v6i42+8eknR8ZsUU6Gd23eUnaOh/jW5a6L8JFPN3JddO8v34tKvE0kaGCjyZ/c3DA2Wf01OTDTw/JI0NDTUyDw52NIEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQgdAEAACQITs02e63fZ3tS0oWBAB1o38BqMNCtjS9RdKmUoUAQEH0LwCLlhWabB8m6SWSzi9bDgDUi/4FoC65W5o+KuntkqYK1gIAJdC/ANRi3tBk+6WS7o2IH82z3hm219tev3PHjtoKBIC9tTf9a9eunQ1VB6DX5Gxp+n1JL7N9m6TPSzrB9menrxQR50XEmohYMzQ8XHOZALBXFty/BgeHmq4RQI+YNzRFxDsj4rCIWC3pNEnfiYjXFq8MABaJ/gWgTnxPEwAAQIaBhawcEVdKunK+9YaGhvSkJ/7WXpaUp7+/v+j4LQ888EDxOe6+++7ic0jS1q1bi8+xa+d4+TmGdhWfQ5JG77u/+BwHrTyk+BySNLLsUUXH72vo9bgYuf1rYGBABx7wmKK19PU18341IorPMT5e/jUvSWPbRovPYbv4HCMjI8XnkKShZcuKz7F58+bic0jSqlWrGpknB1uaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMhCaAAAAMgwUGbS/XytWrCgx9G6lx29Zvnx58Tl27dpVfA5JmpycLD5HX1/5HG67+BySNDY2WnyOe+65p/gckrR169ai409MTBQdv0mDg4M69NBDi86xbNmyouO39Pf3F5/jgQceKD6HJN19993F5yj9OpGk8fHx4nNIzfxd2bZlS/E5JGnlypWNzJODLU0AAAAZCE0AAAAZCE0AAAAZCE0AAAAZCE0AAAAZCE0AAAAZCE0AAAAZCE0AAAAZ5g1NtodtX2P7ets32X5fE4UBwGLRvwDUKecbwXdKOiEittkelPQ929+MiB8Wrg0AFov+BaA284amiAhJ29LZwbREyaIAoA70LwB1ytqnyXa/7Q2S7pV0eURcXbYsAKgH/QtAXbJCU0RMRsSxkg6TdLztp05fx/YZttfbXj86uu03BwGADlho/9q+vfyBmgH0pgX991xEPCjp3yWdNMN150XEmohYMzLSzBG8ASBXbv/ab7+R5osD0BNy/nvuINsr0s/7STpR0i2lCwOAxaJ/AahTzn/PHSrp07b7VYWsL0bEJWXLAoBa0L8A1Cbnv+dukHRcA7UAQK3oXwDqxDeCAwAAZCA0AQAAZCA0AQAAZCA0AQAAZCA0AQAAZCA0AQAAZCA0AQAAZCA0AQAAZMj5RvAFm5yc0rZtZQ/a29fXTN7bsWNH8TkiovgckmS7+Bzj4+PF55jcOlV8Dkmyyz/HSr9OWqamCj9mzTyFG2FLAwNFWuNuK1asKDp+k/MsX768+ByStGvXruJzTE5OFp+jqb9dTfT70dHyfx8lafM99zUyTw62NAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQgNAEAAGQYKDFoX581PDxcYujdSo/fMjY2VnyOHTt2FJ9DkkZHR4vPYbv4HP2DRZ62vzlPf3/xOcbHx4vPIZV/Hk9OTRYdv0lTU1H8Nblt27ai47f09ZV/X9xU/4qI4nM00b+aes1PTU0Vn8P9zWx3aer1kmPee2z7cNv/bvtm2zfZfksThQHAYtG/ANQp5y37hKS/johrbS+X9CPbl0fEzYVrA4DFon8BqM28W5oi4u6IuDb9vFXSJkmrShcGAItF/wJQpwV9IGl7taTjJF1dohgAKIX+BWCxskOT7WWSLpZ0ZkQ8NMP1Z9heb3v96Gj37LQFAAvpXzu2l//nDwC9KSs02R5U1XAujIivzLRORJwXEWsiYs3IyLI6awSAvbbQ/jW839JmCwTQM3L+e86SPilpU0R8uHxJAFAP+heAOuVsafp9SadLOsH2hrS8uHBdAFAH+heA2sz7lQMR8T1J5b/xCwBqRv8CUCcOowIAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJBhoMSgIWliarLE0HvmcNHhdyt9PyRpfGJX8TmkZh6zpSNLi8+x4oBHF59DkpYsHSk+x/adzfzut43tKDp+FB29WbY1ODhYdI7h4eGi4zc5z9jYWPE5JGnHjrLPYUkaHR0tPofdzB+vgYEif94fpr+/v/gckjQ+Pt7IPDnY0gQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJBh3tBk+1O277W9sYmCAKBO9DAAdcnZ0rRO0kmF6wCAUtaJHgagBvOGpoi4StL9DdQCALWjhwGoC/s0AQAAZKgtNNk+w/Z62+vHRrfVNSwAFNfev7ZvH+t0OQC6VG2hKSLOi4g1EbFm6ciyuoYFgOLa+9d++y3tdDkAuhQfzwEAAGTI+cqBz0n6gaQjbN9p+8/LlwUA9aCHAaiLI6L+Qe37JP1ikcMcKGlLDeWU1it1Sr1Ta6/UKfVOraXrfHxEHFRw/MbsY/1L6p1ae6VOqXdq7ZU6pS7pYUVCUx1sr4+INZ2uYz69UqfUO7X2Sp1S79TaK3U+UvTS490rtfZKnVLv1NordUrdUyv7NAEAAGQgNAEAAGTo5tB0XqcLyNQrdUq9U2uv1Cn1Tq29UucjRS893r1Sa6/UKfVOrb1Sp9QltXbtPk0AAADdpJu3NAEAAHSNrg5Ntt9r+y7bG9Ly4k7XNBfbf207bB/Y6VpmY/ts2zekx/My24/tdE0zsX2u7VtSrV+1vaLTNc3G9itt32R7ynbH/7tjOtsn2f6x7Z/afken69lX9Fr/krq/h/VK/5J6p4fRvxamq0NT8pGIODYt3+h0MbOxfbikF0q6vdO1zOPciDgmIo6VdImkd3e6oFlcLumpEXGMpJ9IemeH65nLRkmvkHRVpwuZzna/pH+U9N8kHS3p1baP7mxV+5Se6F9Sz/SwXulfUu/0MPrXAvRCaOoVH5H0dkldvZNYRDzUdnZEXVpvRFwWERPp7A8lHdbJeuYSEZsi4sedrmMWx0v6aUTcGhHjkj4v6eQO14Tu1PU9rFf6l9Q7PYz+tTC9EJrelDZvfsr2oztdzExsnyzproi4vtO15LB9ju07JL1G3f1OreXPJH2z00X0qFWS7mg7f2e6DM3o+v4l9VYP68H+JdHD9lbX9a+BTk4uSbavkHTIDFedJemfJZ2t6t3E2ZI+pOrJ17h56vxbVZu1u8JctUbE/42IsySdZfudkt4k6T2NFpjMV2da5yxJE5IubLK26XJqxb6nV/qX1Ds9rFf6l9Q7PYz+VZ+Oh6aI+MOc9Wx/QtVn2B0xW522nybpCZKuty1Vm2CvtX18RGxusMTdch9TVS/ib6hDTWe+Om2vlfRSSS+IDn83xgIe025zl6TD284fli5DDXqlf0m908N6pX9JvdPD6F/16eqP52wf2nb2FFU7rHWViLgxIlZGxOqIWK1q8+EzOhWY5mP7yW1nT5Z0S6dqmYvtk1TtX/GyiBjrdD097D8lPdn2E2wvkXSapK91uKZ9Qi/0L6m3eliv9C+JHlaTrutfHd/SNI//ZftYVZu3b5P03ztbziPCB20fIWlK1ZHc39jhembzD5KGJF2e3v3+MCK6slbbp0j635IOkvR12xsi4kUdLkuSFBETtt8k6VuS+iV9KiJu6nBZ+wr6V/16pX9JPdLD6F8LwzeCAwAAZOjqj+cAAAC6BaEJAAAgA6EJAAAgA6EJAAAgA6EJaFDuwTFnO0il7RfYvjYdsPR7tp+ULn9uunzC9qnTxnpcOrjpJts3216dLl9n++dtB5Q9dp7aX9O27oZ0H+a8DYBHDvqXpIhgYWEpsEh6nqR10y47StIRkq6UtGaW2/VL+pmkJ0paIul6SUen634i6aj081+0xpe0WtIxkj4j6dRp410p6cT08zJJS9PP66avu4D79jRJP+v0Y8zCwlJmoX/NvLClCWhQ5B0cc66DVIakR6Wf95f0yzTubRFxg6rvr9ktHRF8ICIuT+tti3m+aM/2SDpW2jW2r0vHJZvu1akuAPsI+hcfzwHdaK6DVL5e0jds3ynpdEkfnGesp0h60PZXUgM513Z/2/XnuDqg7EdsD6XLzpL0nYg4XtLzJZ1re2TauK+S9LmF3zUAj3CP6P5FaAJqZvtq2xsknS/pZW2fodfxLbt/JenFEXGYpAskfXie9QckPUfSWyX9jqpN5mvTde+UdGS6/ABJf5Muf6Gkd6T7cKWkYUmPaw1o+3cljUVEVx4WBMDeo3/NXxCAGkXE70qS7edJWhsRaxc4xIwHqbR9kKSnR8TV6fIvSLp0nrHulLQhIm5NNf2rpGdK+mRE3J3W2Wn7AlWNSZIs6Y/m2Ax/mtjKBDwi0b/mxpYmoPvMdpDKByTtb/spab0TJW3KGGtFaliSdIKkm6U9B5S1bUkv154Dyn5L0pvT5bJ9XGsw232S/ljszwRgZo/o/kVoAhpk+5T0ef6zVB0c81vp8sfa/oZUHaRSUusglZskfTEibkqXv0HSxbavV7VPwNvS7X8njftKSR+3fVMaa1LVO7Bv275R1buwT6RyLkyX3SjpQEnvT5efLWlQ0g1pnLPb7sJzJd3ReucHYN9B/+KAvQAAAFnY0gQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJCB0AQAAJDh/wNNPTxDHV7YvAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 720x360 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "with rasterio.open(\"gdal_translate.tif\") as r1:\n",
    "    with rasterio.open(\"rasterio.tif\") as r2:\n",
    "        fig, axs = plt.subplots(1, 2, figsize=(10, 5))\n",
    "        show(r1, ax=axs[0], title=\"gdal_translate\")\n",
    "        show(r2, ax=axs[1], title=\"rasterio\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r1.transform == r2.transform"
   ]
  }
 ],
 "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.6.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}