logo kso logo uni graz

20130821 - Filament data for 09 UTC
Available data for 20130821: »05 UTC   »06 UTC   »07 UTC   »08 UTC   »09 UTC   »10 UTC   »11 UTC   »12 UTC   »13 UTC   »14 UTC   »15 UTC  
Click on image to open sensitive map!
 
  Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  653 1291 (N30E19 L=119)   252    130        0 ( 643/1229- 669/1358) -90.0 E20-E18 N25-N33
      2 1835 1307 (N04W67 L=205)   562     59        0 (1813/1279-1851/1344)  15.9 W62-W71 N02-N05
      3 1876 1111 (S08W65 L=203)  1106    158        0 (1846/1027-1900/1174) -73.3 W61-W67 S13-S03
      4 1785  853 (S21W50 L=188)   566     93        1 (1754/ 814-1811/ 903)  38.7 W46-W53 S22-S18
      5 1612  819 (S18W34 L=172)   561    165        0 (1562/ 788-1668/ 856)   0.0 W30-W38 S19-S14
      6  870  601 (S16E18 L=121)   583    241        0 ( 828/ 517- 955/ 646)  26.6 E21-E11 S20-S13
      7 1406  443 (S39W15 L=153)   428     76        1 (1372/ 419-1437/ 469)  18.4 W11-W17 S40-S37
  
Download Binary data for this image

» All Filament Data of 20130821 ...
IDL code for reading binary data:
    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    
  • each filament has the value of Fil_ID
  • all values are in pixel
to get the area of a filament:
  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 
be careful: y-origin is upper left! you can transpose the data:
  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y