弹性变换和镜像变换

    最近在看U-Net处理医学图像时,用到了弹性变换等,因为变异组织本身就是弹性变换的,可以通过这种方法进行数据扩增。

    对于图像而言,常用的增加训练样本的方法主要有对图像进行旋转、移位等仿射变换,也可以使用镜像变换等等。

(1)弹性变换

    这里介绍弹性变换(Elastic Distortion)。对原图像进行弹性变换的操作扩充样本以后,对于手写体数字的识别效果有明显的提升。

    首先来详细介绍一下算法流程:

        (1)首先需要对图像中的每个像素点(x,y)产生两个-1~1之间的随机数,Δx(x,y)和Δy(x,y),分别表示该像素点的x方向和y方向的移动距离;

        (2)生成一个以0为均值,以σ为标准差的大小为n*n的高斯核k_nn,并用前面的随机数与之做卷积,并将结果作用于原图像。

    作者在这里提出σ的大小与弹性变换的处理结果息息相关,如果σ过小,则生成的结果类似与对图像每个像素进行随机移动,而如果σ过大,则生成的结果与原图基本类似。不同的σ效果如下

n = 5,σ=4,8,16时


n = 21,σ= 4,8,16


n = 105,σ = 4,8,16时

可以看出来,只有在n足够大(与要处理的图像相比),且σ大小合适时才能够到合适的扭曲图像,如图所示,这里n = 105且σ=8时比较合适。

    部分参考:https://blog.csdn.net/lhanchao/article/details/54234490

(2)其他变换

    沿坐标轴缩放

    沿任意方向缩放

    正交投影

    镜像

    切变:切变是坐标系的变换,非均匀的拉伸。切变时候,角度变化,但是面积或体积不变。也可以理解为坐标轴间的角度变化,造成的扭曲。

    以上变换均属于矩阵变换。

(3)图像变换的一些代码

        x = 1150

        y = 15

        w = 150

        h = 30

        region = im.crop((x, y, x+w, y+h))

        region#图像裁剪 


        #图像resize

        import cv2

        im1 = cv2.imread('1.jpg')

        im2 = cv2.resize(im1, (1024,1024), interpolation=cv2.INTER_CUBIC)

        cv2.imwrite('3.jpg', im2)


        #图像批量操作

        import os

        fpath = "D:/Data/CAG/No"

        dpath = "D:/Data/CAG/1"

        for filename in os.listdir(fpath):   

        im1 = cv2.imread(fpath+"/"+filename)

        im2 = cv2.resize(im1, (1024,1024), interpolation=cv2.INTER_AREA)

        filename.replace('.jpg','.bmp')

        cv2.imwrite(dpath+'/'+filename, im2)

        print("1")


        #水平翻转

        import cv2

        import os

        fpath = "D:/Data/CAG/train1"

        dpath = "D:/Data/CAG/h_flip"

        for filename in os.listdir(fpath):   

            im1 = cv2.imread(fpath+"/"+filename)

            h_flip = cv2.flip(im1, 1)

            strlist = filename.split('.')

            m = int(strlist[-2])+600

            filename1 = strlist[0]+'.'+str(m)+'.bmp'

             cv2.imwrite(dpath+'/'+filename1, h_flip)

        print("1")


#垂直翻转

import cv2

import os

fpath = "D:/Data/CAG/train1"

dpath = "D:/Data/CAG/v_flip"

for filename in os.listdir(fpath):   

    im1 = cv2.imread(fpath+"/"+filename)

    h_flip = cv2.flip(im1, 0)

    strlist = filename.split('.')

    m = int(strlist[-2])+1200

    filename1 = strlist[0]+'.'+str(m)+'.bmp'

    cv2.imwrite(dpath+'/'+filename1, h_flip)

print("1")


#水平垂直翻转

import cv2

import os

fpath = "D:/Data/CAG/train1"

dpath = "D:/Data/CAG/hv_flip"

for filename in os.listdir(fpath):   

    im1 = cv2.imread(fpath+"/"+filename)

    h_flip = cv2.flip(im1, -1)

    strlist = filename.split('.')

    m = int(strlist[-2])+1800

    filename1 = strlist[0]+'.'+str(m)+'.bmp'

    cv2.imwrite(dpath+'/'+filename1, h_flip)

print("1")

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 158,560评论 4 361
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,104评论 1 291
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 108,297评论 0 243
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 43,869评论 0 204
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,275评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,563评论 1 216
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,833评论 2 312
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,543评论 0 197
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,245评论 1 241
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,512评论 2 244
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,011评论 1 258
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,359评论 2 253
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,006评论 3 235
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,062评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,825评论 0 194
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,590评论 2 273
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,501评论 2 268

推荐阅读更多精彩内容

  • 深入理解傅里叶变换Mar 12, 2017 这原本是我在知乎上对傅立叶变换、拉普拉斯变换、Z变换的联系?为什么要进...
    价值趋势技术派阅读 5,545评论 2 2
  • 直方图变换 灰度变换 点运算 几何变换 直方图变换 1.灰度直方图 灰度直方图:数字图像中每一灰度级像素出现的频次...
    hyfine阅读 4,455评论 0 0
  • 本章介绍了基于elastix的基本配准概念。 更高级的配准主题将在第6章中讨论。图像配准是医学影像领域的重要工具。...
    peterpan_hai阅读 8,692评论 1 10
  • 1.今日事今日畢,迴向給自己,建立對自己的信任感。為自己點贊。 收拾新家;給新客戶送國學機;老客戶問題國學機直接換...
    李豪亮阅读 144评论 0 0
  • 01 在百度里搜索无趣有以下几种意思: 没意思,没趣;缺乏悠闲、温雅、风韵、活泼、逼真、兴趣或风趣。 百度百科的意...
    辰苓阅读 759评论 0 0