WMTS_Practice_Examples.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\"\n",
    "--------------------------------------------------------------------\n",
    "Author: Aaron Valoroso\n",
    "Date: July 16th, 2018\n",
    "License: MIT\n",
    "Email: valoroso99@gmail.com\n",
    "---------------------------------------------------------------------\n",
    "\"\"\"\n",
    "# This will only work if you assume that the Mercator projection has equal distances between each latitude gratitude.\n",
    "import numpy as np\n",
    "from itertools import product\n",
    "\n",
    "url = 'http://c.tile.openstreetmap.org/{0}/{1}/{2}.png'\n",
    "zoom_level = 0\n",
    "\n",
    "long_1= -89.033203\n",
    "lat_1 = 31.278551\n",
    "\n",
    "long_2 = -72.685547\n",
    "lat_2 = 23.966176\n",
    "\n",
    "# Scrub our points and get them into our format.\n",
    "the_lats = [lat_1, lat_2]\n",
    "the_lats = [np.around((90 - y), decimals=2) for y in the_lats]\n",
    "the_lats.sort()\n",
    "print(the_lats)\n",
    "\n",
    "the_longs = [long_1, long_2]\n",
    "the_longs = [np.around((180 + x), decimals=2) for x in the_longs]\n",
    "the_longs.sort()\n",
    "print(the_longs)\n",
    "\n",
    "# Firgure out our zoom points\n",
    "our_y = (170/(2**zoom_level))\n",
    "our_x = (360/(np.power(2, zoom_level)))\n",
    "print(our_y)\n",
    "print(our_x)\n",
    "\n",
    "# Firgure out our tles\n",
    "tile_lats = [int(y / our_y) for y in the_lats]\n",
    "tile_longs = [int(x / our_x) for x in the_longs]\n",
    "print(tile_lats)\n",
    "print(tile_longs)\n",
    "\n",
    "# Get the range of tiles\n",
    "lat_range = range(tile_lats[0], (tile_lats[1] + 1))\n",
    "long_range = range(tile_longs[0], (tile_longs[1] + 1))\n",
    "print(lat_range)\n",
    "print(long_range)\n",
    "\n",
    "# Get the pictures\n",
    "for x, y in product(long_range, lat_range):\n",
    "    print(x,y)\n",
    "    print(url.format(zoom_level, x, y))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# What we used to get the equations below. https://en.wikipedia.org/wiki/Web_Mercator\n",
    "import numpy as np\n",
    "from itertools import product\n",
    "\n",
    "url = 'http://c.tile.openstreetmap.org/{0}/{1}/{2}.png'\n",
    "zoom_level = 10\n",
    "\n",
    "long_1= -89.033203\n",
    "lat_1 = 31.278551\n",
    "\n",
    "long_2 = -72.685547\n",
    "lat_2 = 23.966176\n",
    "\n",
    "the_lats = [lat_1, lat_2]\n",
    "the_lats = [int(1/(2*np.pi)*2**zoom_level*(np.pi-np.log(np.tan(np.pi/4+np.radians(y)/2)))) for y in the_lats]\n",
    "print(the_lats)\n",
    "\n",
    "# Scrub our points and get them into our format.\n",
    "the_longs = [long_1, long_2]\n",
    "the_longs = [int(2**(zoom_level-1)*(x/180+1)) for x in the_longs]\n",
    "print(the_longs)\n",
    "\n",
    "# Get the range of tiles\n",
    "lat_range = range(the_lats[0], (the_lats[1] + 1))\n",
    "long_range = range(the_longs[0], (the_longs[1] + 1))\n",
    "print(lat_range)\n",
    "print(long_range)\n",
    "\n",
    "# Get the pictures\n",
    "for x, y in product(long_range, lat_range):\n",
    "    print(x,y)\n",
    "    print(url.format(zoom_level, x, y))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import requests\n",
    "def download(destination, url, filename):\n",
    "    temp_path = None\n",
    "    if destination[-1] != '/':\n",
    "        destination += '/'\n",
    "\n",
    "    if os.path.exists(destination):\n",
    "        response = requests.get(url, verify=False)\n",
    "        temp_path = destination + filename\n",
    "        if response.status_code == 200:\n",
    "            # Have to download the files in order, not when they finish. Download each file at a time.\n",
    "            open(temp_path, 'wb').write(response.content)\n",
    "#             chunk_size = 64 * 1024\n",
    "#             with open(temp_path, 'wb') as f:\n",
    "#                 for content in response.iter_content(chunk_size):\n",
    "#                     f.write(content)\n",
    "    else:\n",
    "        raise ValueError(\"Sorry, but the download destination that you requested is not a path.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from PIL import Image\n",
    "def image_stitch(dir_path):\n",
    "    \n",
    "    dir_path += '/'\n",
    "    \n",
    "    #----------------------------------------------------------------------------\n",
    "    list_of_pictures = list()\n",
    "    lists_of_list_of_pictures = list()\n",
    "    list_iterator = 0\n",
    "    stored_number = ''\n",
    "    temp_number = ''\n",
    "    \n",
    "    for iterator, filename in enumerate(os.listdir(dir_path)):\n",
    "        temp_number = filename[0] + filename[1]\n",
    "        if iterator == 0:\n",
    "            stored_number = temp_number\n",
    "            \n",
    "        if temp_number != stored_number:\n",
    "            lists_of_list_of_pictures.append(list_of_pictures)\n",
    "            stored_number = temp_number\n",
    "            list_of_pictures = list()\n",
    "            \n",
    "        list_of_pictures.append(filename)\n",
    "\n",
    "    lists_of_list_of_pictures.append(list_of_pictures)\n",
    "    stored_number = temp_number\n",
    "    #----------------------------------------------------------------------------\n",
    "    total_width_of_list = 0\n",
    "    tile_width = 0\n",
    "    growing_width = 0\n",
    "    \n",
    "    total_height_of_list = 0\n",
    "    tile_height = 0\n",
    "    growing_height = 0\n",
    "\n",
    "    heights = list()\n",
    "    widths = list()\n",
    "    stitched_image = None\n",
    "    new_image = None\n",
    "    creating_a_bigger_picture = None\n",
    "    \n",
    "    for item in lists_of_list_of_pictures:\n",
    "        images = [ Image.open(dir_path + i) for i in item ]\n",
    "\n",
    "        for second_item in images:\n",
    "            (width, height) = second_item.size\n",
    "            widths.append(width)\n",
    "            heights.append(height)\n",
    "        \n",
    "        # We need to make sure that the individual height and width of each time is an int and not a float.\n",
    "        # Hence why I had to do the int() around our division. \n",
    "        total_height_of_list = sum(heights)\n",
    "        tile_height = int(total_height_of_list/len(heights))\n",
    "        total_width_of_list = sum(widths)\n",
    "        tile_width = int(total_width_of_list/len(widths))\n",
    "\n",
    "        # Create a new image for our picture. So our width is not going to change but our height \n",
    "        # will. Our height is the number of tiles on top of one another.        \n",
    "        new_image = Image.new('RGB', (tile_width, total_height_of_list))\n",
    "\n",
    "        # Here I am putting the rows of tiles together. \n",
    "        for iterator2, item in enumerate(images):\n",
    "            if iterator2 == 0:\n",
    "                new_image.paste(im=item, box=(0,0))\n",
    "            else:\n",
    "                growing_height += 256\n",
    "                new_image.paste(im=item, box=(0, growing_height))\n",
    "            \n",
    "#       Put our two stacked images together into one image.\n",
    "        if stitched_image is None:\n",
    "            stitched_image = new_image\n",
    "        else:\n",
    "            growing_width += 256\n",
    "            creating_a_bigger_picture = Image.new('RGB', ((total_width_of_list+256), total_height_of_list))\n",
    "            (width, height) = new_image.size\n",
    "            creating_a_bigger_picture.paste(im=stitched_image, box=(0,0))\n",
    "            creating_a_bigger_picture.paste(im=new_image, box=(growing_width,0))\n",
    "            stitched_image = creating_a_bigger_picture\n",
    "\n",
    "        widths = list()\n",
    "        heights = list()\n",
    "        images = list()\n",
    "        new_image = None\n",
    "        creating_bigger_picture = None\n",
    "        growing_height = 0 \n",
    "\n",
    "    stitched_image.show()\n",
    "    \n",
    "    return"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# What we used to get the equations below. https://en.wikipedia.org/wiki/Web_Mercator\n",
    "import numpy as np\n",
    "import shutil\n",
    "import os\n",
    "from itertools import product\n",
    "\n",
    "switch = False\n",
    "\n",
    "if switch == True:\n",
    "    url = 'http://c.tile.openstreetmap.org/{0}/{1}/{2}.png'\n",
    "    zoom_level = 7\n",
    "\n",
    "    long_1= -89.033203\n",
    "    lat_1 = 31.278551\n",
    "\n",
    "    long_2 = -72.685547\n",
    "    lat_2 = 23.966176\n",
    "\n",
    "    the_lats = [lat_1, lat_2]\n",
    "    the_lats = [int(1/(2*np.pi)*2**zoom_level*(np.pi-np.log(np.tan(np.pi/4+np.radians(y)/2)))) for y in the_lats]\n",
    "\n",
    "    the_longs = [long_1, long_2]\n",
    "    the_longs = [int(2**(zoom_level-1)*(x/180+1)) for x in the_longs]\n",
    "\n",
    "    # Get the range of tiles\n",
    "    lat_range = range(the_lats[0], (the_lats[1] + 1))\n",
    "    long_range = range(the_longs[0], (the_longs[1] + 1))\n",
    "\n",
    "dir_path = os.getcwd()\n",
    "dir_path += '/image_example'\n",
    "\n",
    "if switch == True:\n",
    "    if not os.path.exists(dir_path):\n",
    "        os.makedirs(dir_path)\n",
    "    else:\n",
    "        shutil.rmtree(dir_path)\n",
    "        os.makedirs(dir_path)\n",
    "\n",
    "    # Get the pictures\n",
    "    for x, y in product(long_range, lat_range):\n",
    "        filename = \"{0}-{1}.png\".format(x,y)\n",
    "        download(dir_path, url.format(zoom_level, x, y), filename)\n",
    "\n",
    "image_stitch(dir_path)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Version-1\n",
    "#                 stored_number = int(images[0][0])\n",
    "#                 image1 = Image.open(dir_path + images[0])\n",
    "#                 image2 = Image.open(dir_path + images[1])\n",
    "#                 (width1, height1) = image1.size\n",
    "#                 (width2, height2) = image2.size\n",
    "#                 result_width = width1 + width2\n",
    "#                 result_height = height1 + height2\n",
    "                \n",
    "#                 result = Image.new('RGB', (width1, result_height))\n",
    "#                 result.paste(im=image1, box=(0,0,width1,height1))\n",
    "#                 result.paste(im=image2, box=(0,height1,width1, result_height))\n",
    "#                 result.show()\n",
    "\n",
    "# Version-2\n",
    "#     for iterator, filename in enumerate(os.listdir(dir_path)):\n",
    "#         images.append(filename)\n",
    "#         if iterator % 2 == 0:\n",
    "#             if iterator == 2:\n",
    "#                 imgs = [ Image.open(dir_path + i) for i in images ]\n",
    "#                 for item in imgs:\n",
    "#                     (width, height) = item.size\n",
    "#                     widths.append(width)\n",
    "#                     heights.append(height)\n",
    "                \n",
    "#                 total_height = sum(heights)\n",
    "#                 total_width = sum(widths)\n",
    "#                 print(total_height)\n",
    "#                 result = Image.new('RGB', (int(total_width/len(widths)), total_height))\n",
    "                \n",
    "#                 for iterator, item in enumerate(imgs):\n",
    "#                     the_height += 256\n",
    "#                     if iterator == 0:\n",
    "#                         result.paste(im=item, box=(0,0,256,256))\n",
    "#                     else:\n",
    "#                         result.paste(im=item, box=(0,(the_height - 256),256,the_height))\n",
    "                \n",
    "#                 result.show()\n",
    "                \n",
    "#                 return\n",
    "\n",
    "#Version-3\n",
    "#     images = list()\n",
    "#     heights = list()\n",
    "#     widths = list()\n",
    "    \n",
    "#     stitched_image = None\n",
    "#     stored_number = ''\n",
    "#     temp_number = ''\n",
    "#     dir_path += '/'\n",
    "#     bigger_picture_width = 0\n",
    "#     bigger_picture_height = 0\n",
    "    \n",
    "#     total_width = 0\n",
    "#     total_height = 0\n",
    "#     the_height = 0\n",
    "\n",
    "#     for iterator, filename in enumerate(os.listdir(dir_path)):\n",
    "#         temp_number = filename[0] + filename[1]\n",
    "#         if iterator == 0:\n",
    "#             stored_number = temp_number\n",
    "            \n",
    "#         if temp_number != stored_number:\n",
    "#             imgs = [ Image.open(dir_path + i) for i in images ]\n",
    "#             for item in imgs:\n",
    "#                 (width, height) = item.size\n",
    "#                 widths.append(width)\n",
    "#                 heights.append(height)\n",
    "\n",
    "#             total_height = sum(heights)\n",
    "#             total_width = sum(widths)\n",
    "\n",
    "#             result = Image.new('RGB', (int(total_width/len(widths)), total_height))\n",
    "\n",
    "#             for iterator, item in enumerate(imgs):\n",
    "#                 the_height += 256\n",
    "#                 if iterator == 0:\n",
    "#                     result.paste(im=item, box=(0,0,256,256))\n",
    "#                 else:\n",
    "#                     result.paste(im=item, box=(0,(the_height - 256),256,the_height))\n",
    "\n",
    "#             # Put our two stacked images together\n",
    "#             if stitched_image is None:\n",
    "#                 stitched_image = result\n",
    "#             else:\n",
    "#                 bigger_picture_width += 256\n",
    "#                 new_result = Image.new('RGB', ((total_width+256), total_height))\n",
    "#                 (width, height) = result.size\n",
    "#                 new_result.paste(im=stitched_image, box=(0,0))\n",
    "#                 new_result.paste(im=result, box=(bigger_picture_width,0))\n",
    "#                 stitched_image = new_result\n",
    "\n",
    "#             stored_number = temp_number\n",
    "#             widths[:] = []\n",
    "#             heights[:] = []\n",
    "#             imgs[:] = []\n",
    "#             images[:] = []\n",
    "#             result = None\n",
    "#             the_height = 0 \n",
    "        \n",
    "#         images.append(filename)\n",
    "\n",
    "#     stitched_image.show()\n",
    "\n",
    "#Version-4\n",
    "    \n",
    "#     dir_path += '/'\n",
    "    \n",
    "#     #----------------------------------------------------------------------------\n",
    "#     whole_directory = []\n",
    "#     for filename in os.listdir(dir_path):\n",
    "#         filename2, file_extension = os.path.splitext(filename)\n",
    "#         filename2 += '1'\n",
    "#         whole_directory.append(float(filename2.replace('-', '.')))\n",
    "\n",
    "#     whole_directory.sort()\n",
    "    \n",
    "#     new_directory = list()\n",
    "#     for filename in whole_directory:\n",
    "#         filename = str(filename)\n",
    "#         filename = filename.replace('.', '-')\n",
    "#         filename = filename[:-1]\n",
    "#         filename += '.png'\n",
    "#         new_directory.append(filename)\n",
    "#     #----------------------------------------------------------------------------\n",
    "#     if len(new_directory) > 1:\n",
    "#         list_of_pictures = list()\n",
    "#         lists_of_list_of_pictures = list()\n",
    "#         list_iterator = 0\n",
    "#         stored_number = ''\n",
    "#         temp_number = ''\n",
    "\n",
    "#         for iterator, filename in enumerate(new_directory):\n",
    "#             first_part, second_part = filename.split('-')\n",
    "#             temp_number = first_part\n",
    "#             if iterator == 0:\n",
    "#                 stored_number = temp_number\n",
    "\n",
    "#             if temp_number != stored_number:\n",
    "#                 lists_of_list_of_pictures.append(list_of_pictures)\n",
    "#                 stored_number = temp_number\n",
    "#                 list_of_pictures = list()\n",
    "\n",
    "#             list_of_pictures.append(filename)\n",
    "\n",
    "#         lists_of_list_of_pictures.append(list_of_pictures)\n",
    "#         stored_number = temp_number\n",
    "#         #----------------------------------------------------------------------------\n",
    "#         total_width_of_list = 0\n",
    "#         tile_width = 0\n",
    "#         growing_width = 0\n",
    "\n",
    "#         total_height_of_list = 0\n",
    "#         tile_height = 0\n",
    "#         growing_height = 0\n",
    "\n",
    "#         heights = list()\n",
    "#         widths = list()\n",
    "#         stitched_image = None\n",
    "#         new_image = None\n",
    "#         creating_a_bigger_picture = None\n",
    "\n",
    "#         for item in lists_of_list_of_pictures:\n",
    "#             images = [ Image.open(dir_path + i) for i in item ]\n",
    "\n",
    "#             for second_item in images:\n",
    "#                 (width, height) = second_item.size\n",
    "#                 widths.append(width)\n",
    "#                 heights.append(height)\n",
    "\n",
    "#             # We need to make sure that the individual height and width of each time is an int and not a float.\n",
    "#             # Hence why I had to do the int() around our division. \n",
    "#             total_height_of_list = sum(heights)\n",
    "#             tile_height = int(total_height_of_list/len(heights))\n",
    "#             total_width_of_list = sum(widths)\n",
    "#             tile_width = int(total_width_of_list/len(widths))\n",
    "\n",
    "#             # Create a new image for our picture. So our width is not going to change but our height \n",
    "#             # will. Our height is the number of tiles on top of one another.        \n",
    "#             new_image = Image.new('RGB', (tile_width, total_height_of_list))\n",
    "\n",
    "#             # Here I am putting the rows of tiles together. \n",
    "#             for iterator2, item in enumerate(images):\n",
    "#                 if iterator2 == 0:\n",
    "#                     new_image.paste(im=item, box=(0,0))\n",
    "#                 else:\n",
    "#                     growing_height += 256\n",
    "#                     new_image.paste(im=item, box=(0, growing_height))\n",
    "\n",
    "#             # Put our two stacked row images together into one image.\n",
    "#             if stitched_image is None:\n",
    "#                 stitched_image = new_image\n",
    "#             else:\n",
    "#                 growing_width += 256\n",
    "#                 creating_a_bigger_picture = Image.new('RGB', ((total_width_of_list+256), total_height_of_list))\n",
    "#                 (width, height) = new_image.size\n",
    "#                 creating_a_bigger_picture.paste(im=stitched_image, box=(0,0))\n",
    "#                 creating_a_bigger_picture.paste(im=new_image, box=(growing_width,0))\n",
    "#                 stitched_image = creating_a_bigger_picture\n",
    "\n",
    "#             widths = list()\n",
    "#             heights = list()\n",
    "#             images = list()\n",
    "#             new_image = None\n",
    "#             creating_bigger_picture = None\n",
    "#             growing_height = 0 \n",
    "#         stitched_image.show()\n",
    "#     #----------------------------------------------------------------------------\n",
    "\n",
    "#     return\n",
    "\n",
    "#Version-5\n",
    "# def image_stitch(current_path, max_width, max_height, min_width, min_height, images):\n",
    "#     # Create the new image with the incoming max width and height.\n",
    "#     stitched_image = Image.new('RGB', (max_width, max_height))\n",
    "    \n",
    "#     for file in images:\n",
    "#         incoming_image = Image.open(current_path + file)\n",
    "#         filename, file_extension = os.path.splitext(file)\n",
    "#         x, y = filename.split('-')\n",
    "#         x = (int(x) - min_width) * 256\n",
    "#         y = (int(y) - min_height) * 256\n",
    "#         stitched_image.paste(im=incoming_image, box=(x, y))\n",
    "#         incoming_image.close()\n",
    "#     stitched_image.show()\n",
    "\n",
    "#     return"
   ]
  }
 ],
 "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.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}