Mangepoint source Fortran source

  • Mangeage des petits points dnas une image 
  1.      program mangepoint
  2.      implicit none
  3.      integer source(6000,8000), resul(6000,8000)
  4.      integer ic, il, sigma, iblancg, iblancd, lgr, nbvois, sigmatot
  5.      integer nblig, nbcol, posit, icc, ill, distance2, sigma2
  6.      character*90 ligne
  7.      character*1 bout
  8.      character*2 sigmachar
  9.      character*30 nomfichier, nominput, nomoutput
  1.      sigma = 0
  2.      call getarg(1, nomfichier)
  3.      call getarg(2, sigmachar)
  4.      read (sigmachar, *) sigma
  5.      if (sigma.eq.0) then
  6.        print *, 'Sigma = 0 !'
  7.        stop
  8.      endif
  9.      posit = len_trim(nomfichier) + 1
  10.      nominput = nomfichier(1:posit-1)//'.pbm'
  11.      call system('copy '//nominput//' old'//nominput)
  12.      
  13.      nomoutput= nominput
  14.      sigmatot = (2*sigma + 1)*(2*sigma + 1) ! 
  15.      sigma2 = sigma*(sigma + 1)
  16.      print *, 'sigmatot = ', sigmatot
  1.      call system('convert -compress none '//nominput//nominput)
  2.      open (unit=8, file=nominput)
  3.      open (unit=9, file=nomoutput)
  4.      read (8, '(a)') ligne    ! P1
  5.      write (9, '(a)') ligne
  6.      read (8, '(a)') ligne  
  7.      do while (index(ligne, '#').gt.0) 
  8.        read (8, '(a)') ligne  
  9.      enddo
  10.      write (9, '(a)') ligne
  11.      read (ligne, *) nbcol, nblig
  12.      
  13.      print *, nblig, nbcol
  14.      ic = 1; il =1
  15.      do while ((il.lt.nblig).and.(ic.lt.nbcol))
  16.        read (8, '(a)', end=200) ligne
  17.        lgr = len_trim(ligne)
  18.        iblancg = 1
  19.        iblancd = 0
  20.        do while (iblancd.lt.lgr)
  21.          iblancd = index(ligne, ' ') 
  22.          bout = ligne(iblancg:iblancd - 1) 
  23.          read (bout,*) source(ic,il)
  24.          ic = ic + 1
  25.          ic = mod(ic, nbcol)
  26.          if (ic.eq.0) il = il + 1
  27.          ligne(iblancd:iblancd) = '*' 
  28.          iblancg = iblancd + 1
  29.        enddo
  30.      enddo
  31.      print *, 'Lecture OK'
  32.      goto 300

200 print *, 'Fin de fichier', ic, il 300 do ic = sigma, nbcol - sigma 

  1.        do il = sigma, nblig - sigma
  2.          nbvois = 0
  3.          resul(ic, il) = source(ic, il)
  4.          do icc = ic - sigma, ic + sigma
  5.            do ill = il - sigma, il + sigma
  6.              distance2 = (icc - ic)*(icc - ic) + (ill - il)*(ill - il)
  7.              if ((distance2.le.sigma2).and.(source(icc,
  8.     +        ill).eq.1)) nbvois = nbvois + 1
  9.            enddo
  10.          enddo
  11.          if (nbvois.lt..125*sigmatot) resul(ic, il) = 0
  12.        enddo
  13.      enddo
  14.      print *, 'Fin du traitement'
  1.      do il = 1, nblig 
  2.        ic = 1     
  3.        do while (ic.le.nbcol)
  4.          write (9,*) (resul(icc, il), icc=ic, min(ic+19, nbcol))
  5.          ic = ic + 20
  6.        enddo
  7.      enddo
  8.      close (unit=8)
  9.      close (unit=9)
  1.      end
Retrieved from "http://imslp.org/index.php?title=Mangepoint_source_Fortran_source&oldid=1946797"