radon车牌校正代码 matlab,基于Radon变换车牌矫正（包括垂直矫正和水平矫正）Matlab实现...

clear all

clc

close all

[inputfilename,dirname] = uigetfile('*.*');

inputfilename = [dirname, inputfilename];

im = imread(inputfilename); % For example: 'input.jpg'

grayImage = rgb2gray(im);

%%%%%

%%%%% Edge detection and edge linking....

binaryImage = edge(grayImage,'canny'); % 'Canny' edge detector

binaryImage = bwmorph(binaryImage,'thicken'); % A morphological operation for edge linking

%%%%%

%%%%% Radon transform projections along 180 degrees, from -90 to +89....

% R: Radon transform of the intensity image 'grayImage' for -90:89 degrees.

% In fact, each column of R shows the image profile along corresponding angle.

% xp: a vector containing the radial coordinates corresponding to each row of 'R'.

% Negative angles correspond to clockwise directions, while positive angles

% correspond to counterclockwise directions around the center point (up-left corner).

% R1: A 1x180 vector in which, each element is equal the maximum value of Radon transform along each angle.

% This value reflects the maximum number of pixels along each direction.

% r_max: A 1x180 vector, which includes corresponding radii of 'R1'.

theta = -90:89;

[R,xp] = radon(binaryImage,theta);

imagesc(theta,xp, R); colormap(jet);

xlabel('theta (degrees)');ylabel('x''');

title('theta方向对应的Radon变换R随着x''的变化图');

colorbar

%%%%%

[R1,r_max] = max(R);

theta_max = 90;

while(theta_max > 50 || theta_max

[R2,theta_max] = max(R1); % R2: Maximum Radon transform value over all angles.

% theta_max: Corresponding angle

R1(theta_max) = 0; % Remove element 'R2' from vector 'R1', so that other maximum values can be found.

theta_max = theta_max - 91;

end

correctedImage = imrotate(im,-theta_max); % Rotation correction

correctedImage(correctedImage == 0) = 255; % Converts black resgions to white regions

subplot(1,2,1), subimage(im)

subplot(1,2,2), subimage(correctedImage)