16
$\begingroup$

I want to make a diagram which corresponds to sea level rise in the UK. I'm having trouble with relief plot.

So, corresponding to a 30 metre height increase of the sea, we would get the following diagram

sealevelrise = 30;

ReliefPlot[GeoElevationData[Interpreter["Country"]["UK"], UnitSystem -> "Metric"], DataReversed -> True, PlotRange -> {Full, Full, {sealevelrise, All}}]

Giving,

enter image description here

But I want a colour change at the new sea level instead, e.g from green to blue at the sea line. i.e. we could still see the relief in the white/"underwater". Is there any way to do this?

Clarification: I would like for the relief to be green, and for the sea to be blue and for the relief of both to still be seen. I just would like the point where the sea becomes land to have a colour change.

Side (non-important) question: Is there any way to overlay cities?

$\endgroup$
17
$\begingroup$

I'm not completely sure of what you need, but you might try this (where you can change the color schemes as you prefer):

ImageCompose[
 ColorReplace[
  ReliefPlot[
   GeoElevationData[Interpreter["Country"]["UK"], 
    UnitSystem -> "Metric"], DataReversed -> True, 
   PlotRange -> {Full, Full, {0, All}}
   , ColorFunction -> "SandyTerrain"]
  , White -> Transparent]
 ,
 ColorReplace[
  ReliefPlot[
   GeoElevationData[Interpreter["Country"]["UK"], 
    UnitSystem -> "Metric"], DataReversed -> True, 
   PlotRange -> {Full, Full, {sealevelrise, All}}], 
  White -> Transparent]
 ]

enter image description here

The orange-ish parts in the image are the parts that are now underwater.

$\endgroup$
  • $\begingroup$ @Kuba thanks for the edit, I was in a hurry and forgot to add the plot :) $\endgroup$ – Fraccalo Aug 21 '18 at 11:05

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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