Friday, March 5, 2010

Create LiDAR intensity images using GeoMedia Grid

I tried to create LiDAR intensity images using GeoMedia Grid. It can be done interactively but I found the workflow to be quite tedious especially if you have a large number of tiles to handle. If you have a large number of intensity images to create, it would be more practical to automate the process either by developing a custom GeoMedia Grid program or use alternative software. To create the intensity images, here are the steps to do the job interactively.

Create a Study Area
  1. Start up GeoMedia and open up a geo-workspace and connect to a read/write warehouse with the appropriate coordinate system.
  2. Select Grid | Study Area | Define New.

  3. In the Coordinate Precision Key-in tool bar, type in the precise coordinate of one corner of the tile. Press RETURN.

  4. Type in the precise coordinate of the diagonally opposite corner of the tile. Press RETURN.

    The Define New Study Area dialog box appears.
  5. In the Study Area name field, type in the desired name. In the Cell resolution field, type in the output resolution e.g. 1 and set the units e.g. meters. Press OK.

    The study area is created.
Import a LAS tile file
  1. Follow the steps in my previous post to import a LAS file into the study area created previously.
Triangulate the tile by intensity values

  1. Select Grid | Study Area | List. Choose the previously created study area and imported LAS tile layer.

    The selected study area and layer are active.
  2. Select Grid | Interpolation | Triangulation.


    The Triangulation dialog box appears.
  3. Click Add.

    The Add Source Feature dialog box appears.
  4. In the Feature field, select the imported LAS tile layer in the current study area. Toggle on the Attribute Z-value. In the Attribute drop down list, choose Intensity. In the Units drop down list, choose None.

    The Add Source Feature dialog box may look like this at this point.
  5. Click OK.

    The Triangulation dialog box is updated.
  6. In the Output units field, select None. In the Result layer Name field, type in a meaningful name. Toggle off Place result in map window. Click OK.

    The result layer is created and added to the active study area.
Define the intensity gray scale range
By default, GeoMedia Grid will assign a grey scale color range to the result layer. Normally, the void or null values will be assigned as white. I prefer to color the voids as black. If you have only one tile, you might as well use the default. But if you have a number of tiles, then in order to ensure that all the tiles look the same, you have to create a grey scale xml file and reuse that xml file for all the tiles.
  1. In the Study Area List dialog box, right click on a layer.

    A pop up menu appears.

  2. Choose View Legend.

    The View Legend dialog box appears.

  3. Double click the box on the left of the VOID label as shown above.

    The color picker appears.

  4. Choose black.

    The View Legend(s) dialog box is updated.

  5. Click Apply.

    The color range is applied to the layer.

    Note: Click Save As to save the color range into an xml file. To reuse the xml file, simply click Load and choose the previously saved xml file.
  6. Click Close.

Create the intensity image as a GeoTiff file
  1. In the Study Area List dialog box, click on the newly created layer.
  2. Select Grid | Layer | Export to File.



    The Export to File dialog box appears.
  3. Click OK.

    The Save As dialog box appears.
  4. In the File name field, type in a name. In the Save as type drop down list, choose TIFF. Click Save.

    The Export to GeoTiff dialog box appears.
  5. Choose RGB or Grayscale. Click OK.

    The intensity GeoTiff image file is created.

No comments: