Quero usar pgfplots em látex para plotar dados de um arquivo que foi formatado para uso com a função de índice do gnuplot. O arquivo de dados possui conjuntos de colunas x,y separadas por duas linhas em branco (formato gnuplot).
10 10
20 18
30 36
10 11
20 20
30 41
Com o gnuplot, posso fazer
plot 'data.dat' index 0 u 1:1 with p lw 2 pt 4 ps 1 lc rgb 'red' t 'One, \
'' index 1 u 1:2 with p lw 2 pt 5 ps 1 lc rgb 'blue' t 'Two'
Sem reescrever o arquivo de dados em um formato diferente, posso plotar esse formato de dados com pgfplots? Existe algo no pgfplots que seja equivalente ao índice do gnuplot? O índice pgfplot refere-se aos números das colunas em dados de várias colunas, não ao segundo bloco de dados, pelo que entendi.
Responder1
Primeiro, vamos ver como o gnuplot analisa as coisas. Pelo que entendi index 0 u 1:2
significa: observe o primeiro bloco (indexado em zero) e use a coluna 1 comoxe coluna 2 comosimvalores. Isso pode ser mostrado plotando os dados usando o gnuplottex
pacote:
\documentclass[border=10pt]{standalone}
\usepackage{xcolor, gnuplottex}
\begin{filecontents}[noheader]{data.dat}
10 10
20 18
30 36
10 11
20 20
30 41
\end{filecontents}
\begin{document}
\begin{gnuplot}[terminal=epslatex]
plot 'data.dat' index 0 u 1:2 with p lw 2 pt 4 ps 1 lc rgb 'red' t 'One', \
'' index 1 u 1:2 with p lw 2 pt 4 ps 1 lc rgb 'blue' t 'Two'
\end{gnuplot}
\end{document}
Isso significa: só precisamos anexar a cada linha um índice que identifique o bloco em que estamos atualmente.
Provavelmente é mais fácil fazer isso com algum outro software, mas você pode deixar o TeX converter o arquivo e anexar uma coluna de índice a cada linha. Você pode então filtrar os dados usandoesta abordagem.
O seguinte lê o arquivo e anexa um índice (um número inteiro começando em 1) que é aumentado se duas linhas em branco consecutivas forem encontradas. Os novos dados são gravados no arquivo data_index.dat
que você pode ler com PGFPlots.
\documentclass[border=10pt]{standalone}
\usepackage{pgfplots}
\pgfplotsset{compat=1.18}
\begin{filecontents}[noheader]{data.dat}
10 10
20 18
30 36
10 11
20 20
30 41
\end{filecontents}
\ExplSyntaxOn
\cs_generate_variant:Nn \iow_now:Nn { Ne }
\int_new:N \l_csnl_gnutopgf_index_int
\int_set:Nn \l_csnl_gnutopgf_index_int { 1 }
\bool_new:N \l_csnl_gnutopgf_blank_line_bool
\bool_set_false:N \l_csnl_gnutopgf_blank_line_bool
\ior_new:N \g_csnl_gnutopgf_input_ior
\ior_open:Nn \g_csnl_gnutopgf_input_ior { data.dat }
\iow_new:N \g_csnl_gnutopgf_output_iow
\iow_open:Nn \g_csnl_gnutopgf_output_iow { data_index.dat }
\iow_now:Nn \g_csnl_gnutopgf_output_iow { index ~ x ~ y }
\ior_str_map_inline:Nn \g_csnl_gnutopgf_input_ior {
\bool_if:NT \l_csnl_gnutopgf_blank_line_bool {
\str_if_empty:nT {#1} {
\int_incr:N \l_csnl_gnutopgf_index_int
}
\bool_set_false:N \l_csnl_gnutopgf_blank_line_bool
}
\str_if_empty:nTF {#1} {
\bool_set_true:N \l_csnl_gnutopgf_blank_line_bool
} {
\bool_set_false:N \l_csnl_gnutopgf_blank_line_bool
\iow_now:Ne \g_csnl_gnutopgf_output_iow {
\int_eval:n { \l_csnl_gnutopgf_index_int } ~ #1
}
}
}
\ior_close:N \g_csnl_gnutopgf_input_ior
\iow_close:N \g_csnl_gnutopgf_output_iow
\ExplSyntaxOff
\begin{document}
\pgfplotsset{
only index/.style={
x filter/.code={
\edef\tempa{\thisrow{index}}
\edef\tempb{#1}
\ifx\tempa\tempb\else
\def\pgfmathresult{inf}
\fi
}
}
}
\begin{tikzpicture}
\begin{axis}[legend pos=north west]
\addplot table[x=x, y=y, only index=1]{data_index.dat};
\addlegendentry{One}
\addplot table[x=x, y=y, only index=2]{data_index.dat};
\addlegendentry{Two}
\end{axis}
\end{tikzpicture}
\end{document}