I have two DSM rasters, one is 10km radius photogrammetric data, the other is 2km radius LiDAR data, both centred around the same point.

I need to combine these together, essentially inserting the 2km LiDAR raster into the other and removing that portion of the photogrammetric data, the idea being that the LiDAR is much more detailed in the area closest to the centre and the site that I'm working on.

I've resampled the photogrammetric raster to match the LiDAR so the cell sizes are the same.

What is the best way to go about doing this? I think it's using 'Mosaic to New Raster' but can't get my head around whether to use FIRST or LAST

1 Answer 1


In ArcMap application you could use some of raster merging techniques. I recommend NAGI (No Alteration of Grayscale or Intensity) pansharpening fusion method which is quite useful when combining multispectral (colorful) and one spectral (grayscale) images. Here are detailed steps of it:

Define the functions for the raster datasets by following the steps below:

  1. Add the grayscale hillshade and multispectral RGB layer tint rasters to ArcMap, if they have not already been added.
  2. Open the Image Analysis window by clicking Windows on the top bar menu, then clicking Image Analysis.
  3. In the top section of the Image Analysis window, select both the hillshade and RGB rasters using the Control key and clicking on each raster’s name to highlight it (figure 5).

Overlaying rasters - Image Analysis Window

Figure 5. The Image Analysis window

  1. Click the Pan-Sharpening tool in the Processing section of the Image Analysis window. This will create a new layer, which will be listed as the top layer in the Image Analysis window.
  2. In the Image Analysis window, right-click the newly generated pan-sharpening layer and click Properties.
  3. On the Functions tab, right-click the Pansharpening Function and click Properties.
  4. On the General tab of the Raster Function Properties dialog box, change the Output Pixel Type to 8 Bit Unsigned.
  5. On the Pan Sharpen tab, change the Method to Simple Mean.
  6. Keep the rest of the defaults and click OK.
  7. Right-click Pansharpening Function, click Insert, and click Stretch Function.
  8. Change the Type to Minimum-Maximum.
  9. Check the Use Gamma option.
  10. In the Gamma section of the dialog box, change the Gamma value from 1.0 to 0.5 for each of the three bands.
  11. In the Statistics section of the dialog box, type 5 as the Min and 215 as the Max value for each of the three bands. The final function chain will look like figure 6.

Overlaying rasters - Function Chain

Figure 6. The final function chain

  1. Click OK to check your results.

Above text is a quote from the article where description of the method could be found: https://www.esri.com/arcgis-blog/products/product/imagery/an-alternative-to-overlaying-layer-tints-on-hillshades/

Another publication concerning this technique: http://www.cartogis.org/docs/proceedings/2012/Nagi_AutoCarto2012.pdf

A presentation: https://www.slideshare.net/aileenbuckley/the-nagi-fusion-method

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.