Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1065 1813 (N57E13 L= 94)  1213    307        0 ( 929/1776-1210/1850)  -1.7 E30-W03 N54-N57
      2  370 1631 (N30E69 L= 39)  1796    220        3 ( 308/1536- 445/1686)  70.0 E74-E62 N22-N34
      3 1036 1624 (N39E08 L= 99)  1147    318        1 ( 880/1585-1175/1685)  -4.8 E22-W02 N36-N41
      4  397 1496 (N21E54 L= 54)  1074    242        1 ( 357/1388- 436/1605) -65.0 E58-E49 N14-N29
      5 1614 1372 (N29W38 L=146)  3095    428        7 (1404/1303-1784/1447)  -5.0 W21-W55 N26-N32
      6  336 1077 (S05E47 L= 61)   593    148        0 ( 304/1010- 371/1142)  71.6 E49-E44 S09-N00
      7 1785 1101 (N14W53 L=161)   548    120        0 (1767/1045-1819/1133)  71.6 W51-W57 N10-N16
      8  591 1043 (S04E27 L= 80)   415    134        0 ( 559/ 993- 610/1102)  90.0 E28-E25 S06-N00
      9  856 1035 (S01E10 L= 97)   494    119        2 ( 811/ 993- 893/1086)  51.3 E12-E08 S03-N02
     10  308  914 (S15E48 L= 59)   296     83        0 ( 291/ 870- 319/ 954) -68.2 E50-E47 S17-S12
     11 1151  720 (S16W12 L=119)   763    285        0 (1069/ 668-1206/ 779)  35.5 W07-W15 S19-S12
  

    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

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y