超声弹性成像（ultrasound elastography）代码

ISBI 2011 文章“A STRAIN-BASED ULTRASOUND ELASTOGRAPHY USING PHASE SHIFT WITH PRIOR ESTIMATES AND MESHFREE SHAPE FUNCTION”的代码实现

function displm = midfilterarea（displm）

displm1 = displm;

num_cols=size（displm2）;
num_rows=size（displm1）;

for ix = 2 : num_rows-1
    for iy = 2 : num_cols-1
        displm1（ixiy）=（displm（ixiy-1）+displm（ixiy）+displm（ixiy+1）+...
            displm（ix-1iy）+displm（ix+1iy））/5;
    end
end

for ix = 2 : num_rows-1
    displm1（ix1） = （ displm（ix1）+displm（ix2）+...
        displm（ix-11）+displm（ix+11） ）/4;
    displm1（ixnum_cols） = （ displm（ixnum_cols）+displm（ixnum_cols-1）+...
        displm（ix-1num_cols）+displm（ix+1num_cols） ）/4;
end

for iy = 2 : num_cols-1
    displm1（1iy） = （ displm（1iy）+displm（2iy）+...
        displm（1iy-1）+displm（1iy+1） ）/4;
    displm1（num_rowsiy） = （ displm（num_rowsiy）+displm（num_rows-1iy）+...
        displm（num_rowsiy-1）+displm（num_rowsiy+1） ）/4;
end

displm1（11） = （ displm（11）+displm（12）+displm（21）+displm（22） ）/4;
displm1（num_rows1） = （ displm（num_rows1）+displm（num_rows2）+...
    displm（num_rows-11）+displm（num_rows-12） ）/4;
displm1（1num_cols） = （ displm（1num_cols）+displm（1num_cols-1）+...
    displm（2num_cols）+displm（2num_cols-1） ）/4;
displm1（num_rowsnum_cols） = （ displm（num_rowsnum_cols）+...
    displm（num_rowsnum_cols-1）+displm（num_rows-1num_cols）+...
    displm（num_rows-1num_cols-1） ）/4;


displm = displm1;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1343  2010-10-14 02:46  ultrasound elastography\midfilterarea.m

     文件        531  2010-10-12 06:29  ultrasound elastography\midfilterrow.m

     文件       2062  2015-12-15 22:00  ultrasound elastography\PSPE.m

     文件    1232556  2010-09-11 20:43  ultrasound elastography\simudata_phantom_2per.mat

     目录          0  2015-12-15 21:59  ultrasound elastography

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

              1236492                    5