Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1083 1690 (N34E15 L=307)  1173    180        4 (1001/1634-1160/1731)  -8.7 E21-E08 N32-N36
      2  338 1539 (N09E61 L=261)   584     70        3 ( 303/1509- 366/1566)  31.0 E64-E58 N07-N10
      3 1339 1411 (N23W09 L=331)   254     76        0 (1321/1375-1355/1450) -78.7 W07-W10 N20-N25
      4 1190 1405 (N18W00 L=322)   385     90        0 (1151/1379-1242/1432)   0.0 E02-W03 N16-N19
      5  300 1391 (S01E57 L=265)   886    132        3 ( 242/1373- 369/1420)   0.0 E64-E50 S01-N01
      6 1847 1235 (N28W53 L= 15)   433     65        1 (1828/1198-1871/1260)  45.0 W50-W56 N25-N30
      7 1655 1104 (N13W36 L=358)   629    149        1 (1632/1047-1694/1172) -45.0 W33-W40 N10-N17
      8 1859 1069 (N18W55 L= 17)  1791    173        3 (1818/1002-1923/1152) -18.4 W50-W65 N14-N23
      9 1382 1067 (N04W19 L=341)   574    168        1 (1294/1041-1460/1107)  41.2 W14-W23 N00-N07
     10  504  889 (S26E30 L=292)   281     60        0 ( 473/ 870- 534/ 910)   0.0 E32-E27 S27-S24
     11 1800  765 (S00W58 L= 20)   619     71        2 (1769/ 747-1838/ 792)   0.0 W54-W62 S01-N01
     12 1710  589 (S12W58 L= 20)  1318    236        6 (1676/ 487-1746/ 697) -50.2 W52-W64 S18-S06
     13  289  642 (S46E53 L=269)  2657    304        6 ( 156/ 589- 435/ 684)   1.2 E84-E34 S47-S43
     14  561  561 (S46E22 L=300)   242     64        0 ( 536/ 526- 584/ 586) -31.0 E24-E19 S47-S44
     15 1641  505 (S18W57 L= 19)   425     66        2 (1619/ 466-1658/ 537) -45.0 W54-W61 S20-S16
     16 1493  486 (S24W47 L=  9)   936    198        2 (1413/ 461-1602/ 508)  19.7 W39-W54 S26-S20
     17 1467  330 (S33W57 L= 19)  2194    286        4 (1324/ 300-1596/ 372)  33.7 W46-W67 S38-S26
  

    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