TIPS & TRICKS GEORADAR

(ramac pulsekko gssi data handling)

Puoi vedere i dati georadar Ramac con  You can display gpr Ramac data with

Puoi vedere i dati georadar Pulsekko con  You can display gpr Pulsekko data with

 

Puoi vedere i dati georadar GSSI con  You can see GSSI data

[index][»]

Plotting Ramac georadar data

Ramac data are three files: *.rad (headers), *.rd3 (integer 16bits little-endian) e *.mrk (ASCI marker).

[«][index][»]

(pc little-endian)

Open file for reading

Read and import file

Close file

Plot A

 fid=fopen('file.rd3','r');

A=fread(fid,[samples,traces],'int16');

fclose(fid);

contour(A);

(workstation big-endian)

Open file for reading

Read and import file

Close file

Plot A

 

 fid=fopen('file.rd3','r','ieee-le');

A=fread(fid,[samples,traces],'int16');

fclose(fid); 

contour(A); 

[«][index][»]

(pc little-endian)

Open file

Read and import file

Close file

Plot A

 fid=fopen('file.rd3','r');

A=fread(fid,[samples,traces],'short');

fclose(fid);

mesh(A);

[«][index][»]

(pc little-endian)

 

 

 

You must

install

Numeric

and

Dislin

 import Numeric
import dislin
import struct
from Numeric import *
from struct import *
from dislin import *

fn=open('file.rd3','rb')

s=fn.read()

dislin.metafl('xwin')
dislin.disini()
dislin.pagera()
dislin.hwfont()

setvar('X','samples')
setvar('Y','traces')
setvar('Z','amplitudes')
dislin.titlin('GPRadar section',1)
dislin.title()
A=Numeric.fromstring(s,'s')
shapeA=shape(A)
traces=shapeA[0]/samples
B=reshape(A,(traces,samples))
Xaxis=arrayrange(samples)
Yaxis=arrayrange(traces)
surf3(transpose(B[:,:],Xaxis,Yaxis)

fn.close

[«][index][»]

(pc little-endian)

Open file

Read and import file

Close file

Grid

Plot A

 filen='file.rd3';

fid=mopen(filen,'rb'); 

for j=1:traces,

    for i=1:samples,

            A(i,j)=mget(1,'s',fid);

    end;

end;

mclose();

x=[1:1:samples];y=[1:1:traces];

Sgrayplot(x,y,A); 

[«][index][»]

(pc little-endian)

Convert from 16 integer to 32 float

Transform into *.su

Insert time sample

Plot data

Delete ghost file

 recast < file.rd3 in=short out=float > a.f  

suaddhead < a.f > a.su ns=samples  

sushw < a.su > aa.su key=dt a=samples*100

suximage < aa.su

rm a.f a.su

(workstation big-endian)

In this case we must swap data .... and then like little-endian

swapbytes < file.rd3 in=short > a.sw

recast < a.sw in=short out=float > a.f

...

[«][index][»]

Plotting Pulsekko georadar data

Pulsekko data are two files: *.hd (headers), *.dt1 (headers +integer 16bits little-endian).

 

[«][index][»]

(pc little-endian)

Download four files

Use: [data,head]=dt1readt(file)

plot(data)

 

 dt1read.m

dt1struct.m

headerstruct.m

dimstruct.m
If you have a GSSI radar (like Sir2), you can separate the headers from the data with this my tool.

FREE SOFTWARE!

DSISoft Refraction traveltime tomography GUI.
Mathtools A big portal for your scientific and engineering needs!
Octave Clone Matlab.
Python A big interpreted language.
Scilab Good package supported by French Government. 
Seismic Unix The famous free software for seismic processing. Thank to J.W. Stockwell and CWP!! 
SepLib Another good seismic software from Stanford Exploration Project.
SeisWide Free program for ray tracing and more.

If you know other free software, tell me!

Created 12/08/98

back to my homepage