How do I get the raster cell(s) with the maximum value?
A very simple question, hopefully with a simple answer.
I have a single-band raster file with flood depths, and want to put a point or label at the deepest point ("X marks the spot"). I can get the deepest value using the RasterBandMinMaxExtractor, but how do I get the location of the deepest value.
The solution I've come up with so far is using a RasterCellCoercer converting the entire raster into a point feature class, and then extract the deepest point using a Tester transformer. But this feels to be a bit of a roundabout way, having to vectorise the raster first.
Hopefully/Surely there is a shortcut....?
I have two thoughts, no idea which will end up being more efficient.
1. RasterExpressionEvaluator, deepest value = 1, everything else = NoData, followed by a RasterCellCoercer with extract NoData set to No.
2. Convert the raster to a point cloud (PointCloudCombiner) followed by a PointCloudFilter, and optionally a PointCloudCoercer.
The point cloud tools are quite efficient and there are often gains in converting a raster to a point cloud and processing.
- My original solution = duration 16 min 39.8s, peak memory 881424 kB
- Alternative solution 1 = duration 2 min 28.1s, peak memory 122756 kB
- Alternative solution 2 = duration 1 min 20.5s, peak memory 55836 kB
I think your comment about converting raster to a point cloud may turn out to be a valuable one!
As a follow-up on @jdh I would add:
If you create a point cloud you could use the PointcloudStatisticsCalcultator transformer to get the values.
