What is a C++ implementation or analogue of the local peak detector skimage.feature.peak_local_max?

I need to find local peaks on the depth image to generate the markers for a watershed algorithm in C++.

My Python implementation uses skimage.feature.peak_local_max to generate the markers for a watershed algorithm.

Can't find a C++ implementation or C++ analogue of the skimage.feature.peak_local_max function.

Some algorithm to find a local maximum for a matrix also could be helpful.

distance = scipy.ndimage.distance_transform_edt(image)
local_maxi = skimage.feature.peak_local_max(distance, indices=False, footprint=np.ones((3, 3)), labels=image)
markers = scipy.ndimage.label(local_maxi)[0]
labels = skimage.segmentation.watershed(-distance, markers, mask=image)

Example of the peak finder

Any ideas on C++ implementation of the local peak finder, please? Thank you!