diff --git a/scikits/image/data/multi.fits b/scikits/image/data/multi.fits new file mode 100644 index 00000000..daa20229 --- /dev/null +++ b/scikits/image/data/multi.fits @@ -0,0 +1,988 @@ +SIMPLE = T / Fits standard BITPIX = 16 / Bits per pixel NAXIS = 0 / Number of axes EXTEND = T / File may contain extensions ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2001-09-10' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:33:03' / Time of last modification COMMENT FITS (Flexible Image Transport System) format defined in Astronomy andCOMMENT Astrophysics Supplement Series v44/p363, v44/p371, v73/p359, v73/p365.COMMENT Contact the NASA Science Office of Standards and Technology for the COMMENT FITS Definition document #100 and other FITS information. INSTRUME= 'GMOS-N ' / Instrument used to acquire data OBJECT = 'NGC1068 ' / Object Name OBSTYPE = 'SCI ' / Observation type GEMPRGID= 'GN-ENG20010910' / Gemini programme ID OBSID = 'GN-ENG20010910-15' / Observation ID / Data label DATALAB = 'GN-ENG20010910-15-001' / DHS data label OBSERVER= 'Jorgensen' / Observer OBSERVAT= 'Gemini-North' / Name of telescope (Gemini-North|Gemini-South) TELESCOP= 'Gemini-North' / Gemini-North PARALLAX= 0. / Parallax of Target RADVEL = 0. / Heliocentric Radial Velocity EPOCH = 2000. / Epoch for Target coordinates EQUINOX = 2000. / Equinox of coordinate system TRKEQUIN= 2000. / Tracking equinox SSA = 'Walther ' / SSA RA = 40.66966667 / Right Ascension DEC = -0.01316667 / Declination of Target ELEVATIO= 0. / Current Elevation AZIMUTH = -135.173831944444 / Current Azimuth CRPA = 0. / Current Cass Rotator Position Angle HA = '+01:14:25.87' / Telescope hour angle LT = '05:00:13.7' / Local time at start of observation TRKFRAME= 'FK5 ' / Tracking co-ordinate DECTRACK= 0. / Differential tracking rate Dec TRKEPOCH= 0. / Differential tracking reference epoch RATRACK = 0. / Differential tracking rate RA FRAME = 'FK5 ' / Target coordinate system PMDEC = 0. / Proper Motion in Declination PMRA = 0. / Proper Motion in RA WAVELENG= 4750. / Effective Target Wavelength RAWIQ = '80-percentile' / Raw Image Quality RAWCC = '50-percentile' / Raw Cloud Cover RAWWV = 'Any ' / Raw Water Vapour/Transparency RAWBG = '80-percentile' / Raw Background RAWPIREQ= 'NO ' / PI Requirements Met RAWGEMQA= 'USABLE ' / Gemini Quality Assessment CGUIDMOD= 'Basic ' / Carousel Mode UT = '15:00:14.2' / UT at observation start CONID = 'GMOS v8.0' / Detector controller ID DETECTOR= 'GMOS + Red1' / Detector name DEWAR = 'Cryostat A' / Dewar name DEW-TEMP= -121.220657276995 / Dewar temperature (Celsius) DETSIZE = '[1:6144,1:4644]' / Detector size NCCDS = 3 / Number of CCD chips NAMPS = 1 / Number of amplifiers SHUTTER = 'OPEN ' / Shutter state during observation AMPINTEG= 5000 / Amplifier integration time TIMESYS = 'UT1 ' / Time system used DATE-OBS= '2001-09-10' / UT Date of observation (YYYY-MM-DD) TIME-OBS= '15:01:11.5' / Time of observation GMOSTHT = 'GOOD ' / Shutter health NSUBEXP = 1 / Number of sub exposures UTSTART = '15:01:11.5' / UT at observation start UTEND = '15:02:11.5' / UT at observation end EXPTIME = 60.0029079914093 / Exposure time in seconds ELAPSED = 60.0029079914093 / Elapsed observation time in seconds DARKTIME= 96.9185329675674 / Dark current integration in seconds INPORT = 1 / Number of ISS port where GMOS was located GMOSCC = 'GMOS-N (V2-8)' / GMOS components controller s/w MASKID = -99 / Mask/IFU barcode MASKNAME= 'IFU ' / Mask name MASKTYP = 0 / Mask/IFU type (0=none/-1=IFU/1=mask) MASKLOC = 1 / Mask/IFU location (-1=unknown/0=FP/1=cassette) FILTER1 = 'open1-6 ' / Filter 1 name FILTID1 = 20000018 / Filter 1 barcode FILTER2 = 'g_G0301 ' / Filter 2 name FILTID2 = 20000018 / Filter 2 barcode GRATING = 'MIRROR ' / Grating name GRATID = 30000008 / Grating barcode GRWLEN = 0. / Grating central wavelength (nm) GRORDER = 0 / Grating order GRTILT = 65. / Grating tilt angle (degrees) DTAX = 103.2 / Detector translation X position (microns) DTAY = 56.64 / Detector translation Y position (microns) DTAZ = 2842. / Detector translation Z position (microns) DTAZST = 0. / Focus at observation start (microns) DTAZEN = 0. / Focus at observation end (microns) DTAZME = 0. / Mean focus during observation (microns) DTMODE = 'FOLLOW-XY' / Detector translation stage mode ADCMODE = ' ' / ADC mode GMOSDC = 'GMOS DC ' / GMOS detector controller s/w DETTYPE = 'SDSU II CCD' / Detector array type DETID = 'EEV9273-16-03EEV9273-20-04EEV9273-20-03' / Chip IDs EXPOSURE= 60. / Requested exposure time in seconds ADCUSED = 1 / ADC used? (1=yes/0=no) DETNROI = 1 / No. regions of interest ADCENPST= 0. / Start entrance prism angle of ADC ADCENPEN= 0. / End entrance prism angle of ADC ADCENPME= 0. / Mean entrance prism angle of ADC ADCEXPST= 0. / Start exit prism angle of ADC ADCEXPEN= 0. / End exit prism angle of ADC ADCEXPME= 0. / Mean exit prism angle of ADC ADCWLEN1= 0. / Lower wavelength for ADC calculation ADCWLEN2= 0. / Upper wavelength for ADC calculation DETRO1X = 1 / ROI 1 X start DETRO1XS= 3072 / ROI 1 X size DETRO1Y = 1 / ROI 1 Y start DETRO1YS= 2304 / ROI 1 Y size AIRMASS = 0. / Airmass averaged throughout the exposure AMSTART = 0. / Airmass at start of exposure AMEND = 0. / Airmass at end of exposure NSCIEXT = 3 / Number of science extensions PIXSCALE= 0.1454 / Pixel scale in Y in arcsec/pixel NEXTEND = 3 / Number of extensions OBSMODE = 'IMAGE ' / Observing mode (IMAGE|IFU|MOS|LONGSLIT) GEM-TLM = '2010-07-02T20:32:41' / UT Last modification with GEMINI GPREPARE= '2010-07-02T20:32:41' / UT Time stamp for GPREPARE END XTENSION= 'IMAGE ' / Image extension BITPIX = 16 / Bits per pixel NAXIS = 2 / Number of axes NAXIS1 = 132 / Axis length NAXIS2 = 288 / Axis length PCOUNT = 0 / No 'random' parameters GCOUNT = 1 / Only one group BSCALE = 1.000000E0 / REAL = TAPE*BSCALE + BZERO BZERO = 3.276800E4 / EXTNAME = 'SCI ' / Extension name EXTVER = 1 / Extension version INHERIT = F / Inherits global header ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2010-07-02T19:55:22' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:59:26' / Time of last modification CTYPE1 = 'RA---TAN' / R.A. in tangent plane projection CRPIX1 = 196.375233233945 / Ref pix of axis 1 CRVAL1 = 306.871725219222 / RA at Ref pix in decimal degrees CTYPE2 = 'DEC--TAN' / DEC. in tangent plane projection CRPIX2 = 143.509119025456 / Ref pix of axis 2 CRVAL2 = 37.3733845984976 / DEC at Ref pix in decimal degrees CD1_1 = -2.9266657155837E-4 / WCS matrix element 1 1 CD1_2 = 1.36727521936194E-4 / WCS matrix element 1 2 CD2_1 = 1.36706238220991E-4 / WCS matrix element 2 1 CD2_2 = 2.93036096725254E-4 / WCS matrix element 2 2 RADECSYS= 'FK5 ' / R.A./DEC. coordinate system reference EQUINOX = 2000. / Equinox of coordinate system MJD-OBS = 52162.6261539214 / MJD of start of obseration CCDNAME = 'EEV 9273-16-03' / CCD chip name(s) AMPNAME = 'EEV 9273-16-03, right' / Amplifier name(s) GAIN = 2.1 / Amplifier gain RDNOISE = 3.5 / Readout noise CCDSIZE = '[1:2048,1:4644]' / CCD chip size CCDSUM = '2 2 ' / CCD sum CCDSEC = '[1:2048,1:4608]' / CCD section(s) DETSEC = '[1:2048,1:4608]' / Detector section(s) DATASEC = '[33:1056,1:2304]' / Data section(s) BIASSEC = '[1:32,1:2304]' / Bias section(s) FRMNAME = 'dhsFrame' / Name assigned to a frame of data. FRAMEID = '0 ' / Frame Id DATATYP = ' ' / Type of data frame WCSDIM = 2 LTV1 = 0.4375 LTV2 = 0.4375 LTM1_1 = 0.125 LTM2_2 = 0.125 WAT0_001= 'system=image' WAT1_001= 'wtype=tan axtype=ra' WAT2_001= 'wtype=tan axtype=dec' END <r7^Ņ߅ .8A; #],/kʅ<]yÆņȆˆΆ҆Ն؆چ߆ + !#%())+.0149<=<;=AA: ?b\Ʌ>^z†ĆƆɆ̆Ά҆Ն؆ۆކ ##!#%&)+,+.12249>?@@ADE?#_Ä*><]zĆȆˆ͆φцԆ؆ۆކ + !('&'()+./0256559?CDEEGIF +1~Nd:Zw†ƆˆΆІцӆֆچކ + $))(*+,.02479;;:?BDGIIKMH F!O)VwÆņƆʆΆ҆Ԇֆچ݆ #&()**-1578::;<=@BDEGHJMNH N8 mĆ Gr†ņdžɆ̆Іӆ׆ۆކ + "$')+++-158:;===>BDEFIKMPPJ T˄N')`Ćdžʆ͆Іӆֆۆ߆ + "$'*,.//1369;=@BBCEFGILPSSSL Zلd؅@@tÆdžˆ͆φӆֆچ߆  "$'+.0256899;=@DEFHIJKMPTVUSN ^vU QņˆΆφцՆ؆݆ !$'+.0137:<<<;=ADFGJLMORTVWUTN  bfdž]Ædž̆φІӆ؆܆ +  !!#&),.02469;<;;?BDFHKMPTVXYYXVO f sӆ%eÆƆɆ̆φ҆׆ۆ !$%&')*,.02468;<;=BCEHKLPTVXZ\\]ZP + + + i|چ,iņɆˆ͆Іֆچކ !#&()*++,.0258;=>?ADEGKNPTWWY\]_`]S + + + +!kކ/jÆdž̆ΆІՆچ܆߆ !#%(*+,,-/1259;>@ABDFGJMQTWXY[]_aa_V + + + + +!k/kÆņʆφцՆچ݆ކ  "$'*-.--/2469<>@@@CGHJMPTWXY[]^`ab`X + + + + !k 1l†ĆƆˆцԆ؆܆ކ !"#&)+./.0258;=>ABABEILNQTWXY[]^^`bcbZ + + + !k!3n†ņdžˆцֆنۆކ !&&&()*.0/257;?@ADEEGJNPRUWXXZ]___acee] + + + + + "k"4oĆɆ͆҆ֆنۆ߆ + $*,++*,.0048:=@BEHIILPSTUWYXX\_``acegh_ + + + + + "k"4oņʆφԆֆ؆ۆކ  &,..--./026:=@BDGJKMOSXYXYZYZ]`aabdgii_ + + + + #k"4qÆdž̆цԆ׆ن܆߆  %*.00000036:?DGIJLOPRV\][Z[[\_`bccehig_ + + + + #k#6t†ƆɆ͆цՆن܆ކ +!$).122102359AHLMMPTUVY]]\\\]^`acddfhhfa + #j#8uņɆˆΆ҆ֆچ݆߆ "%'*.12101346;CIPQPRVXZ\^^]]^_`bdeeeghhgd + #k#7sƆʆ̆φӆֆن܆߆ "!!"$'*,-.0212446:?DIPRRSWY\^_____`acfiigghijf + #k$5qƆɆ̆ІԆ׆ن܆߆ "#$&'*/10/145667:>AEIORTVWZ^`abba`abeillkklmmh + #k&7tņȆ̆цԆ׆چ݆ #%&(+0210269999=@CFKOTXZZ]`bdfedcceilmmmnpqok + #k';z†ĆȆ͆Іӆ׆چކ #%'),/00137:;<=?BEIMRX^aba`aeikjgfilmmlmpsrql + #j'?}ÆƆʆφ҆Ն؆ۆ߆  #&)+...03569<>@ACGKQX`enpg`_cjonjilmlllmqssrn + $j'@}Ædž̆цԆ׆چ݆!&*-///24579yĆɆΆ҆׆ۆކ #(,/13447;=?BEHKPUYZ\\]]\]_bfjmopomnqstuvyyr + %m*=yņʆΆ҆؆܆ #(,/1346:=?BEILOSWXY[\\[[]`bfilsxvqqstuxz||u + %m +=zņɆ͆ӆ׆ۆ !%)+.1579<>ADIMPRTVXZ[\\]]_acdfjs{ytqruw{~y + %m ,>|ņʆφԆ׆ۆ + $')*-26:<>@DILPSTUWY[\\]__`cccejqvvsrtx{~~ + &m .@~Ædž̆цԆ؆݆ #&()+/48;>@DHLORUVXYZ[\^_abccbchlpstrtw{} &m /AņʆφӆՆچ߆ !%()+.25:=@DHLNPSVXYYY[\_abddcbejnprstw{} &m 0Cdž̆цֆن݆  #&(*-147;?DILNPSTWYYYZ[\`ceggfgimopqsv{~ &m 0EȆ̆ц؆݆߆ !$''(),/36:=BINOORUVXYZ\]\]bfilmmnopqrsvz~ &m 0EɆΆӆنކ !'+--,,/369BDGILOQSVWXZ]acdegjlmoruvwz~ + 'l1KƆφԆ׆ۆ  ')(*.114;BFKNPQSVWXYZ]`dggikmnopsttv| + &k1KƆІԆ׆ۆ !""&+/26EGIJMSWWUTV[_`_`begilmloruwx{} + &j3OʆԆچ #15/,.26:AFHGIMSVTTTX]a`^_beghjlnqsuwz|~ + &j3N͆׆݆ =moG1159@DFHKMORTVY]_^`cfggimpt{ + + + + %j> _Ɔ؆ +#'*.0368:<>ACEFILORTUY\^_behjklorv~ + + + + %i@ aȆن $(,/147:;@ABCCEINPRUWZ]`bceiosvyz|~ + + + + + + + $i CfɆ܆ #(,/268:=ACCBBBBFLOOQUX[`cedcipvy|}} + + + + + + + + $h Egˆ݆  $(*.367:>CGHFCBDIMOPRUY^adghgkpuz} + + + + + + + $i Fgˆކ "%')-2469CGIJKMPU]a]Z]_aeilpvȈ=ÇȇƂ + + + + + + + %i …Kmφ +%*/24556:=@CEHJKMPTWZYZ^abdhjnv+ԏˋUVLJˇǂ + + + + + + %i ÅLm҆ +  '.465568<@CEFHKLNQTVY[\_bdehkow݉k:6ƇLJ‚ + + + + + + + + + + %i ÅLnԆ  $*398668ADFINSZ`]YWWXX[^afjmpqrtw|ÇȇӇ݇ۇ҇Ȃ + + + + + + + + + + + + + + + + $h ąN!vۆ !&*-257:?CEHJMPUXXVUUVY[_cgknqrtwz‡Ƈ̇ч҇ЇƂ + + + + + + + + + + + + + $i ŅO!wކ "%(,1358=BFJKKMQTVUUUWZ]aehknprvy|ÇLJʇχ҇Їł + + + + + + + + + + + + + $i ŅO!w߆  $&)-1237CGIIIKMQUWWXZ\`efhlnprv{ÇLJʇ̇ЇӇӇ˂ + + + + + + + + + + + + + %i ƅQ$y߆ !&)*+/37<@DGHHIKMQTVX[]^begkoqrux}ÇƇɇ̇χӇԇЂ + + + + + + + %i DžS†&zކ "&(*,/39=@CFGHJLNPRTX\^`cehlprtx{‡ćƇʇ·чՇׇӂ + + + %i ȅTĆ'{ "%'+.15;>@CEFHKNPQRTWY\_beilosvy}ćƇɇ·҇ԇևهԂ + + + %i ɅTĆ(| !"%+036;>ACEFGKOPQSTTVZ^bgklouyz}ÇŇLJ̇҇ՇՇׇهՂ + %h ȅTĆ(}  #&+036:>BCDEFJNPQSTTVZ^cimnqw{|~‡ćŇɇχՇ؇؇ڇ܇ق + %h DžTņ(~ !%)-147:>ACDFHJMOQRSUX\`djnqvz}ÇŇȇ̇҇؇ۇ݇߇܂ + %h DžUƆ)~ "',/269;>@BEIKLNPQRTX[_bfjoty|ćLJ̇чՇ؇݇ + &h ȅUņ)~ + $)-036:=?@CFILNORTTVZ]`cglqw{~ćɇχԇׇهއ &i ȅTņ)} +  &*-147:=@DFGJMNPRUWZ]_`chmsxz~ćɇЇՇ؇܇ &i DžTņ*~  &*.147:=BGGHLOOPRTW[]^`dhnswx|ÇƇʇЇԇه߇ &h DžTņ*  %+/136:>DGHJOPPQRSVYY\aeimrvx{Ňȇ̇ЇՇ܇ &h ƅTƆ*  %+/026;?DFHLOPQRSTVWY]bfimrvxz}Ňʇχׇ҇އ  %g +ƅTƆ+ !'+./38@CDGKNRWZ]`bejnqsuxz{}ćɇ͇Їԇهއ  %e…QÆ) "&,5:88:=@BDFJLOUZ\^acfjmprtuw{~‡Ƈʇ̇Їԇه߇   %eąSņ+ +"%)0579;>BDFJLNQW[^`cehiknrsrtz‡LJɇ͇҇Շڇ  %fȅXɆ.  "$).37:<@BDHLNORVZ_bdehhjnqstv{ÇLJɇ·Ӈ؇އ  &g ΅aц2† !#&*.159<@CEINPRTVY^bdfhhjmoswy|ŇȇˇЇՇۇ   &g Յkۆ:Ɔ  $(+-0368=BEINPRTUY]`cgihikotwy}ćȇ·ԇ؇݇   &g܅vFʆ + %(*-1367;AEIMOQSUX[_cfhijkpuwy}‡ȇчׇڇއ   &hㅁP͆ !$'*.258:>DGILNORUWZ_begiklquwy|‡ʇчևهއ   'i녌XІ  #&*/259=BGJKMNORUX\`cegjknruwy|~ÇʇχӇ؇އ +  'j󅘆`ӆ +"&,0149?DHKLMNPSWZ]`cfhikosuwy{ćLJˇ҇ڇ  (ki؆ "(-136;AFJLNOPRUZafedfgilosuy{}ÇŇƇˇՇ݇  (ks܆ #).26:>CGJMOPRTYglfgjmoruy|ćŇȇχ؇߇   (i%{ #).49<@DGJLMNPT]ӈ{hfilmqtw{‡Ňʇ҇ه߇   &d-  %*069<@CEIJJLOS]/ׇihjlnrux}Ƈ̇҇ه #\셲4Æ !&-47:=@ADGHILPT[rˇwiiklnsw{LJ͇҇؇ !Sۅ9Ȇ  &/7:<>@ACEGINRVZ`mvqjgijlpty}‡ɇχӇڇ  LڄŅ=͆ !(19<=?@ACEGKOTXZ[^`bcehknrw|~ćˇчׇ܇ !Gʄ@҆ + #)18;>??@DFHLPTXYZ]_`adhmqux{}}ҡ҇ه߇  #HÄCֆ !%*069<=?ADFILOSUWZ^abbehntwyz{}Ƈ͇ԇ܇  &L„F؆53!&,1478;>@BFJKMPSVZ^acdfjotwyz{~ȇЇׇ݇ )QDŽG؆/jd."'-3777:=?AEIIKNRW[]_aehlpsvxz|‡ʇЇׇއ  +V΄G؆ *\X-!$(.4877:=?AEGHKOSX[\^aeimpruxz}ćˇЇև݇ +  -ZքD׆')"$).1468;=@BDEIMPTXY[_bfjmoqtvz}LJׇ̇҇݇   /^܄Bֆ"%)-/269<>@BCFJMPTWY\adgjkmpruy~ćɇχӇׇ܇  0b…@Ն"&*-036:=?@ACGKNQUXZ^acfijlpruy}LJˇЇԇ؇݇ !! 0fх<ӆ  #&*.258<>@AADGKNQVZ\_acfhimpsvy}ŇȇˇЇՇڇއ && /k六9φ "$&+/36:>?@ACEHJMQVZ]_bdfgjmqtxz~Ňȇˇчׇ܇އ߇"+) .o6Ɇ  #$'+/37:=>?AEGIKLOUY\_bdehknqux{͇ЇƇ‡Ƈˇчׇ܇އ߇&,$ -q!0ņ "%(*-269<>?ACFIKMPUX[^cefhlorux|هۇɇćɇχև݇ %( +o") !%(),058<@@ABCFJNRVX[_bdehlnqtw|~ʇ̇Çćȇ·ԇۇ!"  )l"y !&(*,/48>CECABEINSWZ^accegjmqux{}‡Ňȇ̇Їׇ    'ks  &)*,/38@JLFBCGJNRV[bhheefhkoswz}ćƇɇ͇Շއ +  &jk݆  &)*-047>JMEBEHJLQV\dllfcdgjlpuz~ÇŇLJ̇ևއ + &jb؆ + $(+.1359ADAAFHILQV\affb`cfghmtyÇŇLJ͇Շ݇ &jYԆ + !&+.2457:=>@DGJMQVZ^__^^beehmsy~ÇƇʇЇև܇  'kꅅPφ !&*/3457:?BFIMPQRUX[_abcglsz}ƇˇЇԇه܇݇ + &jӅdՆ:Ɇ #&),/479;<=>ADGKMPSVXZ]`behlrx|ȇ·҇Շ؇ڇ܇  + &i҅bԆ8Ȇ "%(*.2479:<=?BEHKOSVWY[^cfimsy}ˇՇۇ܇ۇڇ݇ + + &i҅bԆ7Ɇ "&(*-03568:=?AEHLORTWYZ]afkpuz}ɇׇ߇܇  + + + + &j҅bՆ9ˆ $')+-/2446:>?BEIMPQSWZ\^agnswz|~ćЇ܇އ +  + + + + + + + + + + + + + &jхbֆ;̆ !'))+-/0136:=?ADHMPQSXZ\_dipuxyz}ɇ҇ه܇݇  + + + + + + + + + + + + + &jхaԆ;͆ + (.,*,-...159;=?@FLQQTWZ\`ejouwy{~Ƈ̇҇ևه߇ + + + + + + + + + + + + + + &jЅ`ӆ:͆ "-40++,,+,059:<=?EKOQSW[_bfjotvx|Ňʇ·҇և߇ + + + + + + + + + + + + + + + + + &iυ_҆9̆ %-1.*)*++,059:;=AFJMQTX^cfhjnsux{‡Ňʇ·чՇ߇ + + + + + + + + + + + + + + + + + + + + + &iυ^҆8ʆ%(('&')+,-0478:>CGJNQUZ`dgiknqtvzćɇχчՇ߇ + + + + + + + + + + + + + + + + + %iυ_ц6Ȇ!%$"!"%(+-/0357:?CGKNRW]acfhjmpruy}ɇ·чՇ݇ + + + + + + + + + + + + %hυ_ц5Ȇ! #&)-/0369=@BEIOTY^`cfghknquw{LJ̇χԇ܇ + + + + + + + + + + + %g υ_І4Ȇ"$(+-037;>?@BHNTY]`cefgjlptvzŇʇ·Ӈۇ + + + + + + + + + $g ΅_І4dž  "'*,/37;>>?BHMSX\_befgiknruyćȇ͇Ӈۇ  + + + + + + + + +  $g ΅^φ4dž  #&),/37;=>@DHMRVY]adgiijlptx}Çȇ͇ԇ݇  + + + + + + +  $g ͅ]φ4ņ + !$&(,047;<>BEIMRVWZ_dgiijlnrv{‡Ƈׇ̇ + + + + + + + + +  $f ͅ]І5ņ !!#%(-048;=@CEHMQUWY^befhjkmqvzćʇև + + + + + + $f ͅ^ц6ņ !#%)-/27@ADGINRTUV[^_`beimqvz~؈/ + + + + + $g ͅT*† +$)-37:<>@BDEGLQTTVY\^adgilquzԈ   + + + $g ʅO%  &+047:<>@ACEGJNQTUX[^bfgglqty}LJۇއه܇އ߇$ + + + + $g ʅQ( + !',0358;=?@AEHJKPTVWZ]adegkosw{}ȇ͇чևڇ܇އ8: + + + + $g ̅U†, "(,.0258;=>AEIKMQSUWY]`cehlorvy{}ćȇ͇Շڇ܇އ?9 + + + + $g ΅Ydž- + #(+,-0247:=BFIJMOQTWY]_aeiloqtvxzzyÇƇ̇ԇڇۇ܇އ % + + + + + $g҅]Ɇ- #'((+./148CHJMPORXZ[]_adgjkkkmrux{}~LJЇԇԇԇׇ؇ڇ܇߇߂ + $fƅS%zކ !$&()*,159>DHIOSQQUWY\^`cehjkmptxz|}Ƈ͇ЇЇ҇ևׇه܇߇܂ + + $fŅS†&{ކ  &()*+.26;@DFHLONPRV[]^abcegjnswz{|}Ňʇ͇̇ЇӇԇׇڇއق + + + $eŅS†&{߆ + %()*,/38?CCDFHKNPS[ggbaa`beiouxz{{|‡LJˇ̇·Ї҇Շ؇܇ق + + + + $eŅS†%x݆  $'(),/39@BABDGKOQTe|yg`__`chpuxyzz}‡ȇ͇̇·чԇׇڇއ߇ڂ + + + + $eŅT†$x܆ +  %)''+/39=>?ACGKMNScywf_^^_ciosvwxzŇʇ̇·чԇ؇ۇ݇߇ۂ + + + $eŅT†$yކ !),(',05<=<=@CGIIKPYcd_\\]_dkoruwwz‡Ƈʇ͇ЇՇه܇݇އڂ + + + $fŅS${ +  )+''+.3:<::@BEFHJLMQTVX^beinprw}‡ȇˇˇʇˇ·ч˂ + $g ɅXdž)~߆ +!#$'+.0357;=>@CEFHILQSTX]behlosx}ŇLJȇʇ͇·χȂ + #g ʅXƆ(} !$)-/0247:=>@BCEEGKOQTY^cgilpuz}‡ćLJ̇χЇχɂ + #f ʅYƆ(~  %*-/1258;<>@ABCEHKNQUZ^dgjmrw{ÇLJˇ·Їч˂ + "e ̅Ydž)߆ + !&),/1369;<>??@BEHJMRVY]beiotw{~ÇLJʇ̇χч̂ + "e ͅ[Ȇ*݆ "&),/1479;<==>@BEHJNSVX[_ciptw{~ÇƇɇ͇χʂ + "e ΅\Ɇ*~ۆ #&),/1589;==<>@BDGKPUVWY]bjpsv|ҡ·Ȃ + "f ΅]ʆ*}ن #&(+.1467;@@??ABDGLRWXYZ]binrv{桷Ȃ + !e ΅\Ɇ*}؆  $'),0235;AC@?BDFINSWY[\^afkptx|ćɇˇł + !e υ\Ȇ){؆"%'*-.16><=ADIPUX[\`aa`bdhkquvwxyz|‡‡ +  e ΅]ʆ*|׆ +!$')-48:;;=ADGMSVY[^``abfilrvxxxwx{~ +  eυ]Ɇ){׆ "%(,14689BCDGKPRQQSUX[`cfilnpruy{zz|}~ +  fЅ]Ȇ'wӆ +  #&)-39>ACEHKPSRPQRV[_bfijloqswxxz||}~ +  fЅ[Ɔ%w҆ "&*.4:>ABDEGLPPOPRVZ]aeggilnqstvxz{||~ + + fЅ[ņ$vцކ "&*/5;<>@ABDHKMNPTY[]adddfilopqsvxy{z| + eЅ\Ɔ$uІކ  #&)-4889;>@BEIKLPY^\\_aabdgjnopruwyyyz} + + fЅ\Ɔ$tφ݆ !#%'+03568BFGFJPTUVX[_bdgiklmopruvxz{} + + +  gυZĆ"r̆ۆ߆ "#$')+.1347:=ADEEHKORSUY^beijiikmnpruy{|}} + + + +  gυZĆ!rˆنކ !##$'(*,/1248ACDEGLPSW]`begggilnprstxz} + + + + + + fυYÆ odžԆچކ + #$&&&()+-/16:=?@ACGJMRY]^`bbbeimprsttux{}~ + + + + + + f΅Xmņцֆۆ #''&')+-/048;=>?CFGJQWYZ^`_`cglprsttuwz|} + + + + + + f΅XlÆІֆۆ߆  $%&')+-./369;<=ADEIOSTX^a``beimprsuwxyz} + + + + + + fͅWk†Іֆچ݆ +!$&)+,-/3689:ADFJOWbiga^`dgjlnqvxxxz~ + + + + + e̅UjΆӆՆن߆  "&)*+,05989;>@CEIP]msi_]`ceiloruvvxz + + + + + e ʅSj†φӆֆچކ + $(***.3789;>@CFIQ_rwi__abcgkostuvx{ + + + + + e ɅQg†ІՆ؆چކ + "'***-2579;=AEIMQZgkd`bbbbfinrsuwz| + + + + + d ɅReφՆ׆چކ !$(*,034679?CGLRUX`ktupllmmnsxyyz|}~~~ + + + + +eɅPdʆ҆؆ކ #).23348::;=@CGLQWdwunmlknsxyzz{|}~ +eɅOdɆц׆܆ "(,/2224679<>ADHOWgxnkjknrvyyz{||~ +e ȅPcȆІՆن݆ + + $'+/0/02489CLU_rzjfghjmquwxyyxz|}~ +eɅOaɆφ҆ֆچކ "$&'')-2788:=BKSZalunedfgilortwyxwwy||}~ +f˅Oadž͆цԆ׆ۆ߆  "#%%'+0468:=BHQXYZ_``bdfhjloqtwwwvxz{{|}~ + iЅPbƆ̆φ҆Ն؆܆ + !##%),/49<>ADKRSSUY\`bdgijlnqsuvuvyzz{{{|}y + i҅RbƆʆ͆ІԆ׆ۆކ !#%&(,39=?ABFJKMPUY]`bfghiknpqrsuxyyyxy{zs + hυRaņɆ̆Ά҆Նن܆߆  #$$$*29=?@BFGGIMQVZ^adeggilnmnpsuvwwvxzxp + g˅P_†Ɔʆ̆φ҆ֆچ݆߆ +  "##$)08=?@CGIIJKNSX\`bcdfhkkjloqsuutvyyun  + gɅN]ĆȆʆ͆ІԆ؆ۆ݆  "#%)/7=?ADGJKKKMQV[^`bdfiiiilopqssruxxtl  gDžL\ÆƆɆ̆φ҆ֆن܆ !$&(.5;@CEEGIJKLOTX[^adfhhhikmopppqsuusl f ÅH YĆȆ̆φӆֆ؆ۆ߆ +  $&(-39?BCCEHJKKNRUX[_bdeefgikmmlmprqrrk dAW†ƆʆΆӆֆ׆چ݆ "$&+17<>?BEGJJJLPSUY]`bbcdegjkjiloonooh b4SÆdž̆цԆֆن݆ +  "%*0599;@DFHIIJMPRUY\_`abceijhhkmlllmg `$J†ƆʆΆцԆ׆ۆ߆ + "&+0357:>BDFGHIKNPRTWZ]^_adhhgikkjjjjh ]|:~†ƆȆˆφ҆Նن݆ + #'+.0259?@BEGJLLLNQUXZ^dgeccdfggfea RӄgL +QÆņȆ̆φцԆ؆ކ +!$$&*,.0369;<>?AEHJKLMORVZ_dfda``cdefc\ KÄNŅ05jÆdžˆΆцӆ׆܆ !#"$)+,-0357:<>@CFGIKLMQTX\`a`_]]_acec\ +C4tΆMuƆɆ͆Іӆֆۆކ  "&()*-02479?AABGJKLMLMOQQQTZ^_``Z \Ʉ=f7Yy†džɆ̆ΆцԆ؆چ܆ !$(+-/0368:=>?BFIJJJIJLOQQRV[\\\U < +xޅ:7[{Ædžʆ͆І҆Ն؆ۆކ "%')+.1236:;ADFDCDHNTUTSTUWVN/ze̅V6ZyÆɆ̆Άцֆ؆چ݆  !"%')+.26778;>ABCBCFLTYYUSTVUN +=߅!W6XwĆɆ̆ІԆֆن܆ކ !$&)-1455689;=??@BIS[ZTPQSSK F"V 3WxĆȆ͆ІԆ؆چۆ݆  $',/233466579:FMNLKKLMJ + J#T 3WxÆɆΆ҆Ն׆ۆކ !%''(+-.///024446>FJKJIKNK + J!S 2UsņʆΆ҆Նن݆ + """%')*+,-.121028AGIHHLQP + J S 0StĆʆІԆ׆ۆ߆  "$&'()*,.//.05>FGGHLRV + + J!߅Q߆1Vwdžφӆֆن݆  ""#%'(*,--./2;CFEGKQY + K"ޅPކ +3Ww͆҆Ԇ׆چކ  "$')*+-//16<@CEINS + L#ޅN߆ +2VwņΆ҆Ԇ׆چ܆݆ކ +  "$'()-./025;@CFGE + M#݅O߆ +2WyĆΆцԆ׆؆ن܆ކ  !"$'+,,-/28>@BCFILOQOJ + &`Ȅ7*KiDžхمޅ߅߅߅ XTENSION= 'IMAGE ' / Image extension BITPIX = 16 / Bits per pixel NAXIS = 2 / Number of axes NAXIS1 = 132 / Axis length NAXIS2 = 288 / Axis length PCOUNT = 0 / No 'random' parameters GCOUNT = 1 / Only one group BSCALE = 1.000000E0 / REAL = TAPE*BSCALE + BZERO BZERO = 3.276800E4 / EXTNAME = 'SCI ' / Extension name EXTVER = 2 / Extension version INHERIT = F / Inherits global header ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2010-07-02T19:55:22' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:59:26' / Time of last modification CTYPE1 = 'RA---TAN' / R.A. in tangent plane projection CRPIX1 = 65.7568549170406 / Ref pix of axis 1 CRVAL1 = 306.871725219222 / RA at Ref pix in decimal degrees CTYPE2 = 'DEC--TAN' / DEC. in tangent plane projection CRPIX2 = 143.564898774826 / Ref pix of axis 2 CRVAL2 = 37.3733845984976 / DEC at Ref pix in decimal degrees CD1_1 = -2.9278046959101E-4 / WCS matrix element 1 1 CD1_2 = 1.37000324401069E-4 / WCS matrix element 1 2 CD2_1 = 1.36969765153489E-4 / WCS matrix element 2 1 CD2_2 = 2.93069759666565E-4 / WCS matrix element 2 2 RADECSYS= 'FK5 ' / R.A./DEC. coordinate system reference EQUINOX = 2000. / Equinox of coordinate system MJD-OBS = 52162.6261539214 / MJD of start of obseration CCDNAME = 'EEV 9273-20-04' / CCD chip name(s) AMPNAME = 'EEV 9273-20-04, right' / Amplifier name(s) GAIN = 2.337 / Amplifier gain RDNOISE = 3.3 / Readout noise CCDSIZE = '[1:2048,1:4644]' / CCD chip size CCDSUM = '2 2 ' / CCD sum CCDSEC = '[1:2048,1:4608]' / CCD section(s) DETSEC = '[2049:4096,1:4608]' / Detector section(s) DATASEC = '[33:1056,1:2304]' / Data section(s) BIASSEC = '[1:32,1:2304]' / Bias section(s) FRMNAME = 'dhsFrame' / Name assigned to a frame of data. FRAMEID = '1 ' / Frame Id DATATYP = ' ' / Type of data frame WCSDIM = 2 LTV1 = 0.4375 LTV2 = 0.4375 LTM1_1 = 0.125 LTM2_2 = 0.125 WAT0_001= 'system=image' WAT1_001= 'wtype=tan axtype=ra' WAT2_001= 'wtype=tan axtype=dec' END wybiƅDžDžʅЅ؅߅ׅυŅ}ywxywuuwwutvwurdy{pŅ…………ąƅȅɅ˅ͅ΅ЅӅՅ؅܅݅ۅۅۅۅۅۅممڅ܅݅݅݅ۅՅυʅDžDžȅʅ˅̅ͅυЅЅЅх҅҅Ѕ΅ͅͅͅ˅ʅɅȅDžƅŅąÅÅz|{؅܅߅ + +߅߅݅ۅڅمׅօՅԅӅхυ̅˅ʅʅʅɅDžDžɅʅʅDžz|ȅ  !  ""#%()&#!  +߂z|Ն !#(*(&'('')+13.*)*,---.-,,,,+++*+--.142012.*)+,.24321/..//.-+((*)('&%$" z|܆ +  #$#$$&*,--.-,,-4EL=1-.022343221/---.0457:<:8::633468<==;:86775444101000/.-+)('%%$!  z| !"#%'''''')+.10//./9R[E3/024467776410012479;>?=;<=:89:;>ABA@><<<:87876443443110/.-+*'%$"!z| "%&')*+++*)(*-0//0/07GL>3113568:;:9644567889==>?@CFFDCBAA?=;::987679975543110-*)(&$#!z| +!&)++,,----+)*,../00049:6333468:<<<;:99:<<:99<@CB?=?@@@BDFIIGFEDCA?>=<;:98:>?=;:85453/,,.-)&#!!!!"""!  z|  !$)-...../0/.,---./012355555579:<>??ABA>?@?<;;>ACCAABCBBDGIJJJIGFDCBA@?>=;;<:8862.043.*'$$#$&%$$""##"#$#"""!z| + #$$%',//0122243100/./23345678889:;<>ABEINLFCB@???ACEDDEGGEDFILKKKJIGFFEEDCC@>>>@DFDA?><:8425:71-+)''()''&%&)(''('&%$#"z|  "&())+.012356799763224677779;<<===<=@DGJNSQLHEDDCCDGIHHJKKIIJLNNMMLKJJIHHHHECBAACGHEDCB?<977;>:3/.,++++)))*,-+*+,+)(''&"z| #&)*+,-0245689;8679:;<<<>@@AAA@?@CGJLNQROLJIIGGHJLMMOPPNNNOPQQQONMLKJJKJGFEDEGIJIHGDB?===>>;631/.//-,+,/0/--//.-++,+&z|""$'+-...03689::;>BBB@=:;<>?@@ACDDDDDDCCFILNOQRQPPONLKKLOQRSTUSSRRRTUUTROMLKLMLJHGGHJLOOMHEDCCC@>=<:64223210/133101210/..00*z|!%&(,/023347:<<===>ACCCA@??ABCDFFGHGHIHGGILOPRSTSSUUTQPPQSTTUWXWWUUUVWXXWSPMLMONLKJJLNRYWPJFFHGEB?>?=:7676653356666665431143. z|$(*,0356788:=>???@@ACEFFFDCCDEGIJKKKLLKJKMOQRTUVWWXZYVTTUVVWXYYYYXXXYZ[[YWSPNOPPONNNOQX`]SMJJKJFCAAA?=<<<;::989;;;:9998754550%z|'+-035899:;=@AAABCCCFGIIIFDEGIKLMNNMOPNNOQQRSUWYZZ[\\ZYYYYZ[[\[\\\[[[[[ZYXVSQQRTTRQPQSW\ZTQNNOMIFEDCBA@BB@@@???@@?<;;;:976652&z|*/13579::;=?ACDEEEEFIKJIHGGHJLMNOOOPQSSRSTSTUWY[\]^^^^]\\]]^^^^___^]]\[[ZYWUUVW\]YTSTTUUTSRPPPOMJHHGEDDFGFEEDDDEEB@??=<;:8773(z|.3567:<<==>@BEGHHGGJLMKHHIKLLMNPPPPRSTUVVVUVXZ[]_`aabb`^^_``_``abb`_]]]]][XVXZ]bd_YWWVTSSSRPOOOPNLKKHGHKMLKIIIIIHFDDCB?><<;96,z|3679:=?@@??ADFIJJIJLNOMJILNOONOQRQRSTTWXXXWXZ\]_bccdedb``abbaaabccb_^^_`_][Z[^`bda]\[XTQSSRQPOQRPNNMKKMOPONNNMLLKIHGGEB@??=;80z|" 7::;???@BCCCCDEGHIKOPPQTTRQSTSRSTUVWXXXWWY[\ZZ[^aceghhgfdcdeeeedeedca````bbcddeggffdb`_^[XWVWVVWWVTRSSSTVXWUTTSPPOOOPPNKIGECA@>3z|%>CDDDDDDDDEFGHHJMQRRTVWVUUUTUUUVXYZZZYY\^^\[]_bdegjjihgeefghgfffdcbaaabceggfhjkjjgdca_^^][ZXYZZYXWVWWXYZZXWVTQQQPQSTTROLIEDDB8z|&AEFGGGGHGFGHIJKMORTVXYYYXWWWWWVWY[[\\[]_aa`_`beggilmllkigghiihhhffeccddehjjikmoomjgfcbbbb`^\\\\\[ZZ[[[[\[ZYWVUTSRSU[a`XPKGGGE>z|(BEFHIIJKJIIJKMNOPRVZ\[ZZ[[ZZYYXX[\\]^^_`abbcdfhjkklmnopmkiijjjkkjigffghikmmnprsrokjgeddeeda_____^^^__^^^][ZXWXWVUUYdro_SMJIIICz|*AEGHIJMOMLLKLNPQQSX\^]\]^^]\\[ZZ\^^^`abbcdefhjlnnnmopqrqnlkkjknonkjijjkmopqsuvvsoljhffffgfbabbbabbccbba`_^\ZZZYXWW[hvraUOMKLMFz|,AEGIJMPQPNMMNPQRSVZ^``__``^]^]]]]_``bdeefgikmnoqrrssrrrrponmlnqrqnlllmoqrrtuwxxtqmjhhhhhigdccddeghhfedaa`^]]]\[YXX[bhf\UQOOQQJz|-BFHILNPRRPPPPRSTVY\`cbaba_^^__^^^`bbdffghjlopqqstvxxvttsqpppqstsrpnooprssuvwyzyvsoljihijkjgeefhjlmkigdaa`]\___\YXZ[\]]ZVTSTVVOz|-DGHJLNPRSSSRSUVWY\^beeddb``a`````bcdfhiklmorsstuwx{{xvvtsrqsvxwusrqrrsttuwxyzzywuqnkjijmnljihjloqpnlhecaa_^_`a_\\]\[Z[[YWWXYYQz|.GIIKLNRTUUUTUWY[\^adfeeedcddccbbbdefhknppqtuvwxxyz{{zxwvutsvyzywvvvuuvvvvwy{||{yvromklnqqnmmmnruusrolifecba`abb``_]\\\\[ZZ\]]Sz|/IKKLMORTVWVVWZ\]^`bdeefffffeffeddffhkrxyvvxyz{{zz{|}|zxwwwxxz|{zzzyxwyyxvvz}~~|yurpnoqsttrqqrsvyywuromkifdcbcddcba`_`__]\]_dg]z|1IKKMOQSUVWWYZ\^_`bccdefgiigfghgffghjox|z}}}~~}~~|zzz{{{|~}|{zz||{xx|}zvrmlqvwwwvuvvxz{{zxtrpnlifdefggedccccba```bksjz|1HKLOQTVWXXY[]_abddccdhiikkihiiiiijlnsz~}}~}~~}|}}{}}|vnmtyz{{zyzz{}}}|{wtrqpmjhiiihgfeeedddccbdmvlz|2JNOQSVXZ[[\]`bdefeccgjlmmmmmllmmllnrw{|}~~yxz|~}}~~~|yvuutqnnnlkkjiiihgghgfeflpez|3NRRRTWY\]^^_befffeefhjnqqpppppqpooqty}|}}{yyywtssrqqponnmlkllkjijmmaz|4"RUTSUXZ\_`abdgghhhhhikpttsrrrstsqsuw{~~}|||zxwwwzzvttsqpoooomklmk^z|6%UXWWXZ\^acdefghjllkjjlrvvutttuuuuvx{~}{zy|{yxvtrrqqqommlh^z|8'X[[[]_`aceggfgjmoonmmotwvvuuuuvwxy{~~||~}|zwusttsrqonlh^z|9'X^^_accdeghhhiloqrrqqruwvvvvvvwyz{}{wuuvwusqpnlgXz|:(Y_abdfgghhijjkmoruvuttuwvwwwxz{||~{wvvxywtrpmkfOz|:)[acehkkjjjjjjlnpsvwvuuvwwyz{}~}zyyzxvtronmiWz|;+]cehkoolkkklmnprtvuuuuwyz{}}|||yvtqooolbz|=-_dgjnsspmmnopqrsuvuuuvxz|}~}{wsppppofz|>/`egkpturpqrrsttuvxwvwxy{~~|wtrrsrofz|@1bfilqsssssstuvwwxyyxyz{|}zwvvutpez|A3dikmpqqrsstvxz{{zzyyz{|~†}|{yxvrfz|B4fjlnopprstuwy{}}||{{|~~~}zujz|D6glnprqqsuwxxz|~~~~}wmz|G8hmqstsstwzzz{~zpz|H9iortuttvy||{|}~}qz|H:jpsvwwwwy{||||~~pz|Ipuwy{{{{{|}~}jz|L@rwxz{{{}~~~{kQz|NCty{||||~oY=z|PEv{~~~~mY<̂z|RHx}nV8ƅpz|SJ|̆҆~jQ5e8 z|TL~Æ܆̆p[LIP[fv}gJ*䅿`.τz|VO††φԆņtR/2Li|eG#Յ[+ńlDz|YRÆÆÆĆĆĆƆ†††††ÆƆÆÆzX0؅Յ;_yv^A˅sFe@!z|ZT†ņƆdždždžƆƆƆņņņĆĆĆĆņƆƆņƆÆÆÆ}]5ՅrNBJ`}ȅ75;I`}ńz|b^†ĆņņĆÆÆÆĆdžɆɆˆ͆ІӆՆ׆׆ֆԆӆԆԆԆԆԆӆ҆ц҆ӆӆ҆цІІ҆ӆ҆цˆWޅTτ^1 䃻fB+-B_@_oodWOKNWg~Ą<`wz|daƆÆĆņƆdžƆƆƆdžʆ̆͆φ҆Ԇֆن܆܆نֆֆ׆׆׆׆׆ֆԆՆՆՆՆԆԆӆՆ׆׆Ն҆Æf#ۅ<鄙O؃w^E,#5Oy7Q\]_es„:fӅz|ebņɆȆņƆÆņƆdžɆʆɆɆ̆φІ҆Ԇֆ؆چ݆ކ݆چنچۆۆۆۆۆن؆نننن؆؆؆چ܆܆چӆ<셗;܄}!σ[;.**$*Q#IcwƄ6^݅7Dz|fcĆʆ͆̆ɆdžĆÆÆĆdžɆʆ̆Ά͆ΆцԆՆֆ؆ۆ݆߆߆݆܆݆ކކކ߆߆ކކކކ݆݆܆݆ކކކӆnX\ $"!;r%[ 6`օ#BZm{}z|geÆɆ͆φІΆˆȆƆdžʆ̆ΆцӆӆԆ؆ننۆ݆ӆ_(9V܂ȂĂԂ),-+%9qCԅRׅ?\uz|hg†ņȆˆφӆӆцΆ̆͆ΆІ҆Ն׆؆ۆ݆ކ߆ֆW􅄅 !܂Ȃ "-32-$*H̄t̅'}Å>[sz}ih†ĆņƆdžˆІӆՆԆӆ҆ӆԆֆ؆نۆކ܆Z􅁅Ђڂ)13/''7[G sՆ%Wtz}jhĆņdžȆɆˆІԆֆ׆׆ֆ׆نچۆ݆ކ熺iԂς *-+&&=joڅEdz}kiÆĆƆȆʆˆ͆ІՆ؆نچچچ݆߆߆Ն%957Ƃ&)'+G̈́,q؆6~z}ml†ĆƆdžʆ̆͆φцԆ׆نۆކކ߆ 􆲆O߅eee˂؂ +(5>KmZ&@z}nmÆƆdždžɆ̆ΆφцԆֆ׆؆ۆ߆  ֆs"ۂ‚‚ւ,Hax:FDz}mm†ĆņdžɆˆ̆̆͆φцӆՆ׆نۆ܆ + + + (>ȄN˃MЂǂɂԂ=fЃC΅n +Mz}moĆƆȆɆˆ͆Іццц҆Ն׆؆ن܆߆   #ņVt) -Z9\~ۅ`&_ÆÆz}pq†ƆɆʆ̆ΆІӆՆԆԆֆ؆چۆ݆ +   (#DτX僈PAOlÃ+[҄CwM{ĆdžƆ†z}rrĆɆˆ̆ΆІӆֆ؆؆نچۆ܆߆  +  !+0\&[ +׃ƃЃ@rׅ,Qs߆0V{ÆȆʆdžĆz}rtÆņĆņȆˆ͆φцԆ؆ۆ܆݆ކކކ + !$,2"솝E慇-ۄm^czʄ-^օ/Jc{ņʆˆˆɆƆƆz}svĆƆdžɆ̆͆φӆ׆نۆ߆ #&&%$%&&&&%%%%%%$" !#%&'*.)܆?텡]'/X7Umdž̆ΆΆ͆ˆɆȆņÆz}uyÆdžɆʆ͆φІӆ׆ۆ݆߆ +  ! !##$),,--.///.--,+***)''')***+/3.؆P܅:bʆцӆ҆҆ԆՆԆӆ҆φ̆ˆȆņĆÆz}v{†Ɔʆ̆ΆцӆՆ؆ކ + !$&(''(**+,*+034677667766410/00///100005=?/ ކ{S8,&!(@`ц݆߆܆ۆن؆׆Ԇцφ͆ʆdžņÆz}vz††Ædžˆ͆цՆ؆ۆ߆ + !$'*,/11124565349>?@@?>=>?@@=9657766788767;AB;+ֆֆކ܆چ؆Ն҆цΆ̆ʆdžÆz}vy†ÆĆņƆʆ͆φӆنކ + #%),/2469:<<=>?>=<>DIKLKIGFGIJIFC@??>==>?@@?@CDB=6-"ކۆ؆ֆԆцІΆ͆Ɇņ†z}vyÆĆĆņȆˆ͆Іӆ؆݆ #'*.37;=>@CDGHHGGFFGJNRUVVUSQRSSQOMKJHEEFFEFHIJKHC@<83-&"  ݆ۆنՆӆӆ҆φˆȆƆÆz}x|†ĆņƆȆˆΆцԆ؆܆ "!"'-159=ADEFHLOQQOMNOQSWZ\_`aba_^^\ZXVTSQOPQONOQRRQMIFC@;730.,& ߆܆؆׆׆ՆцΆˆdžÆz}|ÆƆȆʆ͆φ҆Ն׆ۆ߆ "$()(+18<@DILNMMPTXYXUSUX[_dghjklmllkigda_^]\[[[YYZ[ZXWSPMKHC>9531+# ߆ۆچنֆ҆φ̆ȆņÆz}ÆņȆˆφ҆Ԇ׆ن܆ކ +$*,-.04:?CHNUXXWWY[_aa^\\`ejnrtvxxwwwurqpmllkifedcdeeca_[WUSNHC>854/(" +߆݆چ׆ԆІ̆ɆȆdžĆz}ƆɆ͆ІӆՆ؆ۆކ +"*/246:>AEHMW]`abbbcgijkkjkosv{}{{{z|~|uomklnoomjgc_]YSNHB<861+'" ߆݆چֆцΆ̆ˆɆņz}ÆɆΆ҆Նֆ؆ۆ߆ $+26:=@DGJMQW^adgjllnqqsy}|zzz|xttvxxxvsokgd_ZUOIB<950+&"  ݆نԆцφ͆ʆdžz}ņ̆цֆنچ݆  &-8>@BFINQRV\addfkorvy{{{xuokga\WQIB>;50+'%" ߆ۆ׆Ԇ҆ІΆʆz}ȆφՆن݆߆ +!'-5@HJLPTX[]`eijkmqux~}xurjb\WQKFA;61-+(&! ކۆ؆ֆԆ҆Ά‚z}ˆ҆ֆ܆ #(+/5;CLQU[`dgikorstvz|~ukc^ZWQIA=9531-'" +ކ܆نֆՆ҆Ȃz}…ΆԆ؆ކ &/469?DLSY`hmrvxxz}~‡·‡yphda\SJEA=;83-($ ܆ن׆ֆ̂z}ÅІ׆܆ +#*2:>@CJRZaelsy‡LJƇƇчه̇pje]TNIEB>93.)# ߆܆چ׆΂z}ąՆۆ +&/7?DGKPW_gnrw}ćЇׇ݇݇χʇ͇Ӈ҇χ҇և̇ʇˇɇɇɇƇLJ׈4?҇smg`WQMJE?:4.(" +݆نЂz}ƅن߆ &/9AHMQW]cjqx}Ƈ̇χ·ʇȇ͇Շ؇߇ ݇҇·Їч͇̇ˇŇƇч؇ڇهևևׇׇԇчχχ·ˇŇɈBKyc{tqk`YURMG@93-(# ߆܆҂z}ƅ܆%/7@GOV]dintzÇׇ̇'&҇Їڇ߇ևʇÇćˇ҇ه߇߇ڇчƇÇ׈Xs{|ztjc^ZTMF?93-(# ߆Ԃz}Dž݆ &08@GOWbjpw~ȇӇ߇0/'%'#ׇهۇ݇ۇчԈ{|tmhb[TMF@93-'  ӂz}Džކ +%/7?GPYdow}Ň̇ԇ߇ +.AQZ^K)):;8:=>3'B5 + هڇć|vrjb[TNF>82+% ӂz}Dž߆ +$.6=EP]it~χۇ߇߇BeֈdKELQUY_c\G53>JJDHS@ + + ۇ҇ɇ{tkc\TLD>70*" +Ղz}Ʌ !-7>EN]n}Çч߇ +Bp։"q͊2_ӈi`glq~{icaZWbwY1   ܇χ}skd[RKD<6/& قz}Ʌ*5?FLWgy·݇#3:04ZƈX܊.Fψow}spbP9   ܇ŇzskcZSKC<4,%! +z}ʅ '2>HNT]j{χ /Rbhlsx|ԉF+|gوxMY׈ڈڈш#w>!*,#)'؇~wnd[RJB:3,)'" z}̅ +%/:FO[hqsżއ.Ga~LjBljKӈr>JΈ&K`XHCN_fW3[:/,,7:& #/1)҇wnd\TKC;51,'! z}̅#-7CNWcvч %FdĈ̈ۈ=zۉʉ6ۈZw-BGEE_܊߉T{I;=?:, 060#чvng_UKB<72,&z}ͅ *4>KV^i}ׇ&CdˈوΈȈ̈وHwWوщ'g}r~!q$߉͉ ^(҈]TF4)  .1* ӇÇzsi_TJC=61+# z}΅‡ &3>FR^fnӇ*Eb}ψ߈ԈˈĈĈۈ11ƈӉ,̉}<"芛Q -:$򉱉k!ӈwX9*%$.,$ۇχʇƇ~sg]SKC<6/'! z}ЅLJ"-;GOWcluχ ,AYm|шΈɈ܈ވΈՉuΊ0- +77{QӊK2@VGa)戦sK1(%!!,51&"$ ݇ˇ}qg]SKD;2*$ z}хˇ +)4BOV\gr~Շ):GWgqԈՈEщ,*!,OxGE8T[&ωv;创jF0(%#,96)$+,! +ӇŇ{qf[RJ@6.'  z}Ӆχ &2>LX`gq|և '05BYp͈шˉKɉ;qċh8ߊ\U>͉2∦vO4%#/1,1<;/*'Շȇɇ̇ˇŇzodYPF;2+$z}Ӆ҇+8GWdlu͇ *NzՈ݈ɈʈB׉1fڋ-siN퉸p-󈻈P/ &/@OMFD=/͇݇̇Ӈׇԇɇyl`UKB80)! z}ԅՇ$1>Pesz‡އ $Jwڈ و҈'_|rr։Dyދ*cwNJ7ωN!vL2!".CSRQUK:+҇·ׇևȇtg\RI@7.&! z}օه#4BVm{҇ !! *Gi͈҈҈؈߈ۈ׈ 1EGH^DhJn􊉊S/Ή?׈|^C1,2@JIKOF92#ۇчه҇ć~pe[QG>4,(%z}؅܇4EWjyȇ$10'*;PkȈ͈ɈȈψшΈ[iamy}s]LB?ETp.GIRQ<]ibF݉( шqQC@@BCA>4-2- هއχÇyodYND:2-)#z}څĆ'9IXixׇ0L\[Yaq҈݈߈݈3}v[$ É݊SGtmt8hzNE!*ՈwbQB?C>4(#.3؇̇wj_ULA7/*$z}څƆ%2?M\m}̇'KvĈ &9N[mYŠŊwF?FLOYaZK0`يފ4DҊNJL)2*"Ȉz_GBE>6.(,2#އՇ̇~od\TI>3+' z}څƆ)6CP_p؇ "=fՈ؈ĈɈHpЊ=^H=VpS\fLIhm$N3'VagTA$爪fPLICDG@76+ ܇ՇʇwkbZRF90+$ z}څƆ,:FRasˇ&&'4Qyڈ B܉܉؉ ҉ԉ5.'z}څdž1?HTdvه)-+3EaĈA4  ̉ϊ_9Ќbы;Dec)͊}9@5]JcW/xN6׈xiehr{u_J2݇χ~qf\OC90) z}ۅȆ4BKWg|̇ ...?Wr܈KNr:"É͊cԋS*׏X3֌Ӌ_B_e@$^0!m'ӉcS>̈oP- ##! +ׇȇznbTH<2*!z}ۅɆ6EO\nه&/3MXfx҇!+B_t(Xz 'ߊ 1976Fo5ߍ؎C'7j<ËYD̊cmFp=Ԋ슟;։ՉQ6=)툸xvvfQ>5 E,!抵fo V̉ъՉtPSJ؈zcE;BC/чʇLJ|l^QE8,z}߅ц@P]l‡"5Qfov{ՉAcN-) RxFʍ+oώЍ׍[ˋD势s@%ӊnɉȊ6ɉzmc>\>@D-؇͇͇ȇscUG9. z}ޅц@O\nɇ#9Qi~# Չĉ݉#! +eŠ̊\4 +]"~\ѓRď`Z\T97B>O#銄+މԊ$DɉpX#ማn_nmME;߇߇܇هˇyhYK=1 z}߅ԆAN[ȯ&&z}؇J\m߇(hщS҉tVٜݏpZJBdЊĊy|`ωN +͈XFLK(ۇ +·reWG2z}؇I[lEV}{pgX;8dO + +#nیQ'a^2H‘2̫2O;D(̌ƌXc#ъ⊳r)74􉵉ljω^㈤hMYgD҇vhXG4z}هL^oHMhg\LA=<9:EQ_tuknĊ fNtc5 ˙NN.K~UƵ]bOLы^#Њҋ +芴^=R8쉴 ߉}%\[d<߇އևxgUE2z~ڇM_oهE2SU@$$JbeanҊk (L?aWy _7p:ߑHdWŋZߊʊ&*Z3+xv<;׈mN?ڇև‡xeRC1z~هK[k҇EӉ 6@%6NSF8@RWYg~,iԊ!%5ZȌL'4]a Ѻ5.L3Y)o9Ïj= ɋCϊ݋;FŊ`kъ]׈sE' ׇćuaOA1z~ׇGVf|͇C7C?609JX]j:mʊ*=9@qDϛ0 l@}1jd5Y⌊ݍ+{GH͊MƉԊ<劔vÈoI.36'ևÇq\M@0z~ՇCSdyʇ@xɈ݈݈؈݈ %43.'3RbY[;k,UIC+t2şjWvLm-,ƭ2:ΎӍ<@印_ `Z-g{XHs􈪈oWNK8.DYVE5(ԇlYLA.z}ՆBRcwLJ߇8b|ψ*3+%$4PZOS/^ȋNE8W΍@񛺝9Ʀuǵ UݷX<,;dJ ׎tŊRce"ʉ&8, ʉV׈~unxseiiZF=1Ӈ~hVMC.z}Ԇ@Qbu‡ׇ&=DQf|78&+;EEUJ{38+HVǐ̖מנMħ( CgǓVˎ}ԋ!.~҉ 'rhmvd?),+ ӇuaRKA-z}ӆ=N_r·և߇(Lt߉$;8 *39@[ފ +1jĊf@ +: X?ަ$鮜sBW!xϒ)+Ŋnˊ:׉Ɋ + Љ҉ ۉhꈈWS]paЇjZPK@+z~߅ц:JZnćʇԇHtш42%6:7;ZӉ ]NJ+ދzțȞWϡ졀}5t'ېT Ίt?Vtω֊܉뉟dJKTlڈ` ˇzfZRK@.z~߅І7FVjÇ̇܇Dh| +&+0BID=PÉ"eŊъڊ9_Je=2KlͣN'b֗ !񊌊k ՉkӉΉn툉[QTZlڈS Ƈvh`WMA/z~ޅΆ4CSf{ч;Td|Έ(!!-=OTKLm%pҊ݊ tښb Ξ&H斟+ALE ->ԋ%kMϊJ ≥h\ۉOڈgaeinrwq?$ +܇ulcYMA.z}݅̆1AQcv؇ ':L\sÈ$!'>TWQ^Ί[|jUmÊȊ>ԝ(u1;M^k^%N0U‹2`Ut8 쉡haljˆmntrmcUMB1(&҇xodXLA1z}܅Ȇ.?O`qׇ5AFRau-ITR]Š6F0'PӋb唳):󚳚[t旦[Ԛ vl鐓Xw[\W^gx.:A@Dgˉ(c֊ Vk!Ϗԓϖ VSMȔL DŘP󎢍ˌ ފ؊ŊvhmX1󉘉.Èy`rxT)ЇćwiYMB0z}ۅƆ*7BQdv҇/MVTW[iԈ 1:>;>aމB}wj۔uɓJ|̔i Ȓ͎vŒL;Y}w7ىu#Ȉ{{[% ݇̇{obUJ@-z}مÆ&2=K^n{̇:HMMKTo +..1;Xډ+'ˊ}r܋*݌m{Џ]Б$xєh D id.pU;x{ QшtmruM݇ʇrcXOG<+z}؅".8EVeqLJއ*:BB=DZ{ш'3Qt׊c҉҉"[vꋂDm3:@ߋ+8c)ꊜK ߉ƉЉЉnЈ2%,Xz݊K*tOSjrbB߇ʇ{pke_ZVQME:)z}Յև!%)/9ER_l{· 9IJCES`ilk.*&5Pi~3 +# 0RreljV ';>;AIVwi.(Y#Έbzk>؇‡qiga\XURMD8&z}Յԇ#*5AO\huՇ&42)+@YjpjqЉ!$/Fiʊ?ߊɊeb%P{k쉬R#2$:#‹>C?vCAQ숎Ȉb/ ˇqkgb[WTQKA4z}ԅӇ )3>JVam{ȇ.Pgne`x&;]zJ܊Ԋ>틈4"A_@lZZ8ЉpJ9',CGYfيۊȋW#}nbщY鈎ɈM" އχ{snic]VRNG=0z}ԅ҇(09CMWbp~ч;Xc]T^ !0E`ijxA􊽊?#׋s8+;'ꊭnȉgv^]-F0쉏)|w\8ЇysplcZRNH@7,z}ӅЇ $+3;DNYguևDZVLKdމ$((;HHXyˉ2{ՊԊvьЋ}Y_}r+؊ɊΊXʼnΉىHC }|zG g,:؈bTQI<*Їwtrk^RLGA91(z}҅̇$,4=GSan{ɇ݇ 4QQH@Hf-'%,=hߊ(OȊЊ͋YB0L~Y NJ?鉸ˉȉ:ۈˆ͉*3ՈiL<51,!|vtpfXMGC=6/$z}хʇ &.7AMYerՇ$݉ƉÉtX=%܈ƈ}z[A5+ }N vi`ZTNJFA8-!z}Ѕ† ")0:DLWgx͇܇#/2$3Ojۈ @lϊW|s=܉։aC)ˈ~rv}rbJ1$ ևƇqe[TOKGC=3(z}υ +#+6?GP^l|‡Ӈ%( 4Nbrֈ܈ڈ߈J͊8A"쉺ĉ߉ˉ]*׈ԈֈӈĈxou}xyug]YQ@*  Ƈwk`UMHEA<6,$z}΅ (29AJUcpȇ؇-DRZf~ȈLjĈΈ 9v؉։vs׉ӉÉb$Ԉ̈ˆuggy~}pXF=9. Їwnf[OGB>:4.&z}ͅ +%-3;CLYesˇڇ  +&=KMN_yψ%Rt}vfRK`jH܈Ljˆticfu}vlfW>)χxohaWLC<73-'! +z}ͅ "(.4;DOZfw͇هއ %:A?GZm|҈ڈ$8<:5''CbtseTB&܈Ljyoga_[XUQJ@1 ߇هχ{umd]UK@71-'"z}̅ +#(-3HR`p|LJ̇χЇ̇LJ·އڇ·!3>J`{xlmĈɈ̈Έňy`SJ@<=7& ڇՇӇч̇ȇćunie`\XSMF<2)  +z}˅ #(08?DPakvȇهۇӇۇ #4BVgkdYUav{\F>80+%ՇЇчЇ͇Ƈrf_[YUQLF?6-$ z}Ʌ !*39=GS]iwʇԇׇԇև +1HVWQHEMZfo}{|uY?0+& ܇ۇۇՇЇ·͇ɇ‡zl`YVUQLGA:2)! ܂z}ȅ $+39AJT^juŇ͇͇̇̇чۇ-=B?:9<@K^xzlmoh^Z[]abZK9)#܇؇ԇӇ҇ч·ɇƇƇŇ|zwoe]XTQLGB;3,%܂z}Dž #)19BKT^iu‡ÇćŇ͇ևه؇݇*010-&,B[f\TXXOE><;974." +؇χȇƇȇˇʇƇ}vroke`[VQJD?:3,&  ۂz}Dž !)08@IS^is|ȇχ҇Їчڇ &'  #'*3>?5,))'# ! އԇ̇ćÇŇÇzuolgb^ZTOIC=61*$ ւz}Dž "(.6@KT]fov}ćƇɇʇȇ·ׇڇ݇   ؇чʇÇ~{vpkgb]XRKEA<5.(# ӂz}ƅކ +"(0;DJQZbipx~Çʇˇɇч܇ڇׇчʇć~~}{unic\VQJC>83-'! ҂z}Ņۆ $,5<@HOW^emtx{~̇ۇ܇؇ۇއӇLJʇχ͇·ччӇׇهևч̇Ňzz||xrkd\TNGC>80+&  ߆҂z}Ņ؆݆ '.38@GLS[binps{͇ՇևևڇهӇЇ҇؇ڇՇ͇ćÇƇŇÇ{xwvtojc[SLE?=91)$  ߆тz}ŅՆچچ݆ !&*18>BHOW]adipw}ćɇʇ̇̇Ƈȇ··Ƈ{vsqokgb^YSLE?;82*" ކςz}Ņ҆ֆ׆چކ߆ $+05:?EMSVZ_cfktLJŇ~|xqligea[VSPKE@<83+$ ܆΂z}ąΆӆֆچ܆݆#',26962-&  ۆ͂z}Å̆҆ֆ؆چ݆ $*/5;@CFHHJP[ht~wvxxvslgc^ZVTROKC;50.*& ߆چ΂z}̆ӆՆՆ׆ۆކ߆߆ $).37:<>@@CLXep{~wtrqnkfa\VRONLJE>6/*'$ +߆݆؆˂z}ʆ҆ӆ҆Ԇ؆ۆ܆݆݆܆݆ +#&*.137;;81+%! ߆݆ۆֆɂz}džΆΆΆцԆ׆چچننۆކ +!%(,1678=DLU^gmqtuvxzzz|~zupje`[XUOJHHEA>;61+%! +݆ۆ؆ӆȂz}ņɆɆˆΆ҆Ն׆ֆՆֆچކ +!&+/159>BHPY`eikmopopqrsw{~~|||||{{zvpkgc_XSOKHFEC>:63.(#  ݆چ׆ՆІĂz}†ņdžˆ͆І҆ӆ҆҆Ն؆܆  $',159<@FNTY]abdefgfefjmosxywxzzy{{yvutsqqpqplfa^\XQKGDBBA<72.+&! +߆ۆ׆Ԇ҆͆z}Ćdžʆ̆ΆφφφцԆֆچކ #'+.27:>CHLRVWY\__\Z\^afklmortrpooonljhfeedc`\XVTPJE@==>;60+(# ݆نՆцΆɆz}ÆƆȆʆˆ̆͆φцӆՆ؆܆ކ +  #%'+/259=AFJLOSVVTSSVartiehkkifdccba`^\[ZYWUQNMKGC>:99961,($ ߆܆׆ӆΆˆƆz}†ĆņdžȆɆ̆І҆ӆՆ؆ۆ܆ކ +  "%(+-/37;>AEILLLLMSgl_``_]\ZYYYXWUSSRPOMIFCA>;855431-(%! ߆܆نՆц͆ɆĆz}†ÆÆņņdžˆІ҆ӆՆنچۆۆކ  #%(+.036:?BDDDFM^z|dWVUSQPPPQRQOMMLJHFEB?;87520/.,*'#  ކ܆چ׆ԆІ̆ȆÆz}†ÆĆƆʆ͆ІӆՆ׆؆؆چ݆ !#%'*-059::;=BJWZSNMMLHGHHIJJHGEDA?>=;85210-*))%"  +܆چچ؆Նӆφ̆ɆÆz}ÆÆÆƆʆ͆ц҆ӆՆֆ؆܆ކ߆  !$'+.01248>DHHGGHGDBA@ABBA?=;9875420-,+(%$$  ݆ن؆׆ՆӆцΆˆȆÆz}†Æ†ĆȆˆΆІц҆ӆՆنۆۆ݆߆ + !$')+-27;?AABBA?>=;::99860,./.--+)'$"  ކچ׆ՆԆӆцΆˆȆņz}}††ĆȆˆ̆ΆφІІӆֆ؆نچ݆ !$&).259::;;9998642220-$#''''%"  ݆ۆ؆ՆӆцІΆˆdžņ†z}zĆȆʆ̆͆͆ΆΆІ҆Ն׆چ܆݆݆݆ކ !$(,.1321122110.,,+)(!"""! +ۆ؆ֆӆцІΆˆȆņ†z}uyÆņȆˆ̆ʆʆˆ̆φ҆Ն׆؆چچچۆކ !%'(**))**)**)''%$$! چֆӆ҆І͆ˆɆdžÆz}ru†ĆƆdžƆdžȆʆ͆ц҆҆Ԇֆ׆׆ن܆݆ކ  !"#####$$#"  ݆نֆ҆φ͆ʆɆȆņ†z}sv†ÆĆņƆȆ͆φΆΆцӆԆՆ؆چۆކ  ߆݆݆݆ۆن׆ԆІ͆ˆɆdžņÆz}ux†ÆĆdžˆ͆̆͆φІ҆Ԇֆֆن݆߆ +  +݆چن؆ֆԆӆ҆φ̆ˆȆņ†z}uy†ņȆɆʆ̆ΆφцӆӆԆ׆چ܆ކ    +߆ۆ؆ֆԆ҆цІφ͆̆ʆdžÆz}ux†Æņdžʆ̆ΆφІцԆֆنچچچۆ݆ކ߆߆߆ +  +ކ܆ن׆ԆцІφΆ̆ˆˆɆņz}ux†ÆƆɆˆ͆ΆІ҆Ԇ׆؆׆ֆ׆؆نچۆ݆߆ ݆چ؆ֆԆцΆ͆͆̆ʆɆɆdžĆz}uw†ĆƆȆʆ̆͆ΆІӆՆՆՆՆԆՆֆ؆چۆ݆߆߆چֆԆ҆цφ̆̆ˆɆȆdždžƆÆz}suÆĆƆɆʆʆˆ̆ΆІ҆ӆӆӆԆԆӆԆ׆چ܆ކ߆ކކ݆ۆ׆ӆІΆΆΆ̆ˆɆƆņņƆƆ†z}rs†ĆdžȆȆȆɆʆ̆φІІцццІцӆֆچ݆߆߆݆ۆن׆ӆІΆ̆̆̆ˆɆdžƆÆÆÆz}qs†ÆĆƆdždžȆɆˆˆˆ͆ΆΆΆφ҆Ԇ׆ۆކ߆߆߆݆چ׆Ն҆φ͆ʆʆʆȆƆĆƆ†z}pr†ņņņņƆdždžȆˆˆ̆͆φцӆֆن܆ކކކ߆߆ކކ݆ۆ؆ՆӆцφˆȆȆȆņ†ÆÆƆz}npÆƆ†ÆÆņdžȆȆʆˆ͆ΆцԆ؆چ܆܆܆ކކކކ݆݆݆݆ކކ݆܆ۆنՆ҆цφ͆ɆdžƆņƆz}ln†ÆņƆƆȆɆʆ͆І҆ֆنچچۆ܆݆ކ߆ކ܆ۆۆچچۆۆچۆކ݆ۆچچن؆ՆцΆ͆ˆȆĆƆ†z}kkÆÆÆņƆdžɆ̆φІԆ׆ننچچچۆކ߆ކ݆݆݆܆چنننننننچۆ݆݆݆ۆ؆׆ֆֆ׆Նφ̆̆ʆdž†z}ii†Æņdžʆ͆φІӆՆֆ׆؆نچۆ܆݆ކކ݆݆ކ߆ކ݆ۆچچن؆ֆֆ؆ن؆׆نۆ܆܆ۆۆن׆Ԇ҆҆ӆц͆̆ˆɆƆ†z}ig†ÆņȆˆ͆φІц҆ԆԆՆ؆چچچچچچچۆ܆݆߆ކ߆ކކ߆߆܆ن؆׆ֆԆՆ׆׆ֆֆ؆ۆچچن׆ֆԆ҆φΆΆ͆ˆʆɆȆņz}gf†ÆdžɆʆˆΆІІІцІцԆֆֆ׆ֆֆ׆׆؆ن܆߆݆ۆ܆݆ކ߆߆ކކ݆چنچކ܆؆ֆԆӆӆԆՆՆԆՆ׆؆؆؆ֆ҆цІφ̆ʆɆȆdžƆņÆz}fdÆņdžȆȆʆ̆ΆφΆ͆̆͆φцц҆҆ԆֆֆՆֆن݆݆ن؆نچچۆ܆܆ۆۆۆۆنֆՆ׆ކنֆԆ҆цц҆҆ӆӆԆՆֆֆՆ҆φΆ͆ˆɆȆdžņĆƆz}fdÆĆĆņƆdžɆˆ̆̆̆ʆɆˆ͆͆͆͆ΆцՆՆԆՆ׆׆׆ֆֆֆֆֆֆֆ׆׆׆׆ֆՆՆӆԆچ݆نՆӆ҆ІІцццц҆ӆӆԆԆ҆φ͆ˆɆdžƆƆņÆƆz}fdņ†ĆȆȆɆʆʆʆɆdždžɆʆʆʆɆˆΆІІ҆ԆՆԆӆԆԆԆӆ҆҆҆ԆԆԆ҆҆ӆӆц҆ԆՆ҆ІІφΆІцццц҆҆ццІφ͆ˆȆdžņƆ†z}dbĆȆφˆɆɆɆɆdžƆņƆdždžȆȆɆɆʆˆΆцццц҆ццІІІІццццццІφІІφΆΆ͆̆͆φІφφІцІφ͆ˆˆˆȆņņņ†z}b`ˆц̆ȆdžȆȆƆņĆĆĆņdždždždždžɆ̆͆ΆΆΆ͆͆ΆφφΆφΆΆІӆ҆ІΆΆ͆̆ˆˆˆˆ̆͆Ά͆͆͆Άφφ͆ʆȆȆdžņ††z|`]†Ȇ̆ɆƆņĆņņņĆĆĆĆņņņņdžɆˆˆˆˆʆʆˆ͆φΆΆ͆ˆˆΆццΆ͆͆ˆʆʆʆʆˆ͆Ά͆͆͆͆͆͆͆ˆɆȆƆÆz|_[ÆņƆƆ†ÆƆÆĆĆĆĆĆĆƆȆɆʆʆɆȆȆɆˆ͆͆̆ˆɆɆˆ͆͆ˆˆ̆ˆˆʆʆʆ̆͆̆̆͆Ά͆ˆʆɆȆȆdžƆz|]YƆ†ÆÆĆĆĆĆƆdžȆɆɆɆȆdždžȆȆȆȆȆdždžɆʆɆdždžɆʆˆʆʆ̆͆̆ˆ̆Ά͆ˆʆɆȆdžņĆƆz|[VÆÆĆņƆƆƆƆƆƆdžņĆņƆÆÆÆņdždžņĆĆņƆȆɆʆ̆̆ˆˆ̆͆ˆɆɆȆdžņ†z|ZT†ÆĆĆÆƆ†ÆƆ††ņņ††Æņdžʆˆʆˆ̆̆ʆȆƆƆņƆz|ZT†ÆƆȆȆɆˆ̆ˆɆdžņĆÆz|[T†ņƆƆȆɆɆȆȆƆņƆz}[SÆņņƆƆņņņĆÆÆz}ZPÆĆĆƆ†z}YM}††z}UJz~Ć}z|PEuz~Ć҆Նʆyz|L@pw{}}Ć/,Æzz|I/`filnquxzz|~{qz|>._egilosvwxz}}zoz|=-^cehjmrttuwyz{|~}znz|<,]ceghknpqsttuwy{}~}~~||{{zoz|;+\bdfghjmopppruwxyyxyz{}~|zyxwwpz|;+\aceeegjklllnqtvutuvwx{}~~~}zyxvtsmz|<-]`bdccdegghhikoqpqsuuvz|~~~~|{ywvtrpgz|=,\`bbaaaacdeeffikloqrruy|~}zxxwvtrqnbz|<+\_``_^^_abcdedegjmnnptx{}~~~}}~~}zvututrpok`z|;)Y\]]]\\]_`bcdddehjklnquz|}~~}|}}|}|zvtsrrqpnmi^z|9&UYYZZYZ\]_`acccdfhiklnquxz}}{zyz}~{z{|{{yurqpopomlkg[z|7$SWWVVVWY[]^^`abcdfgiklnorvxxwwwwz{|~~{xwxxxwuronmmmljjieZz|5!PTTTTTUXZ\\\^`abccdfiiiknqrrtvwwwwxz}~~~~~~~~{xvuttsrpnmlkkigghhdXz|1MQRRSSTWY[[[]^_`abbcddegiknoqtvvvuvvy|~}~~~~~~~~~}|}}|{|}|{yvtsrqomlkkjihfdeffbWz|/KPPQRSTWXZ[[\\\]`aaaa`acdgklmpstsrstvx{||}~~~}|{z||zz|||}~~}}}}}}~~~}||~|zxxxwxwvsqpoljihhhggededc_Tz|.JOOPQRSUWYZZZZ[\^```_^^_`cgiijnoooprtuvwy{}|zzzzyywuwz||}}{yxyz{|}~~||}~~}}||{{{{{|}~}|||{zzxwvvutuvuronliggffffedeeda\Pz|0KNNNOQRSUXXXXZ[\]^__^]\\\`deefhjklmorstuwyyvtuuuutsrsuvxyz{zyxxxz~|zz{}}}}{zyyyxxxxxy|~}{{||ywutsstsrstspmkhfffededcdddb`\Oz|2KNMMMOQQSUUUVY[\]^]]][ZZZ\^_`bdefhjknprtwwvspoonoppopqqrtwyzywuux~|yxxyyyxwuuvwwuuvvxz{{{||zwutrqqqpoopomkhfeedbbbabbba`^[Pz|0IMLKMOQQRSRSUWY[^^[XXXXXXXYZ\]_``dghjmqsttspmkihilmllnnnoruutttuw}zwwvvuutrqrtvustuuvwxy{|zwutsrponlkklljhfedca__`aa_^^][XNz|-FJJILOPQRQQRTVWYZ[XUUUUUUUVVXZ]]^acegkoqppomlkheeghhikllmopoopstvx{}yuuvvutrpopprttssttttvwyxvutrponmkihhiigeddc`^]^`a`]\\YVTIz|*CGGHJKMOQPQRSSTTUUTTUUSQRTVVVXZ[\_abeilmllklopkeccdefhijlmmlmnqrrsssrqrtsrqpnnnnoqrssssrrsvwusuurmjjigeeggfeba`_^]]^`_][ZXVSODz|(?DFGHHIKOQRRQQPPPQQRTTRPQTUUVVVWY\^_aeghhhgm|vgccccdfghjkkklmopponnmmnoonlllmlkklnqrrqrrsttsrstoigffcbcddcb_\[[\\\]][XVUTSPK@z|% ;ADFGGFHLPQPPOMMMMMOPPPPQQSTTTTTVXZ\]`abbbck|hbbbbceffffghjkmlkjjjjjkkjhfgiiiiiilnoopqqrqponnnjeccca`aa`_^[XVWZZZ[ZVSRQPOLH=z|" 8>BDEDCFJNNMMMLLLKJJJLMNNNOQRQQSTVXZZZ[\]]]bqzoda_^`bddccccegiigfgggggggedefedegghjkllnoponlkihgd`^_``__^][YVTTUVVVWTQOONMLIE;z| 6;>?@ABEHIIJJKLKJJHGHIIJJKMNOOPQSTVWWWXXXXYZ^a`_^\\_bdcaabcdffeefffedccb`bild`bddeghijlmljihfdca^[[\]__][YWTRSSRRRRROMMMLKIFB6z|49:;<>@CDDEFGHJIHHIHHGFEGIKMNOPPQRSRRTTUUVWWXXYZZZ\_cgeaabdeecbcddcb`^]]\^gjb^_``bdfghjigeccb`_\[YXX[]][XURPQRRPONNNMKJJJHEC?2z|16779;>???ACDFFEEGGGGFEEFHIKMOOOOOONNPQRUWWWVUUVWY\^dmkb_`cdc`__`__][YYXXY]^]]]]]`bdeefeca_^]\ZXYXWVXYXVSQOOQQOMMKJKMKHHGDB@=1z|,33468:;;=?@ACCBCDDDEEDFFGGIKMNNNMMLLLMORUVWUSSTVZ\]bmka]]_a`^]\[ZYXVVUVVWWXZ[\\]_`aaaba`^\[ZXVUVVTTUUSQONOQQNKKJHGHKJGFC@?>;0z|(.023568:;<<>@??@A@BCCDFFFGHJKLLLLMLKJJKMPSUUTTTW\]Z[ab]ZZ\]\[ZYWUTSRRRSTVVVWXY[\]^^^^]]\[YXVTSRRPOQRRPNLNPPOKIHGECCEEDB@>>=8/z|%*-/01357899:;;;<=>@ABDDDFHHIJKKKLMKIHGGHKPSTTSTW[YTTXYWVVXYXVUTSRPPQPOOQTTSTUWYZZ[[[ZXWVVVTQOOOMKKMNNNMLMMLIFDB@>;;=?>=<;::6.z|#(*+,-/0234667667:<=?ABCCEGGGIJIIJJIGFDDEGLPQQPQSSQNOQRQQRSTRQONMNNOOMLLOPQQQQTVVVVVVUTRQPONMLKJIGGHIJKJIGC?:41.,)),02100//0.&z|#%&&&'(*,.011113579:<=>?ACCCDFEEEDCB@>>@CGKLMLLLKJKKKLLLMNNMLKIHIKLKJJKLLLLMOQRQONLKKJIHGEDDB@??><;;;:73-# + + +  z| "#$%')**,.//036:977898898653457;?ABCDDEFGFEEFFGGFFGFCBBDDCCEEEECCDEDA=:730.-,+)'%# +܅υƅ…ą…ƅʅĂz|˅    %+/246:>AA?>=<;:98777642110010/.,*(" +ۅօх̅ȅƅ…~oc\XSLFCBCIRWVTTUX[Xz|w؅υ҅҅҅҅҅Ӆԅօ؅ۅ݅ޅ$+1442.*&" ߅Ѕytuvqib\VPLGA;6/$݄؄Ԅ΄Ƅ˄Єфτ΄Є҄Մтz|iy………ÅąŅŅƅDžʅ̅ͅЅׅ +ׅυȅwfRB81) ބք΄Ƅ~m\MEB>6.)&%'.7;;98:?D@z{YGmnmllnopqsuwz|}}ɅЅօڅۅօυƅtg\PE<4.*&! ل„}naWRRSQG<3,%  σz{Kn>@?>>@BABDFHLNNNPSUVXYYY[\\]_^\[YWUTSSTW\dlu}|rj`SC0 Մʄ{lX?.&؃σ˃˃˃ƃyqeVG8,&$  !)'XTENSION= 'IMAGE ' / Image extension BITPIX = 16 / Bits per pixel NAXIS = 2 / Number of axes NAXIS1 = 132 / Axis length NAXIS2 = 288 / Axis length PCOUNT = 0 / No 'random' parameters GCOUNT = 1 / Only one group BSCALE = 1.000000E0 / REAL = TAPE*BSCALE + BZERO BZERO = 3.276800E4 / EXTNAME = 'SCI ' / Extension name EXTVER = 3 / Extension version INHERIT = F / Inherits global header ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2010-07-02T19:55:23' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:59:26' / Time of last modification CTYPE1 = 'RA---TAN' / R.A. in tangent plane projection CRPIX1 = -66.6459346493144 / Ref pix of axis 1 CRVAL1 = 306.871725219222 / RA at Ref pix in decimal degrees CTYPE2 = 'DEC--TAN' / DEC. in tangent plane projection CRPIX2 = 143.15396209091 / Ref pix of axis 2 CRVAL2 = 37.3733845984976 / DEC at Ref pix in decimal degrees CD1_1 = -2.9402995676385E-4 / WCS matrix element 1 1 CD1_2 = 1.36706049992071E-4 / WCS matrix element 1 2 CD2_1 = 1.37399803289694E-4 / WCS matrix element 2 1 CD2_2 = 2.93132453649777E-4 / WCS matrix element 2 2 RADECSYS= 'FK5 ' / R.A./DEC. coordinate system reference EQUINOX = 2000. / Equinox of coordinate system MJD-OBS = 52162.6261539214 / MJD of start of obseration CCDNAME = 'EEV 9273-20-03' / CCD chip name(s) AMPNAME = 'EEV 9273-20-03, left' / Amplifier name(s) GAIN = 2.3 / Amplifier gain RDNOISE = 3.0 / Readout noise CCDSIZE = '[1:2048,1:4644]' / CCD chip size CCDSUM = '2 2 ' / CCD sum CCDSEC = '[1:2048,1:4608]' / CCD section(s) DETSEC = '[4097:6144,1:4608]' / Detector section(s) DATASEC = '[1:1024,1:2304]' / Data section(s) BIASSEC = '[1025:1056,1:2304]' / Bias section(s) FRMNAME = 'dhsFrame' / Name assigned to a frame of data. FRAMEID = '2 ' / Frame Id DATATYP = ' ' / Type of data frame WCSDIM = 2 LTV1 = 0.4375 LTV2 = 0.4375 LTM1_1 = 0.125 LTM2_2 = 0.125 WAT0_001= 'system=image' WAT1_001= 'wtype=tan axtype=ra' WAT2_001= 'wtype=tan axtype=dec' END q|wvusolkkiededa^[[[YUU\b^UOKD8)ȃiG"ׂ` ā{xuoeWF0݃wR+ւVǁԄބ܄ڄل؄ׄՄф΄̈́̄ʄƄĄÄ„xeM1Ѓ]3قXɁۄׄ؄ل؄քՄ҄Є̈́ʄDŽńĄ}aC$ރb6 +ۂYɁ ߄ۄلׄքԄф΄̄Ʉ„lN. 惽d8 ۂYɁ'*&" +ބ܄لքԄӄτƄrT3烾e8 +ڂY ʁ,.*'%#! +݄ۄلׄԄτȄwW4郿e7 ق[ ˁ./+*('&# ބڄׄՄ҄ʄxX6郾e7 ۂ] ʁ!00-,,+*'#  + ݄ۄڄ؄Ԅ̈́ÄyY7ꃿd7 +܂[ +ʁ!2310/0.)%#  +߄ބ݄ڄՄτĄzZ8ꃾe8 ܂Z +ʁ$4654320,(%#  ߄ۄ؄фĄ{[8郿e8 ۂZ ˁ'7987642.+(&$"  + ބڄ҄Ä{[7ꃿe8 +ڂ[ ˁ(9<::9641-+)'%#! ߄ۄфÄ{[7郾e7 +ڂ\ ˁ*;>>=;8630.-+(&$#"!  لфĄzZ7郾d6 ڂ] ́,=@@?=:84210.,)(&%#! ݄؄фÄzZ7郿d7 +ۂ]́/?CB@><964320/-+)'$"! +߄܄؄Є„yY7g9 ܂_΁2CGEBA?<:86541/.,(%$#  ݄؄ЄĄyY7g: ܂_΁7HKIFEC@><987420-)'&%"  ܄ׄфĄxY8h: ۂ_ρ<:9641-+*)'$"  ބۄׄЄ„xZ8g9 ۂaρ@NQPONLIDA><:852/.,*)&$! +ބڄՄ̈́zZ7g9 ܂a΁?NRTRPRPJDA><:9631.,*&""!  ܄؄ӄ̄yZ7f9 ۂ_́?PUWUTXXOGC@><;9630-+$ "#"!  +߄ۄׄӄ̄yY6郾a2҂W ȁATYZXWXWPIEB@><:7431.&!$&%$"!  +߄܄؄ӄ˄wX3 +݃}L낼EFX\[ZYWUPKHEB?><86651)#'*)&$$#  +߄܄ׄфʄnF胷T$łd)遷K\_]\ZXUQNKHEB@>;9962-)+.-*('%#!  ݄؄ӄ̄m@߃Oe7́N]__]\ZWSPMKHDB@=<:6310021.,+(&$#" ܄Ԅ˄vO#ǃk=܂|Q(ҁR_aaa^\YVSPMKHECA>;755335420.+(&%%" +܄ӄĄfE ̓tI삺Y.Vceeda^\YUQOMKHFD@=;976665420-*)('$  +؄ʄhI'ヽkA삼W)݁Zikieb`^ZVRPOMKIFB??=;;:86531/-+)'$  ܄˄gH'ჿvP(ׂ}LЁ_ppkfcb`]YUSQONKHDBB@??=;86431/,)&$" + ބτcE$݃qL&؂a5 āeuslgeedb]YURQPMJHFEEC@>=;8642/,(%%$! + τa@كiDт]7Ёfutokihhd`[XUSRPMKJJIFB@?=;9641.+('&$! ҄~]<Ӄb<ƂzT/ ҁfuurnkjhda^ZWTSQPONMLHDB@?@?;741-)()'# ф}\9σ\4 傽nI%΁Ł́ссāfuvtqnkgda_[WTSSSSQOMJGDBEMMD;50,)((% ЄyX5ǃ}V-܂e?ʁŁԁ灺euxvsplhec`\YVVWWVTSPLHEDM^`N=4.+)(%  ΄rQ/ +僾rK$ӂZ5ځǁǁׁ&5@GH: +ȁewzxuroljgb]ZZZ[ZXWVRNIFFO^aN<3/.+&̈́mJ'ރg@ȂwO)ԁāȁׁ1Kcvp+ԁgy|zwusqole`_^^][YXWSOKHHKPPF:42.' DŽhE Ճ^5 傼lE ЁÁɁׁ1LgȂӂׂՂʂHށi{}|zxwvtnifdba_][YVSOMJIHGE?82,!Ąb>̓}U+ق_:΁Łāʁׁ0Khق ς_䁦l{~~}|zyvqmkhedb`\YWTQNKJHE@8/#„[7ÃsJ ͂{S. ܁ρʁ΁؁/Kgۂ0GXac^P6p遧o~~}{yvsplhfec_[YWSOKIF?4%ބyV0 +⃺jAnH$ځҁҁځ-Ieڂ3Qm] {끨r~~{uokhgda^]ZTNHA7'܄tO)ۃ^6䂶b<ہׁ܁0Kdۂ4QnÃ؃ڃŃy쁨ryrmjhd`^]XPF9(քmG!ӃT)ՂX1߁*Gkڂ6TqɃ'-)僽$큨szrmjgc_[VM>*Ԅe>ǃuJȂtN* (D`Ȃ7Ut̓!;SenhX>ă%tyqljfaYN?+Є_6僼h>炻hB!&A^}ڂEbyЃ +'A[srM"ƃ$큩uxqmh_SA)̈́X0܃]1ڂ\9&@\z؂8Wwԃ+Gb{ȄɄzP#ƃ#쁩wtkbVC+ DŽ{R(ԃ|P$͂yS2&?Zxׂ7Wwփ/Lf˄Մڄ݄݄ׄĄ}R$ƃ"쁩yqcS@' +„tK ɃrE낿nK. $#ᄻkA샿e8 ߂dD+([vƄۄ  ߄ʄX*Ƀ"|聧fiW;̈́}Oe6ۂpZLIO]r͂4Wyރ @`˄  ߄ʄX+˃$}聨91ƄqG냸X*ӂq_WXcv͂3Ux"Aaӄ #" ̄Y,̓&~遨ㄿi;烸Oʂvieky˂5Xz&Ffք#(*++*'#  ΄[.΃)遨a4փPǂ}vwɂ 0V{ă&Hi„ۄ(.01110/,*'$  Є^/΃(遨>&̃l@ʂʂ /Sw +-MmȄ &.3678754331/-,(# ф_0̓$|聨Ճc7̂˂ -Qv -Pq΄#-5;>>===<:875310/-)%  +҄_0̓!z聨lS/ тς*Ot -QtԄ,7>CEGGEB@??@><96300.+(%! ӄ_0̓!z聨قʂɂق +Mr -Quׄ$4AIMNNNNMKHECBBBA>96320.*'%! +ׄ`0̃ y灨ނւӂւ0Ru .Rvڄ(:ITZ\[XVTRQPNKIGFEEB>:8520.+(%! ބ_0ʃx灨؂1Sw 1Uxބ/APZaeffd`][XUSQONLJHGFC?<9620.-*'# ᄸ^/ʃ"|聨2Rwà 3Z5JZdkopomkjfba_[WUSRQOMKIGDA>;62/.-,)&! ݄_0˃)ꁩKsƒ7[~Ʉ 7L^ltwyyywtpnkhfea]ZXVTTROLJHFDA<730/.-+'# +ل`1̓(ꁩ9`˄ *BSbox~}zwspnkjhd`][XWWUQMJIHEA=9620//,($  +لb3Ѓ#z灩*\ӄ.G]oz|yvsqomjfb`^\ZYVSOLJIEA>;7410/-($! ڄc4҃!w恩΄7Ocu~{xurqnkgdba_[XVTQMKIEB>;8520.+'$" ۄc4҃!v恩3Sk~{ywtrolifeb^[XVTQNLJHD@<9531.*&$"  +ۄc4҃ u偨Wz|{yvspnljhc^[YWTROLKJGB=9632/+(&$!  ۄd5҃ u偨}{ywtrpqokf`^\XUSPMLKGA=:8530,)(%#"!  ݄d5҃t䁨~{xvusstrmhdb_YVTQPOKFB?<:730-+)'%#!  ݄f6Ӄ t䁨|ywvtssqligfa[WUSQOKFDA=:730.,*)&$" +ބf7Ӄ!u䁨~{yxusroljhfb]YVTQNJGEA>:730/-+*'%" +ބh9Ӄ!u䁨~{zwuspmigec_\XTPMJGDA=:731/-+)'%! j:ԃ!t䁨~|{xuqmifeca^ZUPMKGCA>;9520.,)(%! +„k:փ sと}yvrnjfdca^[WROKGCB?><941/,)(%! l; ׃ rと}yurplhdb`][YUQLGECA@>:630-*'$"  +m= ׃rと}yusqnkgda][ZVRMIGDB@=:751.*&$"! +m? ׃rと…~ywusqnieb_\ZUPMKIFC@=:853/+'$!  „n?؃q⁨DžÅ|{zwtpkheb]YTPNMKHEB>;8651+(%"  +ńp@كpၨɅŅ…}}{wqmjgb^ZURQOMJGD@=9761,)'$  ȄrAڃoၧȅƅÅ|{xsokfa][WTROMKIFB>;972-+'#  ʄsBۃoDžDžą…|zwtplfa_\YUROLKJHD@<:61.+&"! ̄uD܃oʅ˅Dž…{zvsplfb`_[WSPMLKIF@<951.*&$! ΄uDۃo̅υʅÅ~|yvsokfca`\XUROMKIE@<963.*'%" ΄vEۃnͅυ˅ƅÅ}zwtpmieca]YWTQNKHC@=;83.+)&"  ΄wG܃nͅͅʅȅƅą}{yvrpmifc`]ZUQNJFB?><83/-+&$" +΄yH݃n̅̅ͅ˅ɅƅÅ~|{wtrnkhfc`\WSNIEA?>;732/+'%#  τyGރn̅ЅЅ΅˅ɅŅ~{xusoljgeb\WSOIFCA>;7430,(%#"  ЄyGރm҅҅ͅЅ΅˅ȅŅą…~|xvspnkheb]WRNKIFC?<8530,)'&$!҄yH߃mЅӅӅ҅хЅͅ˅ȅŅÅÅ}yvsqnlhe`\WROMJGD@=952/-+*(%! +ӄzI߃l߁ą҅҅҅ӅӅԅ҅΅Ʌƅą…zusqnlhaZYXTPMJGDA>:61..,*(%! ӄzI߃k߁ŅӅԅӅԅօׅԅ΅ɅƅŅ…~zvtqokf[RUWTPKHGEA>;71..-*'$"҄zIރkށŅօ؅օօ؅مՅ΅ʅȅŅÅzwurokdWOTURNJHGD@=;9520-*'%! ҄zH߃jށƅم܅څڅڅمՅυ̅ʅDžą|yurokbVSVTPMKIGC@>=;851.*(%! ӄ{H߃jށͅ߅݅܅܅مօх΅˅DžŅ{wsojb[ZYUQNLKHDBA@>:51.,)%! Ԅ|I߃jށԅޅ݅ۅمׅԅЅ̅ȅƅÅ}xsojfb^[WTQNLHECBA?;620-)&" Մ|I߃kށم݅ۅمׅՅӅυ˅ȅŅ~xrnjhd`]YVTRNJGFCA?;741-)&" Մ}Jjށ܅݅م؅օӅυ˅ɅƅÅ|vqmjgca^[XXWSOLHDA?<840,)%!քLボj݁݅߅ۅمօхͅ˅ʅȅą~yuqmifdb_\ZZZWSOIEB@<840,)$ քM䃝i݁ޅ߅܅څՅх΅̅ʅɅƅÅ~yvsnheca_][ZZXTPKGDB>951-($! ׄN僝h݁߅߅܅مՅ҅Ѕ΅̅Ʌȅą}{xsnidb`_^\ZZXUPLHEC?:62.*&$  ؄P僞h܁߅߅ۅ؅Յ҅хυͅʅDžÅ~{wsokeb`_`_^\YUQMIEB?:62/+(&! +لP僝g܁݅ۅ؅օԅхυͅʅƅ…|xtplgdbabca^ZVROKGC>:62.+)&" ڄQ䃜f܁߅݅مօԅ҅υ̅ɅŅ…~yuplhfdcccb_[WSPNJD>:51.+(%! ۄR䃛eہ߅ۅօԅ҅ЅͅɅƅązupmjgfdccb_[XTPNKF@;61-+'$! +܄Q䃜dہ܅ׅԅӅхͅʅȅŅÅ~ytqomigecba_[WTPNLGB=82.+'$! ܄Q僝dځۅׅՅԅ҅Ѕͅ˅ȅą…|xtssqmiebaa^ZWTRPMIC?940,(%! ۄP惞eہڅׅׅօԅԅ҅υ̅Džą…}ywuusokgdba^[WUTRPJD?:62.*%! ݄Q烞eہۅممׅׅօԅ҅΅ʅDžƅÅ|yxvsoljhea_\XVTSPIC>;841-'! ߄S烝dځޅ܅څمم؅ԅ҅Ѕͅ˅ȅÅ}zvsolkieb_\YVSQMFA>;842.(# T 胞cځ݅ۅۅ܅څօ҅хυͅɅÅ|wspmkhdba^[WSOJEA?;842.*&# U!郟cځ݅܅݅ޅ܅؅Յ҅ЅͅȅÅ}yurnjgedca]XSOJEB>;741/+'$ ᄶT!ꃠcځ߅݅܅څօӅхͅɅƅ…~{xuplheedb]XTPKFB>:8530,'# ㄷU!ꃟcځ܅ۅڅօӅхЅ΅ʅą~{ywtpkgdc`[VSQLGC?<:852-'# 愻V!郟cځ߅݅܅څօӅ҅҅Ӆ΅Dž~zwwwunhdb^ZVSPKFC@?=:62.*&"  +넿Y"ꃠcځޅ݅ۅ؅ԅ҅ӅՅхɅ~zwuutmgca^ZVSOIECA?=:62.,($""" „[$샢cځ߅ޅ݅مԅӅӅӅυʅÅ{wsqnieb`]YVTOHDCA>;:61.+'%##% ]&cځޅ݅څօՅӅυ˅Ʌƅ{vrolhdb_\XWVPHDC@<9852/*&%$%(! ]&cځ߅݅ۅم؅ׅӅ΅ʅʅȅązurokifc_[XVTOIFD@<9641.+(&$&(^(cځ ޅۅ؅؅ׅօ҅υ̅ͅɅŅ}yvrnkjhc_[WUSOKHEA><83/-+)'%#!a-􃦃bځ +߅݅ۅم؅ׅԅххЅͅɅŅ}yvrmjigc`]YWUQMJEA@>:4/-*(%"kC bځ ߅ޅ܅څ؅ԅ҅҅Ѕͅʅƅ…~zwsnkifdb_][XSNJEA?=:50-*&#„\Ճ%cځ + +߅܅مՅӅ҅Ѕ΅˅DžÅ{xtpligec`_]YSOJEA=;:50-+'$Ą-dځ ޅڅׅԅӅ҅Ѕ̅ȅŅ…|xspmjheb`_]YTPMGA=:940.,)%!Ą*dځ ߅ۅ؅օԅӅх΅˅ȅą|wsqnlifa_^]ZVTRLD>:730.,)&# Ƅ|c.ƃ%cځ +܅؅ՅՅԅӅЅ̅Ʌƅ…}xtromkgaZ[][XVTNE?:620-*(&$  ɄwK%cځ ޅمօՅՅԅх̅ȅƅ…}yvspmlhaXY\[XVRKE@<731.*)'$! ΄~QÃ'cځ +߅܅څׅՅՅӅυ˅ȅą|yvspnlje^\\YWVQKEA>9420-+(&# Մ^&ʃ)cف  +߅݅مօՅԅх΅ʅą}yvromllid_\XWVSNHC?;852/-*(%"܄l2у,cف  +݅مօՅԅ҅΅ȅ…~{wsommmlf`\YWWTPJEB?;731.+)&$ńw<؃.cف  ޅمׅՅԅх̅ƅ|zuqponmhb]ZXWTOIEDA=731/,*(%΄Eރ/cف  ܅مօԅх˅ƅ}{xtrqnkhc]YVURLGEDA<8420.+(& ׄM0b؁   ߅ۅׅӅЅ˅ƅ…~|yurqmjgb\XTQNJFEB?<9631/,(&  V1a؁##߅ۅׅ҅υ˅ƅ…{wsqoljfa[VROKHFD@><97630,(% ̄_3bف&&  ߅ۅׅӅυ˅ƅzurpnkhd_ZUROKIGC@=;::840-)&!քh5d݁()# ߅܅؅ԅх̅ȅ…{vrpnieb^ZVRONKHC@><<<951.+(" ބp7i䁲)*%  ݅څօӅ΅Ʌą|wsplfb`^ZVTRPMIFCA@?=;741/.)"儷s9n큹)+'"ޅ܅؅Ӆ΅ɅDžŅ~{wrmheb`^ZXVRNMLJHFCA?><99;=92,$t;t)-*% ܅ׅ҅ͅʅȅŅ}zvojfdb`][XUOLLLMLHECBCEGO\bXI;72*DŽv:Ճ+iꁺ&:;4.(! +܅؅Յ҅΅Ʌƅ|yurnkjhfc_\ZYWTQNMKHECA=951.+$ لh/̓&d⁳':;4.'  ۅօӅЅͅʅƅ|xurnkjhfc_^\ZXVSQOKGDA><93.,)"ф[#ă"aށ(<<4-& ߅مՅх΅̅Ʌą|xvsolkigea_]ZYWTQOKGDB?=82.+' ̄|M_ہ)=<4.' + ݅؅ӅЅͅʅƅ{yvsqnljheb_\ZYVROMJGEB?;5/-)$˄p?]ځ*<;50( +܅؅ԅхͅȅą|yvsqomkheb`]YVTQNKHED@<72.+'!„c0\ف*:94/(  +܅؅օ҅ͅȅÅ}{vspomkheb`]YVSPLHFDB>952/+& ۄW#냝Z؁,::5.' +مօׅօ҅̅DžÅ{wromkjgeb`]ZVSOKGDC@>952/+&!ׄ~Mვ Y؁.=<70' +݅Ѕ҅ͅӅЅ˅ƅ…{wsomjgfeda]ZVRNKGC@?=:51.*'#ՄxFك Xׁ)<<70& +ޅӅυхх΅Ʌƅ}yurnjgggeb]YTPMKFA>=<951-*'#քvA ҃ Xׁ%894-$ ߅مօӅυ˅ȅą}xuqmkigda\WSOMJD@><:751-*&"ׄyAσY؁'75/'! +ۅ؅ׅ҅̅ȅƅ…~zwtqnjfb_ZVSPMIDA?<9741.+%!؄|B΃Yف,94*$߅ׅׅօхɅŅÅ~{xvtpjfb^ZWTQNIEB?<96320+& ׄx?̓Y؁,83*#݅օՅԅͅƅ…~{xvtpkgc_[WTRNIEA>;85320+& ׄt<΃Xׁ,84+$ ߅؅ӅυɅÅ|xvtpkgc`\XURLHEB>;8642/*%!ׄt<ЃWׁ,83,% م҅̅Dž}yvrokfa^[WUQIDDB?<9752-)% Մu>ЃWց,73,% ޅ؅Ӆͅƅ~|xtpmid`^[XVPEACB@=:752-($Ԅu?ЃWց*51*# ܅؅ӅͅDž~|zwsojea^\[YVQHDDB@?<962-&!Մv?ЃWց%2.'! ߅ۅׅ҅̅ƅ~|{wtqnhb_]\[YVQKHEBA@>:51+%!Մu?уXց-+$ +߅څօхʅą~{zyuqolfa_^\[ZWPKIFCA@=950*%!Մt@ ҃ Xׁ*(" ݅څׅӅʅÅ~{yxspnjda`_\ZYVQMJGC?=:850+%!ԄvB ӃXׁ('  ߅څׅօӅʅÅ|zvsrojeba_\YWURPLHC?;8752,% ԄwB ԃWց)' ܅ׅӅ҅υʅą}zvturlgdb`]ZWUSPLGC?;8543,#ՄvA ԃWց)' ޅمՅхυ˅Dž…|ywuwtnida_]ZWVRNJFB?;8431*$ Ԅu@ ԃWՁ&$܅؅ԅЅͅɅŅ{xuttrlgc`^[YXUPLIFB>:7420+&!҄t@ ԃ WՁ! ܅ׅӅЅ˅Džązvsrpnjfb_\ZZXTOLJGC>:7420,'! ӄu? +ӃVՁ ߅ۅׅӅυʅDžązusqnkheb`][[XSMJIGC>:7420+&! ӄu? +ӃVՁ + ޅۅׅӅͅȅŅÅyvtqmhedca_^[WRNJHGC>:631.)# Ԅv@ ӃUՁ݅مׅхʅŅÅzvsqnifedb`_[VROLHEA=9410-'" Մu? ԃTԁ ߅ۅׅԅЅɅą…}wsqpmigda`^[VSPLGC?<720/-(" Ԅt? ԃSԁ݅مօ҅ͅȅą…ysqpnkgb``^YUSQLF@=:51/.-)# Ԅt@ ԃRӁޅۅ؅ՅЅ˅DžÅysqomieb``^XSQPKE?;730.-,)$ ԄvB ԃPҁ܅څׅӅ΅Ʌƅąxsqolhfc`^\WQNMID?9531.-,)#ԄxC ӃNс߅ۅ؅Յυ˅ȅDžƅ…|wtrpmjgc_\ZVQNKGC?9542/-+($" ԄxB ҃LЁޅمօх̅ɅȅȅDžŅ…{xuronkfa][YVSQMGB>9752/,+)&# ҄u@ ҃Jρ܅؅Յх΅̅ɅDžDžŅÅ|ytpnmjd_\ZWTSRME@<9752/-,+($ τs? +ԃH́߅ۅׅԅӅ҅υ˅ȅƅą}xspnkgb^[XTRRPJC>;97410..,)&!τs@كF́݅مօԅԅԅυ˅ɅƅÅ{uromiea]ZURPPMGA><:7410/.,(%"Єs@܃Eˁ݅؅օօօӅ΅ʅɅŅ…~xsqnkhd_\XTQONKFA><;9520/.*'$ ҄t? ׃Dʁ܅؅ׅׅօ҅΅˅ȅą{vspligb^[XTPNKJFA><;962/.-*%"҄u@ ԃCʁޅ܅م؅ׅԅЅυ̅ȅÅ{wrnjhe`][WSOLIGD@=<;:62/.,)$ ҄u@ ӃBɁׅۅمׅՅ҅Ѕυ̅Dž…|wrmifb^\YVRNJHEB>==<9630.+'" ӄu@ ӃBɁ΅ۅ؅ׅՅ҅Ѕͅʅƅ~{vplifb^ZWTRNKHFB>==<8531-)%! ҄uA ҃BɁ҅܅م؅Յ҅υ˅ȅŅ|xtpmjfb^YVUROLIGC@=<:6320,)$  фuA +уBɁԅ߅܅څ؅Յх̅ȅDžązurolieb^ZWURNKHFEA><9520/,)%  +Єv@ ЃBʁ҅ޅۅمօ҅ͅɅƅąyspmjgda]ZXUQNJFEEB?=;730.+)%  фu?уCʁ҅ޅۅׅӅυʅƅą~xrnkheb_\ZWTPLGDDCA>=;851,)'#фt>уCʁ҅߅څՅхͅȅŅ…}xrmjgda^\ZWTOJEBCB?==;851+('#Єs>҃Dˁ҅ޅڅԅЅ̅Džą{wsnjgc`_\ZWSOKHFEDA@><840-*'$τs>҃Eˁԅ߅ۅ؅Ӆυ˅DžÅ}wurmifc`_]ZVROMJHEDCB@=720-*(% ΄r> ӃF́Ӆ݅مօхͅʅƅ…}ytrokhfc`^\YVSPNKGECBA@<61.+('$̈́q> +ԃǴӅޅڅ؅ՅЅ˅ɅŅ{wrpmihfc`^[YVTRNJGECB@>941.*'%"̄r> +ԃǴӅ߅ۅ؅ՅӅυ˅ȅą~{vrnjhhfb_][XWUSNJGFCA?;731.*&# ̈́s> +ԃǴ҅݅څׅՅӅ΅ʅŅ|xurnjhgeb_]ZYXURMIHFC?=:52/,)&" ΄s? +ԃǴӅ݅مօӅυ˅ƅ…{usqnkigeb`][ZXTQMJHEA><951-+)&! ΄s@ ԃF́Յޅ؅Ӆ΅ʅDžąytrolihgdba^\ZWTQMJGC@>;961-+)%  τs@ ԃǴ҅ޅׅЅɅƅą…~xtqmjgfedcb`]YVTPMIFC@=:840-+'" τs@ +ӃH́̅څӅͅDžŅÅ{wuqmigecccc`][XSOKIGD@<:73/-+&! τuA ԃH́ƅׅׅӅ΅ʅDžą{xwtpmjhecbbb`][XSNJHGD@=;72/-*&  ΄uB ԃÍÅ҅ӅЅ̅ʅDžÅ}ywuroljhfeca`^[ZWSOKIHGB<852/,)% ̈́uB ՃÍЅӅЅ΅̅ȅÅ~{xvtspmjigfda_\YWURPLIJJC3+/0.+(# ̈́tA փÍυӅ҅Ѕͅȅą}zxvutrnljhfc`^[XVTQNKHIIA-").,)%#!̈́tA փJ́υ҅хυ˅ƅ…~zwvvuromjgeb`][WVTPLIGEC>3+,-*&$$" ΄s@ փJ΁…Ѕхυ̅DžÅ~zwuutromjgec_\YVTSOKHEA=:61/,)'%%# +τs@ փK΁΅υ˅ȅą|zwutsrqmifeb^[XURQMJFB>;840.,)'&%# +΄s@ ՃK΁̅̅ȅŅÅ~{zwuttsqmifda^[XTQNJGC@=962/.,)'%%# +̄r? ԃL΁ȅȅŅą…}zzxvutspmieba_[XTQMHEB@<741/.+(%$$# ʄr? ӃLρŅƅŅ|{zxwuromid``^YVTPLGDB?;74420+'%##! ʄr@ ԃLρ…ƅą~}|yupmkhda^[WURMJGEC@;66640,(&#! Ȅq? ՃMρÅƅą~~~ztnjhfda]ZWUPKIFDCA=9751.+)'# +DŽp> ՃMρąƅÅ}}~}zslhgec_[YXTOJHEDDDA=94/,*)(# DŽp= +ՃNЁƅƅ…~|zxqkhgeb^ZWWTNIGECCDB>93/,*)(" DŽq> փNρŅŅ~{xsnjhfd`]YVTQLHFEB@@?<62/-)'$ Ƅq? ׃Nρ~zuokigdb^ZXURNJGECA??>:63/+(%! ńp? ؃M΁|xsnjhfc`\XVTQMIGEB@?><974/*&#  ńo? ؃M΁~|zuqnkgda^ZVUROLIGEA?><:853.)%" ńp?كNρ}{wqnljfb_]XUROLJIGEB@>;9741.*&#  +ƄqAڃOρ}zxtoljgc`^[WROMJHHHEB?><9632/+'# ƄrCۃOρ~zxvrnkiea_]ZVRNLJHGGDA?><8422/+'# DŽsCۃNρ}ywuqnkhda_\YVRNMLJIFCA?=:6311.*'# ɄtDۃNρ~{xusomjfda_\YUQONLKIFCA>;7420/,)&" ɄuDكPЁ|zxvspmkheca_]YTQPNKJHEC@<8530.-+(%" ȄtD؃Pс{xvtqnljheca`]YTRQMJIGDB>;841.,+*'$! ȄtB؃Pρ|xusqnkigfcb`\WUSOKJIFDA=:841-+*(&"  DŽsBك P΁|xvtqmjihfdb_[XVRLIIHEC@=9620.+*(%! DŽrBك P΁~zxwtpmjihfc`^\YVQLIGFEC@<830/-+*(%  +DŽrBكPρ|yxwtqnkihea^_^YTPLIGFED@;620.+*)'$ ƄrBڃPρ~{yxvtspligd`^^\WRNKHFEDB?9410.+*)&#  ńrBۃQЁ|{yvtsrnjihea^\YUPLIGFCA@=8410.,)(%#! +ńqAۃ Rс}zyxvtqokihgeb_\YVPKHGFB?=;742/-,)'%$! ƄqAۃ +Sс~zxwuusqmigfeca^]\XQKHFD@=;9752.,+)'&# ńpA܃ Sс~{xwtsrqoljhfda^]\ZWPKHEB><:8762.,)''%  ÄoA܃ RЁ|zxvsqpnmlkigda][ZWSOKGC@=;98752/+(('$ ބpAۃ RЁ~}{yvsqpnmmligea]ZWTPMJGC?<:98531.+*)($  ݄p@ۃ RЁ|xtsrponmkifc`\YVSOKIFB>;:96420.-,*'# ܄p? ڃ SЁ|xusqqqpmkida^ZWURMJGC@=;:85430...)%# ܄o> ؃ TЁ~}{yusqqqoljgb^[XVTQMIFB?=;976430/.+&$# ܄m> ك TЁ}{zywutsssrnkiea\YWUSPLIEB@>:76531/.,($#! ܄k= ڃ TЁ~{xvusrstvuqmjfc_ZWVTQMJGDBA>:76420.,)&$" ܄l= ۃ TЁ|yvtsrrrtuuqkgc_\XUTROKGECA@>;97530-*'%#! ڄn?ڃ Uс}zxvusrqqqrspid`\XUSROLIFCA?><:88630,(&#! +؄n@ڃ Uс{yxwvtspoonnlgb^YVTROLJHEB?>=;98763/*'$" ׄm@ڃ Uс}zxwwvtqnlmmkgda\WTSRNLJFCA@>=<:8652/*&#  քnB܃ Vс~{xwwvurokjkkhda^ZVTSPMKIFB@?>><:8753/*$  ׄvL߃ Vс|xvwwvsqmihhgeb_[WTSQNLJHFC@>>=<98752-(" لĄ^&⃐ Uс~{wwwwurolhgfdb`\YWSQOMJHFEDA?=<;87530+&" ߄Єr,߃RЁ{~zxxxvtqmkhfec`]ZXWTQNKIFEEDB@=<:862/+'$! ܄ʄu"ЃyMρ{|zyywuroligedb^[XVVTQMKHFEDCB@=:974/+&#"  ՄRcG΁||{yxvspnlhdba_]ZWUTRPLJGEDCBA>;9751-($"  ԄZփKڂÁ{~}zxvroljfa^]\[YVTRPMKIFEDCA?<9753/+'$! ҄^(탰u2ɂ9ˁx~~|ywupljgc^\[ZYWUTQNLJGFDCB@=;9631-)&$! ӄl7ȃX0ȁw}|zxurnjgc`][ZYWVUSQMKIGEDCA><:852/+&#"  +҄zK߃r=%ā{~|zyvrolhea^\ZXWVUSROLKJGEDC@=:8530-(%"  фY)􃾃T 傋~}~~}zxwtoljfb][YWWVUTRPMKJIFCBA>;8520-)&$" +ф^5҃g4ǂn|}|||{xvtpljhd_XVWWWTSRPNMKJIEBA?<851.,)&%$! Ԅ_9xE₥J恬r~}{zzyvtqmigeb^YWWWWSPOMLMKIGDA@>:62/+)'%$" քa;郺S!򂿂{ʁk}}{zyyxtpnjfec`^[XWWVSPMLLLKHEB@><73/-*('&#  + ؄c>_-͂D쁳h{}{zywvuqmjfcb`][YWVUTSQNLKJIFB?><950-*(''&%# ؄e@ăf6ׂXÁgz}{yxurrpnkgc`_]ZXXWTSRRQONLJGD?<;852.*)'%%&)'  ؄„gBŃj; ݂`΁gx{ywuqmllkhda_\YXWXVTSRPNNMKIGC=:752/,)(&$"%*)  + ل„hCŃl=߂dҁfvywurokihfdb`]ZWVWXWVTRPNMKHGFB=8420-*('%#"#&# +݄ńiDƃk>eӁdtwurpnkifca`_\YVUVXXXVRONKGECA?<730/,)'$""!  +߄ɄjEŃk>ႫeӁ`rsponljgda_][YXWUUUWYVPMKHEC@<==710.*'$#"! ߄ɄlGǃk>ႪdӁ^oomlkifeca^ZXWWVSRQSUSOLIECA>;;:4/-+'$##! +ʄnI ȃl>₪bҁ]mmkjgdddb_\YVUTSQQPOONLIFC@=;8740-*(%#"" ߄˄nJ!Ƀn@オbс]klkhdbbb`]ZWUTQOPQOMLJHFCA=9642/-*'&%#! +ɄmH!ʃnA䂫bЁYhkjfb``^\ZWVTRONOPMJHFEB@>:631.,)'&&%" ʄnH ˃oA䂬cсTeiiea_][YWUTQNLLLKIGEB@=<:730.+)'&%%#  +̄nH!ʃpBカdсRbfec`\ZYVSRQNKJIGFEEE@;98731.,)'%%$#! ΄oI!̓oA䂫cҁQ_baa^[YWSPOOLKJGDCBDF?96530.-*'%##"  +΄nI"̓rC傭dсM\^^]\ZXUQONMLKJFC@>?@;752/-,+(&$"!  τpJ"̓sE傭eҁHY\[ZYXVSPMLKKJHEB>;97642/-+)(&$"! ЄrM$΃rD悮eӁFWYXWVUTQNLJJIGED@<96321/-+('&%#! +фrM&σsE悭eӁDTWVTSRPMLKJIGDCB=9641/.,*(%$$$#  ҄sM%σtF悮dҁ@PTSRQOMJIHHGDA@>9641/.,)'$   +ӄtN%΃tF炭dс=<:641/.,)'%" +ՄsM%σtF炭dс6EHIJIFB>;:85444320-,+)'%#  +ԄtM%ЃuG肮cЁ /35640)#!&))'&$"! +΄rM%уvG邯dЁބ߄ ӄjF ΃vH邯eЁ{yuv{Ԅ܄؄Ԅ҄τʄ„tW6ǃsFꂰeЁ-@Yu„nV<܃kA邱eЁ`mqokc[X\bkwރ4IZht}yrhZJ5ƃ]7傰fЁȂӂւԂ΂ƂĂׂ̂ +!=]ԃ(9FQWYWUSQPNI?3$˃pO+߂fЁEMOMG?9:ALYi{ʂ ,KgЃ%(*('%#σcC"߂sӁ \ No newline at end of file diff --git a/scikits/image/data/simple.fits b/scikits/image/data/simple.fits new file mode 100644 index 00000000..c54dcb1f Binary files /dev/null and b/scikits/image/data/simple.fits differ diff --git a/scikits/image/io/_plugins/fits_plugin.ini b/scikits/image/io/_plugins/fits_plugin.ini new file mode 100644 index 00000000..10fe672a --- /dev/null +++ b/scikits/image/io/_plugins/fits_plugin.ini @@ -0,0 +1,4 @@ +[fits] +description = FITS image reading via PyFITS +provides = imread, imread_collection + diff --git a/scikits/image/io/_plugins/fits_plugin.py b/scikits/image/io/_plugins/fits_plugin.py new file mode 100644 index 00000000..5c673bb3 --- /dev/null +++ b/scikits/image/io/_plugins/fits_plugin.py @@ -0,0 +1,147 @@ +__all__ = ['imread', 'imread_collection'] + +import numpy as np +import scikits.image.io as io + +try: + import pyfits +except ImportError: + raise ImportError("PyFITS could not be found. Please refer to\n" + "http://www.stsci.edu/resources/software_hardware/pyfits\n" + "for further instructions.") + + +def imread(fname, as_grey=True, dtype=None): + """Load an image from a FITS file. + + Parameters + ---------- + fname : string + Image file name, e.g. ``test.fits``. + as_grey : bool + For FITS images, this is ignored (treated as True). + dtype : dtype, optional + For FITS, this argument is ignored because Stefan is planning on + removing the dtype argument from imread anyway. + + Returns + ------- + img_array : ndarray + Unlike plugins such as PIL, where different colour bands/channels are + stored in the third dimension, FITS images are greyscale-only and can + be N-dimensional, so an array of the native FITS dimensionality is + returned, without colour channels. + + Currently if no image is found in the file, None will be returned + + Notes + ----- + + Currently FITS ``imread()`` always returns the first image extension when + given a Multi-Extension FITS file; use ``imread_collection()`` (which does + lazy loading) to get all the extensions at once. + + """ + + hdulist = pyfits.open(fname) + + # Iterate over FITS image extensions, ignoring any other extension types + # such as binary tables, and get the first image data array: + img_array = None + for hdu in hdulist: + if isinstance(hdu, pyfits.ImageHDU) or \ + isinstance(hdu, pyfits.PrimaryHDU): + if hdu.data is not None: + img_array = hdu.data + break + hdulist.close() + + return img_array + + +def imread_collection(load_pattern, conserve_memory=True): + """Load a collection of images from one or more FITS files + + Parameters + ---------- + load_pattern : str or list + List of extensions to load. Filename globbing is currently + unsupported. + converve_memory : bool + If True, never keep more than one in memory at a specific + time. Otherwise, images will be cached once they are loaded. + + Returns + ------- + + ic : ImageCollection + Collection of images. + + """ + + intype = type(load_pattern) + if intype is not list and intype is not str: + raise TypeError("Input must be a filename or list of filenames") + + # Ensure we have a list, otherwise we'll end up iterating over the string: + if intype is not list: + load_pattern = [load_pattern] + + # Generate a list of filename/extension pairs by opening the list of + # files and finding the image extensions in each one: + ext_list = [] + for filename in load_pattern: + hdulist = pyfits.open(filename) + for n, hdu in zip(range(len(hdulist)), hdulist): + if isinstance(hdu, pyfits.ImageHDU) or \ + isinstance(hdu, pyfits.PrimaryHDU): + # Ignore (primary) header units with no data (use '.size' + # rather than '.data' to avoid actually loading the image): + if hdu.size() > 0: + ext_list.append((filename, n)) + hdulist.close() + + return io.ImageCollection(ext_list, load_func=FITSFactory, + conserve_memory=conserve_memory) + + +def FITSFactory(image_ext): + """Load an image extension from a FITS file and return a NumPy array + + Parameters + ---------- + + image_ext : tuple + FITS extension to load, in the format ``(filename, ext_num)``. + The FITS ``(extname, extver)`` format is unsupported, since this + function is not called directly by the user and + ``imread_collection()`` does the work of figuring out which + extensions need loading. + + """ + + # Expect a length-2 tuple with a filename as the first element: + if not isinstance(image_ext, tuple): + raise TypeError("Expected a tuple") + + if len(image_ext) != 2: + raise ValueError("Expected a tuple of length 2") + + filename = image_ext[0] + extnum = image_ext[1] + + if type(filename) is not str or type(extnum) is not int: + raise ValueError("Expected a (filename, extension) tuple") + + hdulist = pyfits.open(filename) + + data = hdulist[extnum].data + + hdulist.close() + + if data is None: + raise RuntimeError("Extension %d of %s has no data" % + (extnum, filename)) + + return data + diff --git a/scikits/image/io/tests/test_fits.py b/scikits/image/io/tests/test_fits.py new file mode 100644 index 00000000..b9894cb6 --- /dev/null +++ b/scikits/image/io/tests/test_fits.py @@ -0,0 +1,77 @@ +import os.path +import numpy as np +from numpy.testing import run_module_suite +from numpy.testing.decorators import skipif +import scikits.image.io as io +from scikits.image import data_dir + + +pyfits_available = True + +try: + import pyfits +except ImportError: + pyfits_available = False +else: + import scikits.image.io._plugins.fits_plugin as fplug + +def test_fits_plugin_import(): + # Make sure we get an import exception if PyFITS isn't there + # (not sure how useful this is, but it ensures there isn't some other + # error when trying to load the plugin) + try: + io.use_plugin('fits') + except ImportError: + assert pyfits_available == False + else: + assert pyfits_available == True + +@skipif(not pyfits_available) +def test_imread_MEF(): + io.use_plugin('fits') + testfile = os.path.join(data_dir, 'multi.fits') + img = io.imread(testfile) + assert np.all(img==pyfits.getdata(testfile, 1)) + +@skipif(not pyfits_available) +def test_imread_simple(): + io.use_plugin('fits') + testfile = os.path.join(data_dir, 'simple.fits') + img = io.imread(testfile) + assert np.all(img==pyfits.getdata(testfile, 0)) + +@skipif(not pyfits_available) +def test_imread_collection_single_MEF(): + io.use_plugin('fits') + testfile = os.path.join(data_dir, 'multi.fits') + ic1 = io.imread_collection(testfile) + ic2 = io.ImageCollection([(testfile, 1), (testfile, 2), (testfile, 3)], + load_func=fplug.FITSFactory) + assert _same_ImageCollection(ic1, ic2) + +@skipif(not pyfits_available) +def test_imread_collection_MEF_and_simple(): + io.use_plugin('fits') + testfile1 = os.path.join(data_dir, 'multi.fits') + testfile2 = os.path.join(data_dir, 'simple.fits') + ic1 = io.imread_collection([testfile1, testfile2]) + ic2 = io.ImageCollection([(testfile1, 1), (testfile1, 2), + (testfile1, 3), (testfile2, 0)], + load_func=fplug.FITSFactory) + assert _same_ImageCollection(ic1, ic2) + +def _same_ImageCollection(collection1, collection2): + """Ancillary function to compare two ImageCollection objects, checking + that their constituent arrays are equal. + """ + if len(collection1) != len(collection2): + return False + for ext1, ext2 in zip(collection1, collection2): + if not np.all(ext1 == ext2): + return False + return True + + +if __name__ == '__main__': + run_module_suite() +