How to open Tonatiuh output .DAT files in MATLAB?

Hi, I have output files from Tonatiuh ray-tracing software but I am unable to open this file in any other software. I was wondering if there is any way to open it in MATLAB so that I can perform the analysis that I need.

Any help would be appreciate in this regard. This is a binary output file with a .DAT extension. I am unable to uplaod the file here so I am attaching a One-Drive link if anybody wants to have a look.

https://drive.google.com/file/d/1e6J5ix29uK20r1WKoIg3zyze7Es1Ksfz/view?usp=sharing

[Image Analyst]

You have to know the format of the file. Do you know that? Otherwise it'll just be a big long string of bytes with no meaning behind them.

[venu]

The format of the file is .dat

Answer 1

[Walter Roberson]

According to https://code.google.com/archive/p/tonatiuh/

The binary format used to store the values is Real64 with big-endian byte ordering. The

first value stored in the file is the power per photon in watts (W). This real value

is then followed by seven-tuples of real values, one per photon . Every seven-tuple

holds the following values in the specific order given: * PhotonID * x-coordinate *

y-coordinate * z-coordinate * surface side ( 0 if is back side , 1 if is front side ) *

Previous PhotonID (0 if the photon is the first of a given ray) * Next PhotonID (0 is

the photon is the last of a given ray)

So you would fopen() the file with 'ieee-be' and you would fread() '*double* with size [7 inf] . The result would store tuples running down columns.

[venu]

Respected sir,

I am not understanding what is present in binary output file when i was extracted by using matlab code it is containing some numeric values like zeros and some numeric values please help me humble request sir

For the reference i am just attched the output what i have received

[Walter Roberson]

[fid, msg] = fopen(FILENAME, 'ieee-be');
if fid < 0
    error('failed to open file "%s" because "%s"',, FILENAME, msg);
end
A = fread(fid, [7 inf], '*double');
fclose(fid);

Now the columns of A correspond to the seven-tuples described above.

[venu]

I am trying to calculate optical efficieny ,number rays receiving the absorber etc,. by using following code but i am getting error can you please help me here by i am sharing pythoncode and .dat and .txt for the reference

https://docs.google.com/document/d/1INbSxrreAxRnkSgHdjtvgs9mnL0WP3oS/edit?usp=drive_link

https://drive.google.com/file/d/1VKEdJ6-cSnUiGlLj9VStN7xjv2MU-V5U/view?usp=drive_link

https://drive.google.com/file/d/1MGNbn1dnm5Ziz5oSgN3oiDbZ5Byhz7po/view?usp=drive_link

[venu]

Respected sir,

when i run the code i am getting the following error

please help me humbel request

[Walter Roberson]

FILENAME = '/MATLAB Drive/datfile/xt_1.dat';
[fid, msg] = fopen(FILENAME, 'ieee-be');
if fid < 0
    error('failed to open file "%s" because "%s"',, FILENAME, msg);
end
A = fread(fid, [7 inf], '*double');
fclose(fid);

[venu]

RESPETED SIR

CODE IS NOT RUNNING

NOW I AM GETTING THIS ABOVE ATTACHED ERROR

[Walter Roberson]

The data in the file is not in the documented format. My guess is that the format is

1. x coordinate
2. y coordinate
3. z coordinate
4. previous photon ID
5. next photon ID

FILENAME = 'sahasra_1.dat';
[fid, msg] = fopen(FILENAME, 'r', 'ieee-be');
if fid < 0
    error('failed to open file "%s" because "%s"', FILENAME, msg);
end
Power_per_photon_W = fread(fid, 1, '*double');
A = fread(fid, [5 inf], '*double');
fclose(fid);

[venu]

Respected sir,

I am trying to calculate optical efficieny ,number rays receiving the absorber etc,. by using following code but i am getting error can you please help me here by i am sharing pythoncode and .dat and .txt for the reference

Can you please help me i am really very grateful to you

https://docs.google.com/document/d/1INbSxrreAxRnkSgHdjtvgs9mnL0WP3oS/edit?usp=drive_link

https://drive.google.com/file/d/1VKEdJ6-cSnUiGlLj9VStN7xjv2MU-V5U/view?usp=drive_link

https://drive.google.com/file/d/1MGNbn1dnm5Ziz5oSgN3oiDbZ5Byhz7po/view?usp=drive_link

[Walter Roberson]

That is python code. We are not a python resource.

[venu]

Respected sir,

Recently i have found Matlab code for calculating flux and number of rays absorbing the receiveretc.,

But i am unable to find where to give the path for .dat file and .txt file

please kindly help me I am really very very grateful to you sir Sorry for disturbance.

https://docs.google.com/document/d/1-o_gszbtw7IDlEs7km-WFHAv3sV_jJw_/edit?usp=sharing&ouid=110023858447791699573&rtpof=true&sd=true

[Walter Roberson]

See attached

[venu]

Respected sir, Where to give path for.dat and .txt files

[Walter Roberson]

Dinfo = dir('*.dat');
Watts_per_photon = 183.19;   %adjust as appropriate
[First, Second, Third] = raytrace_function(Dinfo, Watts_per_photon);

[venu]

Respected sir I am tried many ways but i am unable to figure it out where i am making mistake. I don't know weather ask this or not.can we have professional discussion on what's app chat this is my number [phone] its my humble request sir. Please accept.

[Walter Roberson]

I do not use What's App at all.

[venu]

Respected sir, I want to calculate optical Efficiency buy using photon mapping .dat for phd thessis your small help really create a lot of confidence in me. Can i send you google invite Link

Answer 2

[Benjamin Thompson]

fid = fopen('yourfilename.dat', 'r');

A = fread(fid, 'uint8=>uint8');

fclose(fid);

Then the data is in A as bytes. If you know more about the file format you can change the arguments to fread in order to make the result in A easier to use.

doc fread

[venu]

Sir can you please help me in binary files

[venu]

Sir i have tried but i am unable to extract data can you please data and send me humble request