Welcome, Guest
Username: Password: Remember me

TOPIC: Extract values along a line/ (or multiple x-sections)

Extract values along a line/ (or multiple x-sections) 4 years 8 months ago #35619

  • Karki
  • Karki's Avatar
Dear all,

I want to extract the bottom elevation of my simulation at different cross-section lines. I imported the x-section lines from ArcGIS to Blue Kenue.

There is the option to extract the time-series along a line [Tools>Extract Time Series>Along a Line....]

This option gives the min/max and average values for the total duration of the simulation.
However, I need the elevation along an x-section at the end of my simulation [lets say at T=600secs].

In this case, how can I extract the time series?

Any suggestions are appreciated.

Thank you.

Saroj
The administrator has disabled public write access.

Extract values along a line/ (or multiple x-sections) 4 years 8 months ago #35623

  • a.barton
  • a.barton's Avatar
  • OFFLINE
  • Moderator
  • Posts: 683
  • Thank you received: 178
With COVID-19 I have been asked to work from home so I don't have all of my usually available things with me here...

Once I am back at work again I can look into the details and provide some guidance. Hopefully in 2 weeks...

Sorry for the delay on this...

Kindly,
Alan
The administrator has disabled public write access.

Extract values along a line/ (or multiple x-sections) 4 years 8 months ago #35624

  • Jose_filho
  • Jose_filho's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 42
  • Thank you received: 2
Hi,

I just wrote this script to access and extract data in a point of the mesh. It is in matlab:

"clc
clear all
close all

filename = 'xxx3D.slf';
filename2D = 'xxx2D.slf';

m = telheadr(filename);
m2 = telheadr(filename2D);
IPOIN =18914;
stp = m.NPOIN/m.NPLAN;

ti = 1;%385 KE
tf = 1;%481
tss =[];
u = [];
v = [];
p = [];
k = 0;
sal=[];
for j = ti:tf %m.NSTEPS
m = telstepr(m,j);
m2 = telstepr(m2,j);

for i = 1:m.NPLAN %
wd(j) = m2.RESULT(IPOIN,3);
fs(j) = m2.RESULT(IPOIN,4);
u(j,i) = m.RESULT(IPOIN+(stp*(i-1)),2);
v(j,i) = m.RESULT(IPOIN+(stp*(i-1)),3);
p(j,i) = m.RESULT(IPOIN+(stp*(i-1)),1);
sal(j,i) = m.RESULT(IPOIN+(stp*(i-1)),5);
rd(j,i) = p(j,i) - fs(j); %depth


% IKLE=m.IKLE;
% ikle=double(IKLE);
% tr = triangulation(ikle,long,lat,m.RESULT(:,5));
end

ts = datetime(m2.IDATE(1,1),m2.IDATE(1,2),m2.IDATE(1,3),m2.IDATE(1,4),m2.IDATE(1,5),m2.IDATE(1,6)+(m2.AT));
tss = [tss,ts];



end"

To extract the results along more than one points, just make an external loop with more points. It is a bit slow, maybe you can improve it.

Hope it can help you,

Best,


José
The administrator has disabled public write access.
The following user(s) said Thank You: Karki

Extract values along a line/ (or multiple x-sections) 4 years 8 months ago #35628

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
Hi
In fact you should use the map object function...
Take care to have a well resampled X-section
And just map the result at the time needed on the line

Regards
Christophe
The administrator has disabled public write access.

Extract values along a line/ (or multiple x-sections) 4 years 8 months ago #35636

  • Karki
  • Karki's Avatar
Thank you Coulet.

In fact, I have used 'map object' options on many occasions for extracting velocity at different x-sections.

I totally forgot about this option. Thank you for making me realize the map object tools.
I guess that will do for my purpose.

Regards,
Saroj
The administrator has disabled public write access.
Moderators: Serrer, a.barton

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.