實(shí)驗(yàn)四-位場(chǎng)邊緣識(shí)別程序設(shè)計(jì)實(shí)驗(yàn)-文庫(kù)吧
2025-05-16 17:36 本頁(yè)面
【正文】 _asm(:,:) real,allocatable:: U(:),V(:),W(:,:) integer N_point,N_line integer m,n,m1,m2,n1,n2 real factor_m real xmin,xmax,ymin,ymax,dx,dy cmdfile=39。39。 call read_cmd(cmdfile,factor_m,input_file,output_file_vdr,output_file_thdr,output_file_asm) call read_grd(input_file,N_point,N_line,Xmin,Xmax,Ymin,Ymax) call calculate_mn(N_point,N_line,m,n,m1,m2,n1,n2,factor_m) allocate(field_re(1:m,1:n),field_im(1:m,1:n)) allocate(Px_re(1:m,1:n),Px_im(1:m,1:n),Py_re(1:m,1:n),Py_im(1:m,1:n),Pz_re(1:m,1:n),Pz_im(1:m,1:n)) allocate(field_vdr(1:m,1:n),field_thdr(1:m,1:n),field_asm(1:m,1:n)) allocate(U(1:m),V(1:n),W(1:m,1:n)) call input_grd(field_re,input_file,m1,m2,n1,n2,m,n) call expand_cos_2D(m1,m2,m,n1,n2,n,field_re,field_im) call FFT2(field_re,field_im,m,n,2) CALL cal_dxdy(xmin,xmax,ymin,ymax,N_POINT,N_LINE,dx,dy) call WAVE2D(m,n,dx,dy,U,V,W) call factor(m,n,field_re,field_im,u,v,w,px_re,px_im,py_re,py_im,pz_re,pz_im) call FFT2(px_re,px_im,m,n,1) call FFT2(py_re,py_im,m,n,1) call FFT2(pz_re,pz_im,m,n,1) call deviration(m,n,px_re,py_re,pz_re,field_vdr,field_thdr,field_asm) call OUTPUT_GRD(field_vdr,output_file_vdr,m1,m2,n1,n2,m,n,eigval,xmin,xmax,ymin,ymax) call OUTPUT_GRD(field_thdr,output_file_thdr,m1,m2,n1,n2,m,n,eigval,xmin,xmax,ymin,ymax) call OUTPUT_GRD(field_asm,output_file_asm,m1,m2,n1,n2,m,n,eigval,xmin,xmax,ymin,ymax) deallocate(field_re,field_im,px_re,px_im,py_re,py_im,pz_re,pz_im,field_vdr,field_thdr,field_asm,u,v,w) end program !**************************************************************************** ! 子程序： read_cmd ! 功能：讀取參數(shù)文件 ! 輸入?yún)?shù)說(shuō)明： ! cmdfile：參數(shù)文件名 ! 輸出參數(shù)說(shuō)明： ! input_file:觀測(cè)面位場(chǎng)數(shù)據(jù)文件 ! output_file_vdr：對(duì)場(chǎng)值作水平導(dǎo)數(shù)處理后的數(shù)據(jù)文件 ! output_file_thdr:對(duì)場(chǎng)值作垂向?qū)?shù)處理后的數(shù)據(jù)文件 ! output_file_asm:對(duì)場(chǎng)值作總導(dǎo)數(shù)處理后的數(shù)據(jù)文件 ! factor_m：擴(kuò)邊因子 !**************************************************************************** Subroutine read_cmd(cmdfile,factor_m,input_file,output_file_vdr,output_file_thdr,output_file_asm) implicit none character*80 str character*(*)cmdfile character*(*)input_file,output_file_vdr,output_file_thdr,output_file_asm real factor_m open(10,file=cmdfile,status=39。old39。) read(10,*)str,input_file read(10,*)str,output_file_vdr read(10,*)str,output_file_thdr read(10,*)str,output_file_asm read(10,*)str,factor_m close(10) end Subroutine read_cmd !*************************************************************************** ! 子程序： read_grd ! 功能：從原始觀測(cè) .grd文件中讀取相關(guān)參數(shù) ! 輸入?yún)?shù)說(shuō)明： ! filename_obser：輸入文件名 ! 輸出參數(shù)說(shuō)明： ! N_point,N_line：點(diǎn)數(shù)、線數(shù) ! x_min,x_max： x的最小值和最大值 ! y_min,y_max： y的最小值和最大值 !*************************************************************************** subroutine read_grd(input_file,N_point,N_line,Xmin,Xmax,Ymin,Ymax) implicit none character*(*)input_file integer N_point,N_line real Xmin,Xmax,Ymin,Ymax open(10,file=input_file,status=39。old39。) Read(10,*) Read(10,*)N_line,N_point Read(10,*)Xmin,Xmax Read(10,*)Ymin,Ymax Close(10) end subroutine read_grd !************************************************************************** ! 子程序： calculate_mn ! 功能：確定擴(kuò)邊數(shù)據(jù)點(diǎn)號(hào)位置 ! 輸入?yún)?shù)說(shuō)明： ! factor_m: 擴(kuò)邊比例因子（ ） ! a,b：點(diǎn)數(shù)、線數(shù) ! 輸出參數(shù)說(shuō)明： ! m1,m2: x方向?qū)嶋H數(shù)據(jù)起點(diǎn)位置和終點(diǎn)位置點(diǎn)號(hào) ! m:擴(kuò)邊后數(shù)據(jù)終點(diǎn)位置點(diǎn)號(hào) (起點(diǎn)位置點(diǎn)號(hào)為 1) ! n1,n2: y方向?qū)嶋H數(shù)據(jù)起點(diǎn)位置和終點(diǎn)位置點(diǎn)號(hào) ! n:擴(kuò)邊后數(shù)據(jù)終點(diǎn)位置點(diǎn)號(hào) (起點(diǎn)位置點(diǎn)號(hào)為 1) !************************************************************************** subroutine calculate_mn(a,b,m,n,m1,m2,n1,n2,factor_m) implicit none integer a,b,m,n,m1,m2,n1,n2 integer mtemp,mu,nu real factor_m mtemp=a DO WHILE ((mod(mtemp,2).).and.()) mtemp=mtemp/2 End do IF () THEN m=a*2 ELSE mu=int(log(float(a))/+factor_m) m=2**mu END IF m1=1+(ma)/2 m2=m1+a1 write(*,*)m,a pause mtemp=b DO WHILE ((mod(mtemp,2).).and.()) mtemp=mtemp/2 End do IF () TH
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1