Overview of Akatsuki data products: definition of data levels, method and accuracy of geometric correction

We provide an overview of data products from observations by the Japanese Venus Climate Orbiter, Akatsuki, and describe the definition and content of each data-processing level. Levels 1 and 2 consist of non-calibrated and calibrated radiance (or brightness temperature), respectively, as well as geometry information (e.g., illumination angles). Level 3 data are global-grid data in the regular longitude–latitude coordinate system, produced from the contents of Level 2. Non-negligible errors in navigational data and instrumental alignment can result in serious errors in the geometry calculations. Such errors cause mismapping of the data and lead to inconsistencies between radiances and illumination angles, along with errors in cloud-motion vectors. Thus, we carefully correct the boresight pointing of each camera by fitting an ellipse to the observed Venusian limb to provide improved longitude–latitude maps for Level 3 products, if possible. The accuracy of the pointing correction is also estimated statistically by simulating observed limb distributions. The results show that our algorithm successfully corrects instrumental pointing and will enable a variety of studies on the Venusian atmosphere using Akatsuki data.


Introduction
The Venus Climate Orbiter, Akatsuki, is a Venus weather satellite launched on May 21, 2010 (Nakamura et al. 2011(Nakamura et al. , 2016. The Akatsuki mission observes cloud distributions at several altitudes on both the dayside and nightside over five wavelength regions selected to obtain cloud distributions at different altitudes. Its main scientific purpose is to contribute to the understanding of the atmospheric circulation and cloud formation processes of Venus. To obtain meteorological information, global three-dimensional maps of the clouds are taken by four cameras at ultraviolet and infrared wavelengths; lightning is detected by a high-speed imager, and vertical atmospheric structures are observed by a radio-occultation technique. Since Akatsuki is in an equatorial elongated orbit with a period of about 11 Earth days, complementary to the polar orbit of Venus Express, which had a period of an Earth day, cloud motions and/or temporal variations at altitudes of about 45-70 km can be well observed continuously over a long period. Systematic imaging observations by Akatsuki are widely expected to enable the three-dimensional visualization of Venusian atmospheric dynamics. In this paper, useful information about the Akatsuki data products is presented. The atmospheric super-rotation of Venus is one of the most remarkable phenomena in planetary science. Venusian rotation is very slow, having a period of about 243 Earth days, such that an axisymmetric subsolarto-antisolar (SS-AS) circulation was expected to be predominant (e.g., Stone 1968;Dickinson 1969). However, a number of observations have shown that a fast zonal wind is predominant in the Venusian atmosphere below ~ 80 km altitude (e.g., Schubert 1983). Almost the entire atmosphere of Venus rotates much faster than Venus itself; the zonal wind increases almost linearly with altitude from the ground and reaches about 100 m s −1 at the cloud's top level (~ 70 km altitude), which is about 60 times faster than Venus's solid globe. To explain the generation and maintenance mechanism of the Venusian atmospheric super-rotation, several mechanisms based on thermal tides, meridional circulation, or SS-AS circulation have been proposed (Schubert and Whitehead 1969;Thompson 1970;Fels and Lindzen 1974;Gierasch 1975;Matsuda 1980). Recent numerical studies suggest that both the thermal tide (e.g., Matsuda 2006, 2007) and meridional circulation mechanisms (e.g., Yamamoto and Takahashi 2004) might be viable in the Venusian atmosphere, in which the vertical and horizontal transport of angular momentum due to thermal tides, eddies, and mean meridional circulation plays an important role. Other kinds of waves, such as the four-and fiveday waves found in UV images of the cloud top (e.g., Del Kouyama et al. 2013Kouyama et al. , 2015 and small-scale disturbances suspected to be related to gravity waves (e.g., Sagdeev et al. 1986;Hinson and Jenkins 1995;Peralta et al. 2008;Garcia et al. 2009) may also play important roles in the establishment and maintenance of the atmosphere's super-rotation. Schubert et al. (1980) suggested that the zonal-mean meridional circulation consists of several vertically stacked cells extending from the equator to poles, based on wind profiles obtained by entry probes. Significant poleward winds at the cloud-top level have been observed in previous observations (e.g., Machado et al. 2017). It is noted, however, that the number of observed regions and/or local times has been limited and that the horizontal winds observed at the cloud top may include those associated with the thermal tide and various waves. We have no observational results as yet for the zonally averaged meridional circulation in the Venusian atmosphere. To elucidate the dynamical effects of meridional circulation, waves, and eddies, their threedimensional structures on both the dayside and nightside must be understood (Peralta et al. 2008).
The clouds, which result in a ~ 0.76 integrated-bond albedo for the planet (Moroz et al. 1985), play a major role in controlling the energy balance of the atmosphere. Sulfuric acid cloud particles are produced through photochemistry near and above the cloud tops and through condensation in updrafts within the clouds (Krasnopolsky and Pollack 1994;Imamura and Hashimoto 2001); however, specific dynamical processes contributing to the transport of key species such as SO 2 and cloud droplets are unclear. Moreover, there is an unidentified ultraviolet absorber in the clouds, which causes dark features in the UV, with related localized decrease in the albedo. This shows spatial and temporal variations, but its cause is unknown (Lee et al. 2015). To better understand the climate system of the cloudy planet, simultaneous global monitoring of the clouds' optical properties and the three-dimensional wind fields in various spatial scales is essential.
Atmospheric motion can be derived from cloud motions using cloud-tracking techniques (e.g., Rossow et al. 1990). It is difficult, however, to obtain the cloudmotion vectors (CMVs) accurately because the Venus disk on the detector array changes with time in the field of view (FOV). Uncertainties in the spacecraft attitude, the alignment of the instruments, and the like are also problematic. To derive the atmospheric motions accurately, it is crucially important to reduce attitudinal errors in the spacecraft and correct geometry coordinate information such as latitude and longitude. Since three-dimensional atmospheric motions can be complex, it is important to easily link and compare data at various altitudes obtained by different cameras and/or at different wavelengths. To realize this kind of data analysis, a data-archive design that enables users to treat data without any gaps between scientific instruments is required. It is also important to employ conventional analytical methods in terrestrial atmospheric science for direct and efficient application to the Akatsuki data. We have considered an appropriate format for the Akatsuki data, keeping in mind that the data would become a common heritage of humankind, just as previous planetary-exploration data have become.
In this paper, we present the data-processing schemes, content, and format of the Akatsuki public data in "Definition of data levels" section. We describe the method used to reduce the errors in spacecraft attitude and correct the geometry information, an improved limb-fitting technique originally described by Ogohara et al. (2012), in "Reliability of pointing correction by limb fitting" section. We also provide initial quality and accuracy estimations, which will be useful in the analysis of Akatsuki data in "Reliability of pointing correction by limb fitting" section. "Summary" section is a summary of this paper.

Definition of data levels
Akatsuki image data obtained by four imaging cameras (IR1, IR2, LIR, and UVI) are processed by a common data-production pipeline (Fig. 1). This pipeline is roughly divided into four levels (Levels 0-3), because returning to raw data is not necessary, for example, when calibration and geometry-correction techniques are improved. Table 1 summarizes the standard and ancillary data products of the four imaging cameras. Levels 1, 2, and 3 and the auxiliary data products (which will be defined in the following subsections) will be open to the public about 1.5 years after the arrival of the image data at the Akatsuki Science Data Archives (http://darts.isas.jaxa. jp/planet/project/akatsuki/). In addition to such image products, documents of the directory tree and the header list will be released to the public through the archive site (e.g., for L2b of UVI, http://darts.isas.jaxa.jp/pub/pds3/ vco-v-uvi-3-cdr-v1.0/vcouvi_1001/document/).

Level 0
Level-0 (L0) image data are the uncalibrated count values (i.e., not directly related to physical quantities) of each imaging camera, produced from raw telemetry. Ancillary data contain meta-data about the image (e.g., time of observation and exposure time duration) that are also produced from raw telemetry and ranging data.

Level 1
Level-1 data also consist of uncalibrated count values and are divided into two datasets. Level-1a (L1a) image data are produced in the FITS format by combining the Level-0 (L0) image data and the ancillary data, some of which are stored in the header of the HDU (Header Data Fig. 1 Schematic diagram of the science data-processing pipeline developed for Akatsuki. Science data telemetry packets are stored into the SIRIUS database in ISAS. Level 0 data is a binary image data that is constructed from the compressed and fragmented packet data. Level 1a data is produced in the FITS format by combining the Level 0 image data and the ancillary data, which are stored in the header of the HDU of the FITS file. Level 1b data is produced from Level 1a with geometry information created by processing navigation data using SPICE kernels, which are provided as two separated FITS files. The orientation of the image is also adjusted in this process. Level 2b data are produced from the Level 1b data, whose image count values are converted into physical quantities (radiance for UVI, IR1, IR2, and brightness temperature for LIR). Level 3 data are global-grid data in the NetCDF format on the regular longitude-latitude map produced from the Level 2b data by correcting the pointing errors of the camera using the limb-fitting technique. The Levels 1b, L2b, and L3 data with the SPICE kernels will be open to the public at the Akatsuki Science Data Archives (http://darts.isas.jaxa.jp/planet/project/akatsuki/) Unit) of the L1a FITS file, e.g., observation-program ID, exposure time, filter name, housekeeping data, and observational information such as time indicators (TIs).
Level-1b (L1b) image data are produced from L1a data as two geometry information files are created by processing navigational data in the SPICE-kernel format using the SPICE toolkit. The SPICE kernels related to Akatsuki are summarized in a meta-kernel file (e.g., vco_v01. tm) provided at the Akatsuki Science Data Archive. The ancillary data including observation date, instrument name, and filter name are stored in the header of the HDU of the L1b FITS file. A complete dictionary of the FITS header keywords stored in the L1b (and also L2b in the following subsection) FITS files is also present in the Akatsuki Science Data Archives. The geometry data are composed of multiple two-dimensional arrays, such as images observed, and include Venus longitude, latitude, local (solar) time, solar-incidence angle, emission angle, azimuthal angle (angle between the projection of the incidence and emission vectors onto the cloud-layer surface), and phase angle. One geometry file (geo file) includes the geometry data at the center of each pixel, and the other file (geo4 file) includes information at each pixel limit. The orientation of the image data is adjusted during this L1a-to-L1b processing by a procedure that combines 90° rotation and vertical flipping of the image. In the L1b-image data, the + Y-axis of the spacecraft (+ Y SC ), which is aligned with the solar-cell paddles, runs from the bottom of the image to the top (Ishii et al. 2004). This procedure differs among the cameras and filters because the alignment of each instrument against the spacecraft differs. The operated orientation adjustment is recorded as a value of the P_FLPROT keyword in the header of the HDU in the FITS file. To allow a quick look, a JPEG image is created for each FITS image file. See the catalog files provided in the public data for more details.

Level 2
Level-2 image data are produced from the L1b data (image data and geometry information), whose image count values are converted into physical quantitiesnamely, the radiation flux (for IR1, IR2, and UVI) and brightness temperature (for LIR), which have been calibrated by laboratory data and/or particular observations conducted for IR1 and UVI. For a quick look, JPEG images are also created for each FITS image file. Hereafter, we refer to Level 2 as Level 2b (L2b) to explicitly show the source data are L1b. It should be noted that Level 2a (L2a), which is L1a calibrated, is not produced because the calibrated L2b data are adequate for scientific research.
L2b image data share a FITS file format with L1b, which is fully compliant with the FITS Standard Version 4.0 [DRAFT] published online on July 22, 2016, by the IAUFWG. The image data are stored in the data part of the secondary HDU as an IMAGE extension. In all data stored in the HDU, + Y SC is upward and the − Z-axis of the spacecraft (− Z SC ) is rightward; when using the common coordinate system with the origin at the lower-left corner of the image, the horizontal axis (the first axis) is oriented from left to right, and the vertical axis (the second axis) is oriented from bottom to top. This coordinate system is commonly used in visualization software for FITS files, such as SAOImage DS9.
Each value in the L1b and L2b data arrays represents a value at the center of a pixel. The following three types of special pixels are defined. The first are missing pixels, which represent missing telemetry due to packet loss or override by other data in the data recorder (DR). The second are dead pixels, representing pixels that have been permanently damaged through such phenomena as X-ray hits to the corresponding detector element. The third are saturated pixels, representing data points that fall outside the nominal measurement range due to too many photons being received at the detector element. These three special pixels are defined and recorded in the header of the HDU in the FITS file using the following keywords P_MPIXV (missing pixel flag value), P_DPIXV (dead pixel flag value), P_SPIXV (saturated pixel flag value), and P_SPIXO (saturated pixel offset value). For completeness,

Table 1 Definitions of the Akatsuki data-processing levels
Definitions of the Akatsuki data-processing levels, which are common in UVI, IR1, IR2, and LIR instruments. The L1b, L2b, and L3 data products will be open to the public at the Akatsuki Science Data Archive (http://darts.isas.jaxa.jp/planet/project/akatsuki/)

L0
Level 0 Raw data

L1a
Level 1a Count values in FITS format

L1b
Level 1b Flipped and rotated version of L1a with summary of geometry information in FITS format, and detailed geometry information in FITS format

L2b
Level 2b Calibrated physical quantity with geometry information in FITS format

L3
Level 3 Calibrated physical quantity on longitude-latitude grids in NetCDF format which is preprocessed using the limb-fitting technique for pointing correction of the camera, including corrected geometry information the following keywords are also defined in the header of HDU: P_MPIXN (the number of missing pixels), P_ DPIXN (the number of dead pixels), and P_SPIXN (the number of saturated pixels). For more details, see the catalog file (dataset.cat) provided in the Akatsuki dataset. The geometry data included in the L1b-and L2b-image data are calculated using the VCO SPICE-kernel dataset, which is available from DARTS for ISAS/JAXA and from the PDS SPICE Archives at JPL/NASA. It should be emphasized that the geometry calculations are performed by assuming specific altitudes: 70 km for the UVI 283-and 365-nm channels, as well as the IR2 2.02-µm channel, 60 km for the IR1 0.90-(dayside), 0.90-(nightside), 0.97-, and 1.01-µm channels, 65 km for the LIR pic channel, 1 and 50 km for the IR2 1.735-, 2.26-, and 2.32µm channels. The assumed altitude is recorded in the FITS file header under the keyword S_CLDALT. Some of the geometric quantities appearing in the labels and headers of the FITS files are in J2000 coordinates. For the Venusian geometry, IAU_VENUS, which is defined as the planetocentric coordinate system for Venus, is used.

Level-3 data
One of the remarkable advantages of Akatsuki is that it allows continuous global monitoring from an equatorial orbit. Using sequential global images, we can visualize the propagation of atmospheric waves and estimate atmospheric motion by measuring the displacement of cloud features on images over a time interval Δt. This traditional method is applicable to all cameras and filters except for LAC. However, it is difficult to track a cloud feature and to measure its displacement in meters because the size and position of the Venus disk temporally vary in the FOV. Data analyses, such as an estimation of the phase speed of atmospheric waves, also suffer from the same difficulty. To promote a variety of scientific studies, we provide Level-3 (L3) data, in which observed variables are mapped onto a regular (equi-spaced) longitude-latitude grid. L3 data are produced from the L2b image data and are distributed in the Network Common Data Form (NetCDF) file format, which is widely used in the atmospheric science community and simplifies the application of conventional analytical methods to the Akatsuki image data. The resolution of the L3 longitudelatitude map is fixed at 0.125° × 0.125° (2880 × 1440 grids for 360° longitude and 180° latitude) for all image data obtained from our four image cameras. Figure 2 shows a sample of longitude-latitude images generated using the SPICE toolkit and kernels (https:// naif.jpl.nasa.gov/naif/toolkit.html). The dark area around 1 A 60-km altitude is assumed for the LIR pic channel in the dataset used in this study (internal release version v20170601).
the Venus limb is an area outside of the disk of Venus (i.e., in space) on the L2 image but is recognized as being inside the Venus disk based on the SPICE kernels. This type of mismapping results from errors in instrumental pointing associated with the attitude of the spacecraft and the alignment of the instruments. Although errors in the alignment of instruments can be reduced by stellar observations, errors in the spacecraft attitude still remain. Since the attitude determination accuracy is around 0.01° (1.75 × 10 −4 rad) or worse (Nakamura et al. 2011), the positions of cloud patterns in L2b images fluctuate by more than a pixel due to errors in the attitude of the spacecraft (Ogohara et al. 2012). Assuming Δt = 2 h, a positional error of 1 pixel leads to errors of 3-10 m s −1 in cloud speed, depending on the distance of the spacecraft from the center of Venus and the instantaneous field of view (IFOV) of the instruments. An error of a few meters per second in speed may be negligible compared to super-rotation; however, considering that 1 pixel is the ideal cloud-position error in an image with high accuracy in attitude determination, it cannot be ignored in comparison with the meridional component and the disturbances of the CMVs.
Such pointing errors in cameras can be corrected by calculating the position vector of Venus's center, − → L V , in instrumental coordinates (Fig. 3). The Venus disk is imaged as an ellipse. Because the spacecraft-Venuscenter distance, the wavelength-dependent cloud altitude, and the IFOV can be assumed to be known, − → L V can be uniquely determined if the positions of D and the semimajor axis of the elliptical Venus disk in Fig. 3 are determined by limb fitting (see Ogohara et al. 2012 for details). Therefore, an accurate determination of the Venusian limb on an image is required to accurately Fig. 2 Sample of a longitude-latitude map without the pointing correction. Sample of a longitude-latitude map derived from a UVI 283-nm image without the pointing-correction algorithm developed in this study (2015-12-07 05:19:58). The dark region is space, which should not be included in the longitude-latitude map convert an L2 image into a longitude-latitude map. Ogohara et al. (2012) showed that such a method of pointing correction offered a sufficiently high accuracy for Venusian atmospheric science using UV images taken by the Venus Monitoring Camera onboard Venus Express (VMC/Venus Express). The pointing-correction accuracy of multi-wavelength Venusian images taken by Akatsuki is estimated in "Reliability of pointing correction by limb fitting" section.
L3 data include observed variables (i.e., radiance or brightness temperature) and geometry data (i.e., incidence, emission, phase, and azimuthal angles) calculated based on corrected pointing information determined by the ellipse-fitting technique (Fig. 4). The contents of the L3 data are summarized in Table 2. The assumed cloud altitudes are identical to those described in the L2 data, which depend upon wavelength. The horizontal resolution is 0.125° × 0.125° (2880 × 1440 grids), as described above. Each grid interval on the equator is approximately 13 km, which is comparable to the spatial resolution around the sub-spacecraft point in L2 images other than LIR, when the distance of the spacecraft from Venus's center is roughly 7 × 10 4 km. To verify the results of limb detection and limb fitting, Level 3x (L3x) NetCDF and FITS files are also created. The L3x FITS file has a similar structure to the geometry information (geo) file generated at the same time as L1b, so the L3x FITS file can be used as a pointing-corrected geometry information file instead of the normal geometry information file.

Reliability of pointing correction by limb fitting
In this section, the accuracies of the limb fitting and pointing correction are evaluated by simulating the Venusian limb observed by the cameras (Ogohara et al. 2012). A few updates on the limb-fitting algorithm from Ogohara et al. (2012) are described in the following subsection. To simulate the Venusian limb, the shape of the Venusian limb observed by cameras was investigated, revealing that the limb could not be approximated by an ellipse plus random noise, as was done by Ogohara et al. (2012) for the VMC-UV images of Venus. Thus, the method used for simulating the Venusian limb observed by Akatsuki, rather than VMC, is reported. The artificial elliptical limbs are then randomly generated, using our method, and the accuracy in the pointing correction is evaluated by fitting ellipses to the artificial limbs.

Method of limb fitting
The ellipse-fitting technique used in this study is an improved version of the method that Ogohara et al. (2012) developed for UV images taken by VMC. There are two points of difference. First, limb points used for ellipse fitting are detected by fitting the hyperbolic tangent (tanh) function to one-dimensional radiance distributions in the X-or Y-directions, whereas the limb points were detected by Ogohara et al. (2012) using a combination of linear and tanh functions (referred to as the xtanh function, hereafter). This change was made because we empirically confirmed that the misfitting rate is lower and the variance of detected limb points is smaller in the tanh cases than in the xtanh cases for the Akatsuki data. The use of tanh fitting for limb-point detection contributes to the increase in limb points, especially in nightside images with low contrast between space and the Venusian disk. Moreover, we also changed the optimization algorithm. Ogohara et al. (2012) adopted a method in which an ellipse was updated step by step using multiple optimization algorithms, starting with the optimization algorithm presented by Taubin (1991) in the first step. However, some VMC-UV images clearly show incorrect elliptical results. In such cases, we adopt the HLS method (Iwamoto et al. 2009), instead of the Taubin method, to improve the ellipse fitting, since we found that the HLS method is more robust than the Taubin method when applied to Akatsuki images. The above two modifications did not decrease the accuracy of the elliptical fitting or pointing correction as long as VMC-UV images were used.

Simulating limb-point distributions for accuracy estimation
Although Ogohara et al. (2012) presented the accuracy of the ellipse fitting and pointing correction for VMC-UV images, the correction accuracy may depend upon cameras and wavelengths. In particular, the accuracy for an image of Venus's nightside taken by Akatsuki's IR2 camera may not be as high as that for VMC-UV images. Hence, we simulate the typical Venusian limb observed by the four cameras onboard Akatsuki, as Ogohara et al. (2012) did, and estimate the accuracy of the pointing correction.
Letting (p xi , p yi ) be the position of the ith limb point (1 ≤ i ≤ M), we define the argument of this point, θ i , as follows: Here, x c , y c is the position of the ellipse's center determined by fitting an ellipse to M limb points. The defined range of θ i is (− 90°, 90°) for p xi > x c and (90°, 270°) for  , which include longitude-latitude maps of radiance, emission, incidence, phase, and azimuthal angles. In the case of LIR, brightness temperature instead of radiance is included in Level 3 data p xi < x c . Furthermore, the radial residual between the ith limb point and the ellipse, r(θ i ), is defined as where (x, y) is a point on the ellipse. Figure 5 shows images observed by the four cameras and the radial residuals, r, obtained by applying ellipse fitting to them. The radial-residual distribution in the 283-nm images (Fig. 5h) appears to have a "3" shape, and in contrast, that in the 2.02-μm images (Fig. 5d) seems to have an "ε" shape with a phase opposite to the 3-shape. The limb points in the LIR image ( Fig. 5f ) are successfully detected at most arguments, and the distribution of their radial residuals from the ellipse is smooth. Although the 3-shape is also seen in the nightside image (1.01 μm, Fig. 5b), the limb points are distributed widely around the 3-shape due to the low contrast between the Venusian disk and space. Ogohara et al. (2012) investigated the accuracy of the ellipse fitting, assuming that the positions of the limb points detected in Venusian images could be simulated by adding Gaussian noise to a true ellipse.
Their assumption was realistic as far as they used VMC-UV images; however, it is obviously unrealistic in the case of images taken by Akatsuki because they did not consider the 3-shape or the ε-shape of the radial residuals shown in Fig. 5. We generate limb-point distributions randomly by adding Gaussian noise to a distribution of r that can be approximated by polynomials such as the 3-shape and ε-shape, fit an ellipse to the limb points, and investigate variations in the sub-spacecraft points in images determined by ellipse fitting. Based on the 1.01-μm, 2.02-μm, LIR, and 283-nm images shown in Fig. 5, r can be expressed as the summation of the noise and a sixth-order polynomial function of argument θ as follows: Here, θ is defined in the range of [− 90, 270), because the limb located on the right-hand side of an image becomes discontinuous. Solid lines shown in Fig. 5 are r(θ ) = aθ 6 + bθ 5 + cθ 4 + dθ 3 + eθ 2 + f θ + g.
functions obtained by fitting Eq. (4) to the radial residuals of limb points from the ellipse and successfully approximating them. Figure 6a shows the functions derived by fitting Eq. (4) to the radial residuals of limb points detected from 32 images (283-nm) for which the spacecraft-Venus-center distance is less than 1.6 × 10 4 km. Most of these functions seem to be 3-shape-symmetric about θ = 180 • , although they differ in amplitude. Figure 6b-g shows correlations between a and the other coefficients of Eq. (4). It is obvious that coefficients b to g can be approximated by a linear function of a. We obtained similar results for all other wavelengths, though correlations between the coefficients became low, especially for nightside images. Our results imply that the six coefficients (b-g) of Eq. (4) can be roughly determined by the other coefficient (a) for all wavelengths of Akatsuki's camera. Therefore, if we obtain the average ā and the standard deviation â of the coefficient a, and the standard deviation σ of the residuals from Eq. (4), which cannot be explained by Eq. (4), from several images, we can statistically simulate radial-residual distributions of detected limb points consisting of a sixth-order function and a noise component, as shown in Fig. 5. Table 3 lists ā , â, and σ for each wavelength obtained by fitting Eq. (4) to the radial residuals of the detected limb points from the ellipse determined by fitting (see Additional file 1 for the source images). L and S indicate the cases where the spacecraft is close to and far from the center of Venus, respectively. The average noise component that remains after fitting Eq. (4) to the radial residuals of the detected limb points is small enough to be negligible relative to σ. Although the shape of the sixth-order function is not easily presumed from the values of ā and â listed in Table 3 alone, the large σ of the nightside wavelengths means that limb points tend to disperse in the radial direction in nightside images.

Accuracy of pointing correction
To estimate the accuracy of the pointing correction, we introduce a virtual Venusian disk expressed by an ellipse E, whose center is located on C x ct , y ct , and an elliptical arc E′ between the angles θ 1 and θ 2 . The ellipticity and inclination of the semimajor axis of E depend upon C x ct , y ct and are determined based on the geometry described by Ogohara et al. (2012). Figure 7 shows samples of the true elliptical arcs E′ used in the next (See figure on next page.) Fig. 5 Observed images, the 3-shape, and the ε-shape. a, c, e, and g Sample images in the cases where the Venus disk is large (1.01 μm, 2.02 μm, LIR, 283 nm) and b, d, f, and h the distributions of radial residuals of limb points from the ellipses. White dots in (a, c, e, and g) indicate detected limb points used for ellipse fitting. The argument of a vector from the ellipse center toward the positive direction of X is 0° and that of each limb point is measured anticlockwise in (b, d, f, and h). Solid lines are sixth-order polynomial functions obtained by fitting Eq. (4) to radial residuals of the limb points from the ellipse subsection to estimate the accuracy of the pointing correction in the Cases L and S, as listed in Table 3. The number of limb points (M) and the range of the elliptical arc [θ 1 , θ 2 ] are set based on observations and are listed in Table 4. The elliptical arc E′ and the value M for the nightside wavelengths are short and small, respectively, because limb points are not detected at high latitudes, which are much darker than the mid-latitudes in these bands. The size of ellipse E is common for all cameras except for LIR, where it is smaller. By contrast, M for LIR is not much smaller than that for the other cameras because LIR captures the whole Venusian disk and E′ is identical to E. Table 4 also shows the true sub-spacecraft point (x SSCt , y SSCt ) assumed for the test and the angle, θ Vt , between the boresight vector and the true − → L V , as calculated from true ellipse E, based on Ogohara et al. (2012).
We add the radial residuals, r, and Gaussian noise, ε, generated from the parameters listed in Table 3 to the true limb points located on the true elliptical arc, E′. We use d 1 , d 2 , . . . , d M to denote the position vectors of M points that divide E′ equally into M − 1 parts and randomly generate a position vector, q j , of the jth pseudolimb point as follows: Here, c is the position vector of C; r θ j is calculated from Eq. (4) if coefficient a is given randomly based on the average and standard deviation listed in Table 3; and ε is also given based on Table 3.
Next, we estimate the sub-spacecraft point (x SSCe , y SSCe ) and the angle, θ Ve , between the boresight vector and − → L V , by fitting an ellipse to M pseudo-limb points expressed by Eqs. (5) and (6) and calculate the differences from their true values, x SSC = x SSCe − x SSCt , y SSC = y SSCe − y SSCt , and �θ V = θ Ve − θ Vt . M pseudolimb points are generated 1000 times independently using the parameters listed in Table 3. The ensemble average (denoted by an overline) and the ensemble standard deviation (denoted by a hat) of �x SSC , �y SSC , and�θ V are listed in Table 5.
The absolute values of �x SSC , �y SSC , and �θ V in Cases L and S are comparable to �x SSC , �y SSC , and �θ V , respectively, for dayside wavelengths. Such biases in the subspacecraft point associated with ellipse fitting were not recognized in results by Ogohara et al. (2012). The same trials as Ogohara et al. (2012) with r(θ ) = 0, in which just the noise component is given randomly, showed that �x SSC , �y SSC , and �θ V were very close to zero. Therefore, these biases result from the 3-or ε-shape of the limb-point distributions and the small standard deviation in the noise component. What may be a problem for practical use are variances in �x SSC , �y SSC , and �θ V , because we can rotate − → L V in addition when calculating the radiance on each longitude-latitude grid, such that these biases vanish. Differences x SSC and y SSC are smaller than 0.1 and 0.01, respectively, in the dayside Cases L and S. Difference x SSC is larger than y SSC , since the limb is assumed to be located on the left-hand side of the images (θ 1 = 115 • , θ 2 = 245 • ). There is not a large difference in �θ V between Cases L and S. However, the impact of the pointing accuracy upon cloud-feature positions in longitude-latitude coordinates is larger in Case S than in Case L, since the distance between the spacecraft and Venus's center in Case S is longer than that in Case L. Variations in cloud-feature positions (1σ) are, for example, 5.8 × 10 −2 km (UVI 365 nm) and 8.6 km (IR1 0.97 μm) when the spacecraft-cloud-layer distance is 10 4 km. When this distance is 10 5 km, the 1σ values are 3.6 × 10 −1 km (UVI 283 nm) and 7.0 km (IR2 2.26 μm). However, it should be noted that the above standard deviations in pointing are estimated under ideal conditions. Users should care about the effects of satellite thermal conditions and uncertainties in the limb altitudes upon the positions of cloud features if, for example, they are attempting to perform cloud tracking using L3 products ("Interpretation of the estimated accuracy" section).
The accuracies of sub-spacecraft-point determination at the nightside wavelengths are markedly worse than those at dayside wavelengths. This is because the length of the limb tends to be short and the detected limb points tend to be dispersed, due to low contrast between space and Venus (especially at high latitudes). There is no IR1 nightside image corresponding to Case S in which the spacecraft is far from Venus. In most nightside images, x SSC , y SSC , x SSC , and y SSC range from O(10 −3 ) to O(10 −1 ), but they are O(10 0 ) at some wavelengths. The ellipse fitting of IR1 0.97-μm images in Case L resulted in completely incorrect ellipses in some of the 1000 time trials conducted. It is likely that x SSC and y SSC seem to be large due to the small number of outliers. Note that the values for IR1 0.9-μm (nightside) in Case L listed in Table 5 are unreliable because only four observed images are available for the derivation of ā,â, and σ.
The accuracy of sub-spacecraft-point determination using images taken by LIR with an arc angle of 360° is much better than those obtained by the other instruments. However, it is not always better than that of the other cameras if it is measured in �θ V , �θ V because of the low spatial resolution of LIR. The difference in LIR

Table 3 Parameters calculated from detected limb points in each wavelength
The average ā and standard deviation â in coefficient a, and the standard deviation σ in residuals from Eq. (4). Coefficient a and residuals from Eq. (4) are derived by fitting Eq. (4) to radial residuals of limb points detected in images listed in Additional file 1 from the ellipses. Cases L and S mean the cases where the distance between the spacecraft and the Venus center is shorter than 1.6 × 10 4 km and longer than 3.5 × 10 4 km, respectively from the other cameras listed in Table 5 shows that | y SSC | is larger than | x SSC |. This is probably because the distributions of detected limb points tend to be longer in the Y-direction than the true ellipse.

Interpretation of the estimated accuracy
It is not clear what gives rise to the 3-shape and the ε-shape, and their reasons are beyond the scope of this study. However, we can determine the sub-spacecraft point and − → L V at the accuracies listed in Table 5, if the limb points follow such radial-residual distributions for any reason. Accuracies x SSC , y SSC and x SSC , y SSC listed in Table 5 range from O(10 −3 ) to O(10 −1 ) at most wavelengths. They are the accuracies of an algorithm that derives an optimal ellipse approximating the distribution of given points. In this study, we detect limb points by fitting a tanh function to the radiance distributions in the X-or Y-directions and assume that these limb points are located at a cloud altitude that depends upon wavelength. However, it is not clear whether one-dimensional distributions of radiance can be expressed by a tanh function, or whether it would be reasonable (or acceptable) to assume constant cloud-top altitudes along local times and latitudes. It remains unclear how the cloud altitudes depend upon latitude and local time at each wavelength. At present, we do not know the typical profile of radiance around the limb or where the profile of the limb should be defined. The positional accuracy of the limb points is reasonably lower than that listed in Table 5. Therefore, the pointing-correction accuracies listed in Table 5 should be understood to be the most optimistic values. Figure 8 shows sequential longitude-latitude maps (nightside) taken at a 2.26-μm wavelength by the IR2 camera. A feature (Feature A) seen around 115°W, 10°S in Fig. 8a is not shifted in Fig. 8b-d (meaning that Feature A is moving westward at 80 m s −1 ). However, Feature A suddenly jumps westward by about 2° in Fig. 8e. Because the time interval is 2 h, this jump corresponds to about Fig. 7 Schematic view of elliptical arcs used for accuracy estimation. Schematic view of elliptical arcs used for accuracy estimation of pointing correction, to which bias and noise calculated based on Eq. (4) and Table 3 are added, in the cases of (a, b) dayside images and (c, d) LIR images. Dotted lines are true ellipses, and solid lines are elliptical arcs simulating the Venus limb. A cross in each panel is the center of the true ellipse  a 30-m s −1 increase in the easterly wind speed. Furthermore, the westward jump is not local but is seen over the whole domain, as shown by the white arrows in Fig. 8. The radius of Venus in the original IR2 images around the time when the jump occurred was approximately 120 pixels, so the jump of 2°, which was observed near the subspacecraft point, corresponds to an erroneous shift in the pointing correction by ~ 4 pixels. This was the worst case, but jumps larger than ~ 20 m s −1 were found in about one-third of all IR2 2.26 μm cases where cloud tracking was successfully performed during the period from 2016-07-11 to 2016-09-06 (see Ikegawa and Horinouchi 2016;Horinouchi et al. 2017a, b for our cloud-tracking method). This indicates that, in the nightside cases, the actual error of the pointing correction can become O(1) pixels, which is much greater than the O(10 −1 )-pixel error indicated in Table 5. We are still investigating the reason for this discrepancy, but we suspect that it is related to the nature of the nightside radiance at ~ 2 μm, which is often very dark and obscure at the limb. Overall, we strongly recommend that L3-product users estimate the mapping accuracy for their particular topics (e.g., cloud tracking) by themselves. One way to do so is to compare cloud-tracking results from multiple combinations of successively obtained images. See the online supplement ("Methods" section) of Horinouchi et al. (2017b) for more details. Note that those authors used a different algorithm to correct pointing and succeeded in reducing the jumps. However, the performance of their algorithm for dayside and IR1-nightside images is unknown, so we did not use it to produce the Level-3 data. We plan to improve the fitting for nightside images in the Level-3 pipeline in the future. Unlike nightside images, discontinuities in cloud-feature positions rarely occur in dayside sequences. Comparing Figs. 2 and 4, the accuracy of longitude-latitude mapping (i.e., the pointing accuracy) has clearly been improved by ellipse fitting. The algorithm whose performance has been investigated by this study enables a variety of studies on the Venusian atmosphere.

Summary
We developed a data-processing pipeline to generate Akatsuki imaging data from the UVI, IR1, IR2, and LIR cameras. This pipeline has multiple levels, which we call Levels 1a, 1b, 2b, and 3. Levels 1 (1a and 1b) and 2 (2b) consist of uncalibrated and calibrated data (radiances for UVI, IR1, and IR2, and brightness temperature for LIR), respectively. Level-1b data are based on the data from Level 1a by adding geometry information such as the phase, incidence, and emission angles. The data for Levels 1a, 1b, and 2b are stored in the FITS format, commonly used for planetary science. Level-3 data are global maps on the regular longitude-latitude grids, which are produced based on Level-2b data. The errors in the navigation data and the instrument alignment, which cannot be neglected when producing Level-3 data, are corrected by an improved limb-fitting technique in which the pointing of the boresight vector is corrected by fitting an ellipse to the Venusian limb. Level-3 data are stored in the NetCDF format commonly used for atmospheric science. The accuracy of pointing correction has been statistically estimated. The result shows that the algorithm we developed and used for Level-3 data works well; the uncertainties in the sub-spacecraft point are about O(10 −2 -10 −1 ) pixels. It is expected that the Level-3 data will be useful for deriving the CMVs accurately, except for some of the nightside images. The pointing-correction Fig. 8 Sequential longitude-latitude maps of radiance (IR2 2.26 μm). Sequential longitude-latitude maps of radiance taken in 2.26 μm by IR2 camera on August 26, 2016 (08:03, 10:03, 12:03, 14:03, 16:03, 18:03, 20:03, and 22:03), which are subtracted images between the original radiance images and the images smoothed by a Gaussian filter with half-width of 2° in both directions. For visibility, the cloud feature shown in (a) is shown being shifted westward with 80 m s −1 in the following panels. Cloud patterns over the whole domain shown in the panels shift westward discontinuously in (e) probably due to mismapping, as indicated by the white arrows. Feature A mentioned in the text is denoted by the black arrows method developed for the Akatsuki data pipeline can be applied to the Venusian images obtained from previous missions. It is noted, however, that the present error estimation of the pointing correction is based on assumed ideal conditions. Therefore, we must pay attention to various factors such as the thermal conditions of the satellite and/or the uncertainty in the assumed cloud altitudes, which could affect the pointing correction. We hope that the Akatsuki public data will be widely used for Venusian atmospheric science.