Adding lidar data to a mosaic dataset

This topic applies to ArcGIS for Desktop Standard and ArcGIS for Desktop Advanced only.

Lidar data, stored as either LAS files, LAS datasets, or terrains, can be added to a mosaic dataset directly, without converting them to a raster before adding them. The output of the mosaic dataset will be a raster.

Why add lidar data to a mosaic dataset

If you're creating a mosaic dataset with your elevation data, you can add the lidar data. This data can be combined with other elevation data, such as SRTM or NED raster data sources.

Lidar data can be managed within a mosaic dataset by adding the LAS files, LAS datasets, or terrain dataset. For example, you can add all your LAS files, from all your different acquisitions, and use the mosaic dataset as the single source for finding the data as it's needed.

Using the mosaic dataset, you can also serve the lidar data as an image service providing access to the interpolated surface, or even to allow users to download the source files.

Learn about managing elevation data in a mosaic dataset

How to add lidar data to a mosaic dataset

Lidar data can be added to mosaic datasets using their associated raster types.

When these types of data are added to a mosaic dataset they are converted into raster data. Each raster type adds a unique function which can be edited but not added to or removed from the item in the mosaic dataset.

The output properties are unique to these raster types. Since the inputs involve some sort of interpolation from points, this can be quite computationally intensive and therefore slow to display. There is an option to create caches at the base pixel size for the inputs to improve performance. Without the cache you may have to wait several minutes for some surfaces to display.

TipTip:

If you do not need to add constraints or make other modifications to the lidar data, then add the LAS files directly to the mosaic dataset using the LAS raster type. Otherwise, create a LAS dataset from the LAS files, add the constraints, then add the LAS dataset to the mosaic dataset.

A LAS dataset can be used to edit the classification used by the LAS files. If the LAS files are modified, they will be updated in the mosaic dataset directly.

LAS file additions to the attribute table

When adding LAS files directly to the mosaic dataset, there are some additional fields that get added to the mosaic dataset's attribute table to support the LAS data.

LAS specific fields

Description

Version

The version of the LAS file. LAS is an industry format created and maintained by the American Society for Photogrammetry and Remote Sensing (ASPRS).

Point Count

Total number of points in the LAS file.

Point Spacing

Average spacing between points in the LAS file in the units of the LAS dataset.

ZMin

Minimum point value in the LAS file.

ZMax

Maximum point value in the LAS file.

Steps to add lidar data to a mosaic dataset

Steps:
  1. First create a mosaic dataset with the Create Mosaic Dataset tool.
  2. In the Catalog window or ArcCatalog, right-click the mosaic dataset and click Add Rasters.
  3. Click the Raster Type drop-down list and click one of the following:
    • LAS
    • LAS Dataset
    • Terrain

    There are properties that need to be set when using either of these raster types, therefore, you cannot run this tool without setting those properties.

  4. Click the Edit Raster Type Properties button Properties.
  5. Click one of the following tabs, depending on the raster type you chose:
    • LAS
    • Las Dataset
    • Terrain

    You will likely want to modify many of these parameters. To learn about each of the parameters on this tab, see either:

  6. You can accept all the defaults, but you must enter a Pixel size.

    The pixel size is the minimum pixel size that will be generated to create the raster. It is better to go with a pixel size that is several times larger than the average point spacing but small enough to identify gaps or voids. Generally, if the pixel size is three or four times greater than the point spacing, the voids in the data should be filled (unless, for example, the voids are due to water). For example, if your data is sampled at 1 meter and your pixel size is 4, you can expect, on average, to get 16 points in a pixel.

  7. If the pixel size is 3-4 times larger than the average point distance, you can safely use binning. If the cell size is smaller than that, you can try binning with void filling turned off. If the resulting raster mainly contains voids and only a few single data cells, binning generally does not produce a meaningful elevation raster. You need to either increase the pixel size or switch to triangulation. If the resulting raster shows enough content with some salt and pepper voids, and maybe a few larger voids, you can use binning with void filling turned on. Click the Void filling drop-down arrow and select either Simple or Plane Fitting/IDW.
  8. If you're adding LAS files, you may want to check Treat each folder as a dataset.
    This will add all the files in a folder as a single item in the mosaic dataset, which is more efficient for the mosaic dataset. You may choose to check this option if all the LAS files in a folder belong together and have the same spatial reference. For example, they may represent a single data collection (project) that are just stored as tiles.

    If the LAS files are unrelated to one another, not part of the same project collection, or don't have the same spatial reference, then don't check this option.

  9. You may want to modify the Cache folder and Number of cached surfaces as they will take some space on disk and need to be moved with the mosaic dataset if it's moved.

    A cache is generated to improve the display time for this data. By default, the cache is stored in a folder next to where the mosaic dataset resides. If you will be modifying the properties of the function to show different views from the same mosaic dataset, then you will want to plan to generate a number of cached surfaces, because one will be created for each modification.

    The cache will be created the first time the dataset is drawn; therefore, it will take some time before the surface can be viewed. Once it's cached the surface will display quickly. Otherwise, you can build the cache by building the overviews for the mosaic dataset, or using the Synchronize Mosaic Dataset tool with the Build Item Cache option checked.

  10. Once you have set all the properties on this tab, you can close the dialog box by clicking OK.
  11. If you won't be adding more data to the mosaic dataset, you may want to check Update Overviews.

    Overviews take time to create and require space for storage, but they will ensure that the statistics area is calculated, the cache is generated, and the mosaic dataset can be viewed at all scales.

  12. If the coordinate system for the data isn't the same as the mosaic dataset, expand the Advanced Options, click the Coordinate System for Input Data button Properties and select a coordinate system.

    You will likely need to specify the coordinate system when adding LAS files.

    The coordinate system specified here is used for all the input data. If some of your data is in a different coordinate system, then add it separately, using this same tool.

  13. Click OK.

Now you can view the mosaic dataset in ArcMap, ArcGlobe, or ArcCatalog. Remember, the cache is created the first time the dataset is drawn; therefore, it will take some time before the surface can be viewed.

If you selected the option to build overviews, then the cache and statistics will be generated. If you did not, then you will have to wait as the cache is generated and since the statistics are likely the entire bit depth, not the range of input values, you should open the Image Analysis window, select the layer, and check DRA to stretch the image based on the values within the extent. If you will not be building overviews, run the Calculate Statistics tool (with a skip factor to save time) to calculate the statistics for the mosaic dataset.

If you will be sharing the mosaic dataset as an image service, it is best to publish the image service, then generate the cache for the mosaic dataset once it's on the server (since the cache is not copied when the mosaic dataset is shared as an image service).

Void filling limitations

If you use the Plane Fitting/IDW void filling method without specifying a Maximum width, then the interpolation will fill all gaps. This includes the gaps that may occur between the points and the extent of the file if the points do not fill the extent of the file. In the example below, the elevation data is shown as a grayscale image and has been rendered as a hillshade with some coloring to help you see the blending. The green lines represent the extent of the LAS file and the red boxes outline areas of blending.

LAS files with interpolation to the file's edges

If you overlay the points from each LAS file you can see that the blue points in the LAS file on the right do not fill the extent of the file, and you can see that the blue and the red LAS files have overlapping extents but the points do not overlap.

Lidar points overlaid on interpolated surface

There are a number of different decisions that can be made to avoid this issue. The first option is to use Simple void filling, or define a value when using Plane Fitting/IDW. This way, the interpolation won't be forced to fill in the gaps along the edge of the file's extent. The examples below use the Simple method.

Void filling using an appropriate option

Another option is to check the Treat each folder as a dataset check box when adding the LAS files to the mosaic dataset. This will add the files as one dataset; therefore, in this case only one footprint would represent these two LAS files. When doing this, there won't be any interpolation along the overlapping extents since the extents are merged. However, void filling will still occur along the edges of the entire dataset.

TipTip:

Adding as a dataset is recommended. This won't impact the download, as a user can still download the individual LAS files from the mosaic dataset (or image service).

If using the Plane Fitting/IDW void filling method without specifying a maximum width or using a significantly large value to fill large voids is necessary, then you can run the Build Footprints tool and recalculate the footprints. First, set the void filling to Simple or None. Then, when rebuilding the footprints, use the By Radiometry computation method and likely take the defaults. This will adjust the footprint to fit the extent of the points. Then you can change the void filling method.

To change the void filling method, if there is only one item in the mosaic dataset's attribute table, then click the Raster field for that item, access the LAS To Raster function, and modify the parameter. If you have many items in the attribute table, you can modify them all at once; see Editing function chains in a mosaic dataset.

Related Topics

3/5/2014