pygmt.Figure.grdview

Figure.grdview(grid, *, region=None, projection=None, zscale=None, zsize=None, frame=None, cmap=None, drapegrid=None, plane=None, surftype=None, contourpen=None, meshpen=None, facadepen=None, shading=None, verbose=None, panel=None, coltypes=None, interpolation=None, perspective=None, transparency=None, **kwargs)

Create 3-D perspective image or surface mesh from a grid.

Reads a 2-D grid file and produces a 3-D perspective plot by drawing a mesh, painting a colored/gray-shaded surface made up of polygons, or by scanline conversion of these polygons to a raster image. Options include draping a data set on top of a surface, plotting of contours on top of the surface, and apply artificial illumination based on intensities provided in a separate grid file.

Full option list at https://docs.generic-mapping-tools.org/6.5/grdview.html

Aliases:

B = frame C = cmap G = drapegrid I = shading J = projection JZ = zsize Jz = zscale

N = plane Q = surftype R = region V = verbose Wc = contourpen Wf = facadepen

Wm = meshpen c = panel f = coltypes n = interpolation p = perspective t = transparency

Parameters:

• grid (str or xarray.DataArray) –

  Name of the input grid file or the grid loaded as a xarray.DataArray object.

  For reading a specific grid file format or applying basic data operations, see https://docs.generic-mapping-tools.org/6.5/gmt.html#grd-inout-full for the available modifiers.

• region (str or list) – xmin/xmax/ymin/ymax[+r][+uunit]. Specify the region of interest. When used with perspective, optionally append /zmin/zmax to indicate the range to use for the 3-D axes [Default is the region in the input grid].

• projection (str) – projcode[projparams/]width|scale. Select map projection.

• zscale/zsize (float or str) – Set z-axis scaling or z-axis size.

• frame (bool, str, or list) – Set map boundary frame and axes attributes.

• cmap (str) – The name of the color palette table to use.

• drapegrid (str or xarray.DataArray) – The file name or a DataArray of the image grid to be draped on top of the relief provided by grid [Default determines colors from grid]. Note that zscale and plane always refer to the grid. The drapegrid only provides the information pertaining to colors, which (if drapegrid is a grid) will be looked-up via the CPT (see cmap).

• plane (float or str) – level[+gfill]. Draw a plane at this z-level. If the optional color is provided via the +g modifier, and the projection is not oblique, the frontal facade between the plane and the data perimeter is colored.

• surftype (str) –

  Specify cover type of the grid. Select one of following settings:

  • m - mesh plot [Default].

  • mx or my - waterfall plots (row or column profiles).

  • s - surface plot, and optionally append m to have mesh lines drawn on top of the surface.

  • i - image plot.

  • c - Same as i but will make nodes with z = NaN transparent.

  For any of these choices, you may force a monochrome image by appending the modifier +m.

• contourpen (str) – Draw contour lines on top of surface or mesh (not image). Append pen attributes used for the contours.

• meshpen (str) – Set the pen attributes used for the mesh. You must also select surftype of m or sm for meshlines to be drawn.

• facadepen (str) – Set the pen attributes used for the facade. You must also select plane for the facade outline to be drawn.

• shading (str) – Provide the name of a grid file with intensities in the (-1,+1) range, or a constant intensity to apply everywhere (affects the ambient light). Alternatively, derive an intensity grid from the input data grid reliefgrid via a call to pygmt.grdgradient; append +aazimuth, +nargs, and +mambient to specify azimuth, intensity, and ambient arguments for that function, or just give +d to select the default arguments [Default is +a-45+nt1+m0].

• verbose (bool or str) –

  Select verbosity level [Default is w], which modulates the messages written to stderr. Choose among 7 levels of verbosity:

  • q - Quiet, not even fatal error messages are produced

  • e - Error messages only

  • w - Warnings [Default]

  • t - Timings (report runtimes for time-intensive algorithms)

  • i - Informational messages (same as verbose=True)

  • c - Compatibility warnings

  • d - Debugging messages

• panel (bool, int, or list) – [row,col|index]. Select a specific subplot panel. Only allowed when in subplot mode. Use panel=True to advance to the next panel in the selected order. Instead of row,col you may also give a scalar value index which depends on the order you set via autolabel when the subplot was defined. Note: row, col, and index all start at 0.

• coltypes (str) – [i|o]colinfo. Specify data types of input and/or output columns (time or geographical data). Full documentation is at https://docs.generic-mapping-tools.org/6.5/gmt.html#f-full.

• interpolation (str) –

  [b|c|l|n][+a][+bBC][+c][+tthreshold]. Select interpolation mode for grids. You can select the type of spline used:

  • b for B-spline

  • c for bicubic [Default]

  • l for bilinear

  • n for nearest-neighbor

• perspective (list or str) – [x|y|z]azim[/elev[/zlevel]][+wlon0/lat0[/z0]][+vx0/y0]. Select perspective view and set the azimuth and elevation angle of the viewpoint [Default is [180, 90]]. Full documentation is at https://docs.generic-mapping-tools.org/6.5/gmt.html#perspective-full.

• transparency (float) – Set transparency level, in [0-100] percent range [Default is 0, i.e., opaque]. Only visible when PDF or raster format output is selected. Only the PNG format selection adds a transparency layer in the image (for further processing).

Example

>>> import pygmt
>>> # load the 30 arc-minutes grid with "gridline" registration
>>> # in a specified region
>>> grid = pygmt.datasets.load_earth_relief(
...     resolution="30m",
...     region=[-92.5, -82.5, -3, 7],
...     registration="gridline",
... )
>>> # create a new figure instance with pygmt.Figure()
>>> fig = pygmt.Figure()
>>> # create the contour plot
>>> fig.grdview(
...     # pass in the grid downloaded above
...     grid=grid,
...     # set the perspective to an azimuth of 130° and an elevation of 30°
...     perspective=[130, 30],
...     # add a frame to the x- and y-axes
...     # specify annotations on the south and east borders of the plot
...     frame=["xa", "ya", "wSnE"],
...     # set the projection of the 2-D map to Mercator with a 10 cm width
...     projection="M10c",
...     # set the vertical scale (z-axis) to 2 cm
...     zsize="2c",
...     # set "surface plot" to color the surface via a CPT
...     surftype="s",
...     # specify CPT to "geo"
...     cmap="geo",
... )
>>> # show the plot
>>> fig.show()

Examples using pygmt.Figure.grdview

Creating a 3-D perspective image

Creating a 3-D perspective image

Draping a dataset on top of a topographic surface

Draping a dataset on top of a topographic surface

Plotting a surface

Plotting a surface