abstract: With the advent of specialized hardware such as Graphics Processing Units (GPUs), large scale image localization, classification and retrieval have seen increased prevalence. Designing scalable software architecture that co-evolves with such specialized hardware is a challenge in the commercial setting. In this paper, we describe one such architecture (*Cortexica*) that leverages scalability of GPUs and sandboxing offered by docker containers. This allows for the flexibility of mixing different computer architectures as well as computational algorithms with the security of a trusted environment. We illustrate the utility of this framework in a commercial setting i.e., searching for multiple products in an image by combining image localisation and retrieval.
author: B Sengupta[^1]; E Vazquez; M Sasdelli; Y Qian; M Peniak; L Netherton; G Delfino
bibliography: large_TMM.bib
title: Large-scale image analysis using docker sandboxing

# Introduction

Large-scale image retrieval has been a mainstay for both academic research and commercial products for several years . A particular emphasis has been put for architectures that are scalable, redundant and most importantly quicker than its predecessors . This has resulted not only in efficient methods for indexing or hashing large databases but also automatic feature extraction using deep neural networks . Tasks that relied on Viola-Jones or SVM based classifiers to produce hand-crafted descriptors have been replaced by deep learning approaches since the ImageNet 2012 competition . The emphasis has lately shifted from "hand-crafted" feature descriptors to "hand-crafted" network architectures. Unlike statistical learning theory only a handful of studies have attempted to move on from this burgeoning resurgence of "deep-learning" architectures to understanding the basic principles on which such networks are based . Nevertheless, deep-neural networks have had immense commercial success.

Whereas efficient algorithms are unequivocally essential to achieve efficient classification, localisation and retrieval, commercial applications introduce new challenges. For commercial applications, a combination of efficient algorithms with a flexible and scaleable architecture is necessary. This paper describes one such architecture (*Cortexica*) for large-scale image localisation and retrieval in a cost-effective and time-efficient manner. Our aim is to present the architecture and tools needed to develop a fast, scalable and flexible system by focusing on a particular commercial application – a multi-product search algorithm. This application is explained in Section . The architecture described has proven to be adequate on a commercial environment (<https://www.cortexica.com/>). An important property of our system is that is not limited to particular algorithms detailed in this paper. Thus, the same architecture can be used with different algorithms for classification, localisation or retrieval.

The rest of the paper is organised as follows. Our approaches for localisation and retrieval are detailed in Sections and , respectively. Subsequently, in Section we explain the combined framework for localisation and retrieval to finally draw conclusions in Section .

# Multi-product search

It is typical in computer vision that a single image can have more than one item that needs to be identified; an additional aim can be to recommend a similar item from a vendor's database. Our approach is tailored for fashion items (clothes, accessories, etc.) yet the same approach can be followed for any type of objects, be it in medical imaging, transportation analytics, amongst many others. For example, we have deployed our framework to analyze streaming video – from walk-through at a fashion show. Our algorithm is able to identify all of the clothing items that a 'fashion model' wears; additionally, returning similar items from a vendor's database for each identified item.

One methodology to achieve such an aim is to first localise the items of interest, i.e. detect the items and find the associated bounding boxes, and then perform retrieval for each detected item. In the following section, we provide an overview of the different methods to perform object localisation as well as details on the architecture and the implementation for achieving good classification accuracy.

The first step of the multi-product search is object localisation. Deep learning has proven to yield high accuracy at relatively low computational cost when using GPU-based computing . In this section, three networks for object localisation are introduced, and compared. The most relevant property of our deployment is that it allows us to hot-swap the neural network used (in a production environment), making it straight forward to introduce neural network modifications or for that matter other non-neural network algorithms.

## Image localisation approach

The state-of-the-art neural networks for object detection that we describe below include: Faster-R-CNN (Faster Region-based Convolutional Neural Network) , SSD (Single Shot Multi-Box Detector) , and R-FCN (Region-based Fully Convolutional Networks) :
**Faster-RCNN**, a quasi-realtime object detection framework, is presented in . The predecessor Fast RCNN illustrated that a convolution neural network (CNN) framework can be successfully applied for object detection, significantly increasing the probability of detecting an object. As a first step, the Visual Geometry Group (VGG) network is trained for image classification , this is then further fine-tuned for object detection. Faster R-CNN includes a region proposal network (RPN), that builds upon Fast R-CNN enabling generation of fast proposals with little overhead. The features of the convolutional part of VGG are shared between two tasks – region proposal and classification. Just like the Spatial Pyramid Pooling (SPP) network , the convolutional part of the network works just with images that have fixed aspect ratios. In our experience, based on prior suggestions , training the classifier on a large number of 'hard' negatives strengthens the classification performance and leads to better precision.
**Single Shot Multi-Box Detector (SSD)** , uses an architecture that is equivalent to many class specific RPNs, each working on different feature maps. This to done to improve the detection of objects at different scales, thereby illustrating the benefits of a multi-scale architecture. SSD uses a fully convolutional layer and an optimized VGG architecture wherein the input aspect ratios are fixed at 500x500 (similar to the classical VGG network). This allows for a competitive execution time, fast enough to be real-time, without compromising on the accuracy of Faster R-CNN. The primary disadvantage of this framework lies in the detection of small objects.
The **Region-based Fully Convolutional Network** has an architecture similar to Faster-RCNN, albeit without the fully connected layer of the network. Not having the fully connected layer allows us to calculate in a single forward pass the loss function for considerably large number of region proposals. These regions are then quickly sorted according to their loss function values, retaining only those regions that supersede a fixed threshold. These are then utilized for gradient computation. This approach is called *online hard example mining (OHEM)* ; such a scheme is generally computationally infeasible to apply on a Faster-RCNN based network. In Faster-RCNN, for each region of interest (ROI), the network calculates a forward pass only when the images propagate to the fully connected layer. Therefore, if there are numerous ROIs, say more than a few hundred, it becomes increasingly time-consuming to train the network.

```latex
\begin{table*}
 \begin{center}
  \begin{tabular}{ l | l l l l l l l l | l }
   Method & time & jacket & dress & skirt & top & trousers & purse & shoe & mAP \\
   \hline Faster R-CNN end2end & 141ms & 83\% & 73\% & 75\% & 69\% & 84\% & 63\% & 85\% & 76.3\% \\
   SSD 300x300 & \textbf{30ms} & 83\% & 73\% & 75\% & 76\% & 85\% & 66\% & 86\% & 77.7\% \\
   SSD 500x500 & 70ms & 84\% & 74\% & 76\% & 76\% & 86\% & 68\% & 87\% & 78.9\% \\
   R-FCN ohem ResNet-50 & 96ms & 88\% & 78\% & 80\% & 77\% & 90\% & 68\% & 90\% & 81.6\% \\
   R-FCN ohem ResNet-101 & 137ms & 89\% & 78\% & 81\% & 78\% & 91\% & 69\% & 90\% & \textbf{82.3\%} \\
    \\
   \end{tabular}
  \end{center}
  \caption{The average precision calculated on each object class and the mean of the APs.
   The bounding boxes are considered only if the intersection-of-union (IOU) is larger than 50\%. The execution times are measured on an NVIDIA Quadro M6000 GPU with 24GB memory.}
   \label{table:timings}
  \end{table*}
```

Removing the fully connected section from Faster RCNN makes the training unstable. Hence, to guaranteee stability the ROIs are subdivided in a 7x7 grid and 49 losses are calculated for each category. Such a scheme allows each part of the grid to recognize features that are typical of that part of the object. Empirically, whilst one does not gain much on classification accuracy by training the network without the fully connected layer, this allows us to use the OHEM method.

To evaluate the performance of the three networks we use a dataset containing 45,000 Street Style images (Figure ; Table ) where seven fashion categories have been manually annotated with a bounding box (bb) around the object. The categories used are: jackets (18k bb), dress (9k bb), skirt (15k bb), tops (30k bb), trousers (13k bb), handbags (23k bb) and shoes (50k bb). The 45k images are split into two sets where 40k images are used for training and 5k images are used as a test set.

To train the networks we used the default parameters chosen in the original papers. Table shows the average precisions (APs) of the results of the models. Figure shows the object detection of the R-FCN ResNet101 model on one image of the validation set. The two newer methods (SSD and R-FCN) both improve over Faster R-CNN. SSD is especially suited when speed is the main concern. The smaller (300x300) achieves realtime performance. The bigger R-FCN is slower, but achieves a higher precision (Table ).

## Architecture/Implementation

The localisation service was designed to be *horizontally scalable* in order to handle a high number of requests per second. Such scalability was made possible by employing a load-balanced micro-service architecture. The micro-service architecture provides a method of developing software applications as a suite of small, modular and independently deployable services in which each service runs a unique process and communicates through a lightweight and well-defined mechanism.

The localisation micro-service is written in C++14 and uses POCO libraries (https://pocoproject.org) to implement a RESTful (Representational State Transfer) API. The service accepts a HTTP multipart request containing a query image and a JSON (JavaScript Object Notation) specifying various parameters. The image is then processed via a pre-specified deep-learning architecture (for example, a R-FCN network), and the localisation results are returned to a client in a JSON containing the classifications, confidences and coordinates of localised items.

A load-balancer is used to distribute the requests among any number of micro-service instances distributed across any number of GPU servers. Currently, we use servers with 4x NVIDIA GeForce GTX 1080, each featuring 2560 cores and 8GB of memory. Each server can host up to 16 instances (4 per GPU) while each instance can process approximately 5 requests/second resulting in a total of 80 requests/second per server. Additional instances can be dynamically deployed on more servers to cope with increasing throughput requirements. The deployment of new instances is straight forward since each micro-service runs from within a Docker container that guarantees that the software will always run the same, regardless of its environment. This is because each Docker container comes with the source-code, runtime executable, system tools and libraries that are needed to run the localiser micro-service.

At a very high-level, all queries are sent to a single load-balanced endpoint. The load-balancer then distributes these requests between a number of docker containers that is further distributed amongst a number of GPU-enabled servers. Each docker container runs a software load-balancing (*haproxy*), which further distributes the requests amongst the number of localiser instances running within a single docker. The maximum number of these instances depend on various constraints, such as GPU compute and memory capacity as well as the size of the model used for localisation. Different docker containers can run on the same or different GPUs as well as run different models for localisation. Various localisation services are made accessible via dedicated ports.

# Image retrieval approach

After having localised a specific object (e.g., coat, jacket, etc.) corresponding to one of the seven categories used for multi-product search, image retrieval is performed against a particular database of inventory items. This retrieval is performed in terms of *visual similarity*. The main factors that make clothes perceptually similar are their colour and texture properties. Such properties need to be extracted in a way that is relevant to human perception. In this context, simple approaches as colour histograms are not adequate. Furthermore, relevant features must be extracted in a compute and memory efficient manner.

Building upon prior work on psycho-physics and human neuroscience, our algorithm encodes the texture and the colour of the image in a similar way as the human brain does. This information is encoded and the matching against a large-scale dataset is performed in an efficient way. Two steps are required to perform retrieval. First of all we extract a signature for every image. Second, we match the signature of a query image to a large dataset of images.

## Architecture for signature extraction

In order to build a signature of an image, the following steps are required: detection of points of interest (key-points) in the image, feature extraction and encoding. The key-point detection aims at identifying the most important locations on the image. Features are extracted from the patches around the key-points, which capture the texture and colour characteristics of the patch. In order to have an efficient protocol to match the descriptors of an image against a large database of images, feature encoding is required. A bag of words model is used to convert the descriptors of all patches to codeword vectors. This gives a unique signature for each image, which is further bit-encoded to save memory space.

A point of interest in an image is a distinctive location in an image that can be robustly localised from a range of viewpoints, rotations, scales, and illuminations. In order to find these key-points in every image we use biologically-inspired non-linear orientation channels, as described in . The key-points are extracted using a pyramid structure, i.e. four scales are used for each image. These scales are used in order to capture information on different level of abstraction – a harmonic (multi-scale) representation for each image is constructed. For each key-point we also construct an associated saliency; this indicates how distinct that point is. A threshold is defined to select the most salient key-points, and a maximum of 512 keypoints is used.

For every key-point, a patch of a fixed size of 32x32 on each scale is selected; features are then calculated for each patch. The algorithm captures texture and colour information for every patch by applying specific filters to the image data. This enables it to identify texture characteristics and colour information. The processing after the extraction of the texture information includes weighting the texture data of each one of the colour channels with intensity data of the image patch. Texture and colour information is combined to obtain a descriptor which allows similar image patches to be identified.

Based on human visual perception, a set of steerable complex wavelet filters are used to extract texture features such as edges and lines/bars in different scales and orientations. If no features are identified, for example, if the image is monochromatic with each pixel intensity being the same, the image will still have a texture, albeit a smooth and uniform one. To capture the colour information, the CIE Lab colour space is used; this mimics the human brain's perception of colour. Particularly, each colour patch is first converted into CIE Lab colour space and normalized to a range of $[0,1]$. 32 band pass filters are applied for every colour channels generating a total of 96 filter responses. Texture weighted colour histograms are subsequently generated. In the end, we obtain a 576 dimensional colour-texture descriptor (3 colour channels x 6 bins x 4 scales x 4 directions x 2 values).

In particular, the texture intensity data describes not only the relative difference between different points in an image patch, but also the absolute intensity of the texture. This enables a descriptor for finding similar images, for example as that perceived by a human viewer, in terms of both texture, intensity and colour. For illustration we present an example scenario in Figure . The key-point and feature extraction is implemented using CUDA for compute time efficiency – the entire calculation takes on an average 220 ms.

In order to encode the feature vectors from the different patches in one signature, a codebook needs to be defined. This is formed from a large number of descriptors of various type of images. The centres of the clusters, which emerge from a $k$-means algorithm, are used to map the vector values. The feature vectors of each image are mapped to a signature using the codebook, which consists of 5000 centres. This results in a signature that is approximately 3.5kB after being bit-encoded. In order to compare the similarity between two images the ${\chi}^2$ distance between the signatures is calculated. The next section discusses how the matching is operationalized in an efficient manner.

## Architecture for matching

After the bag of words (BOWs) have been produced as described in Section , they are stored in a specific data structure that maximises search and retrieval times. In particular, all of the BOW files created with the customers' images are processed and an "Inverted Index" is produced. Here, we explain how this index is created and utilized at query time by a distributed in-memory data-grid.

To understand the structure of the "Inverted Index" we need to first understand the data that needs to be stored. In a single bag of words, for instance for the Colour-Texture descriptor, we usually have 2.5k words in a dictionary of around 5k. This means that if we were to store the BOWs in a big matrix, where on one axis we had the IDs of the images and on the other we had all the words, the matrix would have at least half of its entries being zeroes – a sparse matrix. The logical solution is to store the data not including the zeroes and use an extra array to "index" it. We choose to use an inverted indexing strategy, as is common practice in storing sparse matrices.

More precisely this means that the "Index" array will point to the beginning of each "word" in the data array. The data array is, in the most simple case, just a list of image-IDs, stored continuously for each word. By pointing at the inverted index, we can quickly find images that contain a specific word. Once we have a data array that is indexed by word, it is easy to distribute it across many machines, each of which carry a range of words. In this way, if we have enough instances, using an inverted index gives us the added benefit of making our system robust: in case one of the machines is faulty, the search would still scan through all the customers' images, and it would ignore just a small range of words present on the faulty server, thereby minimally affecting the query results. The size and distribution of the grid is fully customisable depending on the configuration and number of images one needs to store. Each query runs in parallel in all the grid instances, each of which uses multiple threads to access the index.

# Combining localisation and retrieval for multi-product search

In order to achieve a multi-product search based on one image or a video stream, the two methods described above are combined. The multi-product search runs sequentially – localisation runs first based on a deep learning architecture and subsequently retrieval based on colour-texture descriptors is initiated.

The input image gets processed by the micro-service for localisation. The output consists of the identified categories and the associated bounding boxes. There is some extra logic in place after the localisation takes place, which ensures that mutually exclusive items do not appear together. For example, if a dress, a top and a skirt are detected in the same area then the one(s) with the higher confidence are kept. This ensures that there are less false positives, and the results are more visually agreeable to the user.

Each cropped image around every detected bounding box is used as a query against the identified category. The final retrieval results for every category are returned. Using the same snapshot of the video in Figure , the category of each item localized (query object) is returned followed by items that are similar to this query object. Similar items obtained from two queries are shown in Figure .

For each potential category that can be detected by using the trained deep learning model, there is an associated database. Splitting the database according to the object type results in more accurate inference. On the other hand, performing retrieval in a very large database that includes all object types items deteriorates the performance. We use in this case 7 separate databases, one for each product, against which similar items can be found. This means that there are 7 different database IDs and associated inverted files. In our experiment each database includes 100 thousand to 2 million images. Certain categories, like tops, have a vast amount of data available, while for others, like purses, data is limited; this results in databases with varying sizes. The queries against the different databases are performed in parallel. On one hand this keeps the overall timing below 1 $s$ on the other hand this timing is not affected by the number of items detected. In addition, the timing of each query against each database is also approximately similar regardless of its size, since the retrieval is performed under a distributed architecture. Overall the process achieves high accuracy and compute time efficiency without requiring any user intervention.

# Conclusions

In this paper we have described *'Cortexica's multi-product search'*, a framework for large-scale localisation, classification and retrieval. This is being used for a wide variety of commercial applications – from fashion industry, health and safety critical applications to medical imaging. It harnesses software scalability using Docker containers whilst hardware scalability is achieved using GPUs that are deployed on a wide variety of cloud computing providers. Due to the hot-swappable nature of our implementation, one can not only use any number of compute efficient deep-learning framework but also non-parametric Bayesian classification and regression algorithms in near future. This gives us immense flexibility to choose the most efficient statistical and compute efficient algorithm for the application of interest. Our current research aims to have a similar architecture not just for image or video data-streams but also add an added layer of security for sensitive data by employing privacy-preserving machine learning techniques.

[^1]: firstname.lastname@example.com; BS has a dual appointment at Dept. of Bioengineering, Imperial College London.


author: M. Pohlen ; M. Balcells ; R. Lütticke ; R.-J. Dettmar
date: Received 23 December, 2003; accepted 04 March, 2004
subtitle: 3D-photometric thin/thick disk decomposition of eight edge-on S0 galaxies
title:  Thick Disks of Lenticular Galaxies [^1]

# Introduction

```latex
\begin{table*}
\begin{center}
{\normalsize
\begin{tabular}{ l  c c l  r@{.}l  c  c  c r@{.}l c }
\hline
\rule[+0.4cm]{0mm}{0.0cm}
Galaxy
&RA
&DEC
&RC3
&\multicolumn{2}{c}{T}
&Diam.
&$v_{\sun}$
&$v_{\rm vir}$
&\multicolumn{2}{c}{D}
&$M^{\scriptscriptstyle 0}_{\rm B}$ \\[+0.1cm]
%
&\multicolumn{2}{c}{(2000.0)}
&type
&\multicolumn{2}{c}{}
&[\ \arcmin\ ]
&[\ km\,s$^{-1}$]
&[\ km\,s$^{-1}$]
&\multicolumn{2}{c}{[Mpc]}
&[mag]\\
%
\rule[-3mm]{0mm}{5mm}{\scriptsize{\raisebox{-0.7ex}{\it (1)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (2)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (3)}}}
&\hspace*{0.1cm}{\scriptsize{\raisebox{-0.7ex}{\it (4)}}}
&\multicolumn{2}{c}{{\scriptsize{\raisebox{-0.7ex}{\it (5)}}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (6)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (7)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (8)}}}
&\multicolumn{2}{c}{{\scriptsize{\raisebox{-0.7ex}{\it (9)}}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (10)}}} \\
\hline\hline \\[-0.2cm]
NGC\,1596   &04 27 38.1&$-$55 01 40 &.LA..*/ &$-$2&0 &3.9 &1510 &1229   &17&1&$-$19.3\\
NGC\,2310   &06 53 54.0&$-$40 51 45 &.L..../ &$-$2&0 &4.1 &1187 & 943   &13&1&$-$18.6\\
ESO\,311-012&07 47 34.1&$-$41 27 08 &.S..0?/ &   0&0 &3.7 &1131 & 893   &12&4&$-$20.0\\
NGC\,3564   &11 10 36.4&$-$37 32 51 &.L...*/ &$-$2&0 &1.7 &2812 &2639   &36&7&$-$20.2\\
NGC\,3957   &11 54 01.5&$-$19 34 08 &.LA.+*/ &$-$1&0 &3.1 &1687 &1607   &22&3&$-$19.0\\
NGC\,4179   &12 12 52.6&$+$01 17 47 &.L..../ &$-$2&0 &4.1 &1248 &1277   &17&7&$-$19.6\\
NGC\,4521   &12 32 47.7&$+$63 56 21 &.S..0.. &   0&0 &2.6 &2500 &2767   &38&4&$-$20.2\\
NGC\,5047   &13 15 48.5&$-$16 31 08 &.L..../ &$-$2&0 &2.7 &6330 &6292   &87&4&$-$21.7\\
\hline
\end{tabular}
}
\caption[]{Global parameters of the observed lenticular galaxies:
{\scriptsize{\it (1)}} Principal name, {\scriptsize{\it (2)}} right
ascension, {\scriptsize{\it (3)}} declination,  {\scriptsize{\it(4)}}
RC3 coded Hubble-type, and the {\scriptsize{\it(5)}} Hubble parameter T
are taken from \cite{rc3}. The {\scriptsize{\it(6)}} diameter in arcminutes, the
{\scriptsize{\it(7)}} heliocentric radial velocities, and the B-Band absolute magnitude
{\scriptsize{\it(10)}} are taken from  LEDA.
According to the heliocentric radial velocities corrected for the Local Group
infall into the Virgo cluster {\scriptsize{\it(8)}} from LEDA, we estimated the
{\scriptsize{\it (9)}} distances following the Hubble relation with the
Hubble constant from the HST key project
of $H_{0}\!=\!72$ km s$^{-1}$Mpc$^{-1}$ \cite[]{hst_h0}.
\label{sample} }
\end{center}
\end{table*}
```

The knowledge of the detailed distribution of stars in galaxies is of fundamental importance to address the formation and evolution of those systems. To a first approximation, a disk galaxy can be described by a set of distinct stellar entities: a disk population, a bulge component, and a stellar halo. Deep surface photometry of external early-type galaxies and later elaborate measurements in our own Galaxy revealed the need for an additional component of stars. This was called 'thick disk' , since it exhibiting a disk-like distribution with larger scaleheight compared to the inner, dominating 'thin disk'.

There are three distinct families of hypotheses for its creation. The first group considers the thick disk as a separate entity produced in an early phase of enhanced star formation during the initial proto-galactic collapse . Another family of models regards the thick disk as an extension (by dynamical heating) of the thin disk. They assume that after the initial collapse all gas settles down into the galactic plane and starts forming stars. On this thin stellar disk a variety of constant or violent heating mechanisms could act: spiral density waves , encounters with giant molecular clouds , scattering by massive black holes , energy input by accretion of satellite galaxies , or bar bending instabilities . For example, recently used N-body simulations to show that tidal heating in a cluster is sufficient to thicken stellar disks by a factor of 2-3. This kinematic heating and vertical expansion will lead to a significant morphological transformation of a normal spiral galaxy into a lenticular. The third model suggests that thick disks are mostly made of debris material from accreted satellites. Recent cosmological N-body+SPH galaxy formation models of locating the thick disk formation before $z\!\approx\!1$ find that more than half of the thick disk stars are actually tidal debris from disrupted satellites. Therefore the thick disk is not a former thin disk thickened by a minor merger. To decide which of these hypotheses could explain the thick disk phenomenon best we need first a more general and complete statistic of thick disk properties. Naturally, these are rather global ones for external galaxies whereas our particular position in the Milky Way makes it possible to determine much finer details.

Since the work of Tsikoudi and Burstein it appears well known that thick disks are quite common in S0 galaxies. However, none of the more recent detections, except for two galaxies in and a short remark in , quantifying detailed parameters such as the ratio of thick to thin disk scaleheight or scalelength, is actually made in S0 galaxies. All these galaxies are of later Hubble type. In addition, we have not found a detailed 2D thin/thick disk decomposition for any S0 galaxy in the literature. Subsequent numerical decompositions dealing with S0 galaxies after the pioneering work in the early 80's treated the thick disk either as an outer flattened but exponential halo , or as a spheroidal bulge component .

The detections of possible halo or thick disk stellar light in disk galaxies of later Hubble type have been made in a scant handful of mostly nearby edge-on galaxies . Quite recently, suggest the detection of extended, ubiquitous thick disks in a large sample of late-type, edge-on galaxies by means of multi-colour imaging. However, their thick disks are solely detected by vertical colour gradients for which dust extinction complicates the interpretation. In addition, their vertical colour profiles, especially $(R-K)$, typically extend out to only very few vertical disk scaleheights. At those z-heights where they attribute the red colour to an additional component the thin disk may still be dominant over a potential thick disk and even determine the measured colour.

In this paper, we analyse a set of eight edge-on S0 galaxies using the classical approach for the identification of thick disks in external galaxies: the need for an additional disk component when attempting to fit single disk models to the light distribution in a deep image. Thereby we characterise the presence and properties (scalelengths, -heights, and central surface brightnesses) of thick disks by directly measuring their structure.

# Sample Selection & Observations

The data were taken as part of the PhD study on the radial structure on galactic stellar disks by . The galaxies were selected according to the allocated observing time, observatory, and CCD-chip size meeting the following morphological selection criteria. Using images from the Digitized Sky Survey (DSS) we verified that they are edge-on ($i~\rlap{\lower -0.5ex\hbox{$>$}}{\lower 0.5ex\hbox{$\sim\,$}}86$`<!-- -->`{=html} ), undisturbed, and similar to some disk-prototypical cases like NGC 4565 or IC 2531 to make it possible to consistently fit the applied simple disk model (cf. Sect. ). Galaxies with the following characteristics were rejected: spiral arms (indicating towards a lower inclination), a significantly asymmetric or disturbed disk (indicating towards strong interaction), and two-sided or significantly one-sided warped disks. In addition, galaxies which seemed to be dominated by the light of the bulge component, or appeared to be too small ($D_{25}\!~\rlap{\lower -0.5ex\hbox{$>$}}{\lower 0.5ex\hbox{$\sim\,$}}\!2$`<!-- -->`{=html}), or showed only a faint, patchy disk were also excluded. For the observed lenticulars we therefore selected galaxies with a symmetric, smooth disk and a distinct, non-dominant bulge component. It is not unexpected that half of these galaxies show a box- or peanut-shaped (b/p) bulge component, since find that $>\!40$% of all bulges are b/p shaped. Global properties of the finally observed eight S0 galaxies are given in Table . We do not argue that this sample is fully representative of the general population of S0 galaxies but it ensures the best prospects for obtaining consistent models with our 3D modelling technique for all galaxies.

The images (in Johnson R or V filter) were obtained in four observing runs in 1998/1999, three at the Danish 1.54m telescope of the European Southern Observatory (ESO, Chile) and one at the 1.23m telescope on Calar Alto (CAHA, Spain). During all three runs at the ESO the 1.54m Danish telescope was equipped with DFOSC and the C1W7/CCD which is a 2k x 2k LORAL chip providing a field size of $\approx\!13$`<!-- -->`{=html} and a scale of $\approx\!0.39$`<!-- -->`{=html} pixel$^{-1}$. The run at the Calar Alto 1.23 m telescope was done in service mode with the Site#18b chip, a 2048x2048 SITE CCD with 24 $\mu$m pixel size, providing an unvignetted field of $\approx 10$`<!-- -->`{=html} and a scale of $\approx 0.5$`<!-- -->`{=html} pixel$^{-1}$.

The standard CCD reduction techniques —overscan correction; subtraction of remaining large scale gradient in combined, oversan-subtracted, masterbias image; and careful flatfielding– were applied using the IRAF data reduction package. Neither the DFOSC nor the Calar Alto CCD R-band images were affected by fringing. The individual, dithered, reduced short exposures (150 s-600 s) were combined to the final deep image using IRAF's *imcombine* task. These images are rotated to the major axis using the smallest angle of rotation according to their true position on the sky. Table summarises the detailed observational parameters.

```latex
\begin{table*}
\begin{center}
\begin{tabular}{l c c c c c c }
\hline
\rule[+0.3cm]{0mm}{0.0cm}
Galaxy
&Filter
&Date
&Site
&Exp.time
&Seeing
&Coadds \\
 &&[mmyy]&&[min]&[\arcsec]&\#\,x\,[s] \\[+0.1cm]
%
\rule[-3mm]{0mm}{5mm}{
\scriptsize{\raisebox{-0.7ex}{\it (1)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (2)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (3)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (4)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (5)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (6)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (7)}}} \\
\hline\hline \\[-0.2cm]
NGC\,1596   & R & 0198 & ESO & 45.0& 1.3  &8x300, 2x150  \\
NGC\,2310   & R & 0198 & ESO & 51.7& 1.3  &600, 500, 200, 9x200 \\
ESO\,311-012& V & 0599 & ESO & 60.0& 1.6  &12x300 \\
NGC\,3564   & V & 0399 & ESO & 66.3& 1.3  &600, 500, 3x360, 6x300 \\
NGC\,3957   & R & 0198 & ESO & 30.0& 1.3  &3x600 \\
NGC\,4179   & V & 0599 & ESO & 65.0& 1.6  &1x600, 1x480, 1x360, 1x300, 9x240 \\
NGC\,4521   & R & 0699 & CAHA& 60.0& 1.5  &6x600 \\
NGC\,5047   & V & 0599 & ESO & 60.0& 1.5  &6x600 \\
\hline \\
\end{tabular}
\caption{Observing log for the individual combined images with
{\scriptsize{\it (1)}} the galaxy, {\scriptsize{\it (2)}} the
filter, {\scriptsize{\it (3)}} the observing date, {\scriptsize{\it (4)}}
the site, {\scriptsize{\it (5)}} the total coadded on-source
exposure time $t_{\rm int}$, {\scriptsize{\it (5)}} mean seeing conditions
during the observations, {\scriptsize{\it (7)}} the number of individual
images with their $t_{\rm int}$.
}
\end{center}
\label{observ}
\end{table*}
```

During the two ESO observing runs in 1999 several fields, partly enriched with additional stars provided by B. Skiff (priv. comm.), were observed. The standard fields were taken at least three times a night at different airmasses to determine the atmospheric extinction. During the other two observing runs no standard stars were taken. The ESO run in 1998 is calibrated by literature values and for galaxies without catalogued values interpolated according to the measured sky background. Only a rough zero point could be estimated for the Calar Alto run by comparing a galaxy also observed in another calibrated Calar Alto run. For more details about the photometric calibration we refer to .

# Extraction of the Disk Parameters

## 3-Dimensional Disk Model

We have developed a semi-automatic recipe to fit true 3-D single-component luminosity distributions to the 2-D data of edge-on galaxies and determine the galaxy parameters, such as scalelength and scaleheight in a physically meaningful way. Our method is described in detail in and and is only briefly recalled here. The disk model is based upon the fundamental work of . They tried to find a fitting function for the three-dimensional light distribution in disks of edge-on galaxies using the empirically determined exponential radial gradient, $I\!\propto\!\exp(R)$, and adding a description for the vertical distribution, $f(z)$, of the stars. The luminosity density distribution $\hat{L}(R,z)$ can be written as: $$\nonumber
\hat{L}(R,z) = \hat{L}_0 \ \exp{\left(-\frac{R}{h}\right)} \ f_n(z,z_0) \ {\rm H
}(R_{\rm co}-R)
\nonumber
\label{hatl}$$ $\hat{L}$ being the luminosity density in units of \[$L_{\sun}$ pc$^{-3}$\], $\hat{L}_0$ the central luminosity density, $R$ and $z$ are the radial resp. vertical axes in cylinder coordinates, $h$ is the radial scalelength and $z_0$ the scaleheight, and $n$ the index of the vertical distribution function. H$(x_0-x)$ is the Heaviside function and  is the cut-off radius characterising the observed outer radial truncations . To limit the choice of parameters we restrict our models to the three main density laws for the z-distribution ($\exp$, sech, and sech$^2$) following . Due to the choice of our normalised isothermal case $z_0$ is equal to $2 h_z$, where $h_z$ is the usual exponential vertical scale height preferred by many authors: $$\begin{aligned}
\nonumber
f_1(z) &=& 4\ \exp{\left(-2\ \frac{\mid z \mid}{z_0}\right)} \\
\nonumber
f_2(z) &=& 2\ {\rm sech} \left(\frac{2 z}{z_0} \right) \\
\nonumber
f_3(z) &=& {\rm sech}^2{\left(\frac{z}{z_0}\right)}
\end{aligned}$$ For any details about the numerical realisation we refer to . Therefore six free parameters ($i,\ n,\ \hat{L}_0,\ R_{\rm co},\ h,\ z_0$) fit the observed surface intensity on the CCD chip to the model. This model assumes that the vertical distribution is independent of position along the major axis as is known to be true in good approximation . The increase of scaleheight with galactocentric distance as reported by can be described as a combination of two disks each with constant scaleheight but different scalelength (cf. Sect. ). Any deviation from constant scaleheights should be visible in the vertical profiles overplotted by the models, which is not the case (cf. Appendix ). The apparent small vertical shift of some models in the lowest plotted vertical profile is due to the change in radial scalelength, as described in Sect. .

The possible influence on the six free parameters of the dust distribution, which was neglected during the fit, is estimated in . There we have shown that even for a worst case scenario (large optical depth $\tau_R$ and a radially and vertically fairly extended dust lane) our model is able to reproduce the input parameters with an error of a typical 20%. We expect this effect to be even less significant for the present sample of lenticular, dust-depleted galaxies. Deriving individual errors on all parameters is a complex task since the main source of error is not the numerical fitting procedure ($<\!1$%) but the systematic uncertainties in the process of fitting a rather simple, empirical model to real galaxies. An estimation of the errors is given in by changing the applied boundaries of the region used to fit the data to the model . They found differences in $h$ and $z_0$ of about $15\%$ and $\hat{L}_0$ varied about a factor of two. This is in the same range found by when they compared published values of scalelength measurements from different studies.

## Thin/thick disks

noted that the S0 galaxies in his sample are not well described by any combination of a single disk and another spheroidal component, such as a de Vaucouleurs $R^{1/4}$ bulge model . All S0 galaxies reveal a typical, continuous change of slope when one compares the major axis with parallel profiles above/below the plane. The profiles significantly flatten towards cuts higher above the midplane. According to the single component model $\hat{L}(R,z)$ all slopes should be nearly parallel. We infer from the images (cf. Appendix ) that all S0 galaxies show a kind of smooth outer envelope or highly flattened spheroidal component. This deviation from a normal shape cannot be explained by a bulge component. Any large $R^{1/4}$ bulge would have to be apparent on the major axis which is certainly dominated by an exponential (disk) component. It is worth mentioning that in the sample of there are four galaxies, NGC 3390 ($T\!=\!3$), NGC 3717 ($T\!=\!3$), NGC 4696C ($T\!=\!3.4$), and NGC 6504 ($T\!=\!2$)[^2] that showed a similar behaviour but were classified as late-type galaxies.

The outer component could be described as a 'thick disk' according to with a flatter slope than the inner thin disk, equivalent to a larger scalelength. However, the main characteristic of a thick disk is that the observed vertical profiles depart from the simple exponential model. Compared to the model, the outer parts are systematically brighter with increasing distance to the major axis. As mentioned by , the measured increase of the scaleheight with radius ('flaring') for early-type, edge-on disk galaxies can be understood if these galaxies have both thick disk scalelengths and scaleheights larger than for the dominant old disk.

## Fitting Method

To apply our single component fitting method described in the previous section we have to assume that there is a region in the galaxy dominated by only one of the components. If one uses exponential profiles it is obvious from the vertical and radial cuts that there are well separated vertical ranges (around the major axis and the outermost profiles) in which the light of one of the two disks dominates. Fitting thin and thick disk component simultaneously would be the desired approach. However, our single component model has already six free parameters. As shown in , fitting this model to observed data is a non-trivial task. The problem is the application of an idealized model itself, which obviously will not totally accurately describe the measured two-dimensional light distribution. This technique requires continuous human supervision to control the influence and quality of each individual parameter. This would most probably be even more difficult for a parameter set twice as large. We decided therefore to use an iterative fitting routine starting from outside-in since the thick disk clearly dominates the outer parts. The first step is to determine an initial estimate for the outer disk by restricting the region to be fitted and using our single component model. The next step is subtracting the derived full thick disk model from the original image and fit an inner disk to the residual by restricting the fitting range to the inner parts. Then we start from the beginning by subtracting this inner disk model from the original image and fit again the outer disk, now to this residual image. The initially pre-defined fitting regions for the thin and thick disk are sometimes adapted slightly after subtracting one of the components and before starting the second fitting round. As it turned out this process is remarkably stable for most of the galaxies. After one iteration the disk parameters are already the same within the range of one single model. The reason for this is the domination of the thick disk at large vertical z-heights ($~\rlap{\lower -0.5ex\hbox{$>$}}{\lower 0.5ex\hbox{$\sim\,$}}3.5\hbox{$z_0^{\rm n}$}\,\!=\!7h_z^{\rm n}$) measured here with high S/N, and the large radial range used to fit the thick disk component.

To restrict the number of free parameters we decided to prefix the inclination during the fitting process for this sample. For most of the galaxies the symmetric shape of the bulge and disk component indicates an exactly edge-on orientation. Only one of the S0 galaxies (NGC 3957) exhibits a dust-lane making it possible to estimate the inclination at $i\!=\!88\hbox{$^\circ$}\,\pm1\hbox{$^\circ$}\,$ following the method by .

The different density laws for the vertical distribution are similar for large z and only differ around the mid-plane of the luminosity distribution . In the actual choice of one fitting function is done individually for each galaxy depending on the measured profiles. However, near the plane of the galaxy the contribution of a thick disk is much smaller than that of to the thin disk component. First tests with a free choice of the fitting function for the vertical density distribution ($\exp$, ${\rm sech}$, or ${\rm sech}^2$) showed in the case of ESO 311-012 that the iterative fit became unstable for this reason. In the subsequent modelling we have chosen the intermediate fitting function $f_2(z)\!\propto\! {\rm sech}(2 z/z_0)$ in all cases for the thin and thick disk component. There was only one galaxy, NGC 2310, for which the ${\rm sech}$ function did not yield a satisfying convergency of the iterative fitting. However, switching to ${\rm sech}^2$ significantly improved the fit. Finally, only three free parameters are left for each disk ($\hat{L}_0,\ h,\ z_0$).

We want to point out that these decompositions are thus model dependent. An isothermal thick disk contributes less light to the thin/thick disk combination than to an exponential one and the shape of the thin disk will also be different. In addition, vertical scaleheights $z_0$ obtained with an $\exp$ model are systematically larger than in a ${\rm sech}$ model and again larger than in a ${\rm sech}^2$ fit. Depending on the vertical boundaries chosen for the fitting these differences could be more than $30\%$. The choice for the vertical distribution influences also the best-fit scale parameter ratio $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,$. used all combinations of either ${\rm sech}^2$ or $\exp$ models for NGC 4565 and derived $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,$ in the range of $4.3-5.4$.

In contrast with some of the previous 1-D-only fitting methods we do not simply use the deviation from a simple exponential on the minor axis, or parallel profiles; to measure thick disks we use the full 2-dimensional information and are therefore able to fit the radial scalelength of the thick disks.

## Breaks in Radial Profiles

One of the main difficulties while properly fitting our model to the S0 galaxies is their outer disk structure. This will be discussed in detail in and only briefly addressed here. As one can see from the figures in Appendix there are clear breaks in the outer parts of the radial profiles. These are similar to the truncations in more later-type galaxies . For all galaxies (except ESO 311-012) a similar break, slightly less pronounced, is also seen in other radial cuts parallel to the major axis. Our model, according to Eq. (), however, describes only an infinitely sharp truncation  . As shown in this implies a tight coupling between the radial scalelength $h$ and the cut-off radius  , when using our model fitting data with the observed breaks in the profiles. In addition, the sharply truncated model exhibits an intrinsic bending of the profile towards the outer parts . This complicates the visual quality control compared to the more flat infinite exponential model without any truncation. Therefore we decided to use the infinite exponential model, realised within the same fitting program by fixing  to ten times the radial scalelength. For our thick/thin fitting we restricted the fitting region to points inside the observed break radius. However, fitting an infinite exponential model to the intrinsic two-slope profiles is also affected by systematic errors . Depending on the ratio of the inner, shallow slope ( ) up to the break, to the steeper, outer slope ( ) beyond the break radius the determined scalelength will be systematically too small. Assuming a mean ratio of $\hbox{$h_{\rm in}/h_{\rm out}$}\,\!=\!4.4\!\pm\!1.7$ from a large, edge-on sample by makes it possible to quantify the expected offset. The best fitting scalelength  will be about 26% larger compared to the intrinsic radial scalelength .

Although the exact value of the scalelength is thus model dependent, neither $z_0$ nor $\mu_0$ are influenced by this problem. The derived scaleheights are independent of the scalelength and therefore robust, but of course still depending on the chosen density law ($\exp$, sech, or sech$^2$) for the z-distribution.

# Results

For six (75%) out of the eight S0 galaxies we were able to derive consistent thick/thin disk solutions (cf. Sect. ). The resulting best fit parameters are listed in Table and Table for the thick and thin disk, respectively.

```latex
\begin{table*}
\begin{center}
\begin{tabular}{lccc crc c crcc}
\hline
\rule[+0.4cm]{0mm}{0.0cm}
%
%
Galaxy
&Dist.
&Band
&\hbox{$\mu_0^{\rm k}$}\,
&\multicolumn{3}{c}{\hbox{$z_0^{\rm k}$}\,}
&
&\multicolumn{4}{c}{\hbox{$h^{\rm k}$}\,} \\
%
\cline{5-7}\cline{9-12}  \\[-0.15cm]
%
&\raisebox{0.0ex}{[Mpc]}
&\raisebox{0.0ex}{[$\lambda$]}
&\raisebox{0.0ex}{[mag$/$\raisebox{-0.4ex}{\hbox{$\Box^{\prime\prime}$}\,}]}
&\raisebox{0.0ex}{[\,\arcsec]}
&\raisebox{0.0ex}{[kpc]}
&\raisebox{0.0ex}{[\hbox{$z_0^{\rm n}$}\,]}
&
&\raisebox{0.0ex}{[\,\arcsec]}
&\raisebox{0.0ex}{[kpc]}
&\raisebox{0.0ex}{[\hbox{$h^{\rm n}$}\,]}
&\raisebox{0.0ex}{[\hbox{$z_0^{\rm k}$}\,]} \\
%
\rule[-3mm]{0mm}{5mm}{\scriptsize{\raisebox{-0.7ex}{\it (1)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (2)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (3)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (4)}}}
&\multicolumn{3}{c}{\scriptsize{\raisebox{-0.7ex}{\it (5)}}}
&
&\multicolumn{4}{c}{\scriptsize{\raisebox{-0.7ex}{\it (6)}}} \\
%
%
\hline \hline
&&&&&&&&&&& \\[-0.3cm]
NGC\,2310\footnotemark[1]
            &13.1&R&22.0&19.8 & 1.3& 2.4&& 62.8 &4.0$\,\;$& 2.6&3.2 \\
E\,311-012  &12.4&V&22.3&18.8 & 1.1& 2.6&& 41.5 &2.5$\,\;$& 1.9&2.2 \\
NGC\,3564   &36.7&V&21.9&22.0 & 3.9& 5.3&& 20.6 &3.7$\,\;$& 1.8&0.9 \\
NGC\,3957   &22.3&R&22.2&20.8 & 2.3& 3.4&& 33.9 &3.7$\,\;$& 1.6&1.6 \\
NGC\,4179   &17.7&V&22.9&37.3 & 3.2& 3.1&& 51.0 &4.4$\,\;$& 1.9&1.4 \\
NGC\,5047   &87.4&V&22.9&19.7 & 8.3& 3.4&& 39.4 &16.7$\,\;$& 1.7&2.0 \\
%
\hline
%
 mean         &    & &    &    &    &3.4&&     &   &1.9& 1.9\\
\hline
\hline
\end{tabular}
\caption{Results for the thick disk: {\scriptsize{\it (1)}} Galaxy,
{\scriptsize{\it (2)}} distance, {\scriptsize{\it (3)}} filter,
{\scriptsize{\it (4)}} central surface brightness (uncorrected for inclination) of the
thick disk model \hbox{$\mu_0^{\rm k}$}\,, {\scriptsize{\it (5)}} vertical scaleheight \hbox{$z_0^{\rm k}$}\,in
arcsec, parsec, and in units of the thin disk vertical scaleheight \hbox{$z_0^{\rm n}$}\,, and
the {\scriptsize{\it (6)}} radial scalelength \hbox{$h^{\rm k}$}\,in arcsec, kpc, in units
of the thin disk radial scalelength \hbox{$h^{\rm n}$}\,, and in units
of the thick disk vertical scaleheight \hbox{$z_0^{\rm k}$}\,\hspace{0.5cm} $^{\rm a}$ (used sech$^2$ instead of sech)
\label{thickres}
}
\end{center}
\end{table*}
```
```latex
\begin{table*}
\begin{center}
\begin{tabular}{lcc c r@{.}lc c ccc}
\hline
\rule[+0.4cm]{0mm}{0.0cm}
%
Galaxy
&Dist.
&Band
&\hbox{$\mu_0^{\rm n}$}\,
&\multicolumn{3}{c}{\hbox{$z_0^{\rm n}$}\,}
&
&\multicolumn{3}{c}{\hbox{$h^{\rm n}$}\,} \\
%
\cline{5-7}\cline{9-11} \\[-0.15cm]
%
&\raisebox{0.0ex}{[Mpc]}
&\raisebox{0.0ex}{[$\lambda$]}
&\raisebox{0.0ex}{[mag$/$\raisebox{-0.4ex}{\hbox{$\Box^{\prime\prime}$}\,}]}
&\multicolumn{2}{c}{\raisebox{0.0ex}{[\,\arcsec]}}
&\raisebox{0.0ex}{[kpc]}
&
&\raisebox{0.0ex}{[\,\arcsec]}
&\raisebox{0.0ex}{[kpc]}
&\raisebox{0.0ex}{[\hbox{$z_0^{\rm n}$}\,]} \\
%
\rule[-3mm]{0mm}{5mm}{\scriptsize{\raisebox{-0.7ex}{\it (1)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (2)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (3)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (4)}}}
&\multicolumn{3}{c}{\scriptsize{\raisebox{-0.7ex}{\it (5)}}}
&
&\multicolumn{3}{c}{\scriptsize{\raisebox{-0.7ex}{\it (6)}}} \\
%
%
\hline \hline
&&&\multicolumn{2}{}&&&& \\[-0.3cm]
%
NGC\,2310     &13.1&R& 19.6& 8&1 &0.5&& 24.0&  1.5 &  3.0 \\
ESO\,311-012  &12.4&V& 20.5& 7&2 &0.4&& 21.4&  1.3 &  3.0 \\
NGC\,3564     &36.7&V& 19.5& 4&2 &0.7&& 11.3&  2.0 &  2.7 \\
NGC\,3957     &22.3&R& 19.2& 6&0 &0.7&& 21.6&  2.3 &  3.6 \\
NGC\,4179     &17.7&V& 20.2&12&1 &1.0&& 26.2&  2.3 &  2.2 \\
NGC\,5047     &87.4&V& 20.7& 5&8 &2.5&& 22.6&  9.6 &  3.9 \\
%
\hline
 mean             &    & &     &  & &   &&     &     &  3.1\\
%
\hline
\hline
\end{tabular}
\caption{Results for the thin disk: {\scriptsize{\it (1)}} Galaxy,
{\scriptsize{\it (2)}} distance, {\scriptsize{\it (3)}} filter,
{\scriptsize{\it (4)}} central surface brightness (uncorrected for inclination) of the
thin disk model \hbox{$\mu_0^{\rm n}$}\,, {\scriptsize{\it (5)}} vertical scaleheight \hbox{$z_0^{\rm n}$}\,in
arcsec, and parsec, and the {\scriptsize{\it (6)}} radial scalelength \hbox{$h^{\rm n}$}\,
in arcsec, kpc, and in units of the thin disk vertical scaleheight \hbox{$z_0^{\rm n}$}\,
\label{thinres} }
\end{center}
\end{table*}
```

The resulting radial and vertical profiles overplotted by our best fitting model are shown in Appendix . The derived thin disk scalelength and scaleheight values are consistent with those found in studies by and , containing mostly galaxies of Hubble type later than Sa.

The mean ratio of thick to thin scaleheight for the five galaxies with the same vertical model combinations (sech$+$sech) for the thick and thin disk is $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,= 3.6 \pm 1.0$ and for the scalelength $\hbox{$h^{\rm k}/h^{\rm n}$}\,= 1.8 \pm 0.1$. Including also NGC 2310 we find $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,= 3.4 \pm 1.0$ and $\hbox{$h^{\rm k}/h^{\rm n}$}\,= 1.9 \pm 0.4$. We derive central surface brightnesses (uncorrected for inclination) in the range of $22.0~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\hbox{$\mu_0^{\rm k}$}\,~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}22.2$ R-mag$/$`<!-- -->`{=html} and $21.9~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\hbox{$\mu_0^{\rm k}$}\,~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}22.9$ V-mag$/$`<!-- -->`{=html}for our thick disks compared to $19.2~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\hbox{$\mu_0^{\rm n}$}\,~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}19.6$ R-mag$/$`<!-- -->`{=html} and $19.5~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\hbox{$\mu_0^{\rm n}$}\,~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}20.7$ V-mag$/$`<!-- -->`{=html}for the thin disks. This implies that the contribution of the thick disk to the central surface brightness is about 10% of that of the thin disk. The thick disk central luminosity density $\hat{L}_0^{\rm k}$ ranges between 3.5% and about 10% (mean: 5.6%) of the thin disk value. The ratio of the total luminosities of the thick and thin disk is between about one third and one.

The profile on the minor axis in addition to radial cuts high above the major axis (cf. Appendix ) reveals that for four (NGC 2310, NGC 3564, NGC 4179, NGC 5047) out of the six fitted galaxies there is no significant bulge component visible at large vertical height above the disk. Even the apparently bulge dominated galaxy NGC 3564 is well described by a thick disk component. The remaining bulge light after subtraction of the thin/thick disk combination is comparable to that of the other galaxies. This implies that all bulges of our S0 galaxies could not be well described by a traditional de Vaucouleurs $R^{1/4}$ bulge. In addition, any bulge component would be too flat to account for all the light high above the disk at large galactocentric radii.

## Fitting regions

<table>
<caption>Radial and vertical fitting regions for the thick/thin disk components with <span><span><em>(1)</em></span></span> galaxy name, <span><span><em>(2)</em></span></span> beginning and end of thick/thin disk vertical region, and <span><span><em>(3)</em></span></span> beginning and end of thick/thin disk radial region </caption>
<tbody>
<tr>
<td style="text-align: left;">Galaxy</td>
<td colspan="4" style="text-align: center;">vertical region <span class="math inline"><em>z</em></span></td>
<td style="text-align: center;"></td>
<td colspan="4" style="text-align: center;">radial region <span class="math inline"><em>R</em></span></td>
</tr>
<tr>
<td style="text-align: left;"></td>
<td colspan="2" style="text-align: center;">thin</td>
<td colspan="2" style="text-align: center;">thick</td>
<td style="text-align: center;"></td>
<td colspan="2" style="text-align: center;">thin</td>
<td colspan="2" style="text-align: center;">thick</td>
</tr>
<tr>
<td style="text-align: left;"><span>2-5</span></td>
<td style="text-align: right;"></td>
<td style="text-align: left;"></td>
<td style="text-align: right;"></td>
<td style="text-align: left;"></td>
<td style="text-align: center;"></td>
<td style="text-align: right;"></td>
<td style="text-align: left;"></td>
<td style="text-align: right;"></td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;"></td>
<td colspan="4" style="text-align: center;"></td>
<td style="text-align: center;"></td>
<td colspan="4" style="text-align: center;"></td>
</tr>
<tr>
<td style="text-align: left;"><span><span></span></span></td>
<td colspan="4" style="text-align: center;"><span></span></td>
<td style="text-align: center;"></td>
<td colspan="4" style="text-align: center;"><span></span></td>
</tr>
<tr>
<td style="text-align: left;"></td>
<td colspan="4" style="text-align: center;"></td>
<td colspan="4" style="text-align: center;"></td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">NGC 2310</td>
<td style="text-align: right;"><p>0</p></td>
<td style="text-align: left;">12</td>
<td style="text-align: right;">20</td>
<td style="text-align: left;">35</td>
<td style="text-align: center;"></td>
<td style="text-align: right;">32</td>
<td style="text-align: left;">75</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">79</td>
</tr>
<tr>
<td style="text-align: left;">ESO 311-012</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">8</td>
<td style="text-align: right;">8</td>
<td style="text-align: left;">31</td>
<td style="text-align: center;"></td>
<td style="text-align: right;">29</td>
<td style="text-align: left;">88</td>
<td style="text-align: right;">37</td>
<td style="text-align: left;">117</td>
</tr>
<tr>
<td style="text-align: left;">NGC 3564</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">4</td>
<td style="text-align: right;">22</td>
<td style="text-align: left;">32</td>
<td style="text-align: center;"></td>
<td style="text-align: right;">20</td>
<td style="text-align: left;">32</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">52</td>
</tr>
<tr>
<td style="text-align: left;">NGC 3957</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">4</td>
<td style="text-align: right;">16</td>
<td style="text-align: left;">28</td>
<td style="text-align: center;"></td>
<td style="text-align: right;">36</td>
<td style="text-align: left;">56</td>
<td style="text-align: right;">24</td>
<td style="text-align: left;">64</td>
</tr>
<tr>
<td style="text-align: left;">NGC 4179</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">12</td>
<td style="text-align: right;">31</td>
<td style="text-align: left;">55</td>
<td style="text-align: center;"></td>
<td style="text-align: right;">43</td>
<td style="text-align: left;">69</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">88</td>
</tr>
<tr>
<td style="text-align: left;">NGC 5047</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">6</td>
<td style="text-align: right;">18</td>
<td style="text-align: left;">25</td>
<td style="text-align: center;"></td>
<td style="text-align: right;">35</td>
<td style="text-align: left;">47</td>
<td style="text-align: right;">0</td>
<td style="text-align: left;">55</td>
</tr>
</tbody>
</table>

One of the important constraints for fitting empirical, surface photometric models to observed data is the definition of the actual regions which are marked to characterise the individual model components. Therefore we list in Table for each galaxy the radial and vertical ranges where our thick/thin disk fit was applied. In all cases we restricted the two fitting regions to be distinct from each other and outside the dust lane (in the case of NGC 3957) or any bar/ring like feature visible in the radial profile. This is obvious for the first thick disk estimation. However, in principle one could extend the fitting range towards the inner/outer parts for the following thick/thin disk iterations, since the other component in each case has already been subtracted. In addition, we masked by hand bright stars and background galaxies, forcing the program to ignore these regions.

## Comparison with Literature

```latex
\begin{table*}
\begin{center}
\begin{tabular}{l l c c c  c  c  l}
\hline
\rule[+0.4cm]{0mm}{0.0cm}
%
%
Galaxy
&Type
&Model
&\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,
&\hbox{$h^{\rm k}/h^{\rm n}$}\,
&\hbox{$\mu_0^{\rm k}$}\,[mag$/$\raisebox{-0.4ex}{\hbox{$\Box^{\prime\prime}$}\,}]
&\hbox{$\mu_0^{\rm n}$}\,[mag$/$\raisebox{-0.4ex}{\hbox{$\Box^{\prime\prime}$}\,}]
&Reference  \\[+0.1cm]
%
&
&
&
&
&{\scriptsize \raisebox{-0.9ex}{$\stackrel{\rm R}{\rm V}\:$}}-band
&{\scriptsize \raisebox{-0.9ex}{$\stackrel{\rm R}{\rm V}\:$}}-band
& \\
\rule[-3mm]{0mm}{5mm}{\scriptsize{\raisebox{-0.7ex}{\it (1)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (2)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (2)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (3)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (4)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (5)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (6)}}}
&{\scriptsize{\raisebox{-0.7ex}{\it (7)}}} \\
\hline \hline
&&&&&& \\[-0.3cm]
5 galaxies\footnotemark
&S0
&sech$+$sech
&$2.6\leftrightarrow5.3$
&$1.7\leftrightarrow1.9$
& \raisebox{-0.7ex}{$\stackrel{22.2}{\scriptstyle 21.9\leftrightarrow22.9}$}
& \raisebox{-0.7ex}{$\stackrel{19.2}{\scriptstyle 19.5\leftrightarrow20.7}$ }
& this study \\[+0.1cm]
NGC 4710& S0 &sech$+$exp&3.2&&&& DE96 \\[+0.1cm]
NGC 4762& S0 &sech$+$exp&4.6&&&& DE96 \\[+0.1cm]
5 galaxies
&S0
&$-$
&1.8$\leftrightarrow$4.6
&
&
&
& DE97 \\[+0.1cm]
NGC\,6504
&Sab
&$\exp+\exp$
&4.0
&
&\raisebox{-0.7ex}{$\stackrel{23.8}{\scriptstyle - }$ }
&\raisebox{-0.7ex}{$\stackrel{18.2}{\scriptstyle - }$ }
&VA94 \\[+0.1cm]
NGC\,891
&Sb
&$\exp+\exp$
&2.3$\leftrightarrow$6.3
&
&
&
& MO97\\[+0.1cm]
NGC\,891\footnotemark
&Sb
&$R^{1/4}+$sech$^2+\exp$
&3.0
&
&
&
& VA84\\[+0.1cm]
NGC\,4565\footnotemark
&Sb
&sech$^2+$sech$^2+$halo
&2.2
&1.4
&
&
& WU02\\[+0.1cm]
NGC\,4565
&Sb
&sech$^2+$sech$^2$
&4.6
&1.4
&\raisebox{-0.7ex}{$\stackrel{23.7}{\scriptstyle - }$  }
&\raisebox{-0.7ex}{$\stackrel{19.6}{\scriptstyle - }$  }
& SH89 \\[+0.1cm]
NGC\,4565
&Sb
&$\exp+\exp$
&5.4
&1.1
&\raisebox{-0.7ex}{$\stackrel{23.0}{\scriptstyle - }$  }
&\raisebox{-0.7ex}{$\stackrel{18.8}{\scriptstyle - }$  }
& SH89 \\[+0.1cm]
MW\footnotemark (optical)&Sbc
&$\exp+\exp$
&3.0         &1.3         &         &       & LA03\\[+0.1cm]
MW\footnotemark[4] (NIR)   &Sbc
&$\exp+\exp$
&3.3         &1.3         &         &       & OJ01\\[+0.1cm]
ESO\,342-017
&Scd
&$\exp+\exp$
&2.5
&~\rlap{\lower -0.5ex\hbox{$>$}}{\lower 0.5ex\hbox{$\sim\,$}}1.0
&\raisebox{-0.7ex}{$\stackrel{22.1}{\scriptstyle22.5 }$ }
&\raisebox{-0.7ex}{$\stackrel{20.0}{\scriptstyle20.5 }$ }
& NE02\\[+0.1cm]
IC\,5249
&Sd
&$\exp+\exp$
&3.0
&0.6
&
&
& AB99\\
\hline
\end{tabular}
\caption{
%
Comparison with literature:
 {\scriptsize{\it (1)}} Galaxy name
{\scriptsize{\it (2)}} morphological type, {\scriptsize{\it (3)}} ratio
of thick to thin disk scaleheight, {\scriptsize{\it (4)}} ratio
of thick to thin disk scalelength, {\scriptsize{\it (5)}} central surface
brightness of thick disk (if quoted), {\scriptsize{\it (6)}} central surface
brightness of thin disk,
{\scriptsize{\it (7)}} reference DE96:\cite{degrijs1996}, DE97:\cite{degrijs1997}, VA94:\cite{vandokkum1994}, WU02:\cite{wu2002}, MO97:\cite{morrison1997}, VA84:\cite{vdk84}, SH89:\cite{shaw1989}, LA03:\cite{larsen2003}, OJ01:\cite{ojha2001}, NE02:\cite{neeser2002}, AB99:\cite{abe1999} \label{comparison} \newline
{\footnotesize
%
{$^{\rm a}$} NGC\,2310 excluded (only sech$^2$ model for thick disk)
\hspace{0.43cm} {$^{\rm c}$} three component model (thin disk, thick disk, and halo)   \newline
{$^{\rm b}$} three component model (bulge and thin plus thick disk)
{$^{\rm d}$} starcount analysis
}
%
}
\end{center}
\end{table*}
```

only describes the 'thick disk' component qualitatively as being more diffuse than that of to the inner, dominating disk, and possessing a flattened shape. The important observation is that the scaleheight of a fitted exponential to the vertical profiles increases with radial distance $R$. He determines a flattening, defined as the intrinsic z-thickness compared to the diameter $a$, of $z/a\!=\!0.25 - 0.2$ at $\mu_{\rm B}\!=\!25.0$ which is somewhere between thin disks ($z/a~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}0.1$) and E4 ellipticals ($z/a\!=\!0.6$) for five edge-on S0 galaxies. We find values of 0.29 (NGC 2310) or 0.32 (NGC 5047) comparing minor to major axis diameter at the second outermost contour (cf. Fig. ,).

While comparing scale parameter ratios in the literature one has to keep in mind that they are often obtained with very different fitting methods and even different fitting functions for the vertical distribution of thin and thick disk light. As discussed in Sect.  this already implies systematic differences of at least $20\%$. Surprisingly, there are no detailed parameter studies for thick disks of S0 galaxies in the literature that provide both scaleheight and scalelength ratios (cf. Sect. ). The only scaleheight ratios for S0 galaxies (cf. Table ) are given for two galaxies (NGC 4710, NGC 4762) in detail by and for five galaxies with Hubble type $T\!~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\!0.0$ without providing the individual values or applied fitting functions by . They find ratios of $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,\!=\!1.8-4.6$, which falls well within the range of values for our S0 galaxies. There are clearly more studies in the literature that provide detailed parameters (scaleheights and often also scalelengths) for thick disks in galaxies of later morphological type. modeled excess light at large vertical distance one-dimensionally with a thin/thick disk combination for the edge-on Sab galaxy NGC 6504. They found that the ratio of scaleheights is roughly 4, slightly higher than the mean value for our S0 sample. However, their central surface brightness for the thick disk is significantly lower ($\Delta \mu_0\!=\!1.6$ mag) than our value for NGC 3957. This effect is not related to the different fitting functions used since our sech model should tend to result in systematically lower values for $\mu_0$ than in their $\exp$ model. For the Milky Way, recent optical star-count measurements by yield thin disk values of $\hbox{$z_0^{\rm n}$}\,\!=\!300$ pc and $\hbox{$h^{\rm n}$}\,\!=\!3.5$ kpc, and for the thick disk $\hbox{$z_0^{\rm k}$}\,\!=\!900$ pc and $\hbox{$h^{\rm k}$}\,\!=\!4.7$ kpc. These are similar to the infrared 2MASS star counts by with $\hbox{$z_0^{\rm n}$}\,\!=\!260$ pc, $\hbox{$h^{\rm n}$}\,\!=\!2.8$ kpc, $\hbox{$z_0^{\rm k}$}\,\!=\!860$ pc, and $\hbox{$h^{\rm k}$}\,\!=\!3.7$ kpc, keeping in mind that the near-infrared surveys always derive systematically smaller scalelengths than in the optical. The resulting ratio of thick to thin disk scaleheight is very similar to the mean value we find for our S0 galaxies but their ratio of the scalelengths is also smaller. present a list of thick disk local normalizations (relative to the solar neighbourhood) obtained in different studies including their own measurements (cf. their Table 2). The quoted values range from 2% up to 13% with a mean value of about 6.1%, which agrees very well with the mean value of 5.6% for the central luminosity density of our S0 galaxies. used one-dimensional vertical fits to derive the thick and thin disk scaleheights of NGC 891, an edge-on Sb galaxy similar to our Milky Way. They quote a large range of possible $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,$ values which are consistent with our S0 disks. Previous modeling of NGC 891 by yield a ratio of $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,\!=\!3.0$ for a three component (thin plus thick disk, and additional $R^{1/4}$-spheroid) fit which is also in the range of our S0 galaxies. The slightly higher value of the thick disk central luminosity density ($17\%$ of the thin disk) compared to our maximum of $10\%$ could be explained by the additional contribution of the applied $R^{1/4}$-spheroid, especially high above/below the midplane. However, already pointed out that in the case of NGC 891 a combination of thin plus thick disk and a model with only one disk and a $R^{1/4}$-spheroid both give a good description of the data used by . For another prominent edge-on Sb galaxy, NGC 4565, several attempts for a decomposition are available. More recently, presented a thin/thick plus halo decomposition from a deep intermediate-band (6660Å) image deriving values only slightly smaller than for our S0 galaxies (cf. Table ). However, the two component disk models of NGC 4565 by yield significantly higher values for the ratio $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,$. The light put down to a halo is here ascribed solely to a thick disk with a large scaleheight. There is a large variety of model combinations (with different components) possible , which complicates a direct comparison. NGC 5907, another nearby edge-on Sc galaxy, is an even more complex case. It is not yet clear if the multiple detected extended light distribution is a thick disk or a stellar halo component . even concluded that NGC 5907 does not have a faint extended halo at all (but compare with discussion in for this issue). obtained a deep optical image of the edge-on Sd galaxy IC 5249. They detected additional light to that predicted by a single exponential disk and tried to fit a thick disk model. The scaleheight ratio of their best thick disk is well within the range for our S0 galaxies. However, the scalelength of the thick disk is exceptionally smaller compared to the thin disk and their central luminosity density ($\hat{L}_0^{\rm n}/\hat{L}_0^{\rm k}\!=\!7.4$) is noticeable larger than for our sample. Recently, reported for the first time the detection of a thick disk in a low surface brightness galaxy (ESO 342-017 classified as Scd). They find for the scaleheight a ratio of $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,\!=\!2.5$ (close to our values for NGC 2310 and ESO 311-012) with a comparable or somewhat larger scalelength for the thick disk. The scaleheight ratio and unprojected central surface brightness (cf. Table ) are surprisingly similar to the range of values found for our S0 galaxies.

# Discussion

## Thick disks: discrete or continuous?

A key question is: *How to describe a thick disk in general?* With reference to the proposed different formation scenarios described in Sect.  we can assume that the disk component is either characterised as a superposition of two discrete and independent isothermal disk systems (as done here), or built from the contribution of multiple velocity-dispersion components . These two approaches on how to treat a thick disk seem to be incompatible. Therefore any derived parameters for the two-disks model appear useless for the multi-component disk. In addition, the latter seems to be superior since it is consistent with the model of continuous disk heating leading to the well known observable age-velocity dispersion relation .

However, already states in his review that detailed studies of the spatial and kinematical distribution of stars in the Milky Way do not make it possible to decide if the thick disk is a discrete component or a more continuous sequence of stellar populations. Recently, concluded that the latest studies argue for the separate entity picture. This confirms our result that 75% of the chosen S0 galaxies are well described with a distinct, two-component thin/thick disk system. Note that these results exclude all thick disk formation scenarios based solely on heating. Especially the elemental abundance trends found by favour a merger scenario where a satellite galaxy either merges with the parent galaxy or sheds significant amounts of its material to form the thick disk as proposed by .

However, as already noted by it does not seem possible to distinguish between our simple model and more sophisticated ones. Therefore the very good description with the applied discrete two-component system for our galaxies is only a first step towards studying the luminosity distribution of external galaxies. In addition, one has to keep in mind that our model only fits well in a restricted radial range (only out to $2.2$`<!-- -->`{=html} $-3.9$`<!-- -->`{=html} ) where most of the inner part is "hidden" by the central bulge, bar, or ring components. One key to the nature of the thick disk may lie in studying the outer parts where the breaks in the radial profiles of both thin and thick disk should provide an additional constraint. Any vertical colour gradients could provide additional information but the only available colour map for our sample suffers from low S/N and is not conclusive.

## Are thick disks of 0s exceptional ?

The fact that the range of thick disk parameters for all known S0 galaxies are not too different from those of late-type galaxies, even compared to a low surface brightness Scd galaxy, is especially surprising, since at first glance one expects S0 galaxies to posses more prominent thick disks. However, in terms of the thick-to-thin-disk scaleheight ratio our values agree well with those derived for all other galaxies (cf. Table ). Even taking into account the central surface brightness or central luminosity density the comparison yields similar values. Does this point to a general formation process for thick disks independent of bulge-to-disk ratio and Hubble type ? At this stage we are not able to answer this question. Although the numbers of galaxies used for the comparison are the same, all literature values come from different sources, using sometimes very different methods and models to derive their parameters. In the case of NGC 4565 it is obvious that depending on the model, one can find an even larger variety of $\hbox{$z_0^{\rm k}/z_0^{\rm n}$}\,$-values than the full range for our six galaxies. In addition, as pointed out by , and again discussed in detail by , fitting disk-model components to surface brightness data of edge-on galaxies is a delicate business, and comparing the results of different authors can be misleading. To overcome this problem one has to extend this survey of thick disk parameters also to late-type galaxies. To reduce the inevitable influence of their prominent dust lanes this has to be done in the near-infrared. After applying the same fitting method one is able to address the questions if thick disks are really common around all Hubble-types, and if their parameters are really similar, as suggested here. This would entail a common formation scenario for thick disks independent of their normal evolution along the Hubble sequence.

Of our galaxies with successful thick disk models there are only two for which there are rotational velocity measurements in the literature. Therefore any correlation between mass and thick disk parameters could not be derived. Consequently it is still unclear if the mass, as a galaxy characteristic, is related to the thick disk parameters.

## Thick disk scalelengths

We find significantly larger scalelength ratios ($\hbox{$h^{\rm k}/h^{\rm n}$}\,\!=\!1.9$) for the thick disks of our S0 galaxies than in the literature ($\hbox{$h^{\rm k}/h^{\rm n}$}\,\!~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\!1.4$). Is this a possible distinction in thick disk parameters between late-type and our early-type galaxies ? Again, differences in fitting methods, especially for the scalelength in edge-on galaxies, could be responsible for this apparent disagreement. In particular, fitting profiles with clear breaks using a single disk (with an infinitely sharp truncation or infinite exponential as done here) will entail systematic errors (cf. Sect. ), therefore this question could only be addressed unambiguously applying exactly the same method to late-type galaxies.

However, although the exact number is uncertain the scalelength of the thick disk is without doubt systematically larger than that of the thin disk. Except in one case this is also true for all the literature values. Does this larger scalelength imply a different formation process or even restrict the available formation scenarios ? At a first glance the different scalelengths contradict a dynamical heating scenario, since this should only alter the vertical distribution. However, assuming two distinct disk components with radially constant thickness the velocity dispersion scales with radius . It depends therefore on the exact way (radial distribution) the proposed mechanisms (cf. Sect. ) dispenses the energy in the vertical motions of the stars. In the case of heating the disk by satellite accretion, the N-body simulations of show that the scalelength of the disk is nearly unchanged (even slightly smaller than the initial scalelength) for the inner parts (out to $3h_{\rm ini}$). Due to the migration of material outwards in radius, at larger radius (out to $\approx\!6h_{\rm ini}$) the final disk shows a second shallower component (cf. their Fig. 4). However, find in their N-body simulation a global outward transport of disk material leading to a general increase of the disk scalelength of 10% to 60%. Although this increase is below our measured scalelength ratios, the latter seems to coincide in general with observations of larger thick disk scalelengths. One has to keep in mind that in the satellite accretion scenario a pre-existing thin disk gets heated. The thin disk to be observed today must be rebuilt out of gas remaining after the merger process. It is not guaranteed that this thin disk will have the same scalelength as its predecessor since we do not have a unique explanation of how the disk scalelength is determined out of an initial (or new) gas distribution .

Looking from a slightly different angle, one can also try to apply the larger thick disk scalelength as a general argument against an internal heating scenario. Assuming infinitely exponential disks, the larger scalelength implies a larger angular momentum content of the thick disk. Any valid heating mechanism must therefore add angular momentum, which is not the case for the internal heating scenarios. However, note that the disks are not at all infinitely exponential but show clear breaks in their profiles (cf. Sect. ). The disks exhibit the break at roughly similar radial radius (e.g. Fig. ). In this sense the thick disk is truncated "earlier" (in respect to its scalelength) than the thin disk. In addition to a probably lower rotational velocity this could add up again to a similar angular momentum content for thick and thin disk, ruling out the whole argument against the heating scenarios. The key point here is again the origin of these breaks in the radial profiles.

## Residuals

In addition to comparing thick disks across a large range of characteristic parameters of galactic disks, fitting and subtracting a combined disk yields valuable information on the structure and size of galactic bulges. As discussed in Appendix the residual images highlight the deviation of the galaxy from the disk model. Structures only faintly indicated in the profiles become obvious. One important result related to this is that for our sample of S0 galaxies none of the present bulges could be described as a traditional early-type $R^{1/4}$ bulge in agreement with . However, a detailed structural analysis of the remaining central structures has to be done with great care. As discussed in Sect. , choosing the vertical model is not a unique process and therefore any variation of $f(z)$ for the thick and thin disk could alter the shape of the disk profile in the inner bulge parts. Especially for galaxies where the bulge component along the minor axis extends vertically above the disk (ESO 311-012 and NGC 3957) this could be even more unconstrained.

# Comments on Individual Galaxies

**NGC 1596:** (Fig. ) A 'normal' S0 enclosed in a huge, spheroidal, outer envelope possibly triggered by interaction with the nearby physical companion (NGC 1602). The difference from other 0s is obvious from radial cuts. It is the only galaxy in which the profiles bend upwards (not caused by a background gradient) at larger galactocentric distances ($R\!=\!\pm115$`<!-- -->`{=html}). This looks similar to an R$^{1/4}$ bulge component beginning to dominate the profile when extrapolating the curvature of the inner bulge towards the outer extension. However, the overall flatness of such an R$^{1/4}$ bulge would be unusual. In addition, the contours start to show substructure in the very outer parts (in contrast to the smooth contours of ellipticals). The inner light profiles are dominated by bulge and a possible inner disk component. The major axis shows several distinct parts: an inner bulge, a steep straight decline out to $R\!\approx\!\pm80$`<!-- -->`{=html}, an asymmetric slightly shallower part out to $R\!=\!-110$`<!-- -->`{=html}  and $R\!=\!+120$`<!-- -->`{=html}, followed by an upwards bending profile an upbending profile. The SW-side is slightly more extended than the NE-side. After the outer thick disk has been subtracted, the remaining profiles clearly resemble still the thin/thick disk behaviour (changing of slopes towards outer/higher radial cuts), therefore no consistent model could be derived.

**NGC 2310:** (Fig. ) exhibits a box-shaped bulge . The bulge/bar is rotated against the outer disk, therefore we find a slightly tilted residual. The radial profiles show a complex structure with inner, small, very bright nuclear bulge component, and further out a significant dip in the profile (caused by a possible dust ring or mimicked by a subsequent stellar ring), sitting on top of a bar/ring-like structure $R\!\approx\!\pm50$`<!-- -->`{=html}. There is a narrow disk-like structure $R\!\approx\!\pm82$`<!-- -->`{=html}, here assumed to be the main underlying disk, finally showing a straight decline into the noise (NE-side disturbed by nearby star-halo). The best fitting model overestimates the inner disk component due to inner ring-like feature. This feature is responsible for the apparent holes in the residual image, which shows the inner box-shaped bulge in addition to a large ring or bar.

**ESO 311-012:** (Fig. ) Box-shaped bulge . The galaxy position is close to the Galactic plane with many foreground stars overimposed. The bulge component is still visible on the minor axis high above the plane. The radial profile shows a bar-like structure at $R\!=\!\pm20$`<!-- -->`{=html}, followed by an exponential disk component out to a break at $R\!\approx\!\pm83$`<!-- -->`{=html}. The residual is a box-shaped bulge with the indication of a weak bar extending on both sides.

**NGC 3564:** (Fig. ) The galaxy is probably part of small group/cluster with at least two large, nearby physical companions (NGC 3568 and NGC 3557). Taking into account the low surface brightness part it is different from the other 0s. Out to $\mu_{\rm V}\!=\!23.4$ mag$/$`<!-- -->`{=html}it looks similar although with a more dominating bulge. The outer contours, on the other hand, reveal that NGC 3564 is sitting in a huge, roughly spherical stellar envelope, which has probably been built up by a merger with a smaller satellite galaxy still visible in this structure on the NE-side ($(R,z)\!=\!-47\arcsec,+24$`<!-- -->`{=html}, distance unknown). The radial profile shows an inner bulge-dominated part, followed by a bar-like plateau $R\!\approx\!\pm10$`<!-- -->`{=html}, sitting on top of an inner disk structure out to $R\!\approx\!\pm30$`<!-- -->`{=html}which is followed by a straight decline into noise. The residual resembles a rather large, spheroidal inner bulge together with a bar or ring-like structure extending to both sides. It nicely shows the possibly merging satellite. The thick disk component takes most of the light from the envelope and therefore the remaining bulge light is comparable to to that of the other galaxies.

**NGC 3957:** (Fig. ) Galaxy cluster in background. Nearly box-shaped bulge . It is an S0 galaxy with unusual, inner ($R\!~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\!\pm30$`<!-- -->`{=html}) dust lane. The bulge component is still visible on the minor axis high above the plane. The radial profiles show inner bulge, followed by a ring- or bar-like component at $R\!\approx\!\pm30$`<!-- -->`{=html}(visible in cuts above the dust lane), with break at $R\!\approx\!\pm60$`<!-- -->`{=html}  and followed by a straight decline into the noise. The profile probably starts to bend outwards from $R\!\approx\!\pm105$`<!-- -->`{=html} ($\mu_{\rm V}~\rlap{\lower -0.5ex\hbox{$>$}}{\lower 0.5ex\hbox{$\sim\,$}}24.3$ mag$/$`<!-- -->`{=html}) to $R\!\approx\!\pm150$`<!-- -->`{=html}. The residual resembles a nearly box-shaped bulge with an indication of a bar extending to both sides. The apparent holes are caused by absorption in the inner dust lane.

**NGC 4179:** (Fig. ) The inner bar/bulge component is rotated against the outer disk. The major axis shows a bright inner bulge, a disk- or bar-like feature between $R\!\approx\!\pm40$`<!-- -->`{=html}($z\!~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\!\pm5$`<!-- -->`{=html}), a more shoulder-like decline out to $R\!\approx\!\pm80$`<!-- -->`{=html}, here taken to be the main underlying disk, followed by a straight decline into the noise. The residual resembles a spheroidal bulge with an indication of an inner disk or bar extending slightly to both sides.

**NGC 4521:** (Fig. ) The inner component is $\approx\!1$`<!-- -->`{=html} rotated against the outer disk. The radial profile shows a bar- or ring-like structure between $R\!\approx\!\pm 20$`<!-- -->`{=html}out to $z\!~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\!\pm5$`<!-- -->`{=html}; on top of a two-slope disk, inner part out to a break radius of $R\!\approx\!-60$`<!-- -->`{=html}, followed by a steeper straight decline out to $R\!\approx\!-77$`<!-- -->`{=html}and a slightly upwards bending part at the very edge (perhaps flatfield problem). The only galaxy with inconsistent thick/thin disk fitting without obvious reason (cf. NGC 1596). All trials, such as changing vertical fitting function, have failed. When various thick disk models have been subtracted, the resulting thin disks always contributes too much light in the outer/higher z profiles dominated by the thick disk. Because of this the following thick disk is flattened, producing an unstable situation, perhaps caused due to lower image quality, proximity to chip edge, and disturbed by bright star on the NW-side. In principle, the galaxy exhibits the same features in radial/vertical cuts as our other 0s.

**NGC 5047:** (Fig. ) Peanut-shaped bulge . There are many companions around associated with the foreground NGC 5044 group of galaxies ($\Delta v\!\approx\!3600$ km s$^{-1}$). Outer disk has a slightly s-shaped warp, on the SW-side towards N, on the NE-side towards S. The inner disk part is slightly off-centered to the outer component. The radial profiles show a bar-, or ring-like feature between $R\!\approx\!\pm 30$`<!-- -->`{=html}out to $z\!~\rlap{\lower -0.5ex\hbox{$<$}}{\lower 0.5ex\hbox{$\sim\,$}}\!\pm3$`<!-- -->`{=html}, on top of a shoulder-like inner disk out to $R\!\approx\!\pm 50$`<!-- -->`{=html}, followed by a straight decline out to $\mu_{\rm V}~\rlap{\lower -0.5ex\hbox{$>$}}{\lower 0.5ex\hbox{$\sim\,$}}25.0$ mag$/$`<!-- -->`{=html}. The residual shows a nice peanut-shaped bulge with a thin bar-like component. The ratio of the bar length to the maximum of the peanut distortion closely matches simulations explaining the peanut by dynamical processes in the bar potential .

# Isophote maps, radial/vertical profiles, and residual images

The following figures show in the *upper panel* an isophote map of each galaxy rotated to the major-axis. A small arrow in the upper right corner indicates the North. The magnitude of the outer contour ($\mu_{\rm lim}$, defined by a 3$\sigma$ criterion of the background) is indicated in each plot. The consecutive contours are equally spaced by 0.5 mag. The contour lines are plotted with increasing smoothing towards the outer parts. For the inner contours out to where the noise begins to increase no smoothing was applied. The following 2-3 contours are smoothed by replacing each pixel by the mean of $3^2$-pixels "around" the central pixel. For the outer two contours this smoothing is increased to $5^2$-pixels. The *second panel* displays the major-axis surface brightness profile (top) and two parallel radial profiles each above and below the major axis. The exact $z$-positions for the plotted profiles are indicated in the upper left corner of the plot. For a consistent thin/thick disk decomposition the model profiles are overplotted as dashed lines. The *third panel* displays the minor-axis surface brightness profile (top) and two parallel vertical profiles on either side of the minor-axis. Again the exact $R$-positions for each cut are marked in the upper left corner and the best model is overplotted with dashed lines. The *lower panel* is the isophote map of the residual image after subtracting the best fitting thin/thick disk model. The magnitude of the outer contour ($\mu_{\rm lim}$) is again indicated in each plot and the following contours are equally spaced by 0.5 mag.

[^1]: Based on observations obtained at the European Southern Observatory, Chile

[^2]: The classification is done by . It is consistent with that of , but significantly different from that of the UGC catalogue .