Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1199 1892 (N73W29 L=147)   724     77        0 (1178/1858-1231/1920)  14.0 W19-W63 N67-N79
      2 1710 1485 (N34W59 L=176)  1893    367        2 (1575/1357-1838/1573) -33.1 W47-W72 N24-N40
      3 1229 1519 (N36W14 L=131)   296     80        1 (1192/1496-1262/1532) -21.8 W10-W16 N34-N36
      4 1406 1416 (N29W26 L=144)   328    114        2 (1378/1362-1427/1468) -74.1 W24-W27 N25-N32
      5  837 1385 (N25E14 L=104)   399    144        0 ( 776/1350- 899/1408)  26.6 E17-E09 N22-N26
      6 1654 1371 (N26W47 L=165)   353     75        2 (1612/1351-1679/1400)   0.0 W42-W49 N24-N28
      7 1665 1043 (N05W44 L=161)   379    128        0 (1648/ 995-1685/1112)  66.8 W42-W45 N02-N09
      8  217  846 (S12E61 L= 56)  2448    208        5 ( 176/ 754- 266/ 940)  63.4 E68-E56 S18-S06
      9 1749  857 (S06W52 L=170)   488     42        0 (1723/ 835-1774/ 880)  -9.5 W49-W55 S07-S04
     10  343  757 (S17E49 L= 69)   551     70        2 ( 321/ 719- 370/ 788) -76.0 E50-E45 S19-S15
     11  505  690 (S20E35 L= 83)  1142    228        1 ( 439/ 614- 572/ 767)  39.8 E42-E29 S25-S15
     12  774  647 (S22E15 L=102)   848    221        3 ( 691/ 587- 869/ 711) -28.6 E20-E09 S25-S17
     13 1687  630 (S21W52 L=169)   550    138        1 (1638/ 597-1754/ 668)  24.4 W47-W58 S23-S18
     14  308  577 (S30E59 L= 58)   867    121        1 ( 256/ 550- 347/ 624)  14.7 E70-E51 S32-S26
     15 1058  558 (S27W04 L=122)   322     98        0 (1032/ 513-1094/ 600)  66.8 W02-W06 S30-S23
     16 1385  482 (S32W29 L=147)   352    104        2 (1334/ 466-1438/ 498)  -5.7 W24-W34 S32-S30
     17 1136  462 (S34W10 L=128)   568    110        2 (1079/ 431-1173/ 486) -23.2 W06-W13 S36-S31
  

    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