data = fopen('Test.txt');
HKoord = textscan(data, '%f %f %f %f');
HSpalten = cell2mat(HKoord);
xx = HSpalten(:,1);                 % x
yy = HSpalten(:,2)';                % y
zz = HSpalten(:,3)';                % z
C = HSpalten(:,4)';                 % Daten
fclose(data);

cmax=40;                            % Datenbereich auswählen
cmin=20;
C_sel=C;                            % C_sel wird Vektor mit Werten zwischen ..
C_sel(C_sel>cmax)=[];               % .. cmax und ..
C_sel(C_sel<cmin)=[];               % .. cmin.

% Colormap für Werte ais C_sel aufbauen
cmap='jet';                               % cmap die Farbskala "jet" zuordnen
cmap=colormap(cmap);                      % Für colormap den Typ "jet" wählen
yyc=linspace(min(C_sel),max(C_sel),size(cmap,1));     % Farbpalette generieren
cm = spline(yyc,cmap',C_sel);             % Finden der Farbwerte
cm(cm>1)=1;                               % Wert darf nicht >1 sein
cm(cm<0)=0;                               % Wert darf nicht <0 sein

%%%%%%%%%%%%%%%%%%%%%%%%% Plot Ver.1
% k=1;
% for i=1:length(zz)-1
%     if C(i+1)==C_sel(k)
%         colchange=cm(:,k);                    % Colormap anwenden
%         for k=1:length(C_sel)
%         end
%     else
%         colchange=[0.7 0.7 0.7];             % Colormap anwenden
%     end
%     h(i)=line([zz(i) zz(i+1)],[xx(i) xx(i+1)],[yy(i) yy(i+1)],'color',colchange,'LineWidth',3);
% end

%%%%%%%%%%%%%%%%%%%%%%%%% Plot Ver.2
k=2;
for i=1:length(zz)-1
    if C(i+1)==C_sel(k-1)
        colchange=cm(:,k);  % Colormap anwenden
        k=k+1
    else
        colchange=[0.7 0.7 0.7];   % Grau färben wenn außerhalb des Wertebereichs
    end
    h(i)=line([zz(i) zz(i+1)],[xx(i) xx(i+1)],[yy(i) yy(i+1)],'color',colchange,'LineWidth',3);
end

% Colorbar für ausgewählten Wertebereich
view(3)
caxis([min(C_sel),max(C_sel)]);    % Max/Min der Colorbar
colorbar('Ticks',[min(C_sel),(min(C_sel)+max(C_sel))/2,max(C_sel)])