Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  968 1823 (N45E25 L= 36)   416     90        1 ( 935/1794- 992/1849) -23.2 E28-E21 N42-N47
      2 1043 1768 (N42E16 L= 45)   331     58        0 (1007/1756-1077/1779)  31.0 E18-E13 N40-N43
      3  189 1296 (S02E66 L=354)  1271    261        5 ( 156/1207- 213/1364) -77.5 E72-E63 S07-N02
      4 1820 1325 (N29W51 L=111)   510     59        1 (1795/1302-1845/1350)   0.0 W48-W53 N27-N30
      5  142 1165 (S10E68 L=353)   814    197        0 ( 116/1129- 188/1209) -18.4 E74-E61 S12-S07
      6 1683 1100 (N11W40 L=101)   273     86        0 (1650/1058-1707/1149) -45.0 W36-W42 N08-N13
      7 1835  994 (N09W55 L=116)  2021    289        2 (1802/ 848-1854/1137) -58.0 W49-W61 N00-N16
      8 1457  605 (S22W37 L= 97)   674    161        0 (1380/ 572-1533/ 653)  45.0 W32-W40 S25-S17
      9 1217  577 (S29W22 L= 83)   308     71        0 (1182/ 565-1262/ 597)   9.5 W19-W25 S30-S28
     10  986  519 (S38W08 L= 69)   215     46        1 ( 959/ 504-1014/ 538)  26.6 W06-W10 S39-S36
     11 1304  273 (S46W44 L=105)   859    123        3 (1251/ 248-1352/ 304)  23.2 W40-W48 S48-S43
     12  861  209 (S66W12 L= 73)   561     90        2 ( 838/ 184- 891/ 233)   0.0 W08-W16 S69-S63
  

    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