{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "In this demo, we follow [Tyler Sloan](https://quorumetrix.blogspot.com/2018/06/visualizing-lidar-data-with-datashader.html)'s walkthrough where he downloads LiDAR data from the City of Montreal (featuring the 1976 Olympic Stadium) and visualizes it ([original notebook here](https://github.com/tsloan1377/montreal_open_data/blob/master/lidar_datashader_blog.ipynb)).\n", "\n", "Make sure you've worked through the [Montreal data demo](Demo%20using%20Montreal%20LiDAR%20data.ipynb) first." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Imports" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "from mpl_toolkits.mplot3d import Axes3D\n", "\n", "import os\n", "import imageio\n", "from laspy.file import File\n", "import datashader as ds\n", "import datashader.transfer_functions as tf\n", "from matplotlib import cm" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We already downloaded the tiles around the [site of Avebury](https://en.wikipedia.org/wiki/Avebury) from the [data.gov.uk](https://environment.data.gov.uk/ds/survey/index.jsp#/survey?grid=SU17) service. These tiles are in the [aveburytiles](/aveburytiles) folder. Unzip them, and explore! Try changing the various `class` parameters in the visualizations. Reuse code blocks" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [], "source": [ "!LAStools/bin/laszip -i aveburtytiles/SU1477_P_8065_20120113_20120113.laz -o avebury1.las" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Load LIDAR file with LasPy" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [], "source": [ "sample_data = 'avebury1.las'\n", "export_path = 'export//'" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", " | X | \n", "Y | \n", "Z | \n", "class | \n", "
---|---|---|---|---|
0 | \n", "41468781 | \n", "17785233 | \n", "17511 | \n", "2 | \n", "
1 | \n", "41468716 | \n", "17785165 | \n", "17517 | \n", "3 | \n", "
2 | \n", "41468734 | \n", "17785159 | \n", "17522 | \n", "3 | \n", "
3 | \n", "41468793 | \n", "17785220 | \n", "17512 | \n", "1 | \n", "
4 | \n", "41469044 | \n", "17785222 | \n", "17505 | \n", "2 | \n", "
5 | \n", "41468980 | \n", "17785156 | \n", "17515 | \n", "3 | \n", "
6 | \n", "41468912 | \n", "17785083 | \n", "17507 | \n", "2 | \n", "
7 | \n", "41468847 | \n", "17785016 | \n", "17514 | \n", "1 | \n", "
8 | \n", "41468853 | \n", "17784998 | \n", "17511 | \n", "2 | \n", "
9 | \n", "41468914 | \n", "17785062 | \n", "17511 | \n", "2 | \n", "
10 | \n", "41468977 | \n", "17785128 | \n", "17514 | \n", "3 | \n", "
11 | \n", "41469038 | \n", "17785192 | \n", "17514 | \n", "3 | \n", "
12 | \n", "41469312 | \n", "17785213 | \n", "17512 | \n", "1 | \n", "
13 | \n", "41469244 | \n", "17785142 | \n", "17508 | \n", "2 | \n", "
14 | \n", "41469178 | \n", "17785073 | \n", "17508 | \n", "2 | \n", "
15 | \n", "41469113 | \n", "17785004 | \n", "17511 | \n", "2 | \n", "
16 | \n", "41469046 | \n", "17784935 | \n", "17512 | \n", "8 | \n", "
17 | \n", "41468980 | \n", "17784866 | \n", "17513 | \n", "2 | \n", "
18 | \n", "41468984 | \n", "17784849 | \n", "17515 | \n", "2 | \n", "
19 | \n", "41469045 | \n", "17784912 | \n", "17512 | \n", "2 | \n", "
20 | \n", "41469107 | \n", "17784976 | \n", "17511 | \n", "2 | \n", "
21 | \n", "41469167 | \n", "17785039 | \n", "17508 | \n", "2 | \n", "
22 | \n", "41469230 | \n", "17785104 | \n", "17511 | \n", "1 | \n", "
23 | \n", "41469291 | \n", "17785169 | \n", "17509 | \n", "8 | \n", "
24 | \n", "41469354 | \n", "17785235 | \n", "17514 | \n", "3 | \n", "
25 | \n", "41469586 | \n", "17785216 | \n", "17504 | \n", "2 | \n", "
26 | \n", "41469519 | \n", "17785144 | \n", "17499 | \n", "8 | \n", "
27 | \n", "41469456 | \n", "17785080 | \n", "17513 | \n", "3 | \n", "
28 | \n", "41469388 | \n", "17785008 | \n", "17507 | \n", "2 | \n", "
29 | \n", "41469323 | \n", "17784941 | \n", "17514 | \n", "1 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
850032 | \n", "41499857 | \n", "17755556 | \n", "18138 | \n", "2 | \n", "
850033 | \n", "41499914 | \n", "17755615 | \n", "18138 | \n", "2 | \n", "
850034 | \n", "41499972 | \n", "17755674 | \n", "18135 | \n", "2 | \n", "
850035 | \n", "41500000 | \n", "17755624 | \n", "18137 | \n", "2 | \n", "
850036 | \n", "41499942 | \n", "17755565 | \n", "18137 | \n", "2 | \n", "
850037 | \n", "41499884 | \n", "17755507 | \n", "18141 | \n", "2 | \n", "
850038 | \n", "41499827 | \n", "17755447 | \n", "18140 | \n", "2 | \n", "
850039 | \n", "41499769 | \n", "17755388 | \n", "18142 | \n", "2 | \n", "
850040 | \n", "41499710 | \n", "17755328 | \n", "18149 | \n", "3 | \n", "
850041 | \n", "41499653 | \n", "17755270 | \n", "18146 | \n", "2 | \n", "
850042 | \n", "41499595 | \n", "17755211 | \n", "18150 | \n", "1 | \n", "
850043 | \n", "41499538 | \n", "17755152 | \n", "18148 | \n", "2 | \n", "
850044 | \n", "41499480 | \n", "17755093 | \n", "18150 | \n", "2 | \n", "
850045 | \n", "41499423 | \n", "17755035 | \n", "18152 | \n", "2 | \n", "
850046 | \n", "41499632 | \n", "17755061 | \n", "18150 | \n", "2 | \n", "
850047 | \n", "41499689 | \n", "17755119 | \n", "18153 | \n", "8 | \n", "
850048 | \n", "41499748 | \n", "17755178 | \n", "18148 | \n", "2 | \n", "
850049 | \n", "41499805 | \n", "17755237 | \n", "18145 | \n", "2 | \n", "
850050 | \n", "41499862 | \n", "17755295 | \n", "18147 | \n", "1 | \n", "
850051 | \n", "41499921 | \n", "17755355 | \n", "18141 | \n", "2 | \n", "
850052 | \n", "41499980 | \n", "17755413 | \n", "18139 | \n", "2 | \n", "
850053 | \n", "41499992 | \n", "17755346 | \n", "18140 | \n", "2 | \n", "
850054 | \n", "41499934 | \n", "17755287 | \n", "18144 | \n", "2 | \n", "
850055 | \n", "41499875 | \n", "17755228 | \n", "18148 | \n", "3 | \n", "
850056 | \n", "41499819 | \n", "17755170 | \n", "18142 | \n", "8 | \n", "
850057 | \n", "41499761 | \n", "17755110 | \n", "18149 | \n", "2 | \n", "
850058 | \n", "41499702 | \n", "17755051 | \n", "18154 | \n", "8 | \n", "
850059 | \n", "41499906 | \n", "17755054 | \n", "18151 | \n", "8 | \n", "
850060 | \n", "41499965 | \n", "17755114 | \n", "18144 | \n", "8 | \n", "
850061 | \n", "41499987 | \n", "17755065 | \n", "18147 | \n", "2 | \n", "
850062 rows × 4 columns
\n", "