Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1095 1855 (N56W25 L=113)   465     52        1 (1061/1837-1128/1868)   0.0 W21-W29 N53-N57
      2 1467 1686 (N36W46 L=134)  1529    245        3 (1363/1644-1578/1740)   7.9 W34-W63 N33-N38
      3  540 1492 (N34E29 L= 58)   363     75        0 ( 495/1466- 574/1513) -29.7 E33-E26 N32-N36
      4 1650 1397 (N15W47 L=135)   394    110        0 (1598/1360-1686/1426) -33.7 W43-W49 N11-N16
      5  755 1406 (N26E13 L= 75)   280     98        0 ( 711/1379- 803/1425) -33.7 E15-E09 N23-N27
      6  389 1340 (N26E41 L= 47)   830    205        0 ( 296/1277- 451/1382) -26.6 E51-E35 N21-N30
      7  384 1083 (N10E41 L= 47)   555    129        1 ( 354/1025- 407/1142)  71.6 E43-E38 N07-N13
      8  251 1096 (N13E53 L= 35)   638     98        2 ( 217/1053- 285/1132)  63.4 E56-E49 N10-N14
      9  283 1040 (N09E50 L= 38)   887    168        2 ( 237/ 980- 332/1100)  54.5 E54-E45 N05-N12
     10 1844  967 (S13W58 L=146)   600    126        1 (1829/ 901-1862/1020) -90.0 W56-W60 S17-S09
     11  342  915 (N01E46 L= 42)   866    223        2 ( 277/ 864- 414/ 959)   0.0 E51-E39 S02-N04
     12 1384  881 (S13W20 L=108)   436    170        3 (1336/ 832-1427/ 927)  39.8 W17-W23 S15-S10
     13 1710  696 (S29W46 L=134)   413     97        0 (1687/ 652-1737/ 741)  63.4 W44-W48 S31-S25
     14  729  646 (S20E24 L= 64)   256    126        0 ( 695/ 603- 782/ 672)  23.2 E26-E19 S22-S18
     15 1445  607 (S31W24 L=112)   756    170        0 (1367/ 561-1509/ 656)  16.7 W18-W28 S33-S28
     16 1332  623 (S29W15 L=103)   246     55        0 (1308/ 594-1355/ 654)  63.4 W13-W17 S30-S26
     17  742  566 (S25E25 L= 63)   276     85        1 ( 703/ 529- 779/ 589) -53.1 E27-E22 S27-S23
  

    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