3D scanning like a pro

May 8, 2017

Saulius Lukse

We love it when our readers get in touch with us to share their stories. This article was contributed to DIYP by a member of our community. If you would like to contribute an article, please contact us here.

3D scanning like a pro

May 8, 2017

Saulius Lukse

We love it when our readers get in touch with us to share their stories. This article was contributed to DIYP by a member of our community. If you would like to contribute an article, please contact us here.

Join the Discussion

Share on:

Almost any sensor yields more interesting results if mounted on a moving platform. Remember scanning thermometer? It’s time to mount TOF LIDAR (Time Of Flight / Light Imaging, Detection, And Ranging) on two precision rotary stages arranged for pan and tilt operation. Rig provides real-time position data along with distance to an obstacle. Using simple math we can calculate the position in Cartesian coordinate system. Data is collected point-by- point to reconstruct 3D object model.

Used parts

I have been looking for LIDAR for a while and decided to use Garmin as an affordable distance sensor, so bought one for a spin.

Results

Garmin LIDAR outputs data at up to 500 measurements a second which is way faster than the thermal sensor, yet it takes about 15 minutes to scan full sphere yielding true scale 1-2M point cloud. No stitching, no lengthy post processing! And the result is spectacular.

Animated GIF for better depth perception. Calculated point cloud normals and colorized according to their vector direction. Click for better quality on Youtube.

This is another fast scan with higher rotation speed. Result is 300K point cloud.

4.7M point cloud room view from above with colors removed. From this picture, it’s possible to see that distance measurement could be more precise and laser beam angle could be narrower. Density is pretty even which indicates nice rotary stage performance and correct calculation.

Depth map colorized as heat map colors and fixed black pixels. This imaging is useful for fast quality inspection.

Unfortunately, this particular LIDAR is not quite suitable for outdoor mapping. The laser is too weak even in dusk and measurement rate drops, therefore scanning speed has to be reduced and it would take hours to cover a longer range. Click for better quality on Youtube.

Laser barely reaches 4’th floor, therefore making it not suitable for long range outdoor operation.

Software

I have written set of Python scripts to make scans like this. As release is still pending, will briefly describe what they do:

Scan

  1. Establish communication to sensor and motion controller
  2. Initialize motion platform (homing, zeroing)
  3. Start motion at specified speed and tilt incline angle. Rotation count is calculated automatically. Platform rotates at constant angular speed gradually changing tilt angle
  4. Platform position data is being recorded into file along with timestamp
  5. Sensor data is being recorded into separate file along with timestamp
  6. After scan is complete platform moves home and computer beeps to indicate end of scan

Fill dead pixels

  1. Load bitmap picture (generated as intermediate result by aggregate script)
  2. Find black pixels
  3. For each black pixel find neighbouring pixels, average their color and fill black pixel with new value
  4. Save file as new bitmap picture

Finding of black pixels can be done by brute force but it will take ages. Speed up can be done by using numpy functionality:

pixels = np.asarray(img_gray)
coordList = np.where(pixels != 0)
coordinates = zip(coordList[0], coordList[1])

Aggregate

  1. Load previously saved files (position data and sensor values)
  2. For each sensor value interpolate motor position from two nearest points (aligned by timestamp)
  3. Convert depth data to color and coordinates into 2D picture for fast preview (or final picture in case of thermometer)
  4. Convert angle and measured depth to Cartesian coordinate system and save each pixel to *.xyz file.

Pan/tilt/distance conversion to x,y,z is as simple as these formulas:

x = r * sin(th) * cos(tv)
y = r * sin(th) * sin(tv)
z = r * cos(th)

That’s it. Now point cloud file can be opened with preferred 3D manipulation program like MeshLab

Improvements and notes

The main part to be improved is LIDAR. Faster acquisition rate, smaller beam expansion, better precision and greater power is required.

Motion platform operation is quiet, smooth, fast and precise.Though with faster laser fast axis might be limiting factor.

Draft software lacks real-time preview of what’s happening, but that’s easy to update.

One slip-ring is damaged (or bad by design) and sometimes I2C communication is interrupted – added watchdog protection to restore LIDAR operation.

One more major update could be getting rid of the laptop, instead using single board computer with tablet or smartphone connected over WiFi.

About the Author

Saulius Lukšė is an engineer with a demonstrated history of working in the wireless industry. He is a founder of Kurokesu, an e-shop where he sells his designs and solutions like touch probe, Raspberry Pi case, C920 webcam rework kit, high sensitivity USB camera and others. You can check out his products on his website and look him up on Facebook and Twitter. This article was also published here and shared with permission.

Filed Under:

Tagged With:

Find this interesting? Share it with your friends!

DIPY Icon

We love it when our readers get in touch with us to share their stories. This article was contributed to DIYP by a member of our community. If you would like to contribute an article, please contact us here.

Join the Discussion

DIYP Comment Policy
Be nice, be on-topic, no personal information or flames.

Leave a Reply

Your email address will not be published. Required fields are marked *