rec2.m ( File view )

  • By zhengzou 2014-05-09
  • View(s):31
  • Download(s):12
  • Point(s): 2
			

clc;close all; clear all;

f = imread('C:\Documents and Settings\Administrator\桌面\a\060.bmp');
data = double(f);
[M N] = size(f);
figure;imshow(f);

x = [106 117 131 148 168 190 215 247 282 320 360 400 437 472 502 527 549 567 583]-360;
y = 288 - [121 112 100 86   73  60  48  36  27  22  20  22  29  39  50 62  75  86  98];
% x =     [77  85  92  104 118 137 162 195 239 295 359 421 477 520 551 576 595 609 621]-360;
% y = 288-[224 219 214 206 198 189 178 167 155 146 144 146 155 164 174 184 192 200 205];
[MM NN] = size(x);

[R,A,B] = circ(x,y,NN); 


z = R/2-110;
%z=90;
data_then = zeros(M,N);

for i = 1:M
   for j = 361:N
        x = j - 360;  %坐标转换
        y = 288 - i;
       
         u = R*x*( (z*z+x*x+y*y).^0.5-z )/(x*x+y*y);
        v = R*y*( (z*z+x*x+y*y).^0.5-z )/(x*x+y*y);
        %u=(x-x0)*(R.^2-(y-y0)^2)^0.5/R+x0;
       % v=y;
        
        i_index = 288 - v;
        j_index = 360 + u;
        
        i_index = uint32(i_index);
        j_index = uint32(j_index);
        
        if(i_index<1 || i_index>576 ||j_index<1|| j_index>720)
            data_then(i,j) = 0;
        else
            data_then(i,j) = data(i_index,j_index);
        end
   end
end
 

figure;imshow(uint8(data_then));

    
 
 
%x1=47-360;x2=601-360;x3=35-360;
%y1=288-290;y2=288-277;y3=288-70;

x1=361-360;x2=688-360;
...
...
(Please download the complete source code to view)
			
...
Expand> <Close

Want complete source code? Download it here

Point(s): 2

Download
0 lines left, continue to read
Sponsored links

File list

Tips: You can preview the content of files by clicking file names^_^
Name Size Date
060-1.fig270.39 kB01-08-12|10:26
060.bmp406.05 kB18-07-12|09:51
060.fig261.75 kB01-08-12|09:39
080.bmp406.05 kB18-07-12|09:51
090l.bmp406.05 kB18-07-12|09:51
090r.bmp406.05 kB18-07-12|09:51
090r.fig277.17 kB01-08-12|09:34
90.bmp406.05 kB20-07-12|16:14
circ.m533.00 B30-07-12|15:32
correct_image_061g.bmp3.82 MB01-08-12|19:52
fitellipse.m3.13 kB28-07-12|22:38
markpoint.m621.00 B01-08-12|09:57
rec.asv3.29 kB04-08-12|10:40
rec.m3.32 kB04-08-12|11:16
rec2.asv2.02 kB01-08-12|15:54
rec2.m2.02 kB01-08-12|15:55
shiyan5.m3.69 kB27-07-12|15:30
test.asv1.59 kB01-08-12|11:04
test.m1.63 kB01-08-12|16:43
test2.asv3.92 kB01-08-12|10:09
Thumbs.db36.00 kB18-07-12|09:51
trapezial_correction-wrong.m1.50 kB02-08-12|16:59
trapezial_correction.asv3.23 kB02-08-12|17:23
trapezial_correction.m1.66 kB04-08-12|10:09
Untitled.m1.45 kB26-07-12|20:43
Untitled3.m3.98 kB25-07-12|15:44
Untitled5.m3.75 kB23-07-12|08:58
Untitled6.m3.82 kB23-07-12|15:35
截图.jpg125.80 kB29-07-12|20:25
校正0.00 B09-04-13|09:07
...
Sponsored links

rec2.m (3.14 MB)

Need 2 point
Your Point(s)

Your Point isn't enough.

Get point immediately by PayPal

More(Debit card / Credit card / PayPal Credit / Online Banking)

Submit your source codes. Get more point

LOGIN

Don't have an account? Register now
Need any help?
Mail to: support@codeforge.com

切换到中文版?

CodeForge Chinese Version
CodeForge English Version

Where are you going?

^_^"Oops ...

Sorry!This guy is mysterious, its blog hasn't been opened, try another, please!
OK

Warm tip!

CodeForge to FavoriteFavorite by Ctrl+D