• 毛毛爱陈晨
    了解作者
  • matlab
    开发工具
  • 171KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 59
    下载次数
  • 2011-09-13 21:35
    上传日期
医学图像配准中的刚性配准两幅图像,希望大家用得到。
rigid-registration.rar
  • 刚性配准(同尺寸)
  • image.mat
    307.9KB
  • PV.m
    3.5KB
  • Fusion.m
    2KB
  • OneDimSearch.m
    682B
  • GLPF.m
    538B
  • Fx.m
    121B
  • Powell.asv
    1.2KB
  • Registration.m
    1.8KB
  • Fusion.asv
    2.9KB
  • PV.asv
    3.3KB
  • Registration.asv
    1.7KB
  • OneDimSearch.asv
    680B
  • Powell.m
    1.5KB
内容介绍
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %几何变换及互信息计算 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [mi]=PV(x,y,ang,I,J) a=double(I); b=double(J); [M,N]=size(a); hab=zeros(256,256); ha=zeros(1,256); hb=zeros(1,256); if max(max(a))~=min(min(a)) a=(a-min(min(a)))/(max(max(a))-min(min(a))); else a=zeros(M,N); end if max(max(b))~=min(min(b)) b=(b-min(min(b)))/(max(max(b))-min(min(b))); else b=zeros(M,N); end a=double(int16(a*255))+1; b=double(int16(b*255))+1; [width,height]=size(b); u=(width-1)/2; v=(height-1)/2; rad=pi/180*ang; t1=[1 0 0;0 1 0;x y 1]; t2=[1 0 0;0 1 0;-u -v 1]; t3=[cos(rad) -sin(rad) 0;sin(rad) cos(rad) 0;0 0 1]; t4=[1 0 0;0 1 0;u v 1]; T=t2*t3*t4*t1; tform=maketform('affine',T); coordinate_x=zeros(width,height); coordinate_y=zeros(width,height); for i=1:width for j=1:height coordinate_x(i,j)=i; coordinate_y(i,j)=j; end end [w z]=tforminv(tform,coordinate_x,coordinate_y); f_New=uint8(zeros(M,N)); % for i2=1:M % for j2=1:N % source_x=w(i2,j2); % source_y=z(i2,j2); % if (source_x>=width-1||source_y>=height-1||double(uint16(source_x))<=0||double(uint16(source_y))<=0) % f_New(i2,j2)=0; % else % if ((source_x/double(uint16(source_x))==1.0)&(source_y/double(uint16(source_y))==1.0)) % f_New(i2,j2)=b(int16(source_x),int16(source_y)); % else %双线性插值法 % p=double(uint16(source_x)); % q=double(uint16(source_y)); % x11=double(b(p,q)); % x12=double(b(p,q)); % x21=double(b(p+1,q)); % x22=double(b(p+1,q+1)); % f_New(i2,j2)=uint8((q+1-source_y)*((source_x-p)*x21+(p+1-source_x)*x11)+(source_y-q)*((source_x-p)*x22+(p+1-source_x)*x12)); % end % end % end % end % f_New=double(f_New); % a=(a-min(min(a)))/(max(max(a))-min(min(a))); % f_New=(f_New-min(min(f_New)))/(max(max(f_New))-min(min(f_New))); % a=double(int16(a*255))+1; % f_New=double(int16(f_New*255))+1; % for i1=1:M % for j1=1:N % index_x=a(i1,j1); % index_y=f_New(i1,j1); % hab(index_x,index_y)=hab(index_x,index_y)+1; % end % end for i=1:width for j=1:height source_x=w(i,j); source_y=z(i,j); if (source_x>width-1||source_y>height-1||double(uint16(source_x))<=1||double(uint16(source_y))<=1) hab(a(1,1),a(1,1))=hab(a(1,1),a(1,1))+1; else m=fix(source_x); n=fix(source_y); index_b=b(i,j); index_a0=a(m,n); index_a1=a(m+1,n); index_a2=a(m,n+1); index_a3=a(m+1,n+1); dx=source_x-m; dy=source_y-n; hab(index_a0,index_b)=hab(index_a0,index_b)+(1-dx)*(1-dy); hab(index_a1,index_b)=hab(index_a1,index_b)+dx*(1-dy); hab(index_a2,index_b)=hab(index_a2,index_b)+(1-dx)*dy; hab(index_a3,index_b)=hab(index_a3,index_b)+dx*dy; end end end habsum=sum(sum(hab)); index=find(hab~=0); pab=hab/habsum; Hab=sum(sum(-pab(index).*log2(pab(index)))); pa=sum(pab'); index=find(pa~=0); Ha=sum(sum(-pa(index).*log2(pa(index)))); pb=sum(pab); index=find(pb~=0); Hb=sum(sum(-pb(index).*log2(pb(index)))); mi=Ha+Hb-Hab end
评论
    相关推荐
    • 图像配准.zip
      图像配准 mr图像的应用 方便观察病变的位置结构
    • VC图像变换 图像配准 图像分割图像编码等图片处理程序
      VC++功能较多的图像处理程序,实现的功能有:图像变换 图像配准 图像分割、图像编码、图像增强、图像复原、图像特征码提娶图像识别、运动检测、图像显示等,特别是源码里面的一些VC的图像处理类,你可以用到你自己的...
    • 图像配准
      图像配准是许多医学图像分析的基础。 精确图像配准的主要障碍是图像外观的变化。 最近,使用深度神经网络的基于深度学习的注册方法(DLR)的计算效率比传统的基于优化的注册方法(OR)高几个数量级。 但是,DLR的一...
    • sift图像配准
      一个很好的图像配准算法,能够使用同源和异源图像的配准。配准精度高,速度快,具有仿射不变性,旋转不变性,尺度不变性等
    • 图像处理配准图像处理配准
      图像处理配准,图像处理配准,图像处理配准,图像处理配准,图像处理配准,图像处理配准
    • 图像配准 图像拼接技术 sift
      图像配准 图像拼接技术 sift 实现任意图像的配准和拼接 利用c++实现
    • 图像配准技术
      图像处理——图像配准技术 可以应用于图像配准方面 matlab代码
    • 图像配准算法
      LKOFlow 是高斯金字塔图像配准算法,Patrick Vandewalle 含有多种配准算法
    • matlab surf图像配准 图像检索
      matlab surf图像配准 图像检索 图像处理 模式识别 我本人的毕业设计 就是利用的surf来提取的特征,并进行相似度比较后,然后进行检索的。此程序对想利用surf的同学可以得到很大的帮助。surf比sift具有更快的速度,更...
    • matlab surf图像配准 图像检索
      matlab surf图像配准 图像检索 图像处理 模式识别 我本人的毕业设计 就是利用的surf来提取的特征,并进行相似度比较后,然后进行检索的。此程序对想利用surf的同学可以得到很大的帮助。surf比sift具有更快的速度,更...