图像畸变矫正Matlab实现，可以直接运行

图像畸变的矫正，matlab代码，无需配置环境，里面带实验原图和效果图

%%%%%%质心求取程序

%
function [plot_xplot_y]=calculate_zhixin（I）
 I=imread（‘原图.png‘）;%读取原图像  
figure（1）;imshow（I）;%显示原图像  
  
I_gray=rgb2gray（I）;%原图像变为灰度图像  



I=imcomplement（I_gray）;
T=imtophat（Istrel（‘disk‘10））;
I_bw=~im2bw（Tgraythresh（T））;
figure（2）imshow（I_bw）





%level=graythresh（I_gray）;%计算图像I_gray的全局阈值，level为标准化灰度值，其范围为[0 1]  
[heightwidth]=size（I_gray）;%计算灰度图像的长宽  
%I_bw=im2bw（I_graylevel）;%im2bw使用阈值level将灰度图像转换为二值图像．  
%figure（2）;imshow（I_bw）;%显示二值图像（背景为白色）  
  


for i=1:height %%循环中进行反色  
    for j=1:width     
        if I_bw（ij）==1        
            I_bw（ij）=0;    
        else I_bw（ij）=1;   
        end  
    end  
end  
figure（3）;imshow（I_bw）;%显示取反后的二值图像（背景为黑色）  


  
[Lnum]=bwlabel（I_bw8）;%bwlabel标注二值图像I_bw中的目标物体，返回标识矩阵Ｌ和I_bw中目标物体的数量num，８表示连通数．  
plot_x=zeros（num1）;%%zeros（mn）产生m×n的全0矩阵．用于记录质心位置的横坐标  
plot_y=zeros（num1）;%zeros（mn）产生m×n的全0矩阵．用于记录质心位置的纵坐标  
  
for k=1:num  %%num个区域依次统计质心位置      
    sum_x=0;    s

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-04-09 14:25  图像校正\
     文件        1876  2015-03-27 17:23  图像校正\calculate_zhixin.asv
     文件        1874  2015-03-31 15:23  图像校正\calculate_zhixin.m
     文件        1782  2015-03-26 21:01  图像校正\calcultate_zhixin.asv
     文件        1300  2015-03-26 21:01  图像校正\fengezuobiao.asv
     文件        1402  2015-03-26 21:01  图像校正\fengezuobiao.m
     文件        2001  2015-03-27 17:23  图像校正\image_huifu.asv
     文件        1961  2015-03-26 22:48  图像校正\image_huifu.m
     文件        4860  2015-03-26 21:01  图像校正\main.asv
     文件        5083  2015-03-31 15:24  图像校正\main.m
     文件       23183  2015-03-26 21:01  图像校正\原图.png
     文件       26913  2015-03-26 21:01  图像校正\畸变图.png