FE triangle mesh without the mesh node connections visible

2 ビュー (過去 30 日間)
Lara
Lara 2023 年 12 月 14 日
コメント済み: Lara 2023 年 12 月 18 日
I have the following code to generate image results from my .txt file:
It works, but I want to not show the mesh node connections in the plots WITHOUT using the scatter function. The pictures should remain the same without the black line connections. How do i realize this?
clc
close all
clear
% Read the data from the TXT file
fileID = fopen('2D_brain_1_object_src8_det8_simple_inv_inv_final.txt', 'r');
headerLines = 2; % Starting header lines
header = textscan(fileID, '%s', headerLines, 'Delimiter', '\n');
fclose(fileID);
% Extract the value after E=
N_value = str2double(strsplit(header{1}{2}, {'N=', ' '}));
nan_columns = any(isnan(N_value), 1);
N_value(:, nan_columns) = [];
num_coordinates_lines = N_value; % Number of lines for coordinates and properties
% Read the coordinates and properties data
fileID = fopen('2D_brain_1_object_src8_det8_simple_inv_inv_final.txt', 'r');
data = textscan(fileID, '%f%f%f%f%f%f', num_coordinates_lines, 'HeaderLines', headerLines);
fclose(fileID);
X = data{1};
Y = data{2};
Z = data{3};
permittivity = data{5};
conductivity = data{6};
% Read the mesh connections separately using textscan
fileID = fopen('2D_brain_1_object_src8_det8_simple_inv_inv_final.txt', 'r');
meshData = textscan(fileID, '%d%d%d', 'HeaderLines', headerLines + num_coordinates_lines);
fclose(fileID);
% Extract mesh connections
node1 = meshData{1};
node2 = meshData{2};
node3 = meshData{3};
% Create mesh connections matrix
meshConnections = [node1, node2, node3];
% Plotting the mesh based on permittivity levels
figure;
subplot(1, 2, 1);
h1 = trisurf(meshConnections, X, Y, Z, permittivity);
colorbar;
title('Mesh with Permittivity Levels');
xlabel('X [$m$]');
ylabel('Y [$m$]');
zlabel('Z [$m$]');
view(2); % Set view to 2D (X-Y plane)
axis equal
% Plotting the mesh based on conductivity levels
subplot(1, 2, 2);
h2 = trisurf(meshConnections, X, Y, Z, conductivity);
colorbar;
title('Mesh with Conductivity Levels');
xlabel('X [$m$]');
ylabel('Y [$m$]');
zlabel('Z [$m$]');
view(2); % Set view to 2D (X-Y plane)
axis equal
  1 件のコメント
Avni Agrawal
Avni Agrawal 2023 年 12 月 18 日
Hi, can you please attach the imported files for better debugging?

サインインしてコメントする。

採用された回答

KSSV
KSSV 2023 年 12 月 18 日
h1 = trisurf(meshConnections, X, Y, Z, permittivity);
h1.EdgeColor = 'none' ;

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeAnnotations についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by