[Tex/LaTex] addplot3 [area plot] filling the wrong area

3dpgfplots

using some saved data, i use (The data file cane be download from pX.txt where the first column is the x and at each iteration, I plot the other columns with respect to the first column)

    \documentclass{article}
\usepackage{pgfplots,pgfplotstable}

\pgfplotstableread{pX.txt}\datatablePX%          

\begin{document}


\begin{tikzpicture}
\begin{axis}[
    zmin=10^5,
    area plot/.style={
        fill opacity=0.75,
        draw=black!50!black,
        fill=blue!50,
        mark=none,
    },
    xlabel={$x[m]$},ylabel={Time},zlabel={p[bar]}
]

\pgfplotsinvokeforeach{120,119,...,1}{
  \addplot3 [area plot]  table [x index=0, y expr=#1, z index=#1] {\datatablePX};
}


\end{axis}
\end{tikzpicture}

\end{document}

to plot some data in 3d, where I want the area under the curves to be filled. However, the wrong area is being filled as one can notice in the following plot.

enter image description here

It is filling the area to a line plot between the start and the end values. Why doesn't it simply fill with respect to the x-axis?

Best Answer

You will need to cycle the path by hand, which is definitely not the best solution, but it works. So I'm doing the following:

\addplot3[area plot] table[x index=0,y expr=#1,z index=#1] {\datatablePX}
-- (axis cs:\pgfkeysvalueof{/pgfplots/xmax},#1,\pgfkeysvalueof{/pgfplots/zmin})
-- (axis cs:\pgfkeysvalueof{/pgfplots/xmin},#1,\pgfkeysvalueof{/pgfplots/zmin})
-- cycle;

This will look horribly with enlargelimits=true, so it is best to set your limits by hand or using enlargelimits=false.

If you don't want to use zmin for the lower filling boundary, you can specifiy an own key, e.g. zcycle and use the values assigned to this key. See the MWE for details.

Implementation

I also tweaked your code a little, to have a little more speaking axis labels. Use pX.txt, uploaded by OP. I didn't plot all values for better performance.

\documentclass[tikz]{standalone}
\usepackage{pgfplots,pgfplotstable}
\tikzset{
    /pgfplots/zcycle/.initial=,
}
\pgfplotstableread{pX.txt}\datatablePX%          
\begin{document}
\begin{tikzpicture}
    \begin{axis}[
            zmin=0.95*10^5,
        zcycle=10^5,
            area plot/.style={
                fill opacity=0.75,
                draw=black!50!black,
                fill=blue!50,
                mark=none,
            },
            xlabel={Distance $x$ [m]},
            ylabel={Time $t$ [s]},
            zlabel={Pressure $p$ [bar]}
        ]
        \pgfplotsinvokeforeach{120,80,...,1}{
            \addplot3[area plot] table[x index=0,y expr=#1,z index=#1] {\datatablePX}
            -- (axis cs:\pgfkeysvalueof{/pgfplots/xmax},#1,\pgfkeysvalueof{/pgfplots/zcycle})
            -- (axis cs:\pgfkeysvalueof{/pgfplots/xmin},#1,\pgfkeysvalueof{/pgfplots/zcycle})
            -- cycle;
        }
    \end{axis}
\end{tikzpicture}
\end{document}

Output

Related Solutions

[Tex/LaTex] plotting two time series with bounds

This happens because PGFPlots only uses one "stack" per axis: You're stacking the second confidence interval on top of the first. The easiest way to fix this is probably to use the approach described in "Is there an easy way of using line thickness as error indicator in a plot?": After plotting the first confidence interval, stack the upper bound on top again, using stack dir=minus. That way, the stack will be reset to zero, and you can draw the second confidence interval in the same fashion as the first:

\documentclass{standalone}
\usepackage{pgfplots, tikz}

\usepackage{pgfplotstable}

\pgfplotstableread{
temps   y_h y_h__inf    y_h__sup    y_f y_f__inf    y_f__sup    

1   0.237340    0.135170    0.339511    0.237653    0.135482    0.339823    
2   0.561320    0.422007    0.700633    0.165871    0.026558    0.305184    
3   0.694760    0.534205    0.855314    0.074856    -0.085698   0.235411    
4   0.728306    0.560179    0.896432    0.003361    -0.164765   0.171487    
5   0.711710    0.544944    0.878477    -0.044582   -0.211349   0.122184    
6   0.671241    0.511191    0.831291    -0.073347   -0.233397   0.086703    
7   0.621177    0.471219    0.771135    -0.088418   -0.238376   0.061540    
8   0.569354    0.431826    0.706882    -0.094382   -0.231910   0.043146    
9   0.519973    0.396571    0.643376    -0.094619   -0.218022   0.028783    
10  0.475121    0.366990    0.583251    -0.091467   -0.199598   0.016664    
}{\table}

\begin{document}
\begin{tikzpicture}

    \begin{axis}
    % y_h confidence interval
    \addplot [stack plots=y, fill=none, draw=none, forget plot]   table [x=temps, y=y_h__inf]   {\table} \closedcycle;
    \addplot [stack plots=y, fill=gray!50, opacity=0.4, draw opacity=0, area legend]   table [x=temps, y expr=\thisrow{y_h__sup}-\thisrow{y_h__inf}]   {\table} \closedcycle;
    % subtract the upper bound so our stack is back at zero
    \addplot [stack plots=y, stack dir=minus, forget plot, draw=none] table [x=temps, y=y_h__sup] {\table};

    % y_f confidence interval
    \addplot [stack plots=y, fill=none, draw=none, forget plot]   table [x=temps, y=y_f__inf]   {\table} \closedcycle;
    \addplot [stack plots=y, fill=gray!50, opacity=0.4, draw opacity=0, area legend]   table [x=temps, y expr=\thisrow{y_f__sup}-\thisrow{y_f__inf}]   {\table} \closedcycle;

    % the line plots (y_h and y_f)    
    \addplot [stack plots=false, very thick,smooth,blue]  table [x=temps, y=y_h]   {\table};
    \addplot [stack plots=false, very thick,smooth,blue]  table [x=temps, y=y_f]   {\table};
    \end{axis}

\end{tikzpicture}
\end{document}

[Tex/LaTex] Draw box around plot or fill plot area in tikz/pgfplots

You can use the \pgf@pathminx and \pgf@pathmaxy commands, which hold the current path's minimum and maximum x and y coordinates, to draw a rectangle around the plot. Here's a macro \drawplotbox that can be called at the end of the \addplot command. The style of the box is controlled by a key called plot box, the padding of the box can be set using plot box padding:

\documentclass{standalone}
\usepackage{tikz,pgfplots}

\makeatletter
\newcommand{\drawplotbox}{
(\pgf@pathminx,\pgf@pathminy) coordinate (current plot south west)
            (\pgf@pathmaxx,\pgf@pathmaxy) coordinate (current plot north east);
            \draw [plot box] ([shift={(-135:1.41*\pgfkeysvalueof{/tikz/plot box padding})}]current plot south west) rectangle ([shift={(45:1.41*\pgfkeysvalueof{/tikz/plot box padding})}]current plot north east);
}

\tikzset{
    plot box/.style={
        thick, rounded corners, orange
    },
    plot box padding/.initial=5pt
}

\begin{document}
    \begin{tikzpicture}
        \begin{axis}
            \addplot coordinates{
                (0,0)
                (2,1)
                (1,1)
                (2,2)
                (0,2)
            } \drawplotbox;
            \addplot +[plot box/.append style=black, plot box padding=10pt] coordinates{
                (3,3)
                (5,4)
                (4,4)
                (5,5)
                (3,5)
            } \drawplotbox;
        \end{axis}
    \end{tikzpicture}
\end{document}