aboutsummaryrefslogtreecommitdiff
path: root/noao/nproto/ace/xtmaskname.x
blob: 9a55fb29a86428f93f3d81417ac2fadaa2a51976 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
#task test
#procedure test()
#char	fname[SZ_FNAME]
#begin
#	call clgstr ("fname", fname, SZ_FNAME)
#	#call xt_maskname (fname, "im1", READ_ONLY, fname, SZ_FNAME)
#	call xt_maskname (fname, "im1", NEW_IMAGE, fname, SZ_FNAME)
#	call printf ("mname = %s\n")
#	    call pargstr (fname)
#end

# MASKNAME -- Make a mask name.  This creates a FITS mask extension if
# possible, otherwise it creates a pixel list file.  To create a FITS
# extension the filename must explicitly select the FITS kernel or the
# default image type must be a FITS file.  The input and output strings
# may be the same.

procedure xt_maskname (fname, extname, mode, mname, maxchar)

char	fname[ARB]			#I File name
char	extname[ARB]			#I Default pixel mask extension name
int	mode				#I Mode
char	mname[maxchar]			#O Output mask name
int	maxchar				#I Maximum characters in mask name

int	i, fits
pointer	sp, temp

bool	streq()
int	strmatch(), stridxs(), strldxs(), strncmp()
int	envfind(), access(), imaccess()

begin
	call smark (sp)
	call salloc (temp, maxchar, TY_CHAR)

	# Determine whether to use FITS pixel mask extensions.  One may set
	# fits=NO to force use of pl even when FITS mask extensions are
	# supported.
	fits = access ("iraf$sys/imio/iki/fxf/fxfplwrite.x", 0, 0)
	if (fits == YES && envfind ("masktype", Memc[temp], maxchar) > 0) {
	    if (streq (Memc[temp], "pl"))
		fits = NO
	}
	i = strldxs ("]", fname)

	# Check for explicit .pl extension.
	if (strmatch (fname, ".pl$") > 0)
	    call strcpy (fname, mname, maxchar)

	# Check for explicit mask extension.
	else if (strmatch (fname, "type=mask") > 0)
	    call strcpy (fname, mname, maxchar)
	else if (strmatch (fname, "type\\\=mask") > 0)
	    call strcpy (fname, mname, maxchar)

	# Check for kernel section and add mask type.
	else if (i > 0) {
	    if (mode != READ_ONLY) {
		call strcpy (fname[i], Memc[temp], maxchar)
		call sprintf (mname[i], maxchar-i, ",type=mask%s")
		    call pargstr (Memc[temp])
	    }

	# Create output from rootname name.
	} else if (fits == YES) {
	    call strcpy (fname, Memc[temp], SZ_FNAME)
	    if (mode == READ_ONLY) {
		call sprintf (mname, maxchar, "%s[%s]")
		    call pargstr (Memc[temp])
		    call pargstr (extname)
	    } else {
		call sprintf (mname, maxchar, "%s[%s,type=mask]")
		    call pargstr (Memc[temp])
		    call pargstr (extname)
	    }
	} else
	    call strcat (".pl", mname, maxchar)

	# Convert to pl form if required.
	i = stridxs ("[", mname)
	if (i > 0 && mode == READ_ONLY)
	    fits = imaccess (mname, mode)
	if (fits == NO && i > 0) {
	    mname[i] = EOS
	    if (mode == NEW_IMAGE) {
		if (access (mname, 0, 0) == NO) {
		    ifnoerr (call fmkdir (mname))
			mname[i] =  '/'
		    else
			mname[i] = '.'
		} else
		    mname[i] =  '/'
	    } else {
		if (access (mname, 0, 0) == NO)
		    mname[i] = '.'
		else
		    mname[i] =  '/'
	    }
	        
	    if (strncmp (mname[i+1], "type", 4) == 0 ||
	        strncmp (mname[i+1], "append", 6) == 0 ||
	        strncmp (mname[i+1], "inherit", 7) == 0) {
		mname[i+1] = EOS
		call strcat (extname, mname, maxchar)
	    } else {
		i = stridxs (",]", mname)
		mname[i] = EOS
	    }
	    call strcat (".pl", mname, maxchar)
	}

	call sfree (sp)
end