Oberflächendiagramm nur aus Z-Daten

Question

Ich möchte eine vermaschte Oberfläche aus einer Datei zeichnen, die nur Z-Daten enthält, die in einer Matrix angeordnet sind. Die Idee besteht darin, die Z-Daten gegenüber dem X-Index und dem Y-Index darzustellen. Gibt es eine einfache Möglichkeit, dies zu tun, oder muss ich auf die Vorverarbeitung der Datendatei zurückgreifen?

(Ich habe erfolglos im pgfplotsHandbuch und in mehreren StackExchange-Fragen gesucht.)

Beispiel für Dummydaten:

CSV data file `data.txt'
M1,M2,M3,M4
0,0,0,0
3.466,2.058,0,0
4.262,2.976,0.001,0
3.822,3.168,0.006,0.008
2.953,2.936,0.019,0.063
2.065,2.492,0.046,0.265
1.332,1.977,0.092,0.734
0.797,1.478,0.164,1.508

MNWE (für minimales nicht funktionierendes Beispiel ^^):

\documentclass{minimal}
\usepackage{pgfplots,pgfplotstable}  
\begin{document}
\begin{tikzpicture}
   \pgfplotstableread[col sep=comma]{data.txt}\csvdata
   \begin{axis}
      \addplot3[surf] table[x expr=\coordindex,y expr={SOME COLUMN INDEX?}] {\csvdata};
   \end{axis}
\end{tikzpicture}
\end{document}