gdalread

function gdalshade(fname; kwargs…) d = KW(kwargs) band = ((val = find_in_dict(d, [:band], false)[1]) !== nothing) ? string(val) : “1” cmap = find_in_dict(d, CPTaliases)[1] # The color cannot be passed to second call to gdaldem

A = gdaldem(fname, "color-relief", ["-b", band], C=cmap)
B = gdaldem(fname, "hillshade"; d...)
blendimg!(A, B)

end

—————————————————————————————————

““” gdalread(fname::AbstractString, opts=String[]; gdataset=false, kwargs…)

Read a raster or a vector file from a disk file and return the result either as a GMT type (the default) or a GDAL dataset.

  • fname: Input data. It can be a file name, a GMTgrid or GMTimage object or a GDAL dataset
  • opts: List of options. The accepted options are the ones of the gdal_translate utility. This list can be in the form of a vector of strings, or joined in a simgle string.
  • gdataset: If set to true forces the return of a GDAL dataset instead of a GMT type.
  • kwargs: This options accept the GMT region (-R) and increment (-I)

Returns

A GMT grid/image or a GDAL dataset