MATLAB 开发 GPR app记录

1 基本语法问题

1.1 如何 hold on数据

 需要用 plot( axes, "on"); 实现。

0714 今天只完成了B-scan缩放问题。

0715，完成了B-scan相关的缩放工作

考虑通过设定每page显示的GPR 道数，以及 可以缩放2种模式。

还有小bug

1.2 数据格式问题

新的GPR都是32bit，之前的大多为16bit，差别在于读取时字节长度不同，前者为4bytes，后者为2bytes。

处理，判断字节长度，如果为32bit，则数据类型用'long'， 否则用'short'

 datatype = 'long';%'uint32';   
DATA.d=fread(fid,datasize,datatype);
%datasize  :[samples per scan, A-scan number]

1.3 rd7 文件的头文件读取存在问题

rd7的rad头文件较rd3多， 在后续参数读取时报错。可以注释掉。

后面修改。

rd3的 rad文件

 

 rd7的文件

  针对32bit的特点，将mala和dzt 的文件读取问题都解决。写的功能，仅测试了rd3和dzt。 rd7格式，回头处理。

1.4  界面截图

    习惯了LabVIEW的截图，这里也做了一个界面截图。看中部下方的小红点了吧

    后面做一个设置，自动将B-scan生成。

1.5 F-K偏移

  西南交大的熊洪强（他在科学网上开放了代码，主要做信号后处理，做的非常好！！） 对FK函数做了优化。我在此基础上进行功能测试

  其中有一个fft2 为系统的2D fft变换，但之前做ANN建模自己定义了一个fft2的函数，导致一直报错。

 经过调试，程序能正常运行，但是偏移效果没有展示出来。

 dz=(vmig)/2*Range/(ns-1);   这一句未理解。。。

vmig为速度，Range为总时窗，ns为A-scan点数。 这里速度为什么要除2.。。

1如果是计算dz，可以理解为 v/t/2 , 然后通过ns-1给它折算到dz上。

2 其次，做的是局部处理，如局部时段，则f-k偏移对应的时间 也不对。

个人理解为指定时间段的fk处理。。。

------------------------------------------

就是v*t/N ，t=trange/2,N为采样点。

但处理后么有效果。 

找到原因了，Dr. Xiong 对 别人的代码简化了，将深度信息丢失。。。用源码来调试。。

0721 通过测试把F-K使用弄清楚了

voh 速度模型 ，第一列为 速度，第二列为厚度。 这个厚度也是估算的。 在ReflexW中就是用的时间。

第一行为 第一层的速度和厚度； 后面一行是 half base。

2层模型

1层模型  （0 为half base）

 功能虽然简单，但为了方便使用，需要多次测试。

07/25 fk封装好

  做F K 偏移前后，将原来的钢筋点还原了。

1.6 滤波器设计 0727

    参考别人的代码，滤波器设计没有调用matlab的函数，自己做的 。还没有测试。

   把陷波器测试好了，通过选择b-scan 图谱中的点来选择trace，根据fft的结果来选择陷波器频率。

   最近效率太低，

2 特殊问题

2.1 app 无法添加function 了

   原来已经添加的函数也不能用了。。

函数中 多了一个end，导致function功能被屏蔽

2.2 maltab出现很多问题，不知道原因

  比如plot 的figure图不能出来，ginput 卡住，导致滤波器无法选择截止频率。

inputdlg 也不能弹出，无法获得输入。

滤波器调试 花费大量的时间都是在这样的问题上。

2.3 网页中部分汉字不能显示

   各种问题被我遇到了。

--- 网页中有部分汉字不能显示，如 有。 实际上别人是可以看到的。

3 识别算法

3.1 卷积

     杨哲把这个功能理解和落实。我把它封装起来。

    使用过程：

1） GPR信号预处理，

2） 选择测试区域，寻找阈值。目前阈值是人工确定；

3）对预处理后的数据，使用cov进行处理，并根据阈值，判断病害位置。

卷积结果，寻找阈值。

有时候 不能3D可视化，

  识别结果， 还是很不错。

  3 补充

  3.1 切片图

    在做切片图识别 发现比较难确定病害，因此 自己尝试做了切片图

    图中高亮区域为脱空

 

 

如何做不同角度的切片？ 需要建立数组直接的索引关系。。。