C++*********************************************************************
C
C CRCSE1.F
C
C **********************************************************************
C=* FROM: SPIDER - MODULAR IMAGE PROCESSING SYSTEM.   AUTHOR: J.FRANK  *
C=* Copyright (C) 1985-2005  Health Research Inc.                      *
C=*                                                                    *
C=* HEALTH RESEARCH INCORPORATED (HRI),                                *   
C=* ONE UNIVERSITY PLACE, RENSSELAER, NY 12144-3455.                   *
C=*                                                                    *
C=* Email:  spider@wadsworth.org                                       *
C=*                                                                    *
C=* This program is free software; you can redistribute it and/or      *
C=* modify it under the terms of the GNU General Public License as     *
C=* published by the Free Software Foundation; either version 2 of the *
C=* License, or (at your option) any later version.                    *
C=*                                                                    *
C=* This program is distributed in the hope that it will be useful,    *
C=* but WITHOUT ANY WARRANTY; without even the implied warranty of     *
C=* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU  *
C=* General Public License for more details.                           *
C=*                                                                    *
C=* You should have received a copy of the GNU General Public License  *
C=* along with this program; if not, write to the                      *
C=* Free Software Foundation, Inc.,                                    *
C=* 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.      *
C=*                                                                    *
C **********************************************************************
C
C
C
C IMAGE_PROCESSING_ROUTINE
C
C        1         2         3         4         5         6         7
C23456789012345678901234567890123456789012345678901234567890123456789012
C--*********************************************************************

         SUBROUTINE CRCSE1(LUN,LUN2,NSAM,NROW,IR)

         DIMENSION BUF(NSAM),SEC(IR),SNO(IR)
         
         SEC = 0.0
         SNO = 0.0
        
         DO J=1,NROW
            KJ=J-NROW/2-1
            IF(IABS(KJ).LE.IR-1)  THEN
               CALL REDLIN(LUN,BUF,NSAM,J)
               DO I=1,NSAM
                  KI=I-NSAM/2-1
                  R=SQRT(FLOAT(KJ*KJ)+FLOAT(KI*KI))+1.0
                  L=R
                  IF(L.LE.IR-1) THEN
                     XD=R-L
                     SEC(L)=SEC(L)+BUF(I)*(1.0-XD)
                     SEC(L+1)=SEC(L+1)+BUF(I)*XD
                     SNO(L)=SNO(L)+1.0-XD
                     SNO(L+1)=SNO(L+1)+XD
                  ENDIF
               ENDDO
            ENDIF
         ENDDO

         DO    I=1,IR
            SEC(I)=SEC(I)/AMAX1(1.0,SNO(I))
         ENDDO
         CALL WRTLIN(LUN2,SEC,IR,1)

         END