Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  650 1521 (N33E22 L= 11)  4451    699       14 ( 398/1314- 879/1671) -32.0 E54-E06 N17-N47
      2 1355 1615 (N33W30 L= 64)   382    107        0 (1304/1595-1405/1635) -26.6 W26-W33 N31-N35
      3 1146 1491 (N26W12 L= 46)   750    293        0 (1052/1404-1264/1588) -51.8 W07-W18 N19-N33
      4 1673 1516 (N24W54 L= 87)   336     14        1 (1652/1496-1693/1534)  31.0 W51-W56 N22-N25
      5  337 1453 (N32E50 L=344)   384     52        0 ( 311/1436- 370/1473)   9.5 E52-E46 N30-N32
      6  388 1303 (N21E42 L=352)   694    190        1 ( 363/1253- 434/1388) -60.9 E44-E37 N17-N26
      7  232 1241 (N19E56 L=337)   649     86        3 ( 206/1197- 266/1275) -76.0 E60-E52 N16-N20
      8  311 1196 (N15E48 L=346)   550     70        2 ( 274/1173- 354/1221) -20.6 E51-E43 N13-N16
      9  324 1112 (N10E46 L=347)   581    130        1 ( 281/1091- 399/1150)  -5.2 E50-E39 N08-N12
     10  144 1072 (N09E67 L=327)  2038    226        9 (  96/ 978- 226/1125)  31.6 E76-E55 N03-N12
     11  767  750 (S16E18 L= 15)   826    352        0 ( 675/ 623- 863/ 866)  45.0 E26-E11 S22-S09
     12  494  728 (S15E38 L=356)  1264    267        2 ( 468/ 593- 523/ 842)  72.6 E40-E34 S23-S07
     13 1461  728 (S23W26 L= 60)  1095    301        3 (1362/ 677-1578/ 828) -29.4 W19-W35 S26-S15
     14 1654  790 (S20W41 L= 75)   544    154        0 (1585/ 773-1729/ 810)   4.8 W35-W47 S21-S19
     15 1671  480 (S42W52 L= 86)   591     69        0 (1640/ 453-1702/ 513) -21.0 W46-W59 S43-S38
     16  351  459 (S30E64 L=330)   649     73        1 ( 324/ 419- 375/ 486)  21.8 E68-E58 S32-S28
     17  610  437 (S35E38 L=355)   356     72        1 ( 581/ 412- 639/ 466)  23.2 E41-E34 S36-S32
  

    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