Contour plot in the XZ or YZ plane

[jagoosw]

Hi all,

I'm not sure this is currently possible, but I'm stuck trying to plot contours in the XZ (or YZ) planes. For example I have a 3D volume:

using GLMakie

data = randn(10, 10, 10)

fig = Figure()
ax = Axis3(fig[1, 1])
volume!(ax, [1:10;], [1:10;], [1:10;], -[1:10;], data)
contour!(ax, data[:, :, 10])

I'd like to plot similar contours on the XZ face but I can't find a way todo that. Since contour takes x::Vector, y::Vector, z::Matrix I tried just switching which were vectors vs matrix but that doesn't work.

Any suggestion would be great, thanks!

[ffreyer]

You can translate!() and rotate!() a plot to move it around in 3D space, or use transform!(plot, (plane, offset)) like volumeslices:

Makie.jl/src/basic_recipes/volumeslices.jl

Lines 40 to 52 in 714cf7a

	for (ax, p, r, (X, Y)) ∈ zip(axes, (:yz, :xz, :xy), (x, y, z), ((y, z), (x, z), (x, y)))
	plot[Symbol(:heatmap_, p)] = hmap = heatmap!(
	plot, attr, X, Y, zeros(length(X[]), length(Y[]))
	)
	plot[Symbol(:update_, p)] = update = i -> begin
	transform!(hmap, (p, r[][i]))
	indices = ntuple(Val(3)) do j
	axes[j] == ax ? i : (:)
	end
	hmap[3][] = view(volume[], indices...)
	end
	update(1) # trigger once to place heatmaps correctly
	end

[jagoosw]

Thank you!

[jagoosw]

For future reference:

using CairoMakie

fig = Figure()

ax = Axis3(fig[1, 1], limits = (1, 10, 1, 10, 1, 10))

data = [i*k for i=1:10, j=1:10,k = 1:10]
colorrange = (1, 100)

top = contour!(ax, [1:10;], [1:10;], data[:, :, 10]; colorrange)

translate!(top, 0, 0, 10)

side = contour!(ax, [1:10;], [1:10;], data[:, end, :]; colorrange)

rotate!(side, √2/2, 0, 0, √2/2)

bottom = contour!(ax, [1:10;], [1:10;], data[:, :, 1]; colorrange)