Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  971 1914 (N53E39 L=355)   552     69        0 ( 937/1904-1010/1926)  15.9 E42-E35 N51-N55
      2  864 1750 (N35E31 L=  3)  2046    249        3 ( 780/1680- 965/1823) -15.3 E42-E20 N30-N39
      3 1414 1733 (N50W11 L= 45)   638     54        0 (1377/1707-1444/1763)   8.1 W07-W14 N47-N53
      4  530 1659 (N22E51 L=344)  1198    225        3 ( 426/1617- 625/1700)  20.6 E60-E44 N18-N26
      5  566 1583 (N17E43 L=351)   916    256        2 ( 510/1499- 660/1666)  77.0 E47-E40 N11-N24
      6  900 1550 (N21E20 L= 14)   253    104        0 ( 861/1525- 953/1567)  31.0 E21-E16 N19-N22
      7 1691 1399 (N32W38 L= 72)  1560    221        0 (1604/1322-1768/1469) -17.4 W29-W45 N29-N35
      8 1310 1351 (N18W10 L= 44)   389    142        0 (1248/1302-1355/1414) -24.0 W04-W13 N16-N20
      9  848 1128 (S05E12 L= 22)   747    204        2 ( 776/1069- 920/1198)  56.3 E15-E09 S09-N00
     10 1900 1055 (N17W63 L= 97)  2965    388        2 (1858/ 843-1938/1187) -56.3 W56-W71 N05-N24
     11  145  971 (S25E68 L=326)   475     64        2 ( 132/ 931- 154/1008)  68.2 E70-E66 S27-S22
     12  681  878 (S22E18 L= 16)   323     88        1 ( 664/ 832- 698/ 924)  90.0 E19-E16 S25-S19
     13  246  821 (S33E55 L=339)  1061    230        1 ( 185/ 734- 333/ 892)  27.6 E67-E45 S39-S27
     14  492  840 (S28E31 L=  3)   416    123        1 ( 458/ 786- 513/ 888)  66.8 E33-E29 S31-S24
     15  667  589 (S41E15 L= 19)   421     87        0 ( 619/ 576- 713/ 607)  20.6 E18-E10 S42-S38
     16 1516  464 (S26W50 L= 84)   663    120        0 (1469/ 425-1562/ 511)  -4.4 W43-W56 S27-S24
     17 1534  292 (S35W70 L=104)   600     78        1 (1509/ 259-1559/ 321)  -3.2 W62-W80 S35-S33
  

    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