C PROGRAM MEDIAN.FT C ---------------- C C PETER LEMKIN C IMAGE PROCESSING UNIT, DCBD C NATIONAL CANCER INSTITUTE C NATIONAL INSTITUTES OF HEALTH C 9000 ROCKVILLE PIKE C BETHESDA, MD. 20014 C C JAN 16, 1977 C JAN 8, 1977 C JAN 7, 1977 C DEC 14, 1976 C C C C INTRODUCTION C ------------ C MEDIAN.FT IS A CHAINED PROGRAM USED WITH BMON2 C IT RECEIVES ITS ARGUMENTS FROM THE CD AREA AND THE IBM1,IHGH1 C IBM2,IHGH2, JBM, JHGH VARIABLES IN COMMON. COMMON IS RESTORED C FIRST BEFORE THE FUNCTION (TO BE INSERTED INTO THE BODY) IS C EVALUATED. AFTER THE FUNCTION IS PERFORMED, COMMON C IS SAVED AND BMON2 IS CHAINED BACK TO. C ANY COMPUTATIONS ARE DONE, THE ARGUMENTS ARE THEN CHECKED C C MEDIAN COMPUTES THE MEDIAN OF EACH 3X3 NEIGHBORHOOD of AND C STORES THAT VALUE IN . C OPDEFS C ------ S OPDEF TADI 1400 S OPDEF DCAI 3400 C S OPDEF DISP1 6435 S OPDEF DISP2 6436 C C S OPDEF SWAB 7431 S OPDEF SWBA 7447 C S OPDEF MUY 7405 S OPDEF DVI 7407 C S OPDEF DISP2 6436 S OPDEF DISP1 6435 C S OPDEF MQA 7501 S OPDEF MQL 7421 S OPDEF KRS 6034 S OPDEF BSW 7002 C C S OPDEF HPL 6360 S OPDEF HSL 6361 S OPDEF VPL 6362 S OPDEF VSL 6363 C [1] INITIALIZATION WRITE(1,995) 995 FORMAT('MEDIAN 1/16/77 - 4:01PM') ET=TIMER(0) C C C C [2] VERIFY BM SPECS S JMS CKIN S JMS CKOUT C@@S JMS CKIN2 C C C C [3] PROCESS IMAGE DO 200 IY1=KY1,KY2 C S CPAGE 3 S JMS TTYCTL S JMP \998 /DONE C IY=IY1-1 S TAD \IY S DISP2 C C C [3.1] PROCESS A LINE DO 200 JJX1=KX1,KX2 IX=JJX1-1 S TAD \IX S DISP1 C C C C [3.1.1] GET BMI NEIGHBORHOOD MEM=IBM1 IBYTE=IHGH1 CALL GETI1 C C C [3.1.2] COMPUTE: MEDIAN DO 399 MAXTIMES=1,8 C C COMPUTE THE POINTERS S CLA CMA S TAD PI10 /POINTER S DCA 14 /POINTER TO I10(1) S TAD 14 S DCA 15 C S TAD 14 /SET UP POINTER TO I10(2) S IAC S DCA 16 S TAD 16 S DCA 17 C C C DO BUBBLE SORT NOACTION=1 C DO 398 KWC=1,8 C COMPUTE: IX1=I10(KWC) S CPAGE 10 S 6211 S TADI 14 S DCA \IX1 C C COMPUTE: IX2=I10(KWC+1) S TADI 16 S DCA \IX2 C C IF I1(KWC) > I1(KWC+1) C THEN I1(KWC) <==> I1(KWC+1), NOACTION=0; IF(IX1-IX2)302,303,303 C C THEN 302 CONTINUE C COMPUTE: I10(KWC)=IX2 S CPAGE 4 S TAD \IX2 S DCAI 15 C C COMPUTE: I10(KWC+1)=IX1 S CPAGE 4 S TAD \IX1 S DCAI 17 C NOACTION=0 GOTO 398 C C ELSE JUST INCREMENT THE POINTERS S\303, INC 15 S INC 17 S CLA C 398 CONTINUE C C TEST IF DONE, IF NO SWAPS IF(NOACTION)350,399,350 C 399 CONTINUE C C C GET THE MEDIAN FROM THE SORTED LIST 350 IZ=I10(5) C C C C [3.1.3] IZ==>BMJ MEM=JBM IBYTE=JHGH CALL PACK2D 200 CONTINUE C C C C C [999] RETURN GOTO 998 999 WRITE(1,996) 996 FORMAT('BAD BM SPECIFICATION!') 998 ET=TIMER(1) CALL CHAIN('BMON2') C*************************************************** C *PROCEDURE C K O U T C****************************************************** C C C CHECK WHETHER THE OUTPUT BM SPEC IS LEGAL ELSE GOTO 999. S CPAGE 3 SRCKOUT, JMP I CKOUT S CKOUT, 0 /ENTRY C C [1] TEST IF KOUTFILE="BM" S TAD \KOUTFILE S CIA S TAD BMTEXT S SZA CLA S JMP \999 /FAILED C C [2] TEST IF (KOUTFILE(2) LAND '7700)=DIGIT S TAD \KOUTFILE# S BSW S AND (0077 S TAD (-60 /"0" S SPA S JMP \999 /NO, FAILED S TAD (-D8 /"7" TEST S SMA CLA S JMP \999 /FAILED S JMP RCKOUT /OK. C C C*************************************************** C *PROCEDURE C K I N C****************************************************** C C C CHECK WHETHER THE INPUT BM SPEC IS LEGAL ELSE GOTO 999. S CPAGE 3 SRCKIN, JMP I CKIN S CKIN, 0 /ENTRY C C [1] TEST IF BMI1="BM" S TAD \SFILE S CIA S TAD BMTEXT S SZA CLA S JMP \999 /FAILED C C [2] TEST IF (BMI1(2) LAND '7700)=DIGIT S TAD \SFILE# S BSW S AND (0077 S TAD (-60 /"0" S SPA S JMP \999 /NO, FAILED S TAD (-D8 /"7" TEST S SMA CLA S JMP \999 /FAILED S JMP RCKIN /OK. C C C*************************************************** C *PROCEDURE C K I N 2 C****************************************************** C C C CHECK WHETHER THE INPUT BM SPEC IS LEGAL ELSE GOTO 999. S CPAGE 3 SRCKIN2, JMP I CKIN2 S CKIN2, 0 /ENTRY C C [1] CHECK FIRST IBM1 SPEC S JMS CKIN C C C [2] TEST IF BMI2="BM" S TAD \SEXT S CIA S TAD BMTEXT S SZA CLA S JMP \999 /FAILED C C [3] TEST IF (BMI2((6) LAND '7700)=DIGIT S TAD \SEXT# S BSW S AND (0077 S TAD (-60 /"0" S SPA S JMP \999 /NO, FAILED S TAD (-D8 /"7" TEST S SMA CLA S JMP \999 /FAILED S JMP RCKIN2 /OK. C C C ************************************************************ C SUBROUTINE: T T Y C T L (INTERNAL) C ************************************************************ C S CPAGE 3 S RTTYC, JMP I TTYCTL S TTYCTL, 0000 /ENTRY C S KSF /ANYTHING TYPED S JMP NORMAL /NO, RETURN NORMALLY S KRB /GET TYPED CHARACTER S AND (0177 /TAKE CARE OF PARITY PROBLEMS S TAD (-17 /TEST FOR CTRL/O S SNA /SKIP IF NOT CTRL/O S JMP RTTYC /ABORT CALLING ROUTINE (ERROR RETURN) S TAD (-4 /TEST FOR CTRL/S [-17-4=-23(OCTAL)] S SZA CLA /SKIP IF CTRL/S S JMP NORMAL /NOT CTRL/O OR CTRL/S SO RETURN NORMALLY C S SLEEP,KSF /WAIT FOR CTRL/Q S JMP SLEEP /KEEP WAITING S KRB /READ CHARACTER S AND (0177 S TAD (-17 /IS IT A CTRL/O? S SNA /SKIP IF NOT S JMP RTTYC /YES, ABORT S TAD (-2 /TEST FOR CTRL/Q (-17-2=-21 OCTAL) S SZA CLA /SKIP IF SO S JMP SLEEP /NOPE, KEEP SLEEPING C S NORMAL,INC TTYCTL /INCREMENT RETURN ADDRESS FOR NORMAL RETURN S CLA /SAFETY VALVE S JMP RTTYC /RETURN C************** P A R A M E T E R S ************* S BMTEXT, TEXT /BM/ S PI10, \I10 END