Extraplanar emission in isolated edge-on late-type galaxies. I. The H$\alpha$ distribution versus to the old and young stellar discs

Isolated galaxies are the ideal reference sample to study the galaxy structure minimising potential environmental effects. We selected a complete sample of 14 nearby, late-type, highly inclined ($i\geq80^{\circ}$), isolated galaxies from the Catalogue of Isolated Galaxies (CIG) which offers a vertical view of their disc structure. We aim to study extraplanar Diffuse Ionized Gas (eDIG) by comparing the old and young disc components traced by near-infrared (NIR) and Ultraviolet (UV) imaging with the H$\alpha$ emission structure. We obtained H$\alpha$ monochromatic maps from the Fabry-Perot (FP) interferometry, while the old and young discs structures are obtained from the photometric analysis of the 2MASS K$_{s}$-band, and GALEX NUV and FUV images, thereby identifying the stellar disc and whether the eDIG is present. The H$\alpha$ morphology is peculiar in CIG 71, CIG 183, CIG 593 showing clear asymmetries. In general, geometric parameters (isophotal position angle, peak light distribution, inclination) measured from H$\alpha$, UV and NIR show minimal differences (e.g. $\Delta i\leq\pm$10$^{\circ}$), suggesting that interaction does not play a significant role in shaping the morphology, as expected in isolated galaxies. From H$\alpha$ maps, the eDIG was detected vertically in 11 out of 14 galaxies. Although the fraction of eDIG is high, the comparison between our sample and a generic sample of inclined spirals suggests that the phenomenon is uncorrelated to the galaxy environment. As suggested by the extraplanar UV emission found in 13 out of 14 galaxies the star formation extends well beyond the disc defined by the H$\alpha$ map.


INTRODUCTION
Galaxy interactions dramatically impact several galaxy properties, such as their morphology and star-formation history (see e.g.Rampazzo et al. 2016, and references therein).The idea of gathering catalogues of isolated galaxies has been developed at least since the 1970s (e.g. the Catalogue of Isolated Galaxies (CIG), Karachentseva 1973) as a baseline for comparison with galaxies located in denser environments subject to interaction-induced phenomena.In the last decades, catalogues of isolated galaxies, selected from very low galaxy density regions, have been either implemented or revised using redshift measurements (Verdes-Montenegro et al. 2005;Karachentseva et al. 2010;Argudo-Fernández et al. 2013, 2015).A recent H i study of isolated galaxy samples suggests a picture of "nurture free" galaxies (Jones et al. 2018).Other works (see e.g.Rampazzo et al. 2020) have shown that isolated elliptical galaxies suffered past interactions from which they still show the ★ Based on observations obtained at the Observatoire de Haute Provence (OHP, France), operated by the French CNRS.
† E-mail: mminerva@astro.unam.mx'scars', although they have probably not experienced gravitational influences from their close neighbours over the past billion years.
This paper concentrates on the stellar versus H components in isolated spiral galaxies seen nearly edge-on ( ≥ 80 • ).This galactic configuration relative to the sky plane provides more information on the vertical structure of the disc (e.g.Burstein 1979;Fraternali & Binney 2006;Kamphuis et al. 2007;Peters et al. 2017) revealing the extraplanar Diffuse Ionized Gas (eDIG), detected in galaxies with relatively high star formation rate (SFR) surface densities (Rossa & Dettmar 2003a,b).The eDIG is usually traced at optical wavelengths using the H recombination line (e.g Rosado et al. 2013).In addition, as the eDIG might be spatially correlated with the star formation (SF)-gas relations, it has been observed to trace regions of H i dominance within galactic discs, especially those located outside active star-forming regions, making the 21-cm emission line a helpful tool for the eDIG detection (e.g.Zschaechner et al. 2015).Ultraviolet (UV) imaging can also detect the low-intensity outer star formation disc by detecting OB stars located in very low SFR density regions (e.g.Hoopes et al. 2001;Thilker et al. 2002;Gil de Paz et al. 2005;Thilker et al. 2005Thilker et al. , 2007; Gil de Paz et al. 2007b).From the comparison between the H and the UV images of spiral galaxies with low inclination ( ≤ 80 • ), Thilker et al. (2007) showed that the Near Ultraviolet (NUV) emitting component is radially more extended than the emission of the old stellar population in at least 30% of galaxies from the Local Universe (z≈0.02).In high inclined galaxies, in cases where an eDIG layer does exist, the UV haloes tend to form a thick disc having a similar morphology to the UV halo and occurring in about the same place (Hodges-Kluck et al. 2016).Moreover, it has been observed that the H-UV flux ratio is lower in the eDIG than in H ii regions, indicating that the field OB stars, which are situated outside the star-forming regions, are important contributors to the eDIG ionization in most of the galaxies, and even in some galaxies they may be the dominant ionization source (Hoopes et al. 2001;Hodges-Kluck et al. 2016;Jo et al. 2018).
The eDIG morphology shows a wide variety of local morphological features (Rand 1998;Rossa & Dettmar 2003b;Rosado et al. 2013) such as prominent layers of diffuse gas, filamentary structures, or just one or a few patches of extraplanar emission (Rand 1998;Rossa & Dettmar 2003b;Rosado et al. 2013).Some eDIG originating phenomena have been proposed, such as Inter Stellar Medium (ISM) activity in the galaxy disc (e.g.Heald et al. 2006aHeald et al. ,b, 2007;;Rosado et al. 2013;Ho et al. 2016;Jones et al. 2017;Bizyaev et al. 2017), intergalactic medium accretion (e.g.Putman et al. 2012;Putman 2017;Levy et al. 2019;Bizyaev et al. 2022), satellite galaxies interaction (e.g.Walker et al. 1996;Zschaechner et al. 2015) and accretion or intra-cluster interaction (e.g.Tomičić et al. 2021;Sardaneta et al. 2022;Boselli et al. 2022).In parallel, several studies have demonstrated that the UV halo emission is consistent with a reflection nebula produced by dust in the halo (e.g.Hodges-Kluck & Bregman 2014;Shinn & Seon 2015;Hodges-Kluck et al. 2016;Jo et al. 2018;Shinn 2019).Various mechanisms have been suggested as potential sources for diffuse-and-global UV haloes, including galactic radiation, magneto-hydrodynamic phenomena, and dust accretion from the circumgalactic or intergalactic medium (see Shinn & Seon 2015, for a comprehensive examination).The contribution of each process over cosmic time is still unclear as they both produce hot and cold components that occupy roughly the same space (e.g.Hodges-Kluck & Bregman 2014).
Several studies from large surveys such as SAMI1 (Ho et al. 2016), MaNGA2 (Jones et al. 2017;Bizyaev et al. 2017Bizyaev et al. , 2022) ) and CALIFA3 (Levy et al. 2019), and from specific surveys as CHANG-ES4 (Lu et al. 2023) have employed various environmental definitions in the selection of their samples of edge-on galaxies, contributing to varying interpretations of the relationship between extraplanar gas properties and their host galaxies.For instance, they have either discarded major mergers with clear tidal features (Ho et al. 2016;Bizyaev et al. 2022), rejected galaxies that have another object in the field, including stars (Bizyaev et al. 2017) or excluded galaxies with distorted discs in their H images (Lu et al. 2023).Recently, two studies agreed that the eDIG could be a consequence of the accretion of the CircumGalactic Medium (CGM) (Levy et al. 2019;Bizyaev et al. 2022) which has also been proposed as probable UV halo origin (e.g.Hodges-Kluck & Bregman 2014;Shinn & Seon 2015).However, most of the literature still lacks a general agreement on the sources of ionization for eDIG.Since the eDIG may be related to the CGM, Bizyaev et al. (2022) visually examined the outer regions of each galaxy of their sample aiming to locate environmental structures that might be connected to the eDIG.Their findings ranged from small satellites lacking noticeable structure to large satellites even larger than the primary object.Therefore, the consensus or discrepancy in explaining the eDIG origin is likely a consequence of the different definitions of environment adopted.As a reference study case, isolated star-forming galaxies with the possibility of edge-on observations is needed.
In order to study the environmental effects on the eDIG distribution, we selected a sample of isolated nearby late-type, high-inclined galaxies from the CIG catalogue (Karachentseva 1973).Indeed, nearby galaxies are better resolved, late-type galaxies are those which contain the most H emission, and almost edge-on galaxies are those for which it is easier to detect the gas emission at high galactic latitude and for which the column density is higher, even in the disk outskirts.The extraplanar gas is commonly defined as the detectable H emission where the stellar continuum is not detected.This could be done from the continuum-subtracted emission-line image which delimits the stellar disc (e.g.Miller & Veilleux 2003;Rosado et al. 2013;Levy et al. 2019;Tomičić et al. 2021).This is typically what is done when observations with narrow-and broad-band filters around H emission are combined.Alternatively, the old stellar disk in spiral galaxies (Population II stars) could be better defined from redder broad-band in the near-infrared (NIR, e.g.Kamphuis et al. 2007;Ho et al. 2016;Bizyaev et al. 2017Bizyaev et al. , 2022) ) and young stellar population from bluer broad-band, e.g. in the UV, which traces the young stellar populations of ages up to tenths of Myrs (Hoopes et al. 2001;Thilker et al. 2007;Bianchi 2011;Kennicutt & Evans 2012).We used this method in this work.We present the H emission maps of a sample of galaxies.The net H maps (i.e.continuum free) were obtained using the scanning Fabry-Perot spectroimager (FP) GHASP, providing complete two-dimensional coverage of very extended line emission regions, ideal for studying the faint diffuse gas emission.We compared H FP monochromatic maps with NIR, NUV and Far UV (FUV) images available in the literature.We defined a certain threshold in the NIR image under which the stellar density is low and consider that the eDIG is the H emission below this limit.The geometric parameters of the galaxies, such as the isophotal position angle, peak light distribution and inclination, are obtained from the photometric analysis of the UV and NIR images, and will serve as reference when studying the H kinematics in a forthcoming paper.
The paper is structured as follows.In Section 2, we provide details about the sample selection.In Section 3, we describe the data acquisition and reduction processes.In Section 4, we explain how the stellar disc and the extraplanar components are identified.In Section 5, we present the individual results of the photometric analysis of each galaxy in our sample.In Section 6, we discuss our main results.Finally, in Section 7, we provide a summary and the conclusions of our results.In order to compare our results with previous surveys, throughout this work we assumed a Hubble constant of  0 = 70 km s −1 Mpc −1 (e.g.Thilker et al. 2007;Ho et al. 2016;Levy et al. 2019).

THE SAMPLE SELECTION
Isolated galaxies have been the subject of selection and study for many decades.Evidences are accumulating that over at least a few billion years the evolution of isolated galaxies has not been driven by interactions with physically associated companions (see e.g.Verdes-Montenegro et al. 2005;Karachentseva et al. 2010;Karachentsev et al. 2011;Rampazzo et al. 2020, and references therein).In this sense, isolated galaxies are the ideal reference sample to study internal galactic forces and galaxy morphological, dynamical and photometric evolution minimizing possible environmental effects.
The CIG catalogue (Karachentseva 1973) and its more recent revisions (Verdes-Montenegro et al. 2005;Verley et al. 2007b;Argudo-Fernández et al. 2013) are examples of a successful attempt to compile a sample of isolated galaxies.Historically, Karachentseva (1973) obtained her catalogue from the visual inspection of the Digitized Sky Survey (DSS) images of the 27 837 galaxies -and of their surroundings -(at Galactic latitudes || ≥ 20 • ) contained in the CGCG 5 .A list of 1050 galaxies were found to meet the isolation criterion she adopted (see also Karachentseva et al. 2010).
Since we are looking for the incidence of eDIG among isolated galaxies we have selected a sample from the CIG catalogue (Karachentseva 1973) of late-type galaxies (Spirals with morphologycal type > 5, i.e.Sb to be checked) (see Buta et al. 2019) dominated by Population I stars that satisfy the following additional characteristics we derived from the NED data-base: (i) galaxies with high inclination ( ≥ 80 • ) computed with the expression where  and  are the apparent optical major and minor axes of the galaxy respectively; (ii) galaxies with redshift  ≤ 0.02 to ensure that the emission of the redshifted H line fits in the optical wavelength; (iii) galaxies with apparent K  -band magnitude   ≤ 12.The near-infrared is sensitive to nuclear rings and large-scale bars, which might fuel active nuclei (e.g.Kormendy & Illingworth 1982;Eskridge et al. 2002;Jarrett et al. 2003).
A complete sample of 14 nearby isolated late-type edge-on 5 Catalogue of galaxies and of clusters of galaxies (CGCG, Zwicky et al. 1968) https://cdsarc.cds.unistra.fr/viz-bin/cat/VII/190 .galaxies was obtained (see Figure 1).General parameters of our galaxy sample are listed in Table 1.

Isolation degree
To be an isolated galaxy, the CIG catalogue required that no similar size neighbours were found close to the galaxy.This requirement can be described as follows: a galaxy with a standard angular diameter denoted as  1 is classified as isolated when the angular separation  1 between the galaxy and its -th neighbor with angular diameter   𝜖 [0.25, 4] •  1 , meets or exceeds a threshold of  1 ≥ 20   (see equations 1 and 2 in Karachentseva et al. 2010).Since galaxies from the CIG catalogue were collected by visual inspection of DSS images, the companion galaxies are expected to be faint or mostly dwarf galaxies.Recently, the CIG catalogue has been refined using new quantification methods (e.g.Verley et al. 2007b) and modern tools like images with higher resolution than the plates of Palomar (e.g.Argudo-Fernández et al. 2013) and new images from H i emission data (e.g.Jones et al. 2018).In this Section, we discuss how effectively our galaxy sample selected from the CIG catalogue meets these new isolation checks.However, even if a subset of our sample does not meet the strictest isolation criteria, all galaxies in our sample still represent a significant part of the overall isolated galaxy sample.
In the frame of the long-term AMIGA project (Analysis of the interstellar Medium of Isolated GAlaxies6 ), Verley et al. (2007a,b) revised all the positions of the CIG galaxies based upon the digitized photographic plates from POSS-1 and POSS-2 providing a quantification of the degree of isolation of CIG galaxies with the local number density of neighbour galaxies (  ), and the tidal strength (  ) affecting the central galaxy by its neighbourhood.These parameters depend on the position, diameter and mass of the primary galaxy and its -th neighbour.Thus, these parameters Figure 1.Superposition of the respective FP H monochromatic map and GALEX NUV image of the galaxies in our sample.The ellipses fitted to the surface brightness level at 1 (white) and 3 (yellow) of the typical background noise (see e.g.Jarrett et al. 2000) of the respective 2MASS K  -band image are overlain on the respective map.To distinguish the extraplanar material of each galaxy, the old stellar disc plane is traced by the ellipse fitted to the 3 surface brightness level.We have masked the nearest and brightest field stars to each galaxy in the UV-band images.At the distance of each galaxy (see Table 1) the kpc measure equivalent to 30 arcsec is indicated with an arrow.provide a clear picture of the environment around the isolated galaxies: the presence of one single similar size neighbour at a small distance to the CIG galaxy would result in a high value of the tidal strength (  ) estimation, while the local number density (  ) remains quite low as it is averaged over more nearest neighbours (Verley et al. 2007b).In Table 2 we list the number of neighbours and the nearest-neighbour distance to each galaxy in our sample according to Verley et al. (2007a).Then, Argudo-Fernández et al. (2013) re-evaluated the isolation criteria of AMIGA galaxies within a field radius of 1 Mpc using both photometric and spectroscopic data available from the Sloan Digital Sky Survey (SDSS) to refine the   and   parameters.
Figure 2 shows the comparison between the local number density (  ) and tidal strength (  ) of the galaxies in our sample making the distinction between the results of Verley et al. (2007b) and Argudo-Fernández et al. (2013).From Figure 2 and Table 2, we observe that some galaxies classified as isolated in the AMIGA sample have its nearest neighbour at an angular distance of twice or three times their optical diameters, for example CIG 171 and 201.According with Verley et al. (2007a), most of the neighbouring galaxies of CIG galaxies have a diameter ≤ 0.25 1 , implying they are dwarf companions which were not took into account by Karachentseva (1973) when collecting the CIG catalogue but perhaps they can influence the evolution of the main galaxy.Therefore, if CIG 171 and 201 are in region of bona fide isolated galaxies in the AMIGA sample, we may infer that the nearest neighbours are small enough to not exert significant gravitational perturbations on the main galaxy.On the other hand, Verley et al. (2007b) had already determined that the galaxies of our sample, CIG 71 and CIG 936, were outside this region, excluding both galaxies from the AMIGA sample.As Argudo-Fernández et al. (2013) improved the quantification of the isolation degree, another galaxy of our sample, CIG 1003, is also out of the area of isolated galaxies in the   −   plane.In this way, 3 out of 14 galaxies in our sample failed the isolation criteria of the AMIGA sample, probably because they have some faint companions previously undetected by Karachentseva (1973).
Later, Jones et al. (2018) provided integrated H i fluxes of AMIGA galaxies by presenting a catalogue of H i single dish observations.With these new data, Jones et al. (2018) included a review of integral features of AMIGA galaxies allowing to ensure isolation level of the AMIGA catalogue.Since the Argudo-Fernández et al. (2013) work is limited to the SDSS data, Jones et al. (2018) adopted the criteria of Verley et al. (2007b), and excluded all sources with heliocentric velocities below 1500 km s −1 , removed all dwarf galaxies that were in the CIG, and discarded all sources with potentially spurious spectral parameters, leaving a final sample of 544 CIG galaxies.In these terms, 8 out of 14 galaxies in our sample would not be considered well isolated : CIG 71, 159, 171, 201, 416, 922, 936 and 1003.While some members of our sample, selected from the CIG catalogue (Karachentseva 1973), may not fully meet the strictest current isolation criteria, it is important to stress that all galaxies in our sample are still classified as isolated galaxies.The different arguments stated above make us aware of the different degrees of isolation of the galaxies in our sample.However, previous studies indicate that late-type galaxies are usually located in very low-density environments.In this context, we decided not to exclude any galaxy a priori but rather discuss each of them knowing their isolation properties.

H𝛼 observations and data reduction
Observations of the H emission line were performed using the scanning Fabry-Perot (FP) interferometer, GHASP, attached at the Cassegrain focus of the 1.93m telescope at the Observatoire de Haute-Provence (OHP, see Table 3).A GaAs Image Photon Counting System (IPCS) is used to reach low detection levels, offering a total FoV of about 5.8 × 5.8 arcmin 2 with a pixel scale of ∼ 0.68 arcsec pix −1 (Gach et al. 2002).Since the IPCS has no readout noise, exposures of 10 seconds were chosen for each of the 32 scanning steps required to cover the GHASP's free spectral range (FSR).In order to detect the eDIG emission, the scanning sequence was typically repeated ∼35 times, resulting in a total exposure time per galaxy of ∼180 min.We used narrow-band interference filters (mainly FWHM≃15 Å and in a few cases slightly larger up to 24 Å) enabling to select the redshifted H line of ionized hydrogen (6562.78Å).The scanning Fabry-Perot interferometer we used, with  4) distance to the nearest neighbour in arcseconds and kpc respectively given the heliocentric galactic distance listed in Table 1.an interference order of  ≃798 and a finesse F ≃13 at H, allowed GHASP to reach a typical spectral resolution of  ∼ 10 000.Two calibration cubes were obtained, one at the beginning and another at the end of the observation, using a neon lamp with selected narrow line at 6598.95 Å, close to the redshifted nebular wavelength in order to minimize phase shift effect (e.g.Gómez-López et al. 2019).The observation date, filter's central wavelength and FWHM, the total exposure time and the astronomical seeing of each galaxy are listed in Table 4.
For the data reduction, we used the packages based on homemade IDL7 routines, reducWizard and ComputeEverything (see e.g.Epinat et al. 2008).The data reduction technique has been widely reported in Daigle et al. (2006b).In summary, it consists in the following steps: (1) integration with guiding correction of the total data cubes obtained during the observation; (2) application of a Hanning spectral smoothing on the data cube; (3) calibration in wavelength of the integrated data cubes through the computation of a parabolic phase map which is computed from the calibration cube in order to obtain the reference wavelength for the line profile observed inside each pixel, creating a wavelength-sorted data cube by applying the phase map correction to the interferogram data cube; (4) subtraction of the OH sky-lines emission; (5) a Gaussian spatial smoothing with a FWHM selected on the wavelength data cubes; and (6) computation of H monochromatic, continuum, radial velocity and velocity dispersion maps for each wavelength calibrated cube.Finally, astrometric information was attached to the processed files by using the IDL task koords from the KARMA package (Gooch 1996).In Table 4, we also listed the seeing and the final spatial resolution in kpc arcsec −1 after the spatial Gaussian smoothing of the H maps.
We performed two separate processes on the wavelength data cube.First, we applied Gaussian spatial smoothing, and then, independently, we performed adaptive spatial binning using a Voronoi tessellation with a target SNR≃5.To ensure a signal-to-noise ratio SNR≥ 5 in the Gaussian smoothed data cube, we compared the monochromatic H maps derived from both processes and set a threshold for the emission detected in the Gaussian smoothed maps.For a SNR≃5, the GHASP surface brightness detection limit is  = 2.5 × 10 −17 erg sec −1 cm −2 arcsec −2 (see Epinat et al. 2008;Gómez-López et al. 2019;Sardaneta et al. 2022).As no H calibrated image of any galaxy in our sample has been previously published, we followed the calibration procedure outlined by Epinat et al. (2008) to determine the total H flux for the GHASP data (see Appendix B).
The data analysis was made with the program ADHOCw8 , IRAF 9 tasks, the SAO Image DS9 software 10 and our own Python scripts.

Archival data
The DSS is available at the site of the Space Telescope Science Institute (STScI) 11 .The DSS 6.5×6.5 deg 2 plates have been scanned using a modified PDS microdensitometer with a binned pixel scale of about ∼1.0 arcsec pix −1 .In this work, we use the DSS R-band images (5900-7150 Å) of our galaxy sample as a reference to provide a large-band optical view of the target.However, we do not calibrate these images in flux since their zero-point magnitude is currently unknown.The average seeing obtained from the image headers was ∼1.5 arcsec.
The Two Micron All-Sky Survey (2MASS, Jarrett et al. 2000) project uses two highly-automated 1.3 m telescopes, one at Mt. Hopkins, Arizona, and the other at Cerro Tololo, Chile.Each telescope is equipped with a three-channel camera, each channel consisting of a 256 × 256 pix 2 array of HgCdTe detectors, capable of observing the sky simultaneously at J (1.24 m), H (1.66 m) and K  (2.16 m) near-infrared bands with an angular resolution of 2.0 arcsec and a pixel size of 1.0 arcsec pix −1 .Owing to the transparency of interstellar dust within the galaxies and to the dominance of late-type stellar populations in producing galactic near-infrared flux, the 2MASS galaxy images trace the overall stellar mass distribution in these galaxies 12 .To calibrate in flux, the zero-point conversion values are published in the Explanatory Supplement to the 2MASS All Sky Data Release and Extended Mission Products 13 .In this work we use the K  -band images of the nearby isolated late-type edge-on galaxies selected to sketch the stellar Population II.
Galaxy Evolution Explorer (GALEX) was the first ultraviolet (UV) all-sky survey covering a field of view (FoV) of ∼ 1.25 deg using microchannel plate detectors to obtain direct images in the NUV ( eff = 2271 Å) and FUV ( eff = 1528 Å) with a resolution of 4.2/5.3arcsec (FUV/NUV) and a pixel scale of 1.5 arcsec pixel −1 (Bianchi 2011) 14 .Conversion values between GALEX count rates, fluxes, and AB magnitudes are available in the GALEX Guest Investigator Web Site 15 .We use GALEX-NUV and FUV images as tracers of stellar Population I in galaxies in our sample.
Thus, in this work we study H, NIR, and UV images, each one with different resolutions.In the case of the H emission data, the seeing conditions varied across the sample of galaxies, with different galaxies experiencing different levels of atmospheric turbulence due to the diversity of observational conditions during data gathering (see Table 4).For example, the closest galaxy, CIG 71, experienced the worst seeing conditions with a resolution of ∼4.4 arcsec (1.8 kpc), as did the most distant galaxy, CIG 416, with a resolution of ∼4.1 arcsec (0.7 kpc).However, CIG 922, a galaxy located at an average distance, had relatively average seeing conditions with a resolution of ∼2.2 arcsec (0.6 kpc).The seeing for the NIR (2 arcsec) and NUV/FUV (4.2/5.3 arcsec) images is considered constant for the whole sample.However, it is important to be aware that the different seeing disc sizes combined with the different distances of the galaxies blur the data and impact the spatial physical resolution making it more difficult to distinguish fine details such as eDIG filaments or patches.
Finally, since the Infrared Astronomical Satellite (IRAS) satellite did not cover the whole sample of CIG galaxies, Lisenfeld et al. (2007) reprocessed the IRAS MIR/FIR survey data using the ADDSCAN/SCANPI utility for 1 030 out of 1 050 CIG galaxies as part of the AMIGA project.In this work, we use the IR flux densities at 60 and 100  from AMIGA survey (see Table 1).

IDENTIFYING THE STELLAR DISC
To disentangle the galactic disc from the gaseous extraplanar component, we performed a photometric analysis of the 2MASS 2.2m image, tracer of the old stellar population, and of the GALEX NUV and FUV images, connected to the young stellar population in the galaxy.Because of the particular characteristics of the images in each band, to derive the surface brightness profile, we used two different software packages, the ellipse task from the IRAF STSDAS package (Jedrzejewski 1987) and the isophote package from the Python Photutils package 16 (Bradley et al.  2020).Both programs fit an ellipse via an iterative method (devised by Jedrzejewski 1987) as well as quantify the distortion from a perfect ellipse by means of higher-order Fourier harmonics: where  0 is the surface brightness averaged over the ellipse as a function of the azimuthal angle  and,   and   the higher order Fourier coefficients.The first-and second-order coefficients ( 1 ,  1 ,  2 ,  2 ) indicate the errors in the fitting procedure, being zero for a best-fit ellipse, while the third-order coefficients ( 3 ,  3 ) characterize the deviation of the fitted ellipse from the isophote shape.The coefficient  4 measures symmetric distortions from pure ellipticity: when  4 > 0, the isophotes have a 'disc-like' (circular) shape and, when  4 < 0, the isophotes have a 'boxy' (rectangular) shape.The inclination angle () with respect to the sky-plane of the galactic disc can be measured through the ellipticity () of the fitted ellipse.Using diverse tracers of the stellar population serves as a crucial means to differentiate between the particular cases presenting UV extended discs or UV haloes which represent at least one-third of the galaxies in the local universe (Thilker et al. 2007) and, as we will show for our sample in Section 5, approximately two-thirds of the late-type isolated high-inclined galaxies.
The analysis of edge-on galaxies is more difficult than the analysis of galaxies with lower inclination because the high inclination affects the apparent shape of the galaxy along the line of sight (LoS).In photometric analysis of low-inclination galaxies, a smooth bulge-disc transition is typically observed in the resulting plots (e.g.Marino et al. 2010).However, in highly inclined galaxies, the brightness of the inner regions is dominated by the bulge with a small contribution from the thin disc.As the radius increases, the contribution from the disc becomes more important, leading to an abrupt bulge-disc transition in the resulting plots.This provides a noisier behaviour at small radii in the photometric analysis for highly inclined galaxies.
Considering that in late-type galaxies old and young stellar populations are mixed in the disc (e.g.Nersesian et al. 2019), we define the stellar homogeneous disc by using the ellipse delimiting the locus of the stars belonging to the old stellar population since they trace the bulk of the mass and probe the global potential of the disc galaxy.

NIR emission
For the 2MASS 2.2 m image, we applied the method of elliptical isophote fitting using the task ellipse (Jedrzejewski 1987) from the IRAF STSDAS package.As initial parameters we set the galaxy coordinates as the centre, the B-band optical diameter ( 25 ) as the maximum semi-major axis length, and the optical position angle (PA) and ellipticity () published in the NED (see Table 1).We used logarithmic radial sampling with an initial step of ∼0.5 pixel along the semimajor axis in order to derive the surface brightness profile of each galaxy in an iterative process.The centre, PA, and ellipticity of the ellipse were allowed to vary in each iteration.Jarrett (2000) defined the isophotal aperture  20 , derived from the K  -band isophote at   = 20 mag arcsec −2 , as the 2MASS standard aperture, corresponding to roughly 1 of the typical background noise in the -band images.However, because of the high levels of background noise in the 1-2 m atmospheric windows, the 2MASS is not as sensitive to low surface brightness emission from galaxies (Jarrett 2000).In consequence, the K  -band benchmark of   = 20 mag arcsec −2 elliptical isophote aperture might underestimate the total flux between 10% and 20%, depending on the radial distribution according to the Hubble type (Jarrett 2000;Jarrett et al. 2003).On the other hand, discrepancies in the measurements at the   = 20 mag arcsec −2 level between the 2MASS and other K  -band surveys (see e.g.Fingerhut et al. 2010) have suggested potential errors in the 2MASS data reduction process.In fact, Jarrett (2000) explained that systematic components such as H-band airglow variations were not well understood by then and might have induced large errors in the photometry 17 .Consequently, the choice of a brighter isophote to define the galactic disc boundaries, beyond the 1 level at   = 20 mag arcsec −2 , may help to mitigate potential data reduction errors in the 2MASS data set.
To be consistent with previous optical photometric results, Jarrett et al. (2003) derived some galactic properties from the 3 isophote, such as the ellipticity () and position angle (PA), in the 2MASS Large Galaxy Atlas (LGA) and the 2MASS Extended Source Catalogue.These photometric parameters are consistent with our K  -band results as well (see Table 5).Furthermore, the K  elliptical isophotal photometry should include the core, bulge, and disk components, capturing most of the flux from a galaxy (e.g.Jarrett et al. 2000).Therefore, to ensure that a significant portion of the galaxy's light is captured and to avoid empty sky background noise, hereafter we use the ellipse at 3 on the background as reference for the disc traced by the stellar Population II.Although this value may not represent the true galactic radii, it is a practical starting point to define the old population stellar disc boundaries.
The older stellar population in spiral galaxies points out internal structures such as spiral arms, bulges, warps, rings, and bars (Jarrett et al. 2003).Bars play a major role in the secular evolution of galaxies and, through vertical resonances, drive stars above the plane to form peanut/box-shaped pseudo-bulges (Combes & Sanders 1981a).The non-axisymmetry of the bars produces the gas to flow inwards 17   2), ( 3) and ( 4), respectively, and at 3 level in columns ( 5), ( 6) and ( 7).Same parameters determined for the NUV image are listed in columns ( 8), ( 9) and ( 10), and for the FUV image in columns ( 11), ( 12) and ( 13).
fuelling SF in the central regions of their host galaxies (see e.g.Kim et al. 2021, and references therein).Bottom frame of panel (d) of Figure 3, and from Figures E1 to E13, shows the variation of the  4 parameter with respect to the galactic semi-major axis (radius).From these graphics, we found that the only one galaxy in our sample that presents a coefficient  4 < 0 is CIG 329 (see Figure E6) implying it has a 'boxy' bulge.

UV emission
The ellipse task from IRAF works in an iterative manner, if there is contamination such as stars and H ii regions or regions of too low SNR, the program continues generating ellipses but it outputs parameters with undefined values.To overcome this difficulty, other fitting software was used to derive the surface brightness profiles of the NUV and FUV images.The isophote package from Photutils (Bradley et al. 2020), provides tools to fit elliptical isophotes to a galaxy image using an iterative method described by Jedrzejewski (1987) (see equation 2) for each ellipse as well.
The isophote package can also work in an iterative manner by measuring a maximum acceptable relative error in the local radial intensity gradient (typically 0.5).However, when two consecutive isophotes exceed the value specified by the parameter the program prevents ellipses either from prematurely stopping due to the stellar contamination or from growing in low SNR regions.Nonetheless, if the maximum semimajor axis is specified, the intensity gradient is set to 'none' and the algorithm proceeds inwards to the galaxy centre.Although it is not the general case and may not be the best fit, this procedure allows to analyse data with too low SNR in the inner regions of the galaxy, such as in some GALEX NUV and FUV images of galaxies in our sample.This fit algorithm is quite sensitive to the initial guesses.The iteration begins by giving a fixed semimajor axis length selected initially to grow linearly with steps of 0.5 pix and specifying the initial value of the galactic centre, PA and , then the algorithm proceeds to fit isophotes inwards to the galaxy centre.
We initialized the fitting parameters using the galaxy coordinates as the centre, the optical PA and ellipticity obtained from the NED database, and the B-band optical diameter ( 25 ) as the maximum length of the semi-major axis obtained from the NED database (see Table 1).To overcome any discrepancies between the optical and UV radii, we made adjustments to the initial parameters by comparing the UV radius measured using the DS9 image analysis software with the optical radius, and modifying the corresponding initial parameter values if the optical radius was smaller than the UV radius.For most of our targets, we allowed the centre position, ellipticity and PA to vary during the fitting process.Some galaxies (CIG 71, 95, 847, 906 and 1003) present the NUV/FUV brightest knots beyond ∼ 10 arcsec from the peak light distribution of the NIR emission, this last position is assumed the center of the galactic disc and it is probable that the NIR emission is shadowing the actual peak of the UV light.Hence, for these galaxies, we fixed the centre position at the NIR peak light distribution, assuming it indicates the galactic center, and we allowed the ellipticity and PA to vary during the fitting process.Thus, we ensure that the fitted ellipses accurately reflect the morphology of the target galaxies, taking into account any deviations from the ideal elliptical shapes.
We computed the isophotal radii at surface brightness levels of  NUV,FUV = 28 mag arcsec −2 in NUV and FUV, which roughly correspond to the average surface brightness at the optical diameter (see e.g.Gil de Paz et al. 2007a;Marino et al. 2010;Cortese et al. 2012), using the limit AB magnitude in NUV of 20.08 mag and in FUV of 18.82 mag18 .

Photometric and geometric parameters at different wavelengths
Photometric parameters, such as the PA and galaxy inclination, are useful observational constraints to reveal galaxy peculiarities when the galaxy configuration is examined, particularly when different wavelengths are considered.In Table 5 we list the photometric parameters resulting from the ellipse fitting process.Panel (g) of Figure 3 and Figures E1 to E13 shows the surface brightness profiles of GALEX NUV and FUV images, while the resulting ellipses were overlaid on the NUV (panel e) and FUV (panel f) images, respectively (see for example Figure 3).The 1 NUV and FUV isophotes computed after applying a Gaussian spatial smoothing of  ≈4.5 arcsec (3 pix) on the NUV and FUV images were overlaid on the H monochromatic map (panel a) in order to allow a better comparison of the morphology traced at different wavelengths.
Inclination plays a critical role in the analysis of galaxy structure as it is used to correct its main parameters as the brightness (e.g.Stone et al. 2021) and the rotation velocity (e.g.Epinat et al. 2010).Assuming that the image of a spiral galaxy is the projection of a disc with the shape of an oblate spheroid, the inclinations of spiral galaxies can be roughly derived from the ellipticity of apertures used for photometry.If the thickness of a disc is neglected the relation between the inclination () to the celestial plane and the apparent diameter ratio can be described using the relation 1.If the thickness is taken in consideration, the inclination can be estimated using the expression: where  = 1 − / is the ellipticity at the isophotal radius and  0 = / is the flattening, in this case ,  and  are the three spheroid's axes (see e.g.Barbosa et al. 2015).Many works have shown that the flattening is type dependent (see e.g.Haynes & Giovanelli 1984).In this sense, Hall et al. (2012) determined an intrinsic flattening for spiral galaxies,  0 = 0.13, by studying the relationship between the axial ratios in the -band versus the  −  galaxy colour computed at the effective radius for 871 edge-on galaxies from the SDSS.They showed that, although  0 increases with redder, progressively bulge dominated galaxies, the  0 distribution was rather flat for late-type galaxies only with an upturn in the transition between Sa to S0 galaxies.Nonetheless, deriving the inclination from ellipticity may be problematic.For instance, prominent bulges can dominate the axial ratio measurement or in other cases, galaxies may not be axially symmetric owing to tidal effects or high surface brightness bars embedded in low surface brightness discs (Epinat et al. 2008(Epinat et al. , 2010;;Hall et al. 2012;Kourkchi et al. 2019;Stone et al. 2021).All these rather common configurations can lead to large errors.Some solutions have been proposed, for example, in order to avoid contamination due to clumps of star formation, Epinat et al. (2010) suggested that the inclination should be preferentially derived from broad-band images with high resolution, as in the NIR rest frame of the galaxy tracer of the old stellar disc.
In addition to photometric techniques, to compute the galaxy inclination, some studies use visual inspection to evaluate the inclinations.For example, to select the high inclined galaxies from SDSS and CALIFA surveys, Bizyaev et al. (2017) and Levy et al. (2019), respectively, used the visibility of the dust lane in optical images of the galaxies to confirm the edge-on nature of the galaxies, excluding galaxies with dust lanes that were not centred in the midplane or displaying any visible track of spiral arms, bars, or other features.Kourkchi et al. (2019) evaluated this technique by designing an online graphical tool, achieving uncertainties of ±4 • and suggesting to replace the human eye with a machine learning algorithm to compute more accurate inclinations.
In this work, the galactic inclination  was one of the selection paramaters (see Section 2) for our sample.Firstly, the inclination of the CIG galaxies in our sample was computed by the relation 1  4), ( 5), ( 6) and ( 7) inclination computed taking into account the galactic thickness (see the relation 3), using the ellipses fitted to the 1 and 3 level of the NIR image and to the one fitted to the brightness level  NUV = 28 mag arcsec −2 of the NUV image.In column (3), by using the flattening constant  0 = 0.13 determined by Hall et al. (2012), the lower limit of ellipticity 1 −  = 0.13, for which equation 3 is valid, is exceeded by some of the galaxies in our sample in the optical band.
using the galactic apparent major and minor axes from the optical wavelength published in the NED.Next, we computed the inclinations of galaxies in our sample applying the relation 3, for the the 2.22 m and NUV emission images we employed the ellipticity obtained from the ellipse fitting performed to them, and for the optical wavelength image, we used the galactic apparent major and minor axes from the NED.In all cases, we used the constant flattening  0 = 0.13 determined by Hall et al. (2012).In the case of the NIR image we have computed the inclination of both 1 and 3 isophotal levels.
In Table 6 we compare the inclinations obtained.Equation 3 has a lower threshold of effectiveness depending on the flattening value ( 0 ) considered.Since Hall et al. (2012) determined the flattening value  0 in the -band, in the optical band we obtained inconsistent results for some galaxies in our sample having ellipticity 1 −  ≥  0 = 0.13.Therefore, in the following we will discuss the inclination calculated considering the galaxy as infinitesimally thin only in the optical band (see relation 1).Excluding particular cases, the inclination computed with the ellipse fitted to the NIR 3 isophotal level has an averaged difference of ∼2 • with respect to the optical inclination, while the inclination from the NIR 1 isophotal level has an averaged difference of ∼4 • and the mean difference with the inclination from the NUV emission is ∼3 • .Hence, the inclinations computed with the parameters of the ellipse fitted to the NIR image at the 3 isophotal level are the most congruent with the optical inclination, and following Epinat et al. (2010), this may be the best estimation to measure the galaxy inclination.
When computing the inclination, particular cases showing Δ ≳10 • with respect to the optical inclination were detected.Inclinations computed with the NIR 3 isophotal level presented four particular cases: CIG 95, 416, 593 and 936, which have a bright central bulge dominating the NIR emission.In the case of the NUV emission, the particular cases are CIG 201, 416, 906 and 1003, which have the most extended UV haloes and the largest filaments above the disc.In contrast, CIG 936 is the only one particular case obtained for the NIR 1 isophotal level.The galaxy CIG 936 (Figure E12) is an extreme case in our sample displaying its NIR and H emission concentrated in the inner central ∼30 arcsec (∼6.3 kpc) of a disc with optical diameter of  25 () = 185.4arcsec (∼38.9 kpc).This galaxy might be compared with the Sombrero galaxy (NGC 4594), an S0 galaxy whose bulge governs the axial ratio (e.g.Hall et al. 2012), though we cannot consider CIG 936 an early-type system whilst its disc is dominated by the UV emission, tracer of the young stellar population.The results previously described seem to indicate that the apparent discrepancies in the inclination values determined using images in different wavelengths are due to actual asymmetries of the old stellar disc and the presence of an extended UV disc and halo.
In conclusion, inclination measurement is essential for unraveling the complex structures of galactic components, in particular UV extended discs and haloes.By applying different stellar population tracers, this study allowed us to better separate these different components.

RESULTS
In this section we present individual results obtained by comparing properties of the DSS R-band, the 2MASS K  -band, the GALEX NUV and FUV images versus the H monochromatic map from FP GHASP data.For each galaxy we indicate the old stellar population disc using the ellipse fitting to the isophotal level at 3 on the background of the 2MASS K  -band image, and we trace the young stellar population using both the ellipse corresponding to the  NUV,FUV = 28 mag arcsec −2 and the 1 isophote of the GALEX NUV and FUV images.An example is shown in Figure 3.

CIG 71 (UGC 1391):
The H monochromatic map shows a discontinuous disc mainly shifted Southwards from the photometric centre of the NIR-band exceeding the radius of the old stellar disc.A prominent and luminous knot is observed at ∼30 arcsec Southwards from the NIR photometric centre.The K  -band image shows a symmetrical disc within the ellipse fitted to the 3 brightness level, but the disc becomes more extended and warped towards the South.The FUV emission is concentrated near the galactic disc centre and then spreads Southwards at lower brightness.This map also shows a FUV emission cloud in the North.The NUV emission displays a continuous disc extended along the major axis beyond the old stellar population disc.The photometric centre does not match with the H, NUV, and FUV emission maxima.The isolation criteria of Verley et al. (2007b) placed CIG 71 outside the isolation plane (see Figure 2), which was confirmed using H i emission by Jones et al. (2018).Moreover, we observed that CIG 71 has common features observed in interacting galaxies (e.g.NGC 5258 from KPG 389, Fuentes-Carrera et al. 2019).Likely, this galaxy experienced an interaction despite its nearest neighbor is 384 kpc away (see Table 2).
Remaining galaxies are discussed in Appendix D.

DISCUSSION
In the following we discuss our findings at the light of the current literature focusing on three aspects: the stellar and warm gas emission in the disc, the detection of extraplanar warm gas emission and the extension of the UV emission with respect to NIR and H emission.

The emission in the disc
In this work, we have examined the location of both the old and young stellar population in order to determine the boundaries of the galactic disc in highly inclined galaxies to reveal the existence or lack of a gaseous extraplanar component.Since in late-type galaxies the old and young stellar populations are mixed in the disc and arms, and to be consistent with previous photometric results (see e.g.Jarrett et al. 2003), we normalised the stellar disc using the ellipse fitted to the 3 on the background level of the NIR image, which delimits the region occupied by both old and young stellar populations.
Nonetheless, the definition of the galactic disc, which is used to find an extraplanar component, is not immediately obvious.In Figures 3,  and E1 to E13, we plot the DSS R-band image of galaxies in our sample to compare images in different wavelengths with the optical one.We observe that the radius of the R-band image is larger than the radius of the isophotal 20 mag arcsec −2 radius in K  -band.In fact, Jarrett (2000) compared the K  -band to optical-band (B-band) imaging properties (size, total flux, central surface brightness and elongation) in order to understand the stellar light coming from different components.He found that the K-band elliptical isophotal diameter corresponding to   = 20 mag arcsec −2 surface brightness ( 20 ) represents about 70% of the optical radius measured at   = 25 mag arcsec −2 , showing that the NIR sensitivity decreases for the later-type (i.e., 'blue') galaxies, and concluding that as late-type galaxies evolve, their cores gradually shift towards reddening and brightening, resembling early-type disks or spheroids.Certainly, the K  -band predominantly traces the oldest stellar population (stars older than ∼4 Gyr, Fingerhut et al. 2010, and references therein), while optical bands like R or B-band are more sensitive to the younger, bluer stars.Therefore, it is expected that K  -band measurements correspond to a smaller region.Since our main goal is to study the extraplanar ionized gas component in edge-on galaxies, we use 2MASS K  -band imagery looking specifically the location of the old stellar population, tracer of the galactic fundamental structure, with the aim of isolating the oldest stars from the effects of young, blue stars.
On the other hand, the Wide-field Infrared Survey Explorer (WISE) have completed a mid-infrared survey of the entire sky in bands centered at wavelengths of 3.4, 4.6, 12, and 22m.However, its angular resolution of 6.1 arcsec, 6.4 arcsec, 6.5 arcsec, and 12.0 arcsec at these bands respectively (e.g.Wright et al. 2010;Jarrett et al. 2012), is notably inferior compared to the higher angular resolution of 2.0 arcsec achieved by the 2MASS for the K  -band images.Therefore, the selection of the narrow-band K  -band imagery for our analysis offers a clear advantage over the DSS R-band or WISE data, as it allows a more precise distinction of the galactic components we are interested in studying in the galaxies in our sample.This choice allows us to discriminate between galactic disc emission and extraplanar emission, providing us with a more solid basis for our analysis and interpretation of the eDIG in galaxies of our sample.
Considering that the mass of our galaxy sample ranges between 10 9 and 10 10  ⊙ (see Table 7), one might infer that our galaxies are dwarf galaxies with fainter emission in older stellar populations.To explore the classification of dwarf galaxies, we compare our sample with the dwarf irregular (dI) galaxies from Fingerhut et al. ( 2010)'s survey.Firstly, it is worth noting that most of these dI galaxies were too faint to be detected by the 2MASS, implying magnitudes of  ≥ 13.5 mag (Jarrett et al. 2000;Jarrett 2000;Jarrett et al. 2003).Secondly, the mass range of their dI sample falls between 10 7 to 10 9  ⊙ (McCall et al. 2012).In contrast, our sample In panels (f) and (g) the surface brightness profile of the K  -band and NUV/FUV emissions are plotted, respectively.In both panels, yellow vertical lines mark the surface brightness level at 1 (dotted line) and 3 (solid line) of the K  -band image.In panel (g), vertical lines indicate the radii corresponding to the isophotes at a surface brightness level at  NUV,FUV = 28 mag arcsec −2 for the NUV (blue dot-dashed line) and FUV (red dashed line) maps, respectively.Additionally, the vertical black loosely-dashed line at ∼5 arcsec shows the approximate FWHM of the GALEX point spread function (see e.g.Marino et al. 2010).The ellipses fitted to these isophotes are superimposed on the corresponding map with their associated colour.The common elements in in panels (a), (b), (d) and (e) are: a colour-bar indicating the flux in units of erg sec −1 cm −2 , the arrow tracing a scale of 30 arcsec, a black circle embedded in a box indicating the spatial resolution of the image, a black ellipse representing the surface brightness level at 3 on the background of the K  -band image and, a black cross (×) pointing the location of the NIR band peak light distribution.At the bottom of panel (f) we show the relationship between the coefficient  4 and the radius using the K  -band map to determine if the isophotes tend to have 'disc-like' ( 4 < 0) or 'boxy' ( 4 > 0) shape (see Section 4).Finally, in panel (a), on the H monochromatic map were overlaid the 1 isophotes on the background of the NUV and FUV images.consists of galaxies with apparent K  -band magnitude  ≤ 13.5 mag (see selection criteria in Section 2), exhibiting masses one order magnitude higher and having been previously studied in the context of the 2MASS LGA and the 2MASS Extended Source Catalogue.This prior analysis ensures that our galaxies are indeed extended objects with bright NIR photometry.Therefore, in the context of Fingerhut et al. (2010), our sample does not align with the category of dwarf galaxies.Since we are not studying dwarf galaxies, we could claim that the surface brightness isophote at 3 of the K  -band image, selected to trace the stellar disc, encompasses the light of both the bulge and the disc.
Inside the ellipses fitted to the 3 level of the 2.22 m images, we observed regular discs, while the isophote at 1 level shows that the old stars at the borders follow the shape of the halo.For instance, the faintest NIR emission in CIG 71 suggests that the old stars are displaced to the South as well as the H emission, and in CIG 201 the faint old stars trace the same shape of the filamentary structure outlined by the H and UV emission.Regarding the behaviour of the  4 parameter versus the major axis relation (see Section 4), in most galaxies the pattern of the curve changes increasing its dispersion from the radial position of the 3 level.This is more evident in CIG 71, 171, 329, 847 and 936.While a similar break in the UV surface brightness profiles is traced at the radial position of the isophote at 1 level of the K  -band image, as observed in CIG 171, 201 and 906.This break in the surface brightness profiles could define the edge of the stellar disc, however, we still relay in our selection because of the better resolution of the 2MASS than the one of the GALEX images (see Section 3).Moreover, for some galaxies as CIG 329, 922 or 1003, if the stellar disc is defined by the the isophote enclosing   = 20 mag arcsec −2 (1 level), they would appear to be dominated by an old stellar disc and no extraplanar gas would have been detected.In this way, besides to the environmental criteria adopted for selecting a sample, the paradox in selecting the edge of the stellar disc might be another parameter causing discordance in the published results regarding the detection of eDIG.
We have found a variety of morphologies in different wavelengths in galaxies of our sample.With respect to the old stellar population disc, most of the K  -band images (12 out of 14) showed regular discs with a central bulge and some small filaments along the disc that do not break the stellar disc symmetry in general.Two galaxies showed an asymmetric or perturbed stellar disc: CIG 416 (Figure E7) presenting a tail shape drawn in its K  -band image, and CIG 593 (Figure E8) with a truncated old stellar disc that follows the morphology of the disturbed and faint H disc.
On the other hand, the only one galaxy in our sample reliably identified as LINER (Sabater et al. 2012) and whose photometric analysis revealed a coefficient  4 < 0 suggesting a "boxy" bulge is CIG 329 (Figure E6).Although it is difficult to detect a bar in edge-on galaxies, a plausible explanation could be that there is a bar in CIG 329 which may have undergone a vertical bending instability to form this box-shaped bulge.It is well known that two-thirds of star-forming disc galaxies in the local universe are barred (e.g.Eskridge et al. 2000;Whyte et al. 2002;Aguerri et al. 2009;Masters et al. 2011).Numerical simulations of isolated galaxies have shown that a stellar bar can be formed spontaneously (see e.g.Miller et al. 1970;Hohl 1971;Combes & Sanders 1981b;Sellwood & Wilkinson 1993;Athanassoula 2003;Athanassoula et al. 2013;Sellwood 2013;Saha & Elmegreen 2018;Ghosh et al. 2023), but also, that a bar can be formed by a 1:1 prograde flyby interaction and remain several Gyr after the interaction when galaxies are already separated by several Mpc (e.g.Lang et al. 2014).If the stellar bar of CIG 329 was formed by a secular process, why has no bar formed in other galaxies in our sample?It might be that these galaxies are hostile to bar formation or that the bar has been destroyed in the course of their evolutionary pathways.Indeed, the presence of a massive central bulge can suppress the bar instability (e.g., see Saha & Elmegreen 2018;Kataria & Das 2018).Alternatively, bars can be destroyed by a massive central mass concentration and/or inflow of the interstellar gas in the central region (e.g., see Pfenniger & Norman 1990;Shen & Sellwood 2004;Athanassoula et al. 2005;Bournaud et al. 2005;Hozumi & Hernquist 2005;Athanassoula et al. 2013).Furthermore, recent study by Ghosh et al. (2021) revealed that a minor merger can lead to a substantial weakening of bars, and in extreme cases, even a complete destruction of bars in galaxies.
Since most of the galaxies in our sample showed extraplanar material in both UV and H emission, in the next sections we discuss some general results from the sample in order to understand the origin of the eDIG and UV-haloes detected.

The eDIG detection
In general, not all edge-on galaxies show detectable eDIG emission, i.e. a gaseous halo above the layer of H ii regions associated with the disc.The presence of eDIG in late-type spirals is a direct consequence of SFR in the underlying galactic disc (e.g Rossa & Dettmar 2000, 2003a).It has been found that the prominence of gaseous haloes correlates with tracers of star formation in the disc, such as the surface density of FIR emission, the SFR determined by H luminosity, and dust temperatures (e.g.Heald et al. 2006a;Lu et al. 2023).For instance, Rossa & Dettmar (2000) constructed a diagnostic DIG diagram (DDD) that shows the ratio of the flux densities at 60 and 100 m expressed as  60 / 100 versus the ratio of the FIR luminosity (    ) divided by the optical diameter of the 25  ℎ mag arcsec −2 isophote squared ( 2 25 ), with which they demonstrated that galaxies suggesting features associated with starburst activity as  60 / 100 ≥ 0.4 and     / 2 25 ≃ 1×10 40 erg s −1 kpc −2 are candidates to display eDIG layers.
In order to locate our sample in the DDD, we computed the FIR-luminosity using the expression: with   (60) and   (100) in Jy (obtained from Lisenfeld et al. (2007) and presented in Table 1) and  the distance to the galaxy in Mpc (Rossa & Dettmar 2003a).The term  2 25 was expressed in kpc 2 .From our maps, showed in Section 5 and in Appendix E, we detected galaxies showing one or more morphological eDIG features (see e.g.Rossa & Dettmar 2003a;Heald et al. 2006a;Rosado et al. 2013), summarized in Table 7.The brightest H ii regions in Local Group galaxies have a surface brightness detection between magnitude orders of 10 −14 and 10 −15 erg s −1 cm 2 arcsec −2 (e.g.Hodge & Lee 1990;Azimlu et al. 2011;Cedrés et al. 2012), while the estimated mean sensitivity for the eDIG observed in other galaxies, such as NGC 891 or NGC 4565, is on the order of 10 −18 erg s −1 cm 2 arcsec −2 (e.g.Rand 1996;Rossa & Dettmar 2003b;Ho et al. 2016;Jones et al. 2017).For this work, the surface brightness detection limit of GHASP is 2.5 × 10 −17 erg s −1 cm 2 arcsec −2 (Epinat et al. 2008;Gómez-López et al. 2019;Sardaneta et al. 2022), which is shallower than needed for detection of the faintest eDIG, but still within an acceptable range of sensitivity to offer a reliable perspective for our study.
By discriminating the H emission in -direction from the galaxies displaying H emitting gas extended radially, almost all the galaxies in our sample showed H emitting gas extended radially outside from the stellar disc, CIG 1003 being the only exception, whose H emission is concentrated in the central regions drawing a cone shape and ionized gas filaments.In the vertical direction, the H monochromatic maps display eDIG configurations such as filamentary structures in CIG 71, 329, 593 and 1003; layers of diffuse gas in CIG 201, 329, 416; and patches or detached clouds of ionized gas emission in CIG 95, 159, 593, 922, 847 and 1003.In some cases, the ionized gas falls in the outskirts of the stellar disc tracing a warped structure likely from a recent interaction as in CIG 71 and 183.We have not detected any extraplanar H emission in the vertical direction in CIG 171, 906 and 936.
As shown in the top left panel of Figure 4, we compare the position in the DDD of our galaxies with those in the Rossa & Dettmar (2003a) sample which contains 62 spiral galaxies with inclination  > 76 • , an optical diameter of 3 ≤  25 ≤ 12 and radial velocity    ≤ 6000 km s −1 , all of them detected with IRAS in the FIR.We must consider that when defining extraplanar emission, Rossa & Dettmar (2003a) did not distinguish between radially or vertically extended H emission.In this case, all our galaxies would have eDIG.To ensure a more reliable comparison, in Figure 4, we differentiate between galaxies in our sample with and without vertically extended H emission.
Figure 4 shows that 6 out of 14 galaxies (42%) of our sample lie above the threshold of warm galaxies at  60 / 100 ≥ 0.4 (horizontal dashed line).Half of these galaxies, namely CIG 71, 936 and 1003, are located outside the region of bona fide isolated galaxies of the AMIGA sample (see Figure 2).In this context, we explore the possibility that eDIG incidence is correlated with the environment.In fact, Lisenfeld et al. (2007) found that the value of  60 / 100 of the AMIGA sample is lower than the one of interacting samples from the literature indicating that the interaction can increase the dust temperature.The panel (a) of Figure 4 also shows that our sample  (Rossa & Dettmar 2003a) showing the flux densities ratio at 60  and 100  ( 60 / 100 ) versus the ratio of the FIR luminosity (   ) divided by the optical diameter of the 25  ℎ mag arcsec −1 isophote squared ( 2 25 ) in units of 10 40 ergs s −1 kpc −2 .The horizontal dashed line marks the threshold for warm galaxies at  60 / 100 ≥ 0.4.We compare our sample (black circles with the CIG catalogue numbers indicated) with the whole sample (frame a) of Rossa & Dettmar (2003a) (squares) labelled RD03.Then, we use the hierarchy published in databases to distinguish whether the galaxies in the Rossa & Dettmar (2003a) sample are isolated (frame b), in pairs (frame c) or in groups (frame d) (see TableF1).Due to the low number of galactic pairs in the Rossa & Dettmar (2003a) sample, in panel (c) we also compare the mean values and confidence intervals (see Gehrels 1986) of this subsample (square) and that of our sample (circle).In all panels we distinguish between galaxies that have (filled markers) or do not have eDIG (empty markers).and the Rossa & Dettmar (2003a) sample share a similar distribution: galaxies with eDIG from both samples populate the entire plane.
Rossa & Dettmar (2003a) rejected galaxies in compact groups but included widely spaced pairs of galaxies.We associated an environmental classification to their sample using NED and Simbad 19 databases, by distinguishing groups, galaxy pairs and relatively isolated galaxies, whereby galaxies with siblings at a distance at  < 1 • were included in galaxies belonging to a galactic group and, galaxies with siblings at  > 1 • and galaxies with no record were considered isolated (see Table F1).These sub-samples are plotted in the (b), (c) and (d) panels of Figure 4. Rossa & Dettmar (2003a) detected eDIG in 5 out of 9 (56%) paired galaxies and in 21 out of 34 (62%) galaxies in groups: that is, in 26 out of 45 (58%) 19 "The SIMBAD astronomical database", (Wenger et al. 2000).
interacting galaxies, while, for the sub-sample of isolated galaxies, the eDIG was detected in 10 out of 19 (53%) targets.If we do not distinguish between radial and vertically extended emission, we detected H emission extended outside the disc in all galaxies in our sample.However, due to the inaccuracies that may exist in the definition of the stellar disc (see Section 4), we consider having detected eDIG only in 11 out of 14 galaxies (79%) where the diffuse ionized gas extends out of the disc vertically.
From the comparison with the entire sample of Rossa & Dettmar (2003a) we suggest that the environment is not connected with the presence of the eDIG.Confirming that the incidence of eDIG in late-type spirals is a direct consequence of the SFR in the underlying galactic disc as found by Rossa & Dettmar (2003a).The interaction could trigger the SFR for a short period of time.However, we also detect eDIG in a high percentage of galaxies that have been isolated for a long period of time.3) and ( 4), ( 5) and ( 6) and ( 7) and ( 8) averaged maximum radial () and vertical () distance reached by the extra-planar component at NUV, NUV and H emission, respectively; (9) morphological description of the extraplanar warm gas component or eDIG: E: extended emission in radial (r)-or vertical (z)-direction mostly, P: patches, F: filaments, H ii-R: extraplanar H ii region; (10) morphological features in different wavelengths: B: bar, A: asymmetric disc, W: warped disc, K: multiple inner emission knots, Ph: unmatching photometric maxima, UV: disc dominated by the UV halo, L(): Low emission detected in the wavelength band .T: thin disc structure.

The UV-halo
From the comparison between the H and the UV imagery, Thilker et al. (2007) discovered that 30% of low inclined galaxies ( < 80 • ) in the local universe present extended UV emission in the extreme outer disc where the number of H ii regions expected is negligible with respect to the UV emission regions.In the case of high inclined galaxies ( > 80 • ), extended diffuse UV emission up to 15-20 kpc above the disc midplane has been detected (e.g.Hodges-Kluck & Bregman 2014), which always coincides with the presence of the extraplanar H emission, although sometimes the ionized gas component is only detected at lower heights (Hodges-Kluck et al. 2016;Jo et al. 2018).
From the individual maps displayed in Figure 3 and Figures E1 to E13, we found that all galaxies in our sample present radial or vertically extended NUV emission, that all galaxies have radially extended FUV emission, and that 11 out of 14 galaxies show vertically extended FUV emission.These 11 galaxies displaying UV halo agree with the galaxies in our sample showing diffuse ionized gas extended out of the disc vertically (eDIG) (see Table 7).We also observed that in most galaxies in our sample presenting the ionized gas extraplanar component, the H and UV emissions have very similar morphology and the UV halo persists reaching larger heights where the ionized gas is no longer detected, as observed previously by Hoopes et al. (2001) and Hodges-Kluck et al. (2016).And, finally, even in galaxies where we did not observe any extraplanar (−axis direction) H emission (CIG 159, 171 and 936), we still have detected an UV halo in the NUV emission.
We found that 7 out of 14 (50%) galaxies of our sample showing a radially extended UV disc presented at least one feature of what could be perceived as a recent interaction such as disturbed ionized gas disc (CIG 71,95,171,183,416,593,847) and/or mismatch of the gaseous and stellar maxima (CIG 71,95,183,593,847) suggesting that the environment might be a potential cause of these events, even within isolated galaxies.From studies to low-inclined galaxies, it has been proposed that the extended UV disc might be a consequence of young stars associated with low-mass stellar associations located at large galactocentric distances (Gil de Paz et al. 2007b) which were likely triggered by recent/ongoing interaction events or a high specific rate of gas accretion (Thilker et al. 2007).However, it is assumed that isolated galaxies have not experienced gravitational influences from its close neighbours over the past few billion years (e.g.Karachentseva et al. 2010;Verdes-Montenegro et al. 2005;Rampazzo et al. 2016).Moreover, the lack of signals from a short-term event as a past major merger implies that the mechanism that supplies the gas needed to form the UV extension has been working for a prolonged duration to achieve this level of richness (Thilker et al. 2007).Studying UV photometry of cluster galaxies, Cortese et al. (2012) interpreted the extended UV emission as the result of the consumption of their H i reservoir, inferring that this growth could be stopped and even reversed if the atomic hydrogen is removed via some kind of environmental effect.Therefore, the observed UV extension could potentially be the result of scattering dust particles stemming from ongoing regions of star formation (see e.g.Shinn & Seon 2015), potentially introducing attenuation effects that shape the observed UV emission.
We observed that the UV halo (vertically extended UV emission) in galaxies CIG 201, 329 and 1003 does not trace a structured disc, these haloes are more extensive and have complex filamentary structures that match the H and NIR morphology.Moreover, the NIR, H and UV emission detected in CIG 201 and 1003 shows signs of a cone-shape suggesting the presence of a central outflow (see individual galaxy notes in the Appendix D).From the comparison between UV and NIR imagery, it has been shown that UV halo along the major axis relative to the optical radius  25 is unrelated to the prominence of the bulge and, therefore, the radially concentrated UV haloes are unlikely light from the bulge outskirts (Hodges-Kluck et al. 2016).In the case of CIG 201, 329 and 1003, the UV diffuse emission might be expected as a result of scattered UV light that escape from the disc dust (Diffuse extraplanar Dust, eDust) or emission from hot core-helium-burning stars (Hoopes et al. 2005).The simultaneity between the eDust and the extraplanar H emission suggests that a large fraction of the H emission could be originated from the galactic plane and is scattered by the eDust into the galactic halo (Hoopes et al. 2005;Seon et al. 2011;Jo et al. 2018).
On the other hand, the existence of an UV halo does not always mean the existence of eDust (e.g.Shinn 2019).We observed that in 10 out of 14 galaxies (71%) of our sample presenting an UV halo (CIG 71,95,159,171,183,416,593,847,906,922), the UV light reveals an intense brightness and a vertical extension above the galactic centre, which gradually decreases with increasing distance from the galactic centre, rather  2016) observed that the UV haloes are brightest over the parts of the disc with bright H, X-ray, or radio emission, suggesting that they are connected to star formation in the disc, but not necessarily in a way that requires ionizing photons or winds to escape the disc.
Previous studies about the origin of the eDIG have found regions consistent with mixed OB and hot low-mass evolved stars (HOLMES), and OB-shock ionization as the primary driver of eDIG ionization (Flores-Fajardo et al. 2011;Jones et al. 2017;Levy et al. 2019;Rautio et al. 2022).The NUV light is dominated by stars in the turn-off point in the HR-diagram (turn-off stars), while, the FUV emission of old 'normal' stellar populations is dominated by post asymptotic giant branch (PAGB) stars (e.g.Marino et al. 2011b).HOLMES, which include PAGB stars and white dwarfs, are found abundantly in the thick discs and lower haloes of galaxies.Their considerable vertical distribution compared to OB stars and their significant contribution to the UV radiation of galaxies (see e.g.Rautio et al. 2022, and references therein), make them candidates to be a source of eDIG ionization and, consequently, of UV haloes.
However, while late-type galaxies in nearby groups, observed from various inclination angles (ranging from 58°to 90°), display FUV and NUV images practically identical (e.g.Marino et al. 2010), our isolated galaxies, as previously noted in early-type galaxies, reveal that the FUV emission appears more compact and luminous compared to the NUV emission (e.g.Gil de Paz et al. 2007a;Rampazzo et al. 2007;Jeong et al. 2007Jeong et al. , 2009;;Marino et al. 2010Marino et al. , 2011a,b),b).Following Cortese et al. (2012), if the extended UV emission arises from the conversion of atomic hydrogen gas into new stars, it suggests that the stellar population is becoming younger or more metal-poor.As a result, the turn-off stars shift to bluer colours and higher luminosity, increasing the NUV emission.Simultaneously, the contribution to the FUV from PAGB stars decrease because, despite their higher luminosity, the duration of the PAGB phase gets much shorter due to decreased available fuel (Marino et al. 2011b).
In summary, the nature of UV haloes in galaxies are subject to a complex interplay of factors, including SF, ionization processes, and possible CGM or ICM contributions.Further investigation of the specific mechanisms driving UV haloes is needed to fully understand their origin and importance in different galactic environments.

SUMMARY AND CONCLUSIONS
Isolated galaxies constitute a reference model to study the environmental influence on galaxy evolution.In this work we studied the environmental impact on the incidence of eDIG in a sample of 14 nearby ( ≤0.02) isolated late-type edge-on ( ≥ 80 • ) galaxies from the Catalogue of Isolated Galaxies (CIG, Karachentseva 1973).We presented the H emission maps of galaxies in our sample with data obtained from the scanning Fabry-Perot interferometer, GHASP, offering a complete two-dimensional coverage of the emitting line region.We aimed to determine the galactic plane and to examine the distribution of the warm gaseous component with respect to the old and young stellar populations.Thus, we compared monochromatic H images with the NIR 2MASS K  -band and UV from GALEX archive images.
We started by defining the stellar disc using the ellipse fitted to the 3 on the background level of the 2MASS K  -band image, which delimits the region occupied by both old and young stellar populations.However, we are aware of selecting the boundary of the stellar disc might be another parameter causing discordance in the published results regarding the detection of eDIG.Thus, we consider having detected eDIG only in the 11 out of 14 galaxies (79%) where the diffuse ionized gas extends out of the disc vertically.If we do not distinguish between radial and vertically extended emission, we would have detected H emission extended outside the disc in all galaxies in our sample.In this case, we demonstrated that our sample is comparable to the broader sample studied by Rossa & Dettmar (2003a) when studying the incidence of eDIG.
We found that all galaxies in our sample present extraplanar NUV emission in both radial and vertical directions.While 11 out of 14 galaxies (79%) have vertically extended FUV emission, the same galaxies presenting eDIG.Because we showed that at least 79% (∼2/3) of isolated late-type, high inclined galaxies have both extended UV discs and UV haloes, our study reinforces the observations made by Thilker et al. (2007) who found that one-third of the low-inclined galaxies in the local Universe have extended UV discs.Hence, suggesting that the oldest star formation occurring between ∼10 and 100 Myrs extends well beyond the disc defined by the H map which only traces the most recent star formation (younger than 10 Myrs).
We find that the PA obtained with NIR, optical and NUV images varies on average ΔPA≃ ±2.5 • , which is expected for galaxies that have not experienced strong interactions in the last billions of years.However, when computing the inclination, some cases showed Δ ≃ ±10 • with respect to the optical inclination as a consequence of the presence of old stellar population centralized discs and UV extended discs and haloes.
Our study suggests that the eDIG is revealed in a significant fraction of galaxies (79% of the galaxies presented H emission vertically extended) which have been isolated from a long period of time.When compared with the generic sample of highly inclined galaxies by Rossa & Dettmar (2003a) we did not find evidence that the environment plays a role in the presence of the eDIG incidence.This confirms that the presence of eDIG in late-type spirals is a direct consequence of the SFR in the underlying galactic disc as suggested by Rossa & Dettmar (2003a).
We foresee to discuss the H kinematics of the galaxies in our sample in a forthcoming paper.Furthermore, since GALEX images are quite faint we are acquiring Astrosat-UVIT (Tandon et al. 2017) pointed observations of our galaxies in the FUV band.

APPENDIX A: SPATIAL SIGNAL ENHANCEMENT
To enhance the spatial signal-to-noise ratio (SNR), we employed two techniques on the wavelength data cubes.Firstly, we implemented adaptive spatial binning through Voronoi tessellation (Cappellari & Copin 2003;Daigle et al. 2006b,a), aiming for a target SNR of approximately 5 or 7. Additionally, we applied a Gaussian smoothing with a FWHM ≃ 2" (3 pix).To maintain a SNR threshold above 5 in the Gaussian smoothed data cube, we performed a comparison between both H monochromatic maps and we set a limit for the detected emission in the Gaussian smoothed maps.Since we aim to study the morphology of the eDIG, traced by the H line emission, we opted to utilize the H monochromatic map resulting from Gaussian smoothing.This map yields a more comprehensive representation of the morphology of the detected H emitting gas, in comparison to the map obtained through spatial binning via Voronoi tessellation.To illustrate this technique, Figure A1 presents the H monochromatic map for CIG 847 (UGC 11132).

APPENDIX B: CALIBRATION OF THE TOTAL H𝛼 FLUX
As indicated by Epinat et al. (2008), our major scientific goal is not to use the Fabry-Perot technique to make photometric studies but kinematic ones.Therefore, during the observations, we decided not to calibrate our data, saving observing time.On the other hand, there is no H calibrated images of galaxies in our sample in the literature yet to compare our targets which would provide us with a more accurate flux range.However, to obtain the total H flux of our galaxies, we can follow previous total H flux indirect calibrations made by Epinat et al. (2008) and by Gómez-López et al. (2019) for large galaxy samples observed with the same instrument, GHASP.Epinat et al. (2008) calibrated the total H flux of more than 200 galaxies from the GHASP sample by making a linear regression of the comparison between monochromatic FP images of 69 galaxies and images with calibrated fluxes from narrow-band filters including H and [N ii] which had been already published.The coefficient of the resulting line is 0.48±0.06×10−16 W m −2 ph −1 s, the fitting was forced to pass through the origin for obvious physical reasons.Later, Gómez-López et al. (2019) studied the kinematics of a sample of 152 late-type galaxies from the HRS catalogue 23 .In this context, they calibrated indirectly the FP monochromatic H maps, this time using the calibrated fluxes from Boselli et al. (2015), obtaining a similar coefficient of the resulting line fitted and a GHASP surface brightness detection limits of ∼ 2.5 ± 0.2 × 10 −17 ergs s −1 cm −2 arcsec −2 for 2 hours exposure time, resulting in objects with an average lower threshold of SNR≥ 5.This value of surface brightness detection limit was also obtained by Sardaneta et al. (2022) when, studying the Thus, from the results obtained for galaxies observed with the FP GHASP, we estimated the integrated flux of each galaxy in the sample studied in this work.We list these integrated fluxes in Table B1.

APPENDIX C: STELLAR MASS FROM MIR
Wide-field Infrared Survey Explorer (WISE) is a NASA-funded Medium-Class Explorer mission that consists of a 40 cm space infrared telescope, whose science instrumentation includes 1024 × 1024 pixel 2 arrays of Si:As and HgCdTe.WISE mapped the entire sky at 3.4 m ( 1 ), 4.6 m ( 2 ), 12 m ( 3 ), and 22 m ( 4 ).Each band covers a field-of-view of 47 × 47 arcmin 2 with an angular resolution of 6.0 arcsec in the short band-pass and 12.0 arcsec in the longest one and a 1.375 arcsec pixel scale (Jarrett et al. 2012).We infer the stellar mass using the mass-to-light (/) ratio of galaxies with  1 ( ⊙ ) = 10 −0.4(  − ⊙ ) , where  is the absolute magnitude of the source in  1 ,  ⊙ = 3.24 mag, and  1 −  2 is the rest-frame colour of the source.Mid-infrared magnitudes and stellar masses are listed in Table C1.

APPENDIX D: INDIVIDUAL GALAXIES
We provide the individual notes as discussed in Section 5.
CIG 95 (UGC 1733).Figure E1: The H map shows a detached cloud from the disc that exceeds radially the old stellar population disc observed in the NIR-band emission.The brightest H emission region is much wider than the old stellar disc and exceeds it in the z-direction.The 2.22m image reveals a central bulge surrounded by an asymmetric warped old stellar population disc that extends Southward.The FUV emission displays an extended and warped disc.The NUV emission traces a well-bound large disc with a brighter Southeastern tip and a sharpened, fainter Northwestern tip.The H, FUV and NUV emissions match in maxima emission but it is located to the Southeastern edge of the old stellar population disc.From the SDSS optical images, Buta et al. (2019) defined CIG 95 as a Sc warped (see Table 1).Although the galaxy appears to be slightly pulled to the Southeast, indicating a possible interaction, CIG 95 meets all the isolation criteria of the AMIGA sample (see Section 2).
CIG 159 (UGC 3326).Figure E2: The H emission displays a discontinuous disc composed of bright clouds uniformly distributed along the major axis with a size similar to that of the K  -band 1 isophote.CIG 416 (UGC 5642).Figure E7: The H emission map displays an asymmetric filamentary wide and radially extended disc with a tail-shape at the Western tip and a few detached ionized gas clouds.
The K  -band image of this galaxy shows a small disc dominated by a bright central bulge.Beyond the K  -band isophotal perimeter of the 3 level, the old stellar population disc is asymmetric with a filamentary structure and a tail shape pointing Northwards at the West extreme.The phenomenon triggering the Western tail could be the reason for the non-linear relationship the  4 parameter and the radius, however, neither the graphics nor the K  -band image show the presence of a bar.The NUV and FUV images show even more extended asymmetric discs and larger detached clouds than the H emission image.The images of these three bands present bright knots inside the disc.The photometric centre of the H, NIR, NUV and FUV emissions match quite well.The asymmetric disc is also observed at the optical wavelength since Buta et al. (2019) classified CIG 416 as Sdm peculiar using SDSS data.CIG 416 meets the isolation criteria of AMIGA sample in the optical band, nonetheless, it does not accomplish these criteria using radio band data (Jones et al. 2018).Therefore, the asymmetries, removed clouds, filaments and tails may be a signature of a past interaction.
CIG 593 (UGC 8598).Figure E8: Despite of a long exposure of 223 minutes (∼4 hours), we only detected a few faint H clouds tracing an asymmetric and truncated ionized gas distribution only visible on the North side, slightly displaced with respect to the major axis of the old stellar population.It looks an 'integral sign arm' starting from a small bar.Indeed, Buta et al. (2019) classified this galaxy as a strongly barred (SB, see Table 1).The stellar disc shows a central bright bulge.The most diffuse NIR emission displays an asymmetric and truncated disc.The extended NUV emission disc covers the entire galaxy and draws a warped disc as the one traced by the NIR emission.The FUV emission brightest knots at the North are located in the same position as the H clouds.The maxima of the NUV, FUV and NIR emissions are displaced by ∼10 arcsec in the galactic centre.CIG 593 meets all the AMIGA sample isolation criteria (see Section 2).
CIG 847 (UGC 11132).Figure E9: The H emission depicts a radially extended asymmetric and warped disc more prominent and brighter towards the South than to the North along the major axis with the maximum located towards the Southeast from the NIR photometric centre.The galaxy has a ionized gas tail-shape at the Southern tip and a detached cloud at the Northern edge.The 2.22m emission reveals an almost symmetric elongated disc with small warp.The FUV emission only shows large disengaged clouds along the major axis instead of a dense disc giving as a result an inaccurate surface brightness profile.The NUV emission traces a long and warped disc with feathering along the borders of the disc, populating mainly the Northwest and with the maximum immediately to the North of the NIR photometric centre.Both H and NUV emissions trace a tail shape at the Southeastern tip.The maxima of the NUV and NIR emissions match within ∼5 arcsec in the galactic centre, while the FUV and H emission maxima are located ∼10 arcsec and ∼15 arcsec Southwards from them, respectively.Using the H i image of CIG 847, Jones et al. (2018) found that this galaxy fulfils the most strict isolation criteria of AMIGA sample (see Section 2).
CIG 906 (UGC 11723).Figure E10: After an exposure of 186 minutes (∼3 hours) with FP GHASP, only faint H clouds as knots were detected along the galactic major axis.A faint ionized gas cloud is located over the Northern edge of the stellar disc.The K  -band image shows an old stellar population disc elongated and symmetrical with some filaments along the disc and a plume at the Southwest edge and a well-defined dust lane.The NUV and FUV emission display wider and brighter discs than the NIR one, both show filaments all along the edges of the disc and even some detached clouds.The surface brightness profiles of NUV and FUV emission show a discontinuity between the ellipses fitted to the 1 and 3 brightness levels of the K  -band image, this is from ∼30 and 45 arcsec.The maxima of the H, NUV and FUV emissions are displaced ∼10 ′′ Southwards from the NIR one.Despite Buta et al. (2019) classified CIG 906 as Sbc, our analysis to the surface brightness profile of the 2.22 m image showed a 'disc-like' shape.
According with Jones et al. (2018), from H i data, CIG 906 is a well isolated galaxy.
CIG 922 (UGC 11785).Figure E11: The H emission disc is asymmetric and warped, it dominates the Northern galactic region and it is composed by large detached clouds.The ionized gas disc is wider and longer than the old stellar disc.The H cloud furthest North is faint and exceeds radially the old stellar disc tracing an arrow-shape, while to the South there are two H detached bright knots just outside of the NIR disc edge.The old stellar population disc is mostly symmetric and elongated, it has a central bulge.The faintest NIR emission shows a warp at the Northeastern and a plume outside the Southwestern tip of the stellar disc, both peculiarities are pointing to the Southeast.The FUV emission traces a wider and denser disc than the H emission disc but both match in shape: the rounded Southern tip and the Northern arrow-shaped.The NUV emission displays the most extended disc, with large plumes at one side of the disc pointing to the Northwest, and a detached cloud at the Northern tip.The maxima of the four emissions matches within the central ∼5 arcsec.Buta et al. (2019) classified CIG 922 as Scd warped galaxy (see Table 1) and, from H i data, Jones et al. (2018) determined that this galaxy is not strictly isolated: the morphology of CIG 922 might be a result of some interaction.
CIG 936 (UGC 11859).Figure E12: The K  -band and most of H emissions were detected in the central 30 arcsec of a galaxy with an optical diameter of  25 () = 185.4arcsec (see Table 1 and Figure E12a).The H monochromatic image shows most of the warm gas compressed at the centre within the old stellar population disc, drawing a cone-shape.The K  -band image shows an asymmetric disc with a central bulge.The FUV emission draws an extended thin disc (isophotal radius  FUV = 47.6 arcsec, see Table 5) whose tips are warped pointing to the Northwest.The NUV disc emission is wider than the FUV one, with filaments mainly along the Southeastern side of the disc with the largest and broadest filament at the galactic centre.At ∼70 arcsec Southwards from the galactic centre, there is a NUV bright knot on the disc matching with the DSS R-band image and a faint H cloud, which is probably a SF region because it is not detected at the NIR emission.The maxima of the four emissions matches in the galactic centre.Jones et al. (2018) classified this galaxy as H i-not-isolated (see Section 2).
CIG 1003 (UGC 12304).Figure E13: The ionized gas is compressed in the central region of the disc which presents three bright knots, one at the centre and the other two symmetrically placed ∼15 arcsec away.The H emission map displays filaments of ionized gas in all the edges, but the most prominent of them are around the central knot tracing a cone-shape pointing to the Southwest.The K  -band image shows an elongated and symmetric stellar disc with a central bulge.
The diffuse NIR emission displays a cone-shape at the galactic centre pointing Southwest and a hook-shaped extension at the Northern tip pointing Northeast.The NUV emission image shows a prominent filament of ∼15 arcsec length pointing to the Southwest; in addition, the cone-shape in the central region is more noticeable in this emission.The FUV and NUV emission maps reveal very extended asymmetric and warped discs.Also, the brightest knots in these bands, match the Southern H knot.All the photometric centres defined by the ellipse fitting lay on the central H knot agreeing within ∼5 arcsec.CIG 1003 displays thinner discs Northwards and wider Southwards in the NIR, FUV, and NUV emissions.CIG 1003's multiwavelength morphology, position on the region of bona fide isolated galaxies (see Figure 2) and H i classification as not isolated (Jones et al. 2018), suggest that it has a galactic wind, possibly due to a recent interaction.F1) and their subsample of galaxies presenting eDIG compared with the Hubble galaxy type of our galaxy sample and our subsample of galaxies presenting eDIG.

Figure 2 .
Figure 2. Comparison between the local number density (  ) and tidal strength (  ) parameters for the CIG galaxies in our sample.These parameters were obtained first, based on digitized photographic plates from POSS-1 and POSS-2 (Verley et al. 2007b,a) (red triangles) and later, on available photometric and spectroscopic data from SDSS-DR9 (Argudo-Fernández et al. 2013) (blue circles).The CIG catalogue numbers are indicated on top the circles and below the triangles.According to Verley et al. (2007b), the horizontal (  = 2.4) and vertical (  = −2) lines enclose the region of bona fide isolated galaxies in the AMIGA sample.

Figure 3 .
Figure 3. Example of the layout of the graphs and maps of each galaxy.CIG 71 (UGC 1301).Multiwavelength maps of the galaxy: (a) the H monochromatic map from FP-GHASP data; (b) the 2MASS K  -band image; (c) the DSS R-band image; (d) and (e) the GALEX NUV and FUV images, respectively.In panels (f) and (g) the surface brightness profile of the K  -band and NUV/FUV emissions are plotted, respectively.In both panels, yellow vertical lines mark the surface brightness level at 1 (dotted line) and 3 (solid line) of the K  -band image.In panel (g), vertical lines indicate the radii corresponding to the isophotes at a surface brightness level at  NUV,FUV = 28 mag arcsec −2 for the NUV (blue dot-dashed line) and FUV (red dashed line) maps, respectively.Additionally, the vertical black loosely-dashed line at ∼5 arcsec shows the approximate FWHM of the GALEX point spread function (see e.g.Marino et al. 2010).The ellipses fitted to these isophotes are superimposed on the corresponding map with their associated colour.The common elements in in panels (a), (b), (d) and (e) are: a colour-bar indicating the flux in units of erg sec −1 cm −2 , the arrow tracing a scale of 30 arcsec, a black circle embedded in a box indicating the spatial resolution of the image, a black ellipse representing the surface brightness level at 3 on the background of the K  -band image and, a black cross (×) pointing the location of the NIR band peak light distribution.At the bottom of panel (f) we show the relationship between the coefficient  4 and the radius using the K  -band map to determine if the isophotes tend to have 'disc-like' ( 4 < 0) or 'boxy' ( 4 > 0) shape (see Section 4).Finally, in panel (a), on the H monochromatic map were overlaid the 1 isophotes on the background of the NUV and FUV images.

Figure 4 .
Figure 4. Diagnostic DIG diagram(Rossa & Dettmar 2003a) showing the flux densities ratio at 60  and 100  ( 60 / 100 ) versus the ratio of the FIR luminosity (   ) divided by the optical diameter of the 25  ℎ mag arcsec −1 isophote squared ( 2 25 ) in units of 10 40 ergs s −1 kpc −2 .The horizontal dashed line marks the threshold for warm galaxies at  60 / 100 ≥ 0.4.We compare our sample (black circles with the CIG catalogue numbers indicated) with the whole sample (frame a) ofRossa & Dettmar (2003a) (squares) labelled RD03.Then, we use the hierarchy published in databases to distinguish whether the galaxies in theRossa & Dettmar (2003a) sample are isolated (frame b), in pairs (frame c) or in groups (frame d) (see TableF1).Due to the low number of galactic pairs in theRossa & Dettmar (2003a) sample, in panel (c) we also compare the mean values and confidence intervals (see Gehrels 1986) of this subsample (square) and that of our sample (circle).In all panels we distinguish between galaxies that have (filled markers) or do not have eDIG (empty markers).

Figure A1 .
Figure A1.H monochromatic map of CIG 847 (UGC 11132) obtained after being applied a Voronoi tessellation (VT) with signal-to-noise ratio SNR ≃ 5 (top panel) and a spatial Gaussian smoothing with FWHM ≃ 2" (3 pix) (bottom panel).We use the VT binned map to mask the Gaussian smoothed map to ensure that all regions of this map have SNR ≥ 5.

Figure F1 .
Figure F1.Histogram of the Hubble galaxy type distribution of the whole galaxy sample of Rossa & Dettmar (2003a) (RD03) (see TableF1) and their subsample of galaxies presenting eDIG compared with the Hubble galaxy type of our galaxy sample and our subsample of galaxies presenting eDIG.

Table 2 .
Verley et al. (2007a)of galaxies in our sample and the distance to the nearest one measured byVerley et al. (2007a).

Table 4 .
Journal of observations  exp : total exposure time; (6) Seeing in arcs; (7) Final spatial resolution in kpc per arcseconds after having applied a spatial Gaussian smoothing on the H wavelength data cubes (see Section 3.1).

Table 5 .
Explanatory Supplement to the 2MASS All Sky Data Release and Extended Mission Products IV https://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html.Photometric parameters resulting from the ellipse fitting to the 2MASS K  -band and GALEX NUV and FUV images.
https://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec4_5e.htmlCIG galaxy name listed in column (1).Major axis length (), position angle (PA) measured from the North to the East and ellipticity () of the ellipses fitted to the K  -band image at 1 level on the background are listed in columns (

Table 6 .
Inclinations computed with different methods and wavelengths.

Table 7 .
Maximum radial () and vertical () distance reached by the extra-planar component in UV and H emission.

Table C1 .
Stellar mass and luminosity computed with the mid-and far-infrared magnitudes respectively.
Sabater et al. (2012) the galaxy shows a large thin warped disc with a central bulge.The GALEX FUV image shows a few faint clouds mostly at the photometric centre.The NUV emission draws a warped disc brighter on the East side than on the West side.The GALEX NUV and FUV images have low SNR which may yield spurious surface brightness profiles.The location of the H emission is consistent with the NUV emission.The emission maxima of the four emissions presented in FigureE2coincide within ∼5 arcsec.Further spectral data and deeper, higher resolution UV images are required to draw a conclusion about the morphology of this apparently quiescent galaxy.However,Jones et al. (2018)found that in the H i emission CIG 159 is a not an isolated galaxy FigureE3: The H emission is brighter in the South than the North likely shadowed by dust and stars emitting in UV.In the South the ionized gas has a dense distribution while in the North it traces detached clouds.The old stellar population disc is mostly symmetric with feathering along the edges.The NUV and FUV emissions trace a slightly warped disc brighter in the North than in the South of the NIR photometric centre.The UV surface brightness profiles are flat in the inner ∼60 arcsec and fall exponentially outward.There is no H emission over the stellar disc in the z-direction, but the H and UV discs are more extended radially than the old stellar population disc.The maximum emissions of NIR, H, NUV, and FUV agree within a ∼5 arcsec at the centre of the galaxy.Jones et al. (2018)classified CIG 171 as not isolated using H i emission data.FigureE4:The H emission displays an asymmetric disc displaced to the Northeast outwards of the old stellar population disc.It also shows two big bright regions located at the North of the photometric centre and several remnant clouds that draw two tails in both sides of the H disc.The 2.22m emission shows a mostly symmetrical disc with a bulge and two plumes at the ends of the galactic major axis.The NUV and FUV emissions agree with the structure of the H monochromatic map.The maxima of the NUV, FUV and 2.22m emissions match within a circle of ∼5 arcsec radius in the galactic centre.No additional spectral or photometric data are available yet in the literature to allow a conclusion about the morphology of this galaxy.inthecentralregion of the galaxy.A filamentary structure displayed by the ionized gas is consistently observed across all wavelengths.The K  -band image shows a small disc that is symmetric but outwardly warped, tracing a similar shape to the H and UV emission.The NUV, FUV, and H maps display extended discs that exceed more than twice the area occupied by the NIR emission.The location of the H bright knots match with the brightest knots seen in NUV and FUV emissions.The plateau shape of the NUV and FUV surface brightness profiles in the inner ∼40 arcsec may be due to the emission of these bright regions uniformly distributed within a specified area.All the maps of CIG 201, presented in FigureE5, outline the same filamentary structure and display a cone shape at the galactic centre pointing Eastwards reaching different heights above the stellar disc.From SDSS data,Buta et al. (2019)classified CIG 201 as a SAc galaxy and identified an inner pseudoring (see Table1).In addition, from H i emission data,Jones et al. (2018)determined that this galaxy is not interacting.The multi-wavelength morphology of CIG 201 suggests that the galactic disc may host a galactic wind.FigureE6: The H emission is distributed in a nearly symmetric ring structure.Indeed,Buta et al. (2019)classified it as a SAb galaxy identifying an inner pseudoring (see Table1).The 2.22m emission displays a long and symmetric disc with a central bulge and dust lane, with filaments along the edge.Following the photometric analysis of the K  -band image, this galaxy has a 'boxy' bulge ( 4 < 0).The brightest FUV emission regions agree with the H emission, while the faintest ones draw a filamentary warped disc and few detached clouds in the North.The NUV emission reveals a central dense distribution of the young stellar population, the NUV brightest clouds are consistent with the 2.22m emission and the faintest draw filaments along the disc and the Northern cloud traced by the FUV emission.Based on the SDSS spectroscopy,Sabater et al. (2012)classified CIG 329 as a LINER.CIG 329 satisfies all the AMIGA sample isolation criteria (see Section 2).
CIG 201 (UGC 3979).FigureE5:The H emission map reveals a concentrated distribution of ionized gas, forming bright knots predominantly located