add new gnuplot capabilities (surface, scatter3d, wireframe, ...)

This commit is contained in:
t-bltg 2021-07-25 15:21:23 +02:00
parent 30e728213e
commit c0e35f2358
3 changed files with 171 additions and 76 deletions

View File

@ -623,20 +623,29 @@ const _gaston_attr = merge_with_base_supported([
# :camera, # :camera,
# :contour_labels, # :contour_labels,
]) ])
const _gaston_seriestype = [:path,
# :path3d, const _gaston_seriestype = [
:path,
:path3d,
:scatter, :scatter,
:steppre, :steppre,
# :stepmid, :stepmid,
:steppost, :steppost,
# :histogram2d, :histogram2d,
# :ysticks, :xsticks, :ysticks, :xsticks,
# :contour, :contour,
:shape, :shape,
# :straightline, :straightline,
:scatter3d,
:contour3d,
:wireframe,
:heatmap,
:surface,
:mesh3d,
] ]
const _gaston_style = [:auto, const _gaston_style = [
:auto,
:solid, :solid,
:dash, :dash,
:dot, :dot,
@ -644,7 +653,8 @@ const _gaston_style = [:auto,
:dashdotdot :dashdotdot
] ]
const _gaston_marker = [:none, const _gaston_marker = [
:none,
# :auto, # :auto,
:circle, :circle,
:rect, :rect,
@ -656,7 +666,8 @@ const _gaston_marker = [:none,
:+ :+
] #vcat(_allMarkers, Shape) ] #vcat(_allMarkers, Shape)
const _gaston_scale = [:identity, const _gaston_scale = [
:identity,
# :ln, # :ln,
# :log2, # :log2,
:log10, :log10,

View File

@ -22,18 +22,11 @@ end
function _before_layout_calcs(plt::Plot{GastonBackend}) function _before_layout_calcs(plt::Plot{GastonBackend})
# Initialize all the subplots first # Initialize all the subplots first
plt.o.subplots = G.SubPlot[] plt.o.subplots = G.SubPlot[]
nr, nc = plt.o.layout = size(plt.layout)
n = 0 n, sps = gaston_get_subplots(plt, 0, plt.layout)
sps = Array{Any}(undef, nr, nc) @assert n == length(plt.subplots)
for r 1:nr, c 1:nc # NOTE: row major
l = plt.layout.grid[r, c]
sps[r, c] = get(l.attr, :blank, false) ? nothing : plt.subplots[n += 1]
end
for c 1:nc, r 1:nr # NOTE: col major plt.o.layout = gaston_init_subplots(plt, sps)
gaston_init_subplot(plt, sps[r, c])
end
# Then add the series (curves in gaston) # Then add the series (curves in gaston)
for series in plt.series_list for series in plt.series_list
@ -59,8 +52,7 @@ for (mime, term) in (
"text/plain" => "dumb", # NEED fixing TODO "text/plain" => "dumb", # NEED fixing TODO
) )
@eval function _show(io::IO, ::MIME{Symbol($mime)}, plt::Plot{GastonBackend}) @eval function _show(io::IO, ::MIME{Symbol($mime)}, plt::Plot{GastonBackend})
xsize = plt.attr[:size][1] xsize, ysize = plt.attr[:size]
ysize = plt.attr[:size][2]
termopts="""size $xsize,$ysize""" termopts="""size $xsize,$ysize"""
tmpfile = G.tempname() * "." * $term tmpfile = G.tempname() * "." * $term
@ -79,11 +71,9 @@ for (mime, term) in (
end end
end end
function _show(io::IO, mime::MIME{Symbol("image/png")}, function _show(io::IO, mime::MIME{Symbol("image/png")}, plt::Plot{GastonBackend})
plt::Plot{GastonBackend},)
scaling = plt.attr[:dpi] / GNUPLOT_DPI scaling = plt.attr[:dpi] / GNUPLOT_DPI
xsize = plt.attr[:size][1] * scaling xsize, ysize = plt.attr[:size] .* scaling
ysize = plt.attr[:size][2] * scaling
# Scale all plot elements to match Plots.jl DPI standard # Scale all plot elements to match Plots.jl DPI standard
termopts="""size $xsize,$ysize fontscale $scaling lw $scaling dl $scaling ps $scaling""" termopts="""size $xsize,$ysize fontscale $scaling lw $scaling dl $scaling ps $scaling"""
@ -111,14 +101,44 @@ end
# These functions are gaston specific # These functions are gaston specific
# -------------------------------------------- # --------------------------------------------
gaston_get_subplots(plt, n, layout, level) = begin
nr, nc = size(layout)
sps = Array{Any}(undef, nr, nc)
for r 1:nr, c 1:nc # NOTE: col major
l = layout[r, c]
if l isa GridLayout
n, sub = gaston_get_subplots(plt, n, l)
sps[r, c] = size(sub) == (1, 1) ? only(sub) : sub
else
sps[r, c] = get(l.attr, :blank, false) ? nothing : plt.subplots[n += 1]
end
end
n, sps
end
gaston_init_subplots(plt, sps, level) = begin
sz = nr, nc = size(sps)
for c 1:nc, r 1:nr # NOTE: row major
sp = sps[r, c]
if sp isa Subplot || sp === nothing
gaston_init_subplot(plt, sp)
else
gaston_init_subplots(plt, sp, level + 1)
sz = max.(sz, size(sp))
end
end
sz
end
function gaston_init_subplot(plt::Plot{GastonBackend}, sp::Subplot{GastonBackend}) function gaston_init_subplot(plt::Plot{GastonBackend}, sp::Subplot{GastonBackend})
if sp === nothing if sp === nothing
push!(plt.o.subplots, sp) push!(plt.o.subplots, sp)
else else
dims = RecipesPipeline.is3d(sp) ? 3 : 2 sp.o = GastonSubplot(
dims=RecipesPipeline.is3d(sp) ? 3 : 2,
axesconf = gaston_parse_axes_args(plt, sp) # Gnuplot string axesconf=gaston_parse_axes_args(plt, sp), # Gnuplot string
sp.o = GastonSubplot(dims=dims, axesconf=axesconf, curves=[]) curves=[]
)
push!(plt.o.subplots, sp.o) push!(plt.o.subplots, sp.o)
end end
@ -129,45 +149,91 @@ function gaston_add_series(plt::Plot{GastonBackend}, series::Series)
sp = series[:subplot] sp = series[:subplot]
g_sp = sp.o # Gaston subplot object g_sp = sp.o # Gaston subplot object
seriesconf = gaston_parse_series_args(series) # Gnuplot string if series[:seriestype] (:heatmap, :contour) && g_sp.dims == 2
c = G.Curve(series[:x], series[:y], nothing, nothing, seriesconf ) g_sp.dims = 3 # FIXME: this is ugly, we need heatmap to use splot, not plot
end
x = series[:x]
y = series[:y]
z = g_sp.dims == 2 ? nothing : series[:z]
if z isa Surface
z = z.surf
end
seriesconf = gaston_seriesconf!(sp, series) # Gnuplot string
c = G.Curve(x, y, z, nothing, seriesconf)
isfirst = length(g_sp.curves) == 0 ? true : false isfirst = length(g_sp.curves) == 0 ? true : false
push!(g_sp.curves, c) push!(g_sp.curves, c)
G.write_data(c, g_sp.dims, g_sp.datafile, append = isfirst ? false : true) G.write_data(c, g_sp.dims, g_sp.datafile, append = isfirst ? false : true)
end end
function gaston_parse_series_args(series::Series) gaston_lc_ls_lw(series::Series) = (
gaston_color(series[:linecolor], series[:linealpha]),
gaston_linestyle(series[:linestyle]),
series[:linewidth],
)
gaston_mk_ms_mc(series::Series) = (
gaston_marker(series[:markershape]),
series[:markersize] * GASTON_MARKER_SCALING,
gaston_color(series[:markercolor], series[:markeralpha]),
)
gaston_palette(gradient) = begin
palette = String[]
n = -1
for rgba gradient # FIXME: naive conversion, inefficient ?
push!(palette, "$(n += 1) $(rgba.r) $(rgba.g) $(rgba.b)")
end
'(' * join(palette, ", ") * ')'
end
function gaston_seriesconf!(sp, series::Series)
gsp = sp.o
curveconf = String[] curveconf = String[]
st = series[:seriestype] st = series[:seriestype]
if st == :scatter clims = get_clims(sp, series)
pt = gaston_marker(series[:markershape]) if st (:scatter, :scatter3d)
ps = series[:markersize] * GASTON_MARKER_SCALING pt, ps, lc = gaston_mk_ms_mc(series)
lc = gaston_color(series[:markercolor])
# alpha = series[:markeralpha] # TODO merge alpha with rgb color
push!(curveconf, """with points pt $pt ps $ps lc $lc""") push!(curveconf, """with points pt $pt ps $ps lc $lc""")
elseif st == :path elseif st (:path, :straightline, :path3d)
lc = gaston_color(series[:linecolor]) lc, dt, lw = gaston_lc_ls_lw(series)
dt = gaston_linestyle(series[:linestyle])
lw = series[:linewidth]
# alpha = series[:linealpha] # TODO merge alpha with rgb color
if series[:markershape] == :none # simplepath if series[:markershape] == :none # simplepath
push!(curveconf, """with lines lc $lc dt $dt lw $lw""") push!(curveconf, """with lines lc $lc dt $dt lw $lw""")
else else
pt = gaston_marker(series[:markershape]) pt, ps = gaston_mk_ms_mc(series)
ps = series[:markersize] * GASTON_MARKER_SCALING
push!(curveconf, """with lp lc $lc dt $dt lw $lw pt $pt ps $ps""") push!(curveconf, """with lp lc $lc dt $dt lw $lw pt $pt ps $ps""")
end end
elseif st == :shape elseif st == :shape
fc = gaston_color(series[:fillcolor]) fc = gaston_color(series[:fillcolor], series[:fillalpha])
fs = "solid" fs = "solid"
lc = gaston_color(series[:linecolor]) lc, _ = gaston_lc_ls_lw(series)
push!(curveconf, """with filledcurves fc $fc fs $fs border lc $lc""") push!(curveconf, """with filledcurves fc $fc fs $fs border lc $lc""")
elseif st == :steppre elseif st == :steppre
push!(curveconf, """with steps""") push!(curveconf, """with steps""")
elseif st == :steppost elseif st == :steppost
push!(curveconf, """with fsteps""") # Not sure if not the other way push!(curveconf, """with fsteps""") # Not sure if not the other way
elseif st (:contour, :contour3d)
push!(curveconf, """with lines""")
if st == :contour
gsp.axesconf *= """\nset view map\nunset surface"""
end
levels = join(map(string, collect(contour_levels(series, clims))), ", ")
gsp.axesconf *= """\nset contour base\nset cntrparam levels discrete $levels """
elseif st (:surface, :heatmap)
palette = gaston_palette(series[:seriescolor])
gsp.axesconf *= """\nset palette model RGB defined $palette"""
if st == :heatmap
gsp.axesconf *= """\nset view map"""
end
push!(curveconf, """with pm3d""")
elseif st == :wireframe
lc, dt, lw = gaston_lc_ls_lw(series)
push!(curveconf, """with lines lc $lc dt $dt lw $lw""")
else
@warn "Gaston: $st is not implemented yet"
end end
# label # label
@ -215,7 +281,14 @@ function gaston_parse_axes_args(plt::Plot{GastonBackend}, sp::Subplot{GastonBack
ticks = get_ticks(sp, axis_attr) ticks = get_ticks(sp, axis_attr)
gaston_set_ticks!(axesconf, ticks, letter) gaston_set_ticks!(axesconf, ticks, letter)
end end
# set title {"<title-text>"} {offset <offset>} {font "<font>{,<size>}"}{{textcolor | tc} {<colorspec> | default}} {{no}enhanced}1
ratio = get_aspect_ratio(sp)
if ratio != :none
if ratio == :equal
ratio = -1
end
push!(axesconf, """set size ratio $ratio""")
end
end end
gaston_set_legend!(axesconf, sp) # Set legend params gaston_set_legend!(axesconf, sp) # Set legend params
@ -317,7 +390,11 @@ function gaston_marker(marker)
return 1 return 1
end end
gaston_color(color) = """rgb "#$(hex(color, :rrggbb))" """ gaston_color(color, alpha=0.) = begin
col = single_color(color) # in case of gradients
col = alphacolor(col, alpha == nothing ? 0. : alpha) # add a default alpha if non existent
"""rgb "#$(hex(col, :aarrggbb))" """
end
function gaston_linestyle(style) function gaston_linestyle(style)
style == :solid && return "1" style == :solid && return "1"

View File

@ -1242,7 +1242,14 @@ _backend_skips = Dict(
], ],
:inspectdr => [4, 6, 10, 22, 24, 28, 30, 38, 43, 45, 47, 48, 49, 50, 51, 55], :inspectdr => [4, 6, 10, 22, 24, 28, 30, 38, 43, 45, 47, 48, 49, 50, 51, 55],
:unicodeplots => [6, 10, 22, 24, 28, 38, 43, 45, 47, 49, 50, 51, 55], :unicodeplots => [6, 10, 22, 24, 28, 38, 43, 45, 47, 49, 50, 51, 55],
:gaston => [2, 4, 5, 6, 10, 22, 24, 28, 31, 32, 35, 38, 43, 45, 47, 48, 49, 50, 51, 55], :gaston => [
2, 4, 6,
16, # TODO: support nested layouts
27, # TODO: support polar
30, 31, 47, 48,
49, # TODO: support polar
50, 51, 55
],
) )