 focal plane by dividing the images by a master flat-field image. We adopted the flat-fielding method that was used by the Kilo Degree Survey (KiDS, [@deJong2015]): the master flat-field is achieved by first median-combining and normalizing 8 dome flat-fields and 8 twilight flat-fields, and then multiplying the averaged flat-fields with each other. High-frequency spatial Fourier modes are corrected using the dome flats and low frequencies using the twilight flats. This is based on the pre-assumption that the large-scale illumination of the twilight flat-field matches better the observational situation than that obtained from the dome. On the other hand, in the dome flats the [*S/N*]{} is high, which can be used to capture the pixel-by-pixel variations in the pixel sensitivity.

However, even after applying the flat-field correction small systematic flux variations remain across the instrument. These variations can be corrected by applying an illumination correction. We use the correction models made for the KiDS (see [@Verdoes-Kleijn2013] for details). The models are made by mapping the photometric residuals across the CCD array using a set of dithered Landolt Selected Area (SA) field [@Landolt1992] observations, and fitting a linear model to the residuals. The images are multiplied with this illumination correction. The correction is applied after sky subtraction (see the next sub-section), to avoid the sky residuals being amplified by the illumination correction procedure.

Background subtraction and de-fringing
--------------------------------------

The images contain sky background flux composed of direct and scattered atmospheric emission, and scattered light of bright celestial sources. A careful removal of the atmospheric background light is essential when studying LSB objects such as UDGs. When considering a single image at low luminosity levels, we cannot directly tell which part of the light is diffuse light coming from the sources and which is background light. To bypass the problem, we can make an assumption that the pattern of background light stays constant if the telescope pointing direction is not changed by more than a few degrees. Due to the large dithers between the consecutive integrations the objects are not likely to appear twice in the same pixel, which allows us to produce background models by averaging a stack of images. The intensity of the sky changes thorough the night (1–10% between exposures), especially at the beginning and at the end of the night, which forces to scale the images used for the background model before
focal plane by dividing the images by a master bland - sphere image. We adopted the flat - fielding method acting that was used by the Kilo Degree Survey (KiDS, [ @deJong2015 ] ): the master flat - plain is achieved by first median - combining and normalizing 8 attic bland - fields and 8 twilight bland - fields, and then multiplying the average bland - fields with each early. High - frequency spatial Fourier mode are corrected using the dome flats and broken frequencies using the twilight flats. This is based on the pre - premise that the large - scale illumination of the twilight flat - field match better the observational situation than that obtained from the dome. On the early hired hand, in the dome flats the [ * S / N * ] { } is high, which can be used to capture the pixel - by - pixel variations in the pixel sensitivity. 

 However, even after applying the flat - field correction minor systematic flux variation persist across the instrument. These variations can be corrected by practice an illumination correction. We use the correction model made for the KiDS (see [ @Verdoes - Kleijn2013 ] for details). The models are made by mapping the photometric residual across the CCD array using a set of dithered Landolt Selected Area (SA) field [ @Landolt1992 ] observations, and meet a analogue model to the residuals. The images are multiplied with this illumination correction. The correction is applied after sky subtraction (examine the next sub - section), to avoid the sky residual being amplified by the clarification correction procedure. 

 setting subtraction and de - fringe 
 -------------------------------------- 

 The images incorporate sky background flux density composed of direct and scattered atmospheric emission, and scattered lighter of bright celestial sources. A careful removal of the atmospheric background light is essential when studying LSB objects such as UDGs. When considering a single persona at low luminosity levels, we cannot immediately tell which part of the light is diffuse light coming from the sources and which is background lighter. To bypass the problem, we can make an assumption that the pattern of background light stay changeless if the telescope sharpen direction is not changed by more than a few degrees. Due to the large dithers between the consecutive integration the objects are not likely to appear doubly in the same pixel, which admit us to produce background models by average a stack of images. The intensity of the sky changes thorough the night (1–10% between exposures), particularly at the beginning and at the end of the night, which forces to scale the images used for the background mannequin before
 fofal plane by dividing tht images by a masjee flat-hield ijage. We xdopted the flat-fielding metiod rhat qas used by the Kilo Ddgree Surney (KiDS, [@eeJoig2015]): the master flef-field lf acgleved uy first median-gombining atd normalizing 8 aole flat-fields and 8 twilight flat-fiejds, and tjen multiplyind tht aderafvd flat-fields with each other. Hjgh-freqlency spatial Foutier modes are corrected uding the dome flats anf low frequghciqw using the gwilight flats. This is based on the pre-assumption thag the large-scalg jlpomination of the nwilight flat-npeld madches brtter the obsevvatimnao situation than that obtained from the dome. On dhz other hand, in the dimw flajs tha [*S/N*]{} us figg, xhidh can be used to czpture the pixel-by-pixel variatoogw in the pixem sensytyvity.

However, even after applying the fldt-fjeld correction small ststematic flux variatlons remayn across the instrument. These variations can be worrertdd yn xpppying an illumination correction. We use the cjdrtctpon models made fjr the KiDS (sfe [@Derdoes-Kleijn2013] for dzfajls). The models are made bi mappung the protoketric residuals across the CCD array ufung a set of dithexed Landolt Relevted Srea (SA) field [@Landolt1992] oyservafions, and fltting a mknear model to tfe gesiguals. Tht images are multi[lied witi thix illumknatoon cowrection. Tje covsection is applied aftet sky vubtractioj (see the next sub-section), to avoid the sky rexiguans being ampllfied by the iljumination cortection pxocedufe.

Backgrouhd subtcaction and qe-fringing
--------------------------------------

The lmages contamn sky baskgriund flux cuoposed of direvt and sccutered atmispheric emission, end rdattered light if vright celestiak suursed. E carqxul removal mf tfe xymospferic backgvouvd loght is essential whan sfudying LSB objectx fuch as ODGs. When considering s single image at pow lnminosmty lefelf, we cannot directly tell whicg part of thc light is disfusc lidht coming from the sources and which is backgrouid light. To bypass the problem, we can make sn assumptipn thwt the padtern of background oight stays constsnt if the telescope lointitg digection is not changed by more than a few degrees. Due to the large ditherw betwxeg the consedutife innegxatnons thq oboerts are not likelj to appear twice in the same pieel, which dljows us to produce backgrounc oodels by avetaging a stack of images. Tge intenxity of the sky changes thorougn the night (1–10% between exppsures), ewpecizlly an tie beginning and at yhe enb of tye nighr, whlch forces to rcake the imagfs uxeb for the background model before
focal plane by dividing the images by flat-field We adopted flat-fielding method that Degree (KiDS, [@deJong2015]): the flat-field is achieved first median-combining and normalizing 8 dome and 8 twilight flat-fields, and then multiplying the averaged flat-fields with each other. spatial Fourier modes are corrected using the dome flats and low frequencies using twilight This based the pre-assumption that the large-scale illumination of the twilight flat-field matches better the observational situation than obtained from the dome. On the other hand, the dome flats the is high, which can be to the pixel-by-pixel in pixel However, even after the flat-field correction small systematic flux variations remain across the instrument. These variations can be corrected by an illumination use the models for KiDS (see [@Verdoes-Kleijn2013] The models are made by mapping across the CCD array using a set of Landolt Selected (SA) field [@Landolt1992] observations, and fitting linear model to the residuals. The images are with this illumination correction. The correction is applied after sky subtraction (see the next sub-section), the sky residuals being by the illumination procedure. subtraction de-fringing The images sky background flux composed of direct and scattered atmospheric emission, and light of bright celestial sources. A careful removal of the light essential when studying objects such as UDGs. considering single image at low we directly of light diffuse light coming from sources and which is background To bypass the problem, that the pattern of background light stays constant the telescope pointing direction is not changed more than a few degrees. Due to the large dithers between the integrations the not likely to appear twice in the same which allows us to background models by averaging a stack of images. The of sky changes the night (1–10% exposures), especially at beginning and at of the which to used for the background model before
 focal plane by dividing the imAges by a masTer flAt-fIelD iMage. we adOpted the flat-fiELdinG method that was used by thE Kilo deGRee SURvEy (KiDs, [@deJong2015]): THe MASteR fLaT-fiElD Is AchieVed By first Median-combIniNg And normaliziNG 8 dOme flat-fieLds And 8 twilight fLat-Fields, AnD thEN multIplYing tHe averAGed flaT-fields wiTh EAch othER. High-frEQUeNcy sPatial Fourier modeS ArE Corrected using The domE fLAtS ANd lOw fRequencies UsIng thE TwilighT FlATS. thiS Is based on the pRe-assumptioN ThaT the laRgE-scALe illuMinatIoN Of tHe twilight fLat-fIeld matchEs bettER the obsERvationAl situAtiOn tHan tHAt ObTaiNeD FroM ThE doME. On The other HaNd, In the Dome FLATS the [*s/N*]{} iS higH, whicH can be used to cAptUre tHE piXel-by-Pixel VariAtIons iN the piXel seNsItivity.

However, eVen aFter applyIng ThE flAt-Field COrrectIon SmaLl systeMatic flUX vaRiATIOnS remain across the inStRUMeNt. These vAriatiONs CaN Be correcTeD by ApplYINg an iLlumINaTion corrEction. wE uSe The corrEcTion moDeLs mAde For thE kiDS (See [@VerDoes-KleiJn2013] for DEtails). The modelS Are made by mappINg THE pHOtomEtrIc residuals AcroSS the cCD aRRaY usINg a seT of diThEReD landolt Selected Area (sA) Field [@LAndolT1992] observations, And fitting A LINear modeL to tHE rESiduals. The imagEs are Multiplied WIth this iLlumiNation coRrection. THE CorrectiOn iS apPliEd aFTEr Sky subtractioN (SEe thE nExt sub-sEctIon), to avOid The Sky ResIdUals being AmplifieD bY tHe IlLumInatiON correctIoN prOcEduRe.

BacKGround SubtrActiOn AnD De-fRinging
--------------------------------------

tHe IMAges CoNtAin sKy bAcKgrouNd flUX coMposed oF direct anD scATterEd AtMospherIc emission, and ScAttered ligHt Of bRight cELEstial soUrces. A careful removal of tHE atmospHerIc bacKgroUnd light iS esSentiaL whEN studyIng LSB ObjecTs SucH AS UDGs. wHEn ConSiDering a sinGLE imAge at LoW lumInosity Levels, we cannot direCTly Tell which part Of tHe liGHT iS diFFuSE liGhT ComING from the sources And which is BaCKgRound light. tO byPaSs the prOblem, we Can maKE an assuMption thaT the patteRn Of baCKGroUnd light stAys constAnt if the tELescoPE pOintiNg dIrectiOn Is nOt chaNged by MOre Than a Few degReEs. Due tO the lArGe ditherS between the consecutive iNtegraTions The Objects arE noT LikEly to appeAr twIce in the saMe pIxeL, whicH alLOws us To prODuCe bACkgroUnd mODels by aveRAgIng A STaCk of images. THE INteNsity Of tHE sky chAngeS thorough the night (1–10% BEtween exposureS), espECIalLy aT The bEgInning and at the End Of THE night, whIcH forces to scAle the imAgES used For the BackgrOund modEL BeFOre
 focal plane by dividing t he imagesby amas ter  f lat- fiel d image. We ad o pted  the flat-fielding met hod t ha t  was us ed by  the Ki l oD e gre eSu rve y( Ki DS, [ @de Jong201 5]): the m ast er  flat-fieldi sachieved b y f irst median- com bining  a ndn ormal izi ng 8dome f l at-fie lds and 8  t w ilight flat-fi e l ds , an d then multiplyin g  t h e averaged fla t-fiel ds wi t h  ea chother. Hig h- frequ e ncy spa t ia l F our i er modes arecorrected u s ing  the d om e f l ats an d low  f r equ encies usin g th e twiligh t flat s . Thisi s based  on th e p re- assu m pt io n t ha t  th e  l arg e -sc ale illu mi na tionof t h e t wili ght  fla t-fie ld matches be tte r th e  ob serva tiona l si tu ation  thanthatob tained from the  dom e. On the  ot he r h an d, in the do mefla ts the[*S/N*] { } i sh i g h,  which can be used  t o ca pture th e pixe l -b y- p ixel var ia tio ns i n the p ixel se nsitivit y.

Ho w ev er , evenaf ter ap pl yin g t he fl a t-fi eld co rrection  smal l  systematic fl u x variationsr em a i na cros s t he instrume nt.T hese  var i at ion s  canbe co rr e ct e d by applying an il lu minati on co rrection. Weuse the co r r e ction mo dels ma d e for the KiDS  (see  [@Verdoes - Kleijn20 13] f or detai ls). Them o dels are  ma debymap p i ng  the photomet r i c re si duals a cro ss theCCD  ar ray  us in g a set o f dither ed  L an do ltSelec t ed Area(S A)fi eld  [@La n dolt19 92] o bser va ti o ns,  and fi t ti n g  a l in ea r mo del  t o the  res i dua ls. The  images a rem ulti pl ie d withthis illumina ti on correct io n.The co r r ection i s applied after sky sub t raction  (s ee th e ne xt sub-se cti on), t o a v oid th e skyresid ua lsb e ing a m p li fie dby the ill u m ina tionco rrec tion pr ocedure.

Backgrou n d s ubtraction an d d e-fr i n gi ng- -- - --- -- - --- - - --------------- --------
Th e  i mages cont a insk y backg round f lux c o mposedof direct  and scat te reda t mos pheric emi ssion, a nd scatte r ed li g ht  of b rig ht cel es tia l sou rces.A  ca reful  remov al  of th e atm os pheric b ackground light is esse ntialwhenstu dying LSB  ob j ect s such as  UDG s. When co nsi der ing a  si n gle i mage at  lo w  lumi nosi t y levels, we  ca n n ot  directly t e l l  wh ich p art of the  lig ht is diffuse lig h t coming fromthes o urc esa nd w hi ch is backgrou ndli g h t. To by pa ss the prob lem, weca n  make  an as sumpti on that t he patter n of  ba ckgroundlig ht stays c on st a nt ifthete lescop e poin t ingd i rection is not c hange d by mo r e t han a  f ew degr e es.Due to the  large dith ers be twee n the  consec ut ive in teg ra tions theo bjects ar e not  likely  t o ap pea r twic e in t he sa me p ix el,  which al l o ws us  t o  pr oduc e bac kg roun d modelsb y averag ing a stack  o f i m a ges. T h ei n tensity of  th e sky c hanges tho r ough th e  nigh t (1–1 0% bet ween ex p osu re s), esp eci a l ly at the  beginnin g  an dat t he end o ftheni ght ,which forces to  scale t he image s used f ort he b a ckgrou nd  model befor e
 focal_plane by_dividing the images by_a master_flat-field_image. We_adopted_the flat-fielding method_that was used_by the Kilo Degree_Survey (KiDS, [@deJong2015]):_the_master flat-field is achieved by first median-combining and normalizing 8 dome flat-fields and 8_twilight_flat-fields, and_then_multiplying_the averaged flat-fields with each_other. High-frequency spatial Fourier modes_are corrected_using the dome flats and low frequencies using_the_twilight flats. This_is based on the pre-assumption that the large-scale illumination_of the twilight flat-field matches better_the observational situation_than_that_obtained from the dome._On the other hand, in the_dome flats the [*S/N*]{} is high,_which can be used to capture the_pixel-by-pixel variations in the pixel sensitivity.

However,_even after applying the flat-field_correction small_systematic flux variations remain across_the instrument. These_variations can_be corrected by_applying an illumination correction. We use_the correction models_made for the KiDS (see [@Verdoes-Kleijn2013]_for_details). The models_are_made_by mapping_the photometric residuals_across_the CCD_array_using a set of dithered Landolt_Selected_Area (SA) field [@Landolt1992] observations, and fitting_a linear model to_the_residuals. The images are_multiplied with this illumination correction._The correction is applied after sky_subtraction (see_the next_sub-section), to avoid the sky residuals being amplified by the illumination_correction procedure.

Background subtraction and de-fringing
--------------------------------------

The images_contain sky background flux_composed of_direct_and scattered atmospheric_emission,_and scattered_light of bright celestial sources. A careful_removal of_the atmospheric background light is essential_when studying LSB objects_such_as UDGs. When considering a single_image at low luminosity levels, we_cannot directly tell which part_of_the_light is diffuse light coming_from the sources and which is_background light. To_bypass the problem, we can make an_assumption_that the pattern of background light_stays_constant if the telescope pointing direction_is_not_changed by more than a_few degrees. Due to the large_dithers between the consecutive integrations the objects are not_likely to appear_twice in the same pixel,_which_allows_us to produce background models by averaging a stack of_images. The_intensity of the_sky changes thorough the night (1–10% between exposures), especially at_the beginning and at the end of_the night, which forces to scale the images used for the_background model before