Re: GPR radargaram plot

I think you're missing the travel distance vector. I'm assuming you want to
plot the transmitter location starting at 0. You might need to change this
depending on what you want to plot for your horizontal axes (like if you
want the T/R midpoint).
dist = [0:Header4.TxStepX:(Header4.NSteps-1)*Header4.TxStepX];

Then I believe you should be able to display it using:
imagesc(dist,Fields4.t,squeeze(Fields4.hy))

You'll probably need to play with the aspect ratio to get it to be something
more like you're used to seeing (it's square by default).

Hope that helps.

Kim

"Ryan North" <ryan.e.north@usace.army.mil> wrote in message
news:kqvshr$ahn$1@newscl01ah.mathworks.com...
> Any luck solving this issue?
>
> "Cristina Porcel" <cristinaporcel@gmail.com> wrote in message
> <k6klgg$c6g$1@newscl01ah.mathworks.com>...
>> Hi,
>>
>> I'm simulating a ground penetrating response using GPRmax2D package which
>> output is a binary file. This had been read by Matlab using a .m file,
>> and finally the following structures appear.
>> Field=[1x1]structure
>> Header=[1x1]structure
>> [Header4, Fields4] = gprmax( 'output4.in' )
>> Header4 = title: [1x250 char]
>> iterations: 1357
>> dx: 0.0100
>> dy: 0.0100
>> dt: 2.3587e-011
>> NSteps: 50
>> TxStepX: 5
>> TxStepY: 0
>> RxStepX: 5
>> RxStepY: 0
>> ntx: 1
>> nrx: 1
>> nrx_box: 0
>> tx: 350
>> ty: 340
>> source: 'antena '
>> delay: 0
>> removed: 3.2000e-008
>> rx: 305
>> ry: 340
>>
>>
>> Fields4 = t: [1357x1 double]
>> ez: [1357x1x50 double]
>> hx: [1357x1x50 double]
>> hy: [1357x1x50 double]
>> My question is which command to use to plot the hy vx t structure in a
>> B-scan or radargram usually show in a color scale.
>> Each of the 50, vector of Hy matrix shows the EM field amplitud in the
>> 1357 iterations of variable t.
>> Thanks in advance.