A talk I gave at PyCon US 2019
Terrain,Andrew Godwin@andrewgodwinArt,Python& LiDAR
View Slide
Hi, I’mAndrew Godwin• Django contributor (Migrations/Channels)• Principle Engineer at• I always need more lasers
1.Laser-cut profiles
DEM: Digital Elevation Model
10 9 5 3 29 7 4 3 28 7 3 2 46 8 5 3 54 9 7 4 63 8 8 6 52 5 4 3 310 9 5 3 29 7 4 3 28 7 3 2 46 8 5 3 54 9 7 4 63 8 8 6 52 5 4 3 3
Loads the DEM as a CSVPicks one in N rowsDraws a contour using svgwrite
2.3D-printed cities
(1, 3, 4)(6, 5, 11)(3, 7, -2)
POINT CLOUD → DIGITAL ELEVATION MAP (DEM)python-pcllastools
Top surfaceFully sealed 3D model(with the tile base)
Loads the DEMClips height (top and bottom)Smooths rough featuresWrites out an STL file
How do you write STL?struct.pack(b"normal[0],normal[1],normal[2],x1, y1, z1,x2, y2, z2,x3, y3, z3,0,)
Should I have used NumPy? Yes.Did I use NumPy? Err, no.
github.com/andrewgodwin/lidartile
3.CNC-milled National Parks
Take the US National DEMGet the outline of the National ParkUse QGIS to cut out a park-only DEMhttps://viewer.nationalmap.govhttps://catalog.data.gov/dataset/national-park-boundariesf0a4cToolbox > GDAL > Clip By Extent
Irregular shapes!
MAP PROJECTIONS
Things I Won't Work WithUnicodeNamesTimezonesCurrenciesNetworksAddressesDioxygen Difluoride
Things I Won't Work WithUnicodeNamesTimezonesCurrenciesNetworksAddressesDioxygen DifluorideMap Projections
4.The Future
More US National ParksI do each one as I visit it. There are… 59.
Easier Milling8 hours per piece. Really.
Better STL optimisationMillions of polygons isn't great.
Personal LiDARThanks, autonomous vehicles!
github.com/andrewgodwin/lidartilegithub.com/andrewgodwin/gis_tools
Thanks.Andrew Godwin@andrewgodwin aeracode.org