Linear optics parameters are often estimated using turn-by-turn (TbT) data. In-plane beta functions, which are the most common measurement objectives, can be estimated from the amplitudes or phases of TbT data, providing an overall characterization of the linear lattice. In addition to estimating uncoupled Twiss parameters, TbT data can also be utilized for characterizing coupled linear motion. We investigate several methods for constructing a full normalization matrix at each beam position monitor (BPM). BPMs provide information on beam centroid transverse coordinates, but direct observation of transverse momenta is not available. To estimate momenta, one can use a pair of BPMs or fit data obtained from several BPMs. By using both coordinates and momenta, it is possible to fit the one-turn matrix (or its power) at each BPM, allowing for the computation of coupled Twiss parameters. Another approach involves the minimization of side peak amplitudes in the spectrum of normalized complex coordinates. Similarly, the normalization matrix can be estimated by fitting linear coupled invariants. In this study, we derive and utilize a special form of the normalization matrix, which remains non-singular in the zero coupling limit. Coupled Twiss parameters can be obtained from the normalization matrix. The paper presents the results of applying and comparing these methods to both modeled and measured VEPP-4M TbT data, demonstrating effective estimation of coupled Twiss parameters.