I have generated a relief map according to a method shown in Mathematica documentation, viz:
upperleft = GeoPosition[{36.843626, 13.916897}]
lowerright = GeoPosition[{35.710239, 15.37673}]
demdata =
Reverse[QuantityMagnitude[GeoElevationData[{upperleft, lowerright}]]]
ListPlot3D[demdata, AspectRatio -> (Divide @@ Dimensions[demdata]),
Boxed -> False, Mesh -> None, Axes -> False, Background -> Black,
ColorFunction -> ColorData["HypsometricTints"],
ColorFunctionScaling -> False]
I would like to add one or two particular contours at specific depths in the ocean areas, like at -150 feet and -300 feet. Is there a way to do that?
Mesh -> None
withMeshFunctions -> {#3 &}, Mesh -> {{-150, -300}}
. $\endgroup$