Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: object registration is an important task in computer vision that determines the position and the orientation of an object in camera - centered coordinates @xcite . an object of interest that was detected beforehand in a coarse 2d bounding box is fed into a registration system that can superimpose the desired translation and rotation of the object onto the raw camera image . by utilizing such a system , one can propose promising solutions for various problems related to scene understanding , augmented reality , control and navigation of robotics , _ etc_. recent developments on visual depth sensors and their increasing ubiquity have allowed researchers to make use of the information acquired from these devices to facilitate challenging registration scenarios . + iterative closest point ( icp ) algorithm @xcite , point - to - model based methods @xcite and point - to - point techniques @xcite demonstrate good registration results . however , the performance of these approaches is severely degraded in cases of heavy occlusion and clutter , and similar looking distractors . in order to address these challenges , several learning based methods formulate occlusion aware features @xcite , derive patch - based ( local ) descriptors @xcite or encode the contextual information of the objects with simple depth pixels @xcite and integrate them into random forests . most particularly , iterative random forest algorithms such as latent - class hough forest ( lchf ) @xcite and iterative multi - output random forest ( imorf ) @xcite obtain state - of - the - art accuracy on pose estimation . on the other hand , these methods rely on scale - invariant features , while the exploitation of rich discriminative information inherently embedded into the scale - variability is one important point been overlooked . +        [ fig1 ]    unlike the aforementioned learning - based methods , the ones presented by novatnack _ @xcite utilize the detailed information of the scale variation in order to register range images in a coarse - to - fine fashion . although promising , they extract and match conventional salient 3d key points . however , real depth sensors have several imperfections such as missing depth values , noisy measurements , foreground fattening , _ etc_. as a result , salient feature points used in @xcite tend to be located on these deficient parts of the depth images , and hence , they are rather unstable @xcite . in such a scenario , 3d reconstruction methods that provide more reliable shape information can be utilized @xcite . implicit b - splines ( ibs ) @xcite are techniques that can provide shape descriptors through their zero - sets and reconstruct surfaces . they are based on locally controlled functions that when combined with their control points produce a very rich part - based object representation . + our architecture is originated from these observations . we integrate the coarse - to - fine registration approach presented in @xcite into the random forest framework @xcite using histogram of control points ( hocp ) features that we adapt from recently introduced ibss @xcite . we train our forest only from positive samples and learn the detailed information of the scale - variability during training . we normalize every training point cloud into a unit cube and then generate a set of scale - space images , each of which is separated by a constant factor . the parts extracted from the images in this set are represented with the scale - variant hocp features . during inference , the parts centered on the pixels that belong to the background and foreground clutters are removed iteratively using the most confident hypotheses and the test image is updated . since this removal process decreases the standard deviation of the test point cloud , subsequent normalization applied to the updated test image increases the relative scale of the object ( foreground pixels ) in the unit cube . more discriminative control point descriptors are computed at higher scales and this ensures the refinement of the object pose . in our prior work @xcite , we have evaluated the registration performance of the proposed architecture by only using fixed - size parts . we extend the work engineering an automatic variable size part extraction framework in such a way that we can further exploit the discriminative information provided by the hocp features . this framework first roughly aligns the object of interest by extracting coarsest parts , the ones occupying the largest area in image pixels , and then iteratively refines its alignment based on finer ( smaller ) parts that are represented with more discriminative control point descriptors . note that we employ a compositional approach , _ i.e. _ , we concurrently detect the object in the target region and estimate its pose by aligning the parts in order to increase robustness across clutter . figure [ fig1 ] depicts a sample result of our architecture . to summarize , our main contributions are as follows :    * to the best of our knowledge this is the first time an implicit object representation , implicit b - spline , is adapted into a `` scale - variant '' part descriptor and is associated with the random forests . * we introduce a novel iterative algorithm for the hough forests : it finds out an initial hypothesis and improves its confidence iteratively by extracting more discriminative `` scale - variant '' descriptors due to the elimination of the background / foreground clutter . * we engineer an automatic variable size part extraction framework for the random forests : it first roughly aligns the object of interest by extracting coarsest parts and iteratively improves its confidence based on finer ( smaller ) ones . the rest of the paper is organized as follows : in sect . [ related work ] , we present a review on the object registration . section [ approach ] demonstrates the computation procedure of the hocp features as a scale - variant part representation , their combination with the iterative hough forest ( ihf ) , and the registration process . experimental results are provided in sect . [ experiments ] and finally , the paper is concluded in sect . [ conclusion ] with several remarks , and discussions . a large number of studies have been proposed for the object registration , ranging from the point - wise correspondence based methods to the learning based approaches . iterative closest point ( icp ) algorithm , originally presented in @xcite , requires a good initialization in order not to be trapped in a local minimum during fine tuning . this issue is addressed in @xcite providing globally optimal registration by the integration of a global branch - and - bound ( bnb ) optimization into the local icp . the point - wise correspondence problem is converted into a point - to - model registration in @xcite . the object model is represented with implicit polynomials ( ip ) and the distance between the test point set and the object model is minimized via the levenberg - marquardt algorithm ( lma ) . @xcite propose a 6 dof pose estimation technique utilizing 3d ips on ultrasound images . in the off - line phase , object model is represented with 3d ips and by utilizing its gradient flow , 2d ultrasound image is registered in the on - line process . in @xcite , a coarse - to - fine fast ip - driven registration method is presented . a rough pose estimation is quickly acquired with a coarse ip model ( low degree curve fitting ) and finer models refine the parameters of this rough estimation ( high degree curve fitting ) . hinterstoisser et al . @xcite extract holistic templates from 3d models of the objects and match to the scene at test time . these studies have demonstrated good registration results on the target point sets that are occlusion - free and/or are subjected to the artificial gaussian noises and outliers . + unlike the abovementioned methods , more realistic registration scenarios have been addressed by the point - to - point techniques that build point - pair features for sparse representations of the test and the model point sets @xcite . _ align two noisy point clouds of real scenes by finding correct point - to - point correspondences between the point feature histograms ( pfh ) and feed this alignment to an icp algorithm for fine tuning @xcite . the cluttered and partially occluded objects poses are hypothesized by accumulating the votes of the matching features in @xcite . _ @xcite propose point - pair features for both rgb and depth channels that are conducted in a voting scheme to hypothesize the rotation and translation parameters of the objects in the cluttered scenes . the features proposed in @xcite make use of the visibility context of the scene to tackle the registration . despite achieving good registration results , these techniques underperform when the scenes are under heavy occlusion and clutter , and the objects geometry are indistinguishable from background . + learning - based methods have good generalization across severe occlusion and clutter @xcite . the method presented in @xcite formulates the recognition problem globally and derives occlusion aware features . a set of principal curvature ratios are computed for all pixels in depth images to extract the edgelets . in @xcite , the contextual information of the objects is encoded with simple depth and rgb pixels . this technique improves the confidence of a pose hypothesis using a ransac - like algorithm . @xcite back project the parts inside the initially found coarse bounding box to the image and pass down the forest again . the parts with the lowest contributions are penalized in such a way that finer registration is produced . + the state - of - the - art accuracy on registration is acquired by the iterative random forest algorithms . the part - based strategy , latent - class hough forest @xcite , refines the initially hypothesized object pose by iteratively updating the object class distributions in the leaf nodes during testing . iterative multi output random forest ( imorf ) @xcite jointly predicts the head pose , the facial expression and the landmark positions . the relations between these tasks are modelled so that their performances are iteratively improved with the extraction of more informative features . the ideas , iterative pose refinement during testing and iterative extraction of more discriminative features , form a basis for our iterative hough forest ( ihf ) architecture : during training , we encode discriminative shape information of the hocp features into the forest . despite that the skeleton of our training procedure is similar to the methodology in @xcite , our forest learns the discriminative shape information that will be iteratively exploited at test time . in the course of inference , unlike @xcite , we update the test image itself and the hypotheses confidence by a noise removal process that allows us to extract more informative features from the test images . whilst these approaches @xcite rely on the scale - invariant features to improve the confidence of a pose hypothesis , novatnack _ et al . _ @xcite introduce a framework that registers the range images in a coarse - to - fine fashion by utilizing the detailed information provided by the scale variation . the shape descriptors with the coarsest scale are matched initially and a rough alignment is achieved since fewer features are extracted in coarser scales . the descriptor matching at higher scales produces improved predictions of the pose . inspired by @xcite , we design a `` scale - variant '' implicit volumetric part description , `` histogram of control points ( hocp ) '' , and associate it with the random forest framework . + selecting the part size is important since larger parts tend to match the disadvantages of a holistic template while smaller ones are prone to noise @xcite . in heavily occluded and cluttered scenes , relatively smaller parts perform well whilst the larger ones are more convenient in occlusion / clutter - free scenarios . discriminative information encoded into small sized parts might not be fully exploited by larger parts , most particularly when the object representation is scale - variant . hence , this application - specific / task - dependent part size selection degrades the generalization and it is one of the remaining challenges that should be addressed , apart from occlusion , clutter and/or similar looking distractors . beyond object pose estimation @xcite , there are several part - based solutions proposed for different tasks such as human pose recognition @xcite , 3d face analysis @xcite , or hand pose estimation @xcite to name a few . they experience different part sizes and select the one that performs best , however , none of these solutions investigate how extracting variable size parts can be utilized in a single framework . in this study , we investigate the effect of this _ size variation _ and show that the simultaneous utilization of the parts of varying size can improve 6d object pose estimation , especially in heavily occluded and cluttered depth maps , supplying a rich source of discriminative information . in this section we detail our registration approach by firstly describing the computation procedure of the hocp features as a scale - variant part representation . we then present how to encode the discriminative information of these scale - variant features into the forest . finally , we demonstrate how to exploit the learnt shape information in a coarse - to - fine fashion to refine the pose hypotheses .      given a positive depth image we initially normalize it into a unit cube and then new point clouds at different scales are sampled as follows : @xmath0 with @xmath1 where @xmath2 $ ] is the world coordinate vector of the original foreground point cloud , @xmath3 is the mean of @xmath4 , @xmath5 $ ] is the normalized foreground pixels , @xmath6 is the number of the scales , @xmath7 is the scale factor and @xmath8 is the scale . the constant @xmath9 takes real numbers to generate the point clouds at different scales , starting from @xmath10 that corresponds to the initial normalization . a training image and its samples at different scales are shown in fig . [ fig2 ] ( a ) . + once we generate a set of scale - space images , we represent these point clouds first globally with the control points of implicit b - splines ( ibs ) . ibs is defined through the combination of b - spline tensor products :        @xmath11 where @xmath12 are the coefficients defining a control lattice of size @xmath13 and @xmath14 , @xmath15 , @xmath16 are the spline basis functions . this definition can be reformulated as the following inner product : @xmath17 where the coefficient vector @xmath18 includes the control values @xmath12 and the basis vector @xmath19 depends on the given data point since it sorts the spline basis function products @xmath20 . the basis vectors in eq . [ eq4 ] are computed for the whole point cloud and the coefficient vector @xmath18 is calculated based on the 3l algorithm @xcite . rouhani et al . @xcite construct the spline basis functions @xmath14 , @xmath15 , @xmath16 through the following blending functions : @xmath21 and reformulate eq . [ eq3 ] in order to determine the control point vector @xmath18 of the point clouds that are normalized into the unit cube @xmath22 ^ 3 $ ] : @xmath23 where @xmath24 thus , the unit cube is split into an @xmath25 voxel grid where @xmath26 is the ibs resolution . each control point in @xmath18 is defined with an index - weight pair : the index number indicates the vertex of this grid at which the related control point is located . the weight informs the descriptor significance about the control of the geometry to be represented . the scale - space images in fig . [ fig2 ] ( a ) are globally represented in fig . [ fig2 ] ( b ) with the control point descriptors . we use all control points to represent the structures , but one can sort these descriptors based on their weights and utilize the ones higher than a threshold . in fig . [ fig2 ] ( b ) , the point clouds are lastly reconstructed by using the control points to show the increase in their discriminative power as the scale gets higher . note that ibs resolution @xmath26 determines the complexity of the representation ( level of detail ) in a unit cube , whilst the scale indicates the relative size of the object with respect to the unit cube dimensions . in our architecture , despite sampling point clouds at different scales @xmath27 , ( _ e.g. _ @xmath28 ) , we fix the complexity of the representation , ( _ e.g. _ @xmath29 ) . + ibs is the combination of the locally controlled functions and allows one to propose effective part - based solutions for object registration . we benefit from such a property and partition the globally represented scale - space depth maps into parts . we express the part size @xmath30 in image pixels which also depicts the ratio between the sizes of the extracted part and the bounding box of the global point cloud . in our prior work @xcite , we have extracted and represented the parts that have the same size , that is , the parts growing around every individual pixel at each scale occupy the same area in image pixels . we now extend the work extracting the parts that are different in size . a 3d bounding box defined in metric coordinates is traversed in the unit cube of each scale - space image and the parts are extracted around non - zero pixels . the total number of the data points in this 3d bounding box varies for the point clouds at different scales , and consequently , the size of the extracted parts differs . figure [ fig2 ] ( c ) shows an example of variable size part extraction process in which the red parts are grown around the same data point . note how the part size decreases when the scale of the normalized object point cloud gets higher , since less number of data points are deployed in the 3d bounding box . each part has its own implicit volumetric representation , formed by the closest control points to the part center , the ones lying inside the 3d bounding box along depth direction . such a part description characterizes the locality in a cascaded fashion , growing regions with different characteristics around a point . we encode this information into histograms in spherical coordinates . each of the part centers is coincided with the center of a sphere . the control points of the part are described by the log of the radius @xmath31 , the cosine of the inclination @xmath32 and the azimuth @xmath33 . then , the sphere is divided into the bins and the relation between the bin numbers @xmath34 and the histogram coordinates @xmath35 is given as follows @xcite : @xmath36 where @xmath37 and @xmath38 are the radii of the nested spheres with the minimum and the maximum volumes , @xmath39 are the cartesian coordinates of each descriptor with radius @xmath40 . @xmath38 equals to the distance between the patch center and the farthest descriptor of the related patch . the numbers of the control points in each bin are counted and stored in a @xmath41 dimensional feature vector @xmath42 . figure [ fig2 ] ( d ) illustrates the hocp representations of the parts extracted in fig . [ fig2 ] ( c ) . note that the control points computed at higher scales capture more detailed part geometry . the proposed ihf is the combination of randomized binary decision trees . it is trained only on foreground synthetically rendered depth images of the object of interest . we generate a set of scale - space images from each training point cloud and sample a set of parts @xmath43 as explained in subsection [ varpart ] and annotate those as follows : @xmath44 where @xmath45 is the part center in pixels , @xmath46 is the 3d offset between the centers of the part and the object , @xmath47 is the rotation parameters of the point cloud from which the part @xmath48 is extracted and @xmath49 is the depth map of the part . + each tree is constructed by using a subset @xmath50 of the annotated training parts @xmath51 . we randomly select a template patch @xmath52 from @xmath50 and assign it to the root node . we measure the similarity between @xmath52 and each patch @xmath53 in @xmath50 as follows :    * * depth check : * the depth values of the descriptors @xmath54 and @xmath55 that represent the parts @xmath53 and @xmath52 are checked , and the spatially inconsistent ones in @xmath54 are removed as in @xcite , generating @xmath56 that includes the spatially consistent descriptors of the patch @xmath53 . * * similarity measure : * using @xmath56 , the feature vector @xmath57 is generated and the @xmath58 norm between this vector and @xmath59 is measured : @xmath60 * * similarity score comparison : * each patch is passed either to the left or the right child nodes according to the split function that compares the score of the similarity measure @xmath61 and a randomly chosen threshold @xmath62 . the main reason why we apply a depth check to the patches is to remove the structural perturbations , due to occlusion and clutter @xcite . these perturbations most likely occur on patches extracted along depth discontinuities such as the contours of the objects . they force a test patch to diverge ( occluded / cluttered ) from its positive correspondence by changing its representation , @xmath38 of the sphere , and the histogram coordinates consequently . + a group of candidate split functions are produced at each node by using a set of randomly assigned patches @xmath63 and thresholds @xmath64 . the one that best optimize the offset and pose regression entropy @xcite is selected as the split function . each tree is grown by repeating this process recursively until the forest termination criteria are satisfied . when the termination conditions are met , the leaf nodes are formed and they store votes for both the object center @xmath65 and the object rotation @xmath66 . + depending on the part extraction approach , all parts in @xmath67 ( see eq . [ eq18 ] ) can either be of the same size or of the variable size . from now on , we will refer to the forests trained on variable size parts as the _ ihf - variable size _ , and to the ones learnt by using fixed size parts as the _ ihf - fixed size_.      once we encode the discriminative information of the scale - variant hocp features into the forest , we next demonstrate 6d pose estimation of objects considering that the learnt forest is ihf - variable size . + the proposed architecture registers objects in two steps : the _ initial registration _ and the _ iterative pose refinement_. the _ initial registration _ roughly aligns the test object and this alignment is further improved by the _ iterative pose refinement_. + consider an object that was detected by a coarse bounding box , @xmath68 , as shown in the leftmost image of fig . [ fig3 ] ( a ) . at an iteration instant @xmath69 , the following quantities are defined :    * @xmath70 : the history of the object position predictions . * @xmath71 : the history of the object rotation estimations . * @xmath72 : the history of the inputs ( noise removals ) applied to the test image . * @xmath73 : the history of the set of the feature vectors where @xmath74 . * @xmath75 : the object scale ( the scale of the foreground pixels ) in the unit cube at iteration @xmath69 ( see eq . [ eq2 ] ) . * @xmath76 : the size of the parts extracted at iteration @xmath69 . we formulate the _ initial registration _ as follows : @xmath77        [ fig3 ]    we find the best parameters that maximize the joint posterior density of the initial object position @xmath78 and the initial object rotation @xmath79 . the initial registration process is illustrated in fig . [ fig3 ] ( a ) . the test image is firstly normalized into a unit cube . unlike training , this is a `` single '' scale normalization that corresponds to @xmath10 ( see eq . [ eq1 ] ) . the patches extracted from the globally represented point cloud are described with the hocp features and are passed down all the trees . at this stage , we extract the coarsest patches from the test image , _ i.e. _ , the ones occupying the largest area in image pixels . we determine the effect that all patches have on the object pose by accumulating the votes stored in the leaf nodes as in @xcite and approximate the initial registration given in eq . once the initial hypothesis @xmath80 is obtained , the set of pixels that belong to the background / foreground clutter @xmath81 are removed from @xmath68 according to the following criterion : @xmath82 with @xmath83 where @xmath84 and @xmath85 are the depth maps of the hypothesis @xmath86 at iteration @xmath69 , and of the @xmath68 , @xmath87 and @xmath88 are the scaling coefficients . the efficacy of @xmath89 is illustrated in fig . [ fig3 ] . in the rightmost image of fig . [ fig3 ] ( a ) , the test image and the initial hypothesis are superimposed . this hypothesis is exploited and the test image is updated by @xmath90 as in eq . the updated image is shown in fig . [ fig3 ] ( b ) and assigned as input for the @xmath91 iteration . it is normalized and represented globally . the object `` scale '' ( @xmath92 ) in the unit cube is relatively increased ( compare with the initial registration ) and more discriminative control point descriptors @xmath18 are computed . this is mainly because of the fact that the standard deviation of the input image decreased since we removed foreground / background clutter . as a follow up step , we traverse the 3d bounding box in the unit cube during part extraction , while the increase in the normalized object scale gives rise to extract patches whose size are smaller ( finer ) than the ones extracted during the initial registration . the resulted hypothesis of the @xmath91 iteration is shown on the right . the extraction of finer parts represented with more discriminative control point descriptors along with the noise removal process result in more accurate and confident hypothesis . this pose refinement process is iteratively performed until the maximum iteration is reached ( see fig . [ fig3 ] ( c ) ) : @xmath93 we approximate the registration hypothesis at each iteration by using the stored information in the leaf nodes as we do in the initial registration . if we would demonstrate the 6d object pose estimation considering that the learnt forest is the ihf - fixed size , the only difference in the formulation would be the part extraction viewpoint . instead of traversing 3d bounding box in the unit cube , we would extract the parts with a predefined size in pixels , and at an iteration instant @xmath69 , @xmath76 would remain the same as @xmath94 . in the next section , we will compare the registration performances of the forests that are separately trained on fixed and variable size parts . there are several publicly available datasets @xcite , @xcite to test the performances of the object registration methods . for each object type in these datasets , a set of rgb - d test images are provided with ground truth object pose parameters . we have analysed these images and have found that `` coffee cup '' , `` camera '' and `` shampoo '' ( included in the icvl dataset @xcite ) are some of the best demonstrable objects to test and to compare our registration architecture with the state - of - the - art methods since they are located in highly occluded and cluttered scenes . we further process the test images of these objects to generate a new test dataset according to the following criteria :        * since the hocp features are scale - variant , the depth values of the training and the test images should be close to each other up to a certain degree . in this study , we train the forests at a single depth value , @xmath95 mm , and test with the images at the range of @xmath96 $ ] mm . * the test object instances located at the range of @xmath96 $ ] mm are assumed as detected by coarse bounding boxes ( see fig . [ fig4 ] ) . the image regions included in these bounding boxes are cropped and the new test dataset is generated . the generated dataset includes @xmath97 `` coffee cup '' , @xmath98 `` camera '' and @xmath99 `` shampoo '' images each of which is at the range of @xmath100 $ ] mm since we train the forests used in all experiments with the positive samples at @xmath101 mm depth . the maximum depth is @xmath102 and the number of the maximum samples at each leaf node is @xmath103 for each tree . every forest is the ensemble of @xmath104 trees with these termination criteria . +     and feature dimensions @xmath105 . ( c ) shows the effect of the iteration number . for the corresponding f1 scores see table [ table_params].,height=528 ]    our experiments are two folds : _ parameter optimization _ and _ comparative study_. the architecture parameters have an important impact upon the registration and include the size of the parts extracted during the initial registration @xmath94 , the ibs resolution @xmath26 , the hocp feature dimension @xmath105 ( the number of the bins or quantization parameter ) and the iteration number . once the best parameters are acquired , we compare the performance of our architecture with the state - of - the - art methods in the comparative study . + both experiments use the metric proposed in @xcite to determine whether a registration hypothesis is correct . this metric outputs a score @xmath106 that calculates the distance between the ground truth and estimated poses of the test object . the registration hypothesis that ensures the following inequality is considered as correct : @xmath107 where @xmath108 is the diameter of the 3d model of the test object and @xmath109 is a constant that determines the coarseness of an hypothesis that is assigned as correct . we set @xmath109 to @xmath110 when we refine the parameters of the proposed architecture and the effect of various @xmath109 values is separately examined in the comparative study . .f1 scores of the initial registrations determined for different part sizes ( g ) , ibs resolution ( n ) & feature dimension ( d ) and number of iteration [ cols="^,^,^,^ " , ]     [ table_interclass_aver ]    the pr curves of the coffee cup dataset are depicted in fig . [ fig7 ] for several @xmath111 values and their corresponding f1 scores are presented in table [ table_interclass_cup ] . a short analysis on the images of fig . [ fig7 ] reveals that the increase in @xmath111 value generates higher f1 scores for each approach . according to table [ table_interclass_cup ] , the lchf trained on the color gradient + surface normal features underperforms the lchf trained merely on the surface normals . the main reason of this underperformance is the distortion along the object borders arising from the occlusion and the clutter , that is , the distortion of the color gradient information in the test process . the initial registration of the ihf - fixed size performs better than any versions of the lchf . this shows that the hocp features represent the parts better than the surface normals and color gradients providing robustness across occlusion , clutter and missing depth pixels . when this initial registration ( see the @xmath112 column of table [ table_interclass_cup ] ) acquired from the ihf - fixed size is iteratively refined , more accurate registrations are resulted ( see the @xmath113 column of table [ table_interclass_cup ] ) . because , the _ iterative pose refinement _ module of the ihf - fixed size reduces the amount of the noise in the test depth maps removing foreground / background clutter . this removal process also enables ihf - fixed size to compute more discriminative control points for better shape representation . the ihf - variable size with initial registration + iterative pose refinement outperforms other approaches . the main reason of this high performance is the utilization of the parts that are different in size . the cascaded representation of the locality increases the robustness across clutter , occlusion , missing depth pixels and/or similar looking distractors . the object of interest is first roughly aligned by extracting the coarsest parts . it is highly possible that these initially extracted parts include the portions belonging to the background / foreground clutter since they are the coarsest and are close to a holistic template . despite the fact that we apply a depth check in order to get rid of those noise during testing , it is highly naive . on the other hand , the proposed framework can get rid of those structural perturbations by growing smaller regions as the iteration progresses . apart from that , the control point descriptors computed at later iterations allows the complete framework to represent the shapes in a more discriminative manner . + regarding the camera dataset , we show its pr curves in fig . [ fig8 ] for several @xmath111 values and the corresponding f1 scores in table [ table_interclass_camera ] . the approaches under comparison perform worse on this dataset with respect to the coffee cup . unlike the results obtained from the coffee cup dataset , we see the positive impact of the color gradients when they are utilized along with the surface normals at most of the @xmath111 values ( compare the @xmath91 and the @xmath114 columns of table [ table_interclass_camera ] ) . ihf - fixed size register objects more accurate than any versions of the lchf thanks to the utilization of the discriminative information embedded into the scale - variant hocp features and the iterative refinement of the test depth maps . ihf - variable size significantly outperforms other approaches demonstrating the importance of the simultaneous utilization of variable size parts . the hocp representations of the cascaded regions grown around the same data points allow the algorithm to be aware of occlusion , clutter and missing depth values . more confident registrations are hypothesized as the iteration progresses based on more discriminative representations of the smaller parts . the registration performance of the proposed architecture is shown in figure [ fig8 ] for the shampoo object and the corresponding f1 scores are demonstrated in table [ table_interclass_shampoo ] for varying @xmath111 values . in cases of registering at lower @xmath111 values , our approach shows better performance than the lchfs , however , when we accept relatively rough estimations as correct , _ i.e. _ , higher @xmath111s , our approach underperforms . + since we address the registration problem rather than individual object detection or pose estimation , we integrate the effect of the different error ratios into our comparisons . we average the f1 scores that are computed at each @xmath111 in the range of [ 0.05 - 0.15 ] and report in table [ table_interclass_aver ] . figure [ fig9 ] illustrates several accurate registrations hypothesized by the proposed framework on the camera and the coffee cup objects . we further evaluate the performance of the globally optimized icp algorithm proposed in @xcite on our test dataset . we use the software kindly provided by the authors and set the default parameters . while accurate registration results are obtained on the clean dataset , it diverges in the case of highly occluded and cluttered point clouds ( see fig . [ fig10 ] ) . in this study , we have proposed a novel architecture , _ iterative hough forest with histogram of control points _ , addressing 6d object registration rather than individually estimating either the object s location in a 2d/3d bounding box , or the object s orientation ( roll , pitch , yaw ) . any off - the - shelf detector can accurately provide a coarse 2d or 3d bounding box for the object of interest . various object orientation predictors are also available , however , they depend on clearly segmented target objects . our architecture fundamentally targets to eliminate the shortcomings of these individual detectors and orientation predictors estimating occluded and cluttered objects 6d pose given a candidate 2d bounding box . our ihf is learnt using parts extracted only from the positive samples . these parts are represented with scale - variant hocp features , which we derive from recently introduced implicit b - splines ( ibs ) . + at test time , we apply two different strategies regarding the parts used to train the forest : the first strategy we apply roughly aligns the object and iteratively refines this alignment based on more discriminative control point descriptors that are computed due to the elimination of background / foreground clutter . the part size is fixed and is empirically predefined . on the other hand , the predefined part size might not be generalizable enough across different objects , degrading the registration performance of the proposed study on one object while working well on another one . besides , discriminative information encoded into small sized parts might not be fully exploited by larger parts , most particularly when the object representation is scale - variant . inspiring by these observations , we use variable size parts in the second strategy . an automatic variable size part extraction framework iteratively refines the object s initial pose that is roughly aligned due to the extraction of coarsest parts , the ones occupying the largest area in image pixels . the iterative refinement is accomplished based on finer ( smaller ) parts that are represented with more discriminative control point descriptors by using our _ iterative hough forest_. the experimental results report that our approach show better registration performance than the state - of - the - art methods . caner sahin is funded by the turkish ministry of national education . b.  zheng , r.  ishikawa , j.  takamatsu , t.  oishi , and k.  ikeuchi , `` a coarse - to - fine ip - driven registration for pose estimation from single ultrasound image '' , _ computer vision and image understanding _ , 2013 , vol.117 , no.12 , ( pp . 1647 - 1658 ) .            c. r.  cabrera , r. l.  sastre , and t.  tuytelaars , `` all together now : simultaneous object detection and continuous pose estimation using a hough forest with probabilistic locally enhanced voting '' , _ proceedings bmvc _ , 2014 , ( pp . 1 - 12 ) . s.  hinterstoisser , v.  lepetit , s.  ilic , s.  holzer , g.  bradski , k.  konolige , and n.  navab , `` model based training , detection and pose estimation of texture - less 3d objects in heavily cluttered scenes '' , _ accv _ , 2012 .
Abstract:

Answer: state - of - the - art techniques for 6d object pose recovery depend on occlusion - free point clouds to accurately register objects in 3d space . to deal with this shortcoming , we introduce a novel architecture called _ iterative hough forest with histogram of control points _ that is capable of estimating the 6d pose of occluded and cluttered objects given a candidate 2d bounding box . our _ iterative hough forest ( ihf ) _ is learnt using parts extracted only from the positive samples . these parts are represented with _ histogram of control points ( hocp ) _ , a `` scale - variant '' implicit volumetric description , which we derive from recently introduced implicit b - splines ( ibs ) . the rich discriminative information provided by the scale - variant hocp features is leveraged during inference . an automatic variable size part extraction framework iteratively refines the object s initial pose that is roughly aligned due to the extraction of coarsest parts , the ones occupying the largest area in image pixels . the iterative refinement is accomplished based on finer ( smaller ) parts that are represented with more discriminative control point descriptors by using our _ iterative hough forest_. experiments conducted on a publicly available dataset report that our approach show better registration performance than the state - of - the - art methods . +    object registration , 6 dof pose estimation , scale - variant hocp features , one class training , random forest , iterative refinement .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the cataclysmic variables with very low mass transfer rates ( 10@xmath0 - 10@xmath1 m@xmath2 yr@xmath3 ) and highly magnetic ( @xmath410 mg ) white dwarfs have been termed larps ( low accretion rate polars ; schwope et al . since the low mass transfer means the accretion luminosity is low , the individual stellar components can be viewed . thus , these systems provide a unique comparison to widely separated binaries for the study of how angular momentum losses and heating of the white dwarf affect the binary evolutionary scenario . while normal polars can reach mass accretion rates close to those of larps when they undergo long periods of low mass transfer , the main difference between larps and normal polars in low states appears to be the temperature of the white dwarf . ultraviolet observations of normal polars have shown temperatures of 11 - 14,000k for their white dwarfs ( araujo - betancor et al . in contrast , the white dwarfs in the larps are less than 10,000k ( schmidt et al . 2005 , 2007 , hereafter s05 , s07 ; vogel et al . 2007 , hereafter v07 ; schwope et al . 2009 , hereafter sw09 ) . these low temperatures imply that compressional heating has not taken place and so these objects have been suggested to be pre - polars ( s05 ) or preps ( sw09 ) . figure 5 of sw09 shows a nice plot of the temperatures of 9 preps compared to normal polars and accreting non - magnetic white dwarfs that illustrates this difference . however , the terminology and behavior is not always so cleanly separated . the system sdss121209.31 + 013627.7 ( schmidt et al . 2005b ; burleigh et al . 2006 ; hereafter b06 ) also has a low accretion rate and a low white dwarf temperature but has a brown dwarf secondary . its evolutionary path has been suggested to be either a polar in a low state or a prep . in addition , the recent observation of enhanced activity levels in the larp sdssj204827.9 + 005008.9 ( honeycutt et al . 2010 ) suggest that it may not be a larp / prep or that even preps can have different levels of activity . a secondary criteria for a larp / prep is the secondary underfilling its roche lobe so that roche lobe overflow does not occur ( s05 ; s07 ; sw09 ) . this is generally determined from the system having an spectral type too late to fill the lobe if it has a normal size for its type , as well as the absence of narrow components in the emission lines that can be traced to a stream . however , the honeycutt et al . ( 2010 ) study found evidence for a stream in sdssj204827.9 + 005008.9 , although they could not determine if the system had entered an increased state of mass transfer at the time . while the low accretion rates and sizes of the secondaries imply no roche lobe overflow , there could likely still be some accretion from the stellar wind of the secondary that is funnelled by the high field of the white dwarf ( li et al . the work of webbink & wickramasinghe ( 2005 ) using considerations of energy , magnetic fields , separations and lifetimes supports a wind model for the larps . in this study , we are not so concerned with the correct classification of our objects as larps , preps or normal polars in low states , or how the accretion might occur , but rather with the effects of low levels of accretion on the white dwarf . these effects are visible through time - resolved uv observations . accretion spots with temperatures of 30,000 - 70,000k are easily visible on normal polars even during their sporadic low states of accretion with @xmath510@xmath6m@xmath2 yr@xmath3 ( araujo - betancor et al . 2005 ; gnsicke et al . 2006 ) as large amplitude modulations of the uv light . however , it was somewhat surprising to find that similar modulations were apparent for ef eri ( szkody et al . 2006 ) which had been in a low state for seven years with @xmath7 of 10@xmath0 m@xmath2 yr@xmath3 and the larp mq dra with an even lower @xmath7 of 10@xmath1 m@xmath2 yr@xmath3 ( szkody et al . 2008 ) . while the uv light curves could be approximately matched with either hot spots or cyclotron components ( campbell et al . 2008 ) , each interpretation has its problems : the hot spots require different sizes and geometries while the cyclotron origin requires much higher magnetic fields in the uv than are apparent in the optical . in order to further study the effects of this low level accretion , we used _ galex _ to obtain uv observations of two additional larps ( wx lmi and sdssj103100.55 + 202832.2 ) and the system sdsssj121209.31 + 013627.7 . for convenience , we will refer to the sdss objects as sdss1031 and sdss1212 . the known parameters for the three systems are listed in table 1 . while all three of these objects have comparable low mass accretion rates , they present differences in magnetic field strength , white dwarf temperature , spectral type of secondary and orbital period . the observational data include ground - based apo @xmath8 light curves and fuv and nuv _ galex _ light curves . the _ galex _ observations took place on 2008 feb 16 , 23 and april 5 ( table 2 provides a summary of the times ) . images in both the near ultraviolet ( nuv ) detector ( 1750 - 2800 ) and the far ultraviolet ( fuv ) detector ( 1350 - 1750 ) ( martin et al . 2005 ) were obtained for each source . the time - tag data were calibrated in 240s intervals ( 120s for the nuv of sdssj1212 ) and phased according to the spectroscopic ephemeris given in v07 for wx lmi , and to the photometric ephemeris given in b06 for sdss1212 . the phasing for sdss1212 is accurate to @xmath90.14 phase at the time of the _ galex _ obseravtaions . for sdss1031 , the ephemeris is not known with enough accuracy to phase to prior data so arbitrary phasing was used with the period of s07 . the phased datapoints were then binned into 10 , 15 and 20 phase bins to optimize the s / n vs time resolution . the magnitudes for each binning were then measured using a 9 pixel radius aperture with the iraf routine _ qphot _ , using an average sky over the entire phase interval . flux conversions for the modeling were done using the values given in the _ galex _ online documentation ( fuv m@xmath10=18.82=1.40@xmath1110@xmath12ergs @xmath13 s@xmath14 @xmath3 and nuv m@xmath10=20.08 = 2.06@xmath1110@xmath15 ergs @xmath13 s@xmath3 @xmath16 ) . optical photometry in 2008 april and may ( table 2 ) was obtained for wx lmi and sdssj1031 using the nmsu 1 m telescope at apache point observatory . differential light curves in @xmath17 filters were made from nearby comparison stars on the same ccd frames . for wx lmi , one of the comparison stars used was the calibrated star a from schwarz et al . ( 2001 ) so that the differential magnitudes could be transformed to actual @xmath17 values . as in our previous spot models of ef eri and mq dra , the binsyn program suite ( linnell & hubeny 1996 ) was used to calculate synthetic light curves that could fit the observations . recent papers ( linnell et al . 2007 ; 2008 ) illustrate details of its application to cataclysmic variables , and hoard et al . ( 2005 ) describe its use in simulating hot spots on the magnetic white dwarf yy dra . while a hot spot may be multi - temperature and have a complicated shape as well as radiation influenced by the strong magnetic field , we attempted to minimize the free parameters by using a first approximation to the spot parameters with a circular , isothermal hot spot ( with a white dwarf spectral energy distribution ) . as our intent is to determine if a hot spot alone can account for the fuv and nuv light curves , we ignore cyclotron effects . the harmonics will affect the optical light curves to various degrees ( depending on where they appear in the various filters ) and can add to or dominate the observed amplitude of variation from the spot .    in general , if the objective of the analysis is only to fit the light variations ( the absolute fluxes are unknown ) , the light curves are normalized to a light maxima of 1.0 , and then binsyn can be used to calculate monochromatic light curves . in this mode , the program represents the stellar components as a sum of black body contributions from local photospheric values of @xmath18 determined by adopted values of bolometric albedos and gravity darkening exponents . no synthetic spectra are needed .    however , the wx lmi light curves required special considerations , since those light curves are absolute flux measurements . their simulation , via synthetic photometry from synthetic spectra , produced calculated absolute flux values outside the earth s atmosphere , using transmission profiles of the filters and the quantum efficiency of the ccd , and the distance of the system . in previous papers the synthetic spectra input to binsyn were calculated using tlusty ( hubeny 1988 ) and synspec ( hubeny et al . 1994 ) . however , the hubeny tlusty wd models do not extend to cooler than 8000k . fortunately , cool wd models have been calculated by bergeron et al . ( 1991 ; 1995 ) . bergeron kindly provided a set of 40 wd synthetic spectra with log  @xmath19=8.0 and @xmath18 ranging from 17,000k to 1500k . each synthetic spectrum includes eddington flux values ( @xmath20 ) at 1788 unequally - spaced wavelengths between 43  and 95,000 . synthetic spectra similar to these have been used in the synthetic wd photometry of holberg & bergeron ( 2006 ) . we produced corresponding synthetic spectra in wavelength units ( @xmath21 ) . we input wd synthetic spectra with @xmath18 values of 6000k , 6500k , 7000k , 7500k , and 8000k to binsyn ; the program then used an adopted wd @xmath18 , within the above range , and interpolated among the input spectra to produce a synthetic wd spectrum of the assigned @xmath18 . as representation of the hot spots on the wd ( the impact regions of the incoming mass from the secondary ) requires additional synthetic wd spectra , we added input wd synthetic spectra with @xmath18 values of 10,000k , 11,000k , and 12,000k . synthetic spectra to represent the secondary star present a separate problem . if the secondary fills its roche lobe , @xmath18 and log  @xmath19 vary over the photosphere as determined by the values of gravity darkening and bolometric albedo ( temperature variation due to irradiation by the wd is negligible ) . as binsyn needs synthetic spectra for cool photospheres over a range of @xmath18 and log  @xmath19 , we downloaded a database of marcs   synthetic spectra for this purpose . each marcs synthetic spectrum tabulates physical flux values ( @xmath21 ) for 101,125 wavelengths unequally spaced between 900  and 200,007  and we then prepared corresponding spectra in eddington flux units . to permit accurate interpolation to the local photospheric values of @xmath18 and log  @xmath19 on the secondary star we provided binsyn with 24 marcs spectra for @xmath18 values of 2500k , 2700k , 2900k , 3100k , 3300k , and 3400k , and for each of these , log  @xmath19 values of 3.0 , 4.0 , 4.5 , and 5.0 . the single secondary component point nearest the l1 point has an associated photospheric segment whose flux contribution is negligible . we include an input black body spectrum of 1500k which , together with the other input spectra for the secondary star , brackets this point . marcs spectra span all other grid segments covering the secondary star . the binsyn output includes a synthetic spectrum for the system , the wd , and the secondary star , all for an adopted orbital longitude and orbital inclination , @xmath22 . wx lmi was found in the hamburg quasar survey by reimers et al . ( 1999 ) who determined the orbital period , the spectral type of the secondary and the distance . from the spectrum and photometry , they concluded there were two accretion spots , with magnetic fields of 60 and 68 mg and the accretion rate was about 100 times less than normal polars . further long term photometry was obtained by schwarz et al . ( 2001 ) and modeled to determine the colatitudes ( 90 - 130@xmath23 ) and azimuths ( -70 to -80 , and + 90@xmath23 ) of the spots . these locations placed the spots slightly below the equator of the white dwarf and further from the line - of - centers between the stars than in normal high accretion rate polars . v07 used _ xmm - newton _ to obtain x - ray data as well as some near uv photometry using the optical monitor ( om ) uv filters at 2910 and 3440 . they also used optical spectra to establish an upper limit to the temperature for the white dwarf of 8000k and to constrain the cyclotron model and accretion rate . they concluded that the low accretion rate was consistent with accretion resulting from wind loss from the secondary . their model found small ( 40 km radius ) accretion spots at similar locations to that of schwarz et al . ( 2001 ) but the x - ray vs uv modulations and the hardness ratios were not consistent with heated spots . the x - ray emission was assumed to originate from both coronal activity on the secondary and from the accretion regions . the _ galex _ data explore the regime between the x - ray and om filters . the uv light curves ( with 15 phase bins ) are shown in figure 1 along with the optical data obtained about 2 months later . the double - humped light curve variation that has been taken to be evidence of the accretion areas undergoes an interesting variation with wavelength . figure 1 shows the progression from a double - humped curve in the i filter to a very low level of variability in the optical @xmath24 and then an increasing variation as the fuv is reached . the phasing here is spectroscopic ( phase 0 is inferior conjunction of the secondary ) so the primary pole is evident between phase 0.8 and 0.9 , the secondary pole near phase 0.3 and both poles are self - eclipsed near phase 0.6 . the fuv shows slightly higher flux for the secondary pole than the primary , which corroborates what v07 regarded as a puzzle i.e. that the secondary spot seems hotter than the primary one . the dips at phase 0.3 in @xmath25 and @xmath24 are related to the cyclotron features at 5950  and 5220  which appear in these filters and change strength during the orbit due to the changing viewing angle . v07 determined or adopted some of the parameters needed for a simulation of the wx lmi system including the orbital period , @xmath26 ; @xmath27(adopted ) ; @xmath28 ; @xmath29k ; and @xmath30k , appropriate to an m4.5 spectral type . they also determined some of the spot parameters associated with the magnetic field , including @xmath31 , the spot colatitudes , and @xmath32 , the spot azimuths . these parameters are listed in table 3 . the mass transfer rate ( @xmath7 ) in wx lmi ( @xmath33 ; v07 ) , is listed as a system parameter in table 3 but is not needed in the analysis . @xmath28 , listed in v07 ( table  2 ) for a secondary star of spectral type m4.5 , differs slightly from the calibration of knigge ( 2006 ; 2007 ) , @xmath34 at the wx lmi orbital period , but we preserve the value @xmath35 for consistency with the other v07 parameters . we have assigned error bars in table  3 from the range of @xmath36 values in table  2 of v07 . we adopt a secondary component model which fills its roche lobe ( @xmath37 , table  3 ) . v07 describe wx lmi as a pre - polar candidate , appropriate to the very small @xmath7 , and underfilling its roche lobe . it is convenient for our analysis to assume a secondary star which fills its roche lobe ; after obtaining a solution we tested the sensitivity of our assumption by performing a simulation in which the secondary underfills its roche lobe . the adopted value of @xmath27 has no assigned error bars . the radius of a zero temperature homogeneous hamada - salpeter carbon @xmath38 wd is @xmath39 ( panei et al . 2000 ) ; we adopt this value with no change for the very small correction to a 8000k wd . the wd roche potential , @xmath40 , table  3 , produces the required @xmath41 . the table  3 @xmath29k was determined by v07 but they state that the true value is less than 8000k ; they were unable to pursue this topic because of their lack of lower temperature models . v07 show that @xmath22 is less than @xmath42 ; the presence of detectable ellipticity in the light curve suggests a fairly large @xmath22 ; we adopt @xmath43 . figure  2 shows the final system synthetic spectrum at orbital phase 0.0 . this synthetic spectrum can be compared with the observed spectrum ( figure  1 of v07 ) . we calculated and stored system synthetic spectra for 33 orbital longitudes for each test simulation . each synthetic spectrum extended from 800  to 10,000  and included eddington flux values at 1  spacing , produced by interpolation within the input spectra to binsyn and integrated , within binsyn , over the photospheric segments visible to the observer , and with an adopted limb darkening of 0.6 . our initial simulation neglected any hot spots on the wd as we expected possible hot spots only to make a negligible contribution to the @xmath8 light curves . we subsequently verified this expectation . the apo filter transmission curves and the quantum efficiency curve of the ccd used for the ground - based photometry were digitized by hand . the ccd quantum efficiency drops to 0.0 at 10,000   and this sets the upper wavelength limit for the simulation . we also produced comparable transmission curves for the fuv and nuv filters using the _ galex _ calibrations . program synphot , within the binsyn package , produces the integrated physical flux for each filter resulting from a system synthetic spectrum incident on the apo filter set and ccd , given by the expression    @xmath45    where @xmath46 is the photometric system response function and @xmath47 is the synthetic spectrum incident flux . numerical integration requires a value for the filter transmission and ccd response at each tabular wavelength of the system synthetic spectrum . we used the numerical recipes ( press et al . 1986 ) program * spline * to fit each digitized filter transmission curve as well as the ccd response curve , and the numerical recipes routine * splint * to produce the required transmission or response at the synthetic spectrum wavelengths . the quantity to compare with an observed light curve is the integrated flux divided by @xmath48 , where @xmath49 is the distance to the system in cgs units . we initially adopted the v07 determination of 100  pc . the initial test showed rough agreement between simulation and observed light curves with variable cyclotron contributions . the calculated light curves were systematically too bright in @xmath50 , and @xmath51 , with a strong cyclotron contribution in @xmath25 . we varied the assumed distance to obtain an optimum overall fit by eye and determined a revised distance of 107  pc . this value is in very close agreement with the v07 value of 104  pc in their table  2 for the m4.5 spectral type adopted for the secondary star in this study . the final light curves , compared with the observations , are shown in figure  3 ( note the change in ordinate scale in figure  3d ) . cyclotron emission dominates the @xmath25 band light curve . this is in agreement with the v07 cyclotron model , their figure  7 . the same figure  7 indicates some cyclotron contribution to the @xmath50 , and @xmath51 light curves . figure  3d shows that the synthetic light curve has slightly too much flux for the adopted distance of 107  pc . as figure  2 shows , the secondary component spectrum rises rapidly to a peak near 10,000 . a small departure of the actual wx lmi secondary component from the adopted model could explain this residual . the inclination can not be increased more than a few degrees without producing eclipses ; a reduction in the inclination decreases the light curve amplitude and increases the absolute flux level . we regard the adopted @xmath52 as optimum , considering the light curve contributions by cyclotron emission . figure  7 of v07 makes no prediction concerning cyclotron emission in the fuv or nuv , and the contribution of the secondary component becomes negligible shortward of about 3500  so it does not contribute in the nuv or fuv . we adopted circular models for the two hot spots and experimented with their @xmath18 values and angular radii . initial simulations suggested a wd @xmath53k . initial choice of spot @xmath18 values of 12,000k produced a fuv light curve amplitude that was much too large . lower temperature spots require larger radii to preserve the light curve amplitude , particularly in the nuv . eventually we found an acceptable fit with spot radii of @xmath54 and @xmath18 values of 10,000k . a final correction of the wd @xmath18 to 7900k produced the best fit , by eye ( figure  2 has a wd with @xmath18 of 7900k ) . the fuv and nuv light curves are shown in figure  4 and the final spot parameters are in table  4 . the spots are large ; a view of the wd at orbital phase 0.0 is in figure  5 . we suggest that this result is not inconsistent with the small @xmath7 and impacting stream from the secondary component wind rather than a stream emerging through the l1 point . while our spot sizes are much larger than those of v07 , ( their estimate of the accretion area was about 40 km in radius , whereas our 32@xmath23 spots are about 4700 km in radius ) , their estimate was based on the x - ray and beamed cyclotron regions . these areas are typically much smaller spots than the heated areas visible in the uv ( gnsicke et al . 2006 ) . time - resolved uv spectropolarimetry will be required to obtain further knowledge of cyclotron effects in the uv . for the field strength of 68 mg , the n=6 - 8 harmonics will be in the nuv band . this could account for some of the deviations of the fit in this band .    as a check , we recalculated the @xmath55 , and @xmath51 light curves with and without the spots included in the model . there was a barely detectable spot contribution to those light curves , far too small to require a revision of the model . our model does not prove that hot spots are the cause of the nuv and fuv light variation , it only demonstrates that our adopted parameters can approximate the observed light variations . the deviations likely indicate that the spots have more complicated geometries than our simple assumption of a circular shape . our model assumes the secondary component fills its roche lobe . as a test we changed the value of @xmath56 ( table  3 ) , so that the volume of the secondary was 95% of the initial model value which filled the roche lobe . recalculation of the @xmath44 light curves showed that the amplitudes are almost unaffected . a slight amplitude reduction of the @xmath51 light curve is barely detectable . further reducing the volume to 90% produces the changes shown in figure 6 . thus , the secondary can underfill the roche lobe by a small amount without reducing the theoretical light curve amplitudes sufficiently to be in disagreement with the observed light curve . however , if it underfills by a large factor , the light curve amplitude , especially in @xmath51 , would be too small due to the reduced ellipticity of the secondary . if there is some contribution of the n=2 harmonics in the @xmath51 filter , the cyclotron would dilute the ellipsoidal variation from the secondary , which would necessitate an even larger filling factor , although the variation from the cyclotron would then add to the variablility as well . as with the majority of the larps , sdss1031 was discovered by searches through the sdss spectroscopic survey ( s07 ) . the variations in the polarimetry and photometry determined a short orbital period and the spacing of the cyclotron harmonics revealed the magnetic field . the spectral fluxes between the cyclotron harmonics placed limits on the white dwarf temperature and the late spectral type of the secondary ( table 1 summarizes these values ) . its very short period makes it unusual among the larps ( figure 5 of sw09 ) , and may present problems for its evolution as a prep . however , its low accretion rate , cool white dwarf and late type secondary which underfills its roche lobe ( s07 ) all match the primary characteristics of larps . sdss1031 shares many characteristics of the short period normal polar ef eri which has mostly been is a low state since 1997 but does have occasional high states ( howell et al . 2006 ) . the uv through optical light curves of sdss1031 are shown in figure 7 , where the arbitrary phase in uv and optical is set so that maximum light occurs roughly at phase 0 . while sdss1031 is the faintest of our three sources in the optical , the fuv magnitudes are comparable to wx lmi . the @xmath57 filter shows the cyclotron variation found by s07 with some evidence of a double - humped shape indicative of beaming in an optically thin high cyclotron harmonic . the @xmath25 filter shows the least orbital variation while @xmath24 has an increased amplitude . the amplitude of the uv orbital modulation is much lower than in wx lmi and more asymmetical in shape , showing a gradual rise to maximum followed by a steeper decline . if the variation is due to a hot spot , the similar amplitudes in fuv and nuv mean that the spots must be different sizes and simple circular isothermal spots will not provide a good fit to the observations . to get a crude idea of what kind of spot would be needed , the binsyn model was computed for the parameters shown in table 4 and fit to the _ galex _ data binned by 10 in phase ( figure 8) . the best fit obtained was with a 13,000k spot at -65@xmath23 latitude . the spot radius was 10@xmath23 for the fuv and 23@xmath23 for nuv . however , as expected , the shape of the variation is not fit well . like sdss1031 , this object was found in searching through the sdss spectra ( schmidt et al . 2005b ) because the balmer lines showed zeeman splitting from a magnetic white dwarf as well as weak h@xmath58 emission . spectroscopy by these authors revealed a period near 90 min but no secondary star was detected and an ir measurement at @xmath59 implied a brown dwarf object of spectral type l5 or later . further ir observations by debes et al . ( 2006 ) refined the orbital period and showed a variation in the @xmath60 band that could be explained with cyclotron emission from a small spot on the white dwarf . farihi , burleigh & hoard ( 2008 ) modeled further ir data with an l8-t1 substellar dwarf and a white dwarf with a 7 mg field . b06 obtained x - ray observations with _ swift _ which detected the source but with too low of a count rate to ascertain variability . the flux was best fit with a thermal plasma with kt=1.9 kev and a luminosity of 10@xmath61 erg s@xmath3 . while this luminosity is comparable to that detected in the larp mq dra ( szkody et al . 2004 ) , and the low state of the polar ar uma ( szkody et al . 1999 ) , it is far larger than the x - ray luminosities of l dwarfs ( stelzer et al . the @xmath62 photometry of b06 shows a single hump that increases in amplitude in the blue , with some slight indication of a second small peak half a phase later that was pointed out by koen & maxted ( 2006 ) . the origin of the second peak is not clear ; koen & maxted ( 2006 ) suggested a reflection effect from the secondary ( since the feature seemed to increase at longer wavelengths ) , while b06 mention a possible second pole ( but this is not evident in the ir photometry of debes et al . ( 2006 ) .    b06 model the primary optical hump as a 14,000k hot spot on the white dwarf . however , as they point out , it is not clear what provides the accretion necessary to produce the x - rays and a hot spot . a brown dwarf that underfills its roche lobe should not produce an accretion stream nor a wind outflow . recent photometry by howell et al . ( 2008 ) showed that the strength of h@xmath58 is variable , indicating that the secondary has some activity . it may be that the secondary stars in close binaries have higher activity levels than their single counterparts due to their faster rotations in a binary system . our _ galex _ data on sdss1212 ( figure 9 ) shows a similar single - peaked modulation with increased amplitude in fuv over the nuv and a possible hint of a second small peak near phase 0.5 . the nuv magnitude is the brightest of the three objects . trial models with binsyn using the b06 parameters for the @xmath22 , masses , and temperature of the white dwarf revealed that the width of the hump is a clear function of the spot latitude , with greater width when the spot is in the same hemisphere as the visible pole . a good fit ( figure 10 bottom ) to the nuv light curve was obtained with a 14,000k spot ( the same temperature found by b06 from the optical light curves ) with a radius of 9@xmath23 , but this spot produced a fuv amplitude that was too large . thus , as in mq dra ( szkody et al . 2008 ) , a single spot size and temperature can not fit both the nuv and fuv . if the spot temperature is kept at 14,000k , a smaller spot of 7@xmath23 fits the fuv reasonably well ( figure 10 top ) . the observed light curve shows some scatter and there is a hint that the light curve is asymmetric with a slower rise than decline . the spot sizes are 1% of the white dwarf surface in fuv and 1.7% in nuv , as compared to the 5% found by b06 from their optical modeling . these numbers follow the trend of decreasing spot sizes with smaller wavelengths .    increased time - resolution and s / n will be needed to further understand the reality and nature of the possible peak near phase 0.5 . it is intriguing that this feature is evident from uv through the @xmath57 band . these phases need to be targeted to determine if there is accretion related to a second pole . our uv observations of three polars with extremely low accretion rates enforce the results found from observations of ef eri and mq dra i.e. that all the white dwarfs have areas of enhanced emission even with these low rates , indicating some accretion is still occuring . despite a large range in magnetic field ( 7 - 70 mg ) , the field is apparently strong enough in all three systems to funnel the accreting material to the magnetic pole(s ) of the white dwarf . this appears to happen even in the case of sdss1212 , which likely has a brown dwarf secondary . thus , if the accretion is via a stellar wind from the secondary that is trapped by the field of the white dwarf , it is difficult to provide a wind of this type for sdss1212 . time - resolved spectra that would enable doppler tomography might provide a resolution of this issue . during high states of accretion , the stream flow is usually visible as a narrow component in the balmer emission lines which is mapped to the stream . in normal polars during their low states , high velocity components in the lines have been interpreted as structures similar to prominences close to the secondary star ( kafka et al . 2007 , 2008 ) . however , the balmer emission in larps is weak to non - existent due to the low the accretion rates low so it would be difficult to construct the map . if the increased emission evident at some phases is due to hotter temperature , simple spots of 10 - 000 - 14,000k covering a few percent of the white dwarf surface can approximate the uv light curves , although the geometries of the spots require a more complex model shape than simply circular . alternatively , if there is a large range in field strengths in the white dwarfs , there could be contributions to the uv from higher cyclotron harmonics . improved uv cyclotron models will be needed to test this possibility . since these objects are too faint for _ galex _ spectra or high s / n time - resolved uv light curves that would merit detailed models of the shape or of pursuing cyclotron radiation , further work will require uv spectra and polarimetry with larger telescopes . support for this research was provided by nasa _ galex _ grant nnx08am07 g . we gratefully acknowledge pierre bergeron for providing cool white dwarf models , the marcs database for access to their spectra , eric bullock for work with some of the data files , and mark klaene for providing optical filter curves . araujo - betancor , s. et al . 2005 , , 622 , 589 bergeron , p. , wesemael , f. & fontaine , g. 1991 , , 367 , 253 bergeron , p. , saumon , d. & wesemael , f. 1995 , , 443 , 764 burleigh , m. r. et al . 2006 , , 373 , 1416 ( b06 ) campbell , r. k. , harrison , t. e. , schwope , a. d. & howell , s. b. , 672 , 531 debes , j. h. , lopez - morales , m. , bonanos , a. z. & weinberger , a. j. 2006 , , 647 , l147 farihi , j. , burleigh , m. r. & hoard , d. w. 2008 , , 674 , 421 gnsicke , b. t. , long , k. s. , barstow , m. a. & hubeny , i. 2006 , , 639 , 1039 hoard , d. w. , linnell , a. p. , szkody , p. & sion , e. m. 2005 , , 130 , 214 holberg , j. b. & bergeron , p. 2006 , , 132 , 1221 honeycutt , r. k. , kafka , s. & tappert , c. 2010 , baas , 215 , 416.19 howell , s. b. , walter , f. m. , harrison , t. e. , huber , m. e. , becker , r. h. & white , r. l. 2006 , , 652 , 709 howell , s. b. , harrison , t. e. , huber , m. e. , szkody , p. , walter , f. m. & harbeck , d. 2008 , , 136 , 2541 hubeny , i. 1988 , comp . comm . , 52 , 103 hubeny , i. , lanz , t. , & jeffery , c.  s. 1994 , in newsletter on analysis of astronomical spectra no . s. jeffery ( ccp7;st . andrews : st . andrews univ . ) , 30 kafka , s. , howell , s. b. , honeycutt , r. k. & robertson , j. w. 2007 , , 133 , 1645 kafka , s. ribeiro , t. , baptista , r. , honeycutt , r. k. & robertson , j. w. 2008 , , 688 , 1302 koen , c. & maxted , p. f. l. 2006 , , 371 , 1675 knigge , c. 2006 , , 373 , 484 knigge , c. 2007 , , 382 , 1982 li , j. , wu , k. w. & wickramasinghe , d. t. 1994 , , 268 , 61 linnell , a. p. & hubeny , i. 1996 , , 471 , 958 linnell , a. p. , godon , p. , hubeny , i. , sion , e. m. & szkody , p. 2007 , , 662 , 1204 linnell , a. p. , godon , p. , hubeny , i. , sion , e. m. , szkody , p. & barrett , p. e. 2008 , , 676 , 1226 martin , d. c. et al . 2005 , , 619 , l1 panei , j. a. , althaus , l. g. & benvenuto , o. g. 2000 , , 353 , 970 press , w. h. , flannery , b. p. , teukolsky , s. a. & vetterling , w. t. 1986 , ( cambridge ; university press , cambridge ) reimers , d. , hagen , h .- j . & hopp , u. 1999 , , 343 , 157 schmidt , g. d. et al . 2005 , , 630 , 1037 ( s05 ) schmidt , g. d. et al . 2005b , , 630 , l173 schmidt , g. d. et al . 2007 , , 654 , 521 ( s07 ) schwarz , r. , schwope , a. d. & staude , a. 2001 , , 374 , 189 schwope , a. d. et al , 2002 , in asp conf . vol 261 , the physics of cataclysmic variables and related objects , ed . b. t. gnsicke , k. beuermann , & k. reinsch ( san francisco : asp ) , 102 schwope , a. d. , gomez - moran , a. n. , schreiver , m. r. & gnsicke , b. t. 2009 , , 500 , 867 ( sw09 ) stelzer , b. , micela , g. , flaccomio , e. r. n. & jayawardhana , r. 2006 , , 448 , 293 szkody , p. et al . 1999 , , 520 , 841 szkody , p. et al . 2004 , , 128 , 2443 szkody , p. et al . 2006 , , 646 , l147 szkody , p. et al . 2008 , , 683 , 967 vogel , j. , schwope , a. d. & gnsicke , b. t. 2007 , , 464 , 647 ( v07 ) webbink , r. f. & wickramasinghe , d. t. 2005 , in asp conf . 330 , the astrophysics of cataclysmic variables and related objects , ed . j. m. hameury & j. p. lasota ( san francisco : asp ) , 137    lccc p@xmath63 ( hr ) & 2.78 & 1.37 & 1.47 + optical mag & @xmath25=16.97 & @xmath19=18.26 & @xmath19=17.99 + t@xmath64 ( k ) & @xmath658000 & 9000 & 9500 + b ( mg ) & 61,70 & 42 & 7 + sec type & dm4.5 & @xmath66dm6 & l8-t1 + d ( pc ) & 100 & 270 - 380 & 120 + @xmath7 ( @xmath1110@xmath0m@xmath2 yr@xmath3 ) & 1.5 & 1.5 - 3 & 1 + ref & 1 , 2 & 3 & 4 , 5 , 6 +    lccl wx lmi & feb 16 & 10:23 - 20:42 & _ galex _ nuv , fuv 7 visits + wx lmi & apr 29 & 03:57 - 08:56 & @xmath17 39 exposures + sdss1031 & feb 23 & 10:01 - 20:19 & _ galex _ nuv , fuv 7 visits + sdss1031 & may 5 & 05:04 - 08:15 & @xmath67 12 exposures + sdss1212 & apr 5 & 11:14 - 23:01 & _ galex _ nuv , fuv 8 visits +    llll @xmath68 & @xmath69 & @xmath70(pole ) & @xmath71k + @xmath72 & @xmath73 & @xmath70(point ) & 1640k + @xmath74 & @xmath75 & @xmath70(side ) & 3259k + p & 0.11592364 days & @xmath70(back ) & 3133k + @xmath49 & @xmath76 & @xmath77(pole ) & @xmath78 + @xmath79 & 76.65 & @xmath77(point ) & @xmath80 + @xmath81 & 2.4625438 & @xmath77(side ) & @xmath82 + _ i _ & @xmath83 & @xmath77back ) & @xmath84 + @xmath85 & @xmath86k & log @xmath87(pole ) & 4.95 + @xmath88 & @xmath39 & log @xmath87(point ) & 1.16 + log @xmath89 & 8.05 & log @xmath87(side ) & 4.89 + @xmath90 & 1.0 & log @xmath87(back ) & 4.67 + @xmath91 & 0.6 & @xmath92(1,2 ) & @xmath93 + @xmath94 & 0.25 & @xmath95(1,2 ) & @xmath96 + @xmath97 & 0.08 +    lccc wd t@xmath98 & 7900k & 9500k & 9500 + wd mass ( m@xmath2 ) & 0.6 & 0.6 & 0.6 + wd radius ( r@xmath2 ) & 0.012 & 0.012 & 0.012 + wd log g & 8.0 & 8.0 & 8.0 + sec mass ( m@xmath2 ) & 0.18 & 0.08 & 0.08 + i ( @xmath23 ) & 70 & 70 & 70 + spot t@xmath98 & 10,000 ; 10,000k & 13,000k & 14,000k + spot latitude & -55 ; -45@xmath23 & -65@xmath23 & -15@xmath23 + spot ang rad fuv & 32@xmath23 ; 32@xmath23 & 10@xmath23 & 7@xmath23 + spot ang rad nuv & 32@xmath23 ; 32@xmath23 & 23@xmath23 & 9@xmath23 +
Abstract:

Answer: _ galex _ near ultraviolet ( nuv ) and far - ultraviolet ( fuv ) light curves of three extremely low accretion rate polars show distinct modulations in their uv light curves . while these three systems have a range of magnetic fields from 13 to 70 mg , and of late type secondaries ( including a likely brown dwarf in sdssj121209.31 + 013627.7 ) , the accretion rates are similar , and the uv observations imply some mechanism is operating to create enhanced emission zones on the white dwarf . the uv variations match in phase to the two magnetic poles viewed in the optical in wx lmi and to the single poles evident in the optical in sdssj1212109.31 + 013627.7 and sdssj103100.55 + 202832.2 . simple spot models of the uv light curves show that if hot spots are responsible for the uv variations , the temperatures are on the order of 10,000 - 14,000k . for the single pole systems , the size of the fuv spot must be smaller than the nuv and in all cases , the geometry is likely more complicated than a simple circular spot .