0
$\begingroup$

In an effort to understand how to use GeologicalPeriodData and extending an example provided in the instructions, I created four GeoGraphics plots of bedrock of varying ages. Each sizes nicely, as can be seen by the following code:

bedrockPaleocene = GeologicalPeriodData["PaleoceneEpoch", "BedrockPolygon"];

paleocene = 
 GeoGraphics[{Line @@@ 
 EntityValue[
  EntityClass["AdministrativeDivision", 
 "ContiguousUSStates"], 
  EntityProperty["AdministrativeDivision", "Polygon"]], 
 GeoStyling[Red], bedrockPaleocene}, 
 GeoCenter -> 
 Entity["City", {"NewOrleans", "Louisiana", 
 "UnitedStates"}], 
  GeoRange -> Quantity[1500, "Miles"], 
 GeoBackground -> "ReliefMap"];

bedrockEocene = 
 GeologicalPeriodData["EoceneEpoch", "BedrockPolygon"];
eocene = GeoGraphics[{Line @@@ 
 EntityValue[
  EntityClass["AdministrativeDivision", 
"ContiguousUSStates"], 
  EntityProperty["AdministrativeDivision", "Polygon"]], 
  GeoStyling[Yellow], bedrockEocene}, 
 GeoCenter -> 
 Entity["City", {"NewOrleans", "Louisiana", 
 "UnitedStates"}], 
 GeoRange -> Quantity[1500, "Miles"], 
 GeoBackground -> "ReliefMap"];

bedrockMiocene = 
 GeologicalPeriodData["MioceneEpoch", "BedrockPolygon"];
miocene = 
 GeoGraphics[{Line @@@ 
 EntityValue[
  EntityClass["AdministrativeDivision", 
 "ContiguousUSStates"], 
  EntityProperty["AdministrativeDivision", "Polygon"]], 
  GeoStyling[Green], bedrockMiocene}, 
 GeoCenter -> 
 Entity["City", {"NewOrleans", "Louisiana", 
"UnitedStates"}], 
 GeoRange -> Quantity[1500, "Miles"], 
 GeoBackground -> "ReliefMap"];

bedrockPliocene = 
 GeologicalPeriodData["PlioceneEpoch", "BedrockPolygon"];
pliocene = 
 GeoGraphics[{Line @@@ 
 EntityValue[
  EntityClass["AdministrativeDivision", 
 "ContiguousUSStates"], 
  EntityProperty["AdministrativeDivision", "Polygon"]], 
 GeoStyling[Blue], bedrockPliocene}, 
 GeoCenter -> 
 Entity["City", {"NewOrleans", "Louisiana", 
 "UnitedStates"}], 
 GeoRange -> Quantity[1500, "Miles"], GeoBackground -> 
"ReliefMap"];

 

However, when one uses Show to combine them, the GeoRange is altered:

Show[paleocene, eocene, miocene, pliocene]

Is there a way to constrict subsequent graphics plots to retain the same GeoRange values as the first plot in the combined plot?

$\endgroup$
3
  • $\begingroup$ Can you include an explantation of how the georange is altered and what you expected instead? $\endgroup$
    – MarcoB
    Commented Jul 31, 2023 at 0:57
  • 1
    $\begingroup$ Legended[Show[paleocene, eocene, miocene, pliocene, GeoRange -> Quantity[1500, "Miles"]], SwatchLegend[{Red, Yellow, Green, Blue}, {"Paleocene", "Eocene", "Miocene", "Pliocene"}, LegendLabel -> "bedrock polygon"]] $\endgroup$
    – Bob Hanlon
    Commented Jul 31, 2023 at 1:23
  • $\begingroup$ Apologies for not getting back to this sooner. Bob Hanlon's suggestion worked giving me what I had expected. Thank you. Moving your comment to the answer section. $\endgroup$ Commented Dec 5, 2023 at 4:53

1 Answer 1

0
$\begingroup$

Following Bob Hanlon's comment, this issue is resolved and nice legends added, which is a bonus. Any points go to him.

Legended[Show[paleocene, eocene, miocene, pliocene,
  GeoRange -> Quantity[1500, "Miles"]],
  SwatchLegend[{Red, Yellow, Green, Blue},
  {"Paleocene", "Eocene", "Miocene", "Pliocene"},
  LegendLabel -> "bedrock polygon"]]

– Bob Hanlon Jul 31 at 1:23

$\endgroup$
1
  • $\begingroup$ To my chagrin, it seems that data for other parts of the world besides the USA appear to be lacking, making this less useful than one might hope. Mathematica-oriented geologists are needed to come to the rescue. $\endgroup$ Commented Dec 10, 2023 at 0:16

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.