Added example of setting projection using EPSG, and modified rgb_example to incorporate this

examples/gallery/images/rgb_image.py

@@ -34,10 +34,19 @@
 # Plot the RGB imagery:
 fig = pygmt.Figure()
 with pygmt.config(FONT_TITLE="Times-Roman"):  # Set title font to Times-Roman
-    fig.grdimage(
-        grid=image,
-        # Use a map scale where 1 cm on the map equals 1 km on the ground
-        projection="x1:100000",
-        frame=[r"WSne+tL@!a¯hain@!a¯, Hawai`i on 9 Aug 2023", "af"],
+    (
+        fig.grdimage(
+            grid=image,
+            # Use a map scale where 1 cm on the map equals 1 km on the ground
+            projection="x1:100000",
+            frame=[r"WSne+tL@!a¯hain@!a¯, Hawai`i on 9 Aug 2023", "af"],
+        ),
     )
+    fig.coast(
+        region=[738000, 755000, 2300000, 2318000],
+        projection="EPSG:32604",
+        frame=False,
+        shorelines="1p",
+    )
+
 fig.show()

examples/projections/misc/misc_epsg.py

@@ -0,0 +1,32 @@
+r"""
+Using EPSG codes to specify projection
+===============================
+
+The EPSG registry is a public registry of geodetic datums, spatial reference systems,
+Earth ellipsoids, coordinate transformations and related units of measurement.
+
+These are used frequently in GIS software and data. If data are plotted using a linear
+projection, additional features (e.g. coast) can be made congruent by specifying EPSG.
+See Gallery Example "RGB Image"
+
+Some common EPSG codes are:
+
+EPSG:4326 - WGS 84 - latitude/longitude coordinate system based on the Earth's centre
+of mass, used by the Global Positioning System among others.
+
+EPSG:3857 - Web Mercator - used for display by many web-based mapping tools, including
+Google Maps and OpenStreetMap.
+
+
+**[EPSG:XXXX]**
+
+"""
+
+# %%
+import pygmt
+
+fig = pygmt.Figure()
+fig.coast(
+    region=[-180, 180, -80, 80], projection="EPSG:3857", frame=True, shorelines="1p"
+)
+fig.show()