Setting the z range in bar3 plot.

[LeeoBianchi]

Hi! Amazing package!

I am trying to plot a 3-d histogram through bar3 with the z-axis starting from 0.0. Is this done through the keyword setting base = 0.0? It seem to get ignored!

G = GMT.mat2grid(randn(5, 6)) cmap = grd2cpt(G); bar3(G, lw=:thinnest, color=cmap, base=0.0, figname="prova.png")

Thank you!

[joa-quim]

Hi! Amazing package!

Thanks, I agree 😃

Humm, it works for me

G = GMT.mat2grid(randn(5, 6));
cmap = grd2cpt(G);
bar3(G, lw=:thinnest, color=cmap, base=0.0, figname="prova.png")

[joa-quim]

Sorry, I overlooked the base=0.0 part.

[joa-quim]

https://forum.generic-mapping-tools.org/t/is-this-a-plot3d-s0-bug/3748

[joa-quim]

Fixed in #1138

[LeeoBianchi]

@joa-quim now the base works! thanks a lot!!

Just a note: in order to adjust the frame accordingly, after having changed base, I had to set both region = [5 65 0.5 5.5 0 300] AND limits = (5, 65, 0.5, 5.5, 0, 300). I don't know it is expected but for me it was not obvious to figure out, maybe something worth mentioning in the docs :)

Cheers,
Leo

[joa-quim]

No, that is clearly a bug. region and limits are synonyms.

[LeeoBianchi]

No, that is clearly a bug. region and limits are synonyms.

Hmm okay, then definitely a bug!

Moreover, when I am specifying both region and limits, the base keyword becomes unnecessary.
Here is the code to reproduce it:

G = GMT.mat2grid(randn(5, 6).+300, x=[10 20 30 40 50 60], y=[1 2 3 4 5])
cmap = grd2cpt(G);
bar3(G,
     lw=:thinnest,
     color=cmap,
     base=0.0,
     region = [5 65 0.5 5.5 0 600],
     limits = (5, 65, 0.5, 5.5, 0, 600),
     xaxis=(annot=10, ticks=:auto, label="X"),
     yaxis=(annot=1, ticks=:auto, label="Y"),
     zaxis=(annot=:auto, ticks=:auto, label="Z"),
     show = true,
     figsize = (10.5, 9),
     figname="test.png"
     )

[joa-quim]

Fixed limits/region issue in #1139

Yes, limiting with region to positive only is another way of setting base=0.0

[LeeoBianchi]

Fixed limits/region issue in #1139

Amazing!

Yes, limiting with region to positive only is another way of setting base=0.0

Yes, that makes sense.

One more thing: what if I wanted to put custom annotations on one of the axis, but leaving the ticks where they are? I haven't found a way to make this work.

In particular in my case I would like a 2 in place of the 0. But if I try to use the xticks keyword I get all the ticks compressed around 0 and the Z axis disappears. Am I doing something wrong?

G = GMT.mat2grid(randn(5, 7).+300, x=[0 10 20 30 40 50 60], y=[1 2 3 4 5])
cmap = grd2cpt(G);
bar3(G,
     lw=:thinnest,
     color=cmap,
     region = [-5 65 0.5 5.5 0 600],
     limits = (-5, 65, 0.5, 5.5, 0, 600),
     xticks=(2, 10, 20, 30, 40, 50, 60),
     show = true,
     figsize = (10.5, 9),
     figname="prova.png"
     )

Cheers,
Leo

[joa-quim]

You can do this but there must be a better (cleaner to the reader) way

bar3(G, lw=:thinnest, color=cmap, region = [-5 65 0.5 5.5 0 600], xaxis=(custom=(pos=[2, 10, 20, 30, 40, 50, 60], type=["a","a","a","a","a","a","a"]),), yaxis=:a, show=1)

[joa-quim]

Simpler this way

bar3(G, lw=:thinnest, color=cmap, region = [-5 65 0.5 5.5 0 600], xticks=([2, 10, 20, 30, 40, 50, 60], ("2","10","20","30","40","50","60")), show=1)

[LeeoBianchi]

Amazing, it worked! Thank you so much.

Last question, I promise: is there a way to show a grid (or simply thin horizontal lines for each tick) in the background, on the zy plane, for instance?

Leo

[joa-quim]

No problem with questions, only that I would prefer that they are made in the GMT forum.

The answer is off course yes but needs to dig from the (frame) docs to find the right option.

[joa-quim]

Can we close this?

[LeeoBianchi]

Sure, thanks for all the useful tips!