Giáo trình xử lý ảnh y tế Tập 4 P6

Chia sẻ: Cinny Cinny | Ngày: | Loại File: PDF | Số trang:6

Thêm vào BST

Báo xấu

94
lượt xem 5
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Các vùng mà nhận ít ánh sáng hơn thì có tính dẫn kém hơn so với các vùng nhận được nhiều ánh sáng hơn và một cảnh được chuyển thành sự phân bố độ dẫn trong vật liệu truyền dẫn. Sự phân bố của các độ dẫn này có thể chuyển thành sự phân bố điện tích bằng cách cho một dòng điện không đổi đi qua vật liệu.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Giáo trình xử lý ảnh y tế Tập 4 P6

( N  i )  ( j  k )  j N  (15.18) Cik C Nki k  i  k  i  k = j! ai j i=0 k=0  m m m ở đây C n    và C n  0 n  Cho tất cả m < n. Cân bằng các biểu thức (15.17) và (15.18) cho ta mối quan hệ giữa các hệ số bj và ai: j N  ( N i ) ( j  k ) ai b j    Cki C jNi k  i k  j k (15.19) k i  0 k 0 Nếu vector (N + 1) cột b = [b0, b1, ..., bN]T và vector (N + 1) cột a = [ a0, a1, ...,aN]T biểu diễn các hệ số chuyển đổi và hệ số gốc thì b được tính theo a qua phép nhân ma trận b = Qa (15.20) Ở đây qịj là các phần tử của ma trận Q, cho như sau: j qij   Cki C jNi k  i  k  i k  N i  j  k (15.21) k k 0 15.4 Chuyển một bộ lọc số 2-D thành bộ lọc 2-D tương tự Đầu tiên thiết kế một bộ lọc số 2-D dùng kỹ thuật mô tả trong các chương trước, chúng ta có thể chiếu bộ lọc này thành bộ lọc 2-D tương tự dùng biến đổi song tuyến tính cho trong biểu thức (15.11). Nếu hàm truyền đạt của bộ lọc số 2-D cho bởi: N N m n a zz mn 1 2 m 0 n 0 (15.22) H ( z1 , z 2 )  NN m n  b zz mn 1 2 m 0 n 0 Và hàm truyền đạt của bộ lọc 2-D tương tự được cho bởi N N   cmn s m z n m0n 0 (15.22) H ( s, z )  NN   d mn s m z n m0n  0 Các hệ số c phụ thuộc vào các hệ số a theo công thức _ _ (15.24) cn  Q a n
_ _ Ở đây c n  c 0n , c1n , c 2n ,..., c Nn T và a n  a 0n , a1n , a 2n ,..., a Nn T , và tương tự cho các hệ số d và a. Q là một ma trận (N + 1)  (N + 1) toàn bộ các phần tử cho bởi biểu thức (15.21). Chu kỳ lấy mẫu T dùng cho biến đổi song tuyến tính được chọn bằng chu kỳ quét dòng của tín hiệu truyền hình, và bằng 63.555 micro giây cho tiêu chuẩn U.S. NTSC. Lý do của việc này là làm cho bộ lọc tuần hoàn đối xứng ở trong miền chữ nhật bao bởi  a = {-/T,/T} và  u = {-/T, /T}, ở đây a là tần số tương tự theo chiều ngang và d là tần số số theo chiều dọc. Chương trình C sau đây sẽ chiếu hàm truyền đạt IIR 2-D số vào hàm truyền đạt tương tự. Chương trình 15.1 “ANA2DIIR.C” Chuyển bộ lọc 2-D số vào bộ lọc 2-D tương tự. /* This program is for mapping a 2-D IIR filter to a 2-D analog filter. It also provides you with data files for Plotting magnitude and phase response. */ #define T 63.555e-6 #define pi 3.14159 #define Tt (T) #define M2 (M) #define wa_start (-pi/T) #include #include #include #include #include #include void Transf_matrix(float *,int); main() { float **a,**b,**temp,*Q,c,s,NR,NI,DR,DI,D,ph,wa,wd,dwa,dwd; float wame,wamo,was,sum1,sum2,t1,t2,ph1,ph2,**phase; int i,j,m,n,sign,N,M,k; FILE *fptr,*fptro1,*fptro2; char file_name[14]; clrscr();
printf("\nEnter file name containing IIR filter coefficients-->"); scanf("%s",file_name); fptr=fopen(file_name,"r"); if(fptr==NULL) { printf("\nThis file does not exist."); exit(1); } fscanf(fptr,"%d",&N); a=(float **)malloc((N+1)*sizeof(float *)); for(i=0;i
temp[i][n]=0.0; for(j=0;j
printf("\nEnter file name to store phase response-->"); scanf("%s",file_name); fptro2=fopen(file_name,"w+"); printf("\nEnter the number of points to be generated (e.g.32x32)-->"); scanf("%dx%d",&M,&M); fprintf(fptro1,"%d %d\n",M2,M2); fprintf(fptro2,"%d %d\n",M2,M2); fprintf(fptro1,"%e",0.0); fprintf(fptro2,"%e",0.0); wd=- pi; dwd=2.0*pi/M; for(i=0;i
t1=sign*a[m][n]*wame; t2=sign*b[m][n]*wame; sign*=-1; wame*=was; sum1+=t1; sum2+=t2; } NR+=sum1*c; NI+=sum1*s; DR+=sum2*c; DI+=sum2*s; } /*n-loop.*/ for(n=0;n