Usage
FAQs
Give Feedback
Affiliates
Login
Application List
coding
code translation
subroutine get_entropy(nx,ny,image,sig,med,r,entropy) implicit none integer nx,ny,r real image(nx,ny),sig,med,entropy(nx,ny) real a20,a40,a60,a80,p(5),a(5),val parameter (a20=0.25) parameter (a40=0.52) parameter (a60=0.85) parameter (a80=1.28) integer i,j,ii,jj,ix,iy,k a(1)=a20*sig a(2)=a40*sig a(3)=a60*sig a(4)=a80*sig a(5)=100.*sig do i=1,nx do j=1,ny do k=1,5 p(k)=0. enddo do ii=i-r,i+r ix=mod(ii-1,nx)+1 do jj=j-r,j+r iy=mod(jj-1,ny)+1 val=abs(image(ix,iy)-med) k=1 do while (val.gt.a(k).and.k.le.4) k=k+1 enddo p(k)=p(k)+1. enddo enddo entropy(i,j)=0. do k=1,5 if (p(k).gt.0) entropy(i,j)=entropy(i,j)+p(k)*log(p(k)) enddo enddo enddo return end subroutine
Run