function [selx,sely] = compute_quadrant_selection(j,q, options)% compute_quadrant_selection - compute the indices for selecting subband of a wavelet transform.%%   [selx,sely] = compute_quadrant_selection(j,q, options);%%   q=1 is for phi/psi (on X / on Y)%   q=2 is for psi/phi (on X / on Y)%   q=3 is for psi/psi (on X / on Y)%%   % Suppose that MW contains a wavelet transform %   j = 4; q = 2;%   [selx,sely] = compute_quadrant_selection(j,q);%   MWs = MW(selx,sely);    % MWs contains the quadrant.%%   Depending on the transform you are using, you should %   set options.transform to %       'isotropic2d' (default): 2D standard wavelet basis.%       'quincunx': for quincunx transform (perform_quincunx_wavelet_transform).%       %%   Copyright (c) 2004 Gabriel Peyr