Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  493 1566 (N53E36 L= 73)   566     73        0 ( 477/1523- 516/1604) -60.3 E39-E32 N49-N56
      2  740 1567 (N47E09 L=100)   722     97        2 ( 690/1536- 796/1591) -21.8 E13-E03 N44-N50
      3 1756 1482 (N15W68 L=176)   558     80        2 (1732/1445-1787/1518) -90.0 W66-W71 N12-N17
      4 1725 1462 (N15W63 L=171)   529     67        0 (1706/1431-1743/1493)  78.7 W60-W65 N12-N17
      5  710 1295 (N29E15 L= 93)   375    140        0 ( 671/1245- 767/1339)  35.5 E18-E10 N26-N31
      6  196 1117 (N26E63 L= 45)  1701    357        1 ( 125/1001- 285/1198)  25.6 E75-E52 N18-N30
      7  341 1121 (N24E46 L= 62)   686    159        0 ( 296/1057- 384/1184) -45.0 E51-E42 N19-N28
      8  812  829 (S00E16 L= 92)   706    209        1 ( 754/ 736- 868/ 904) -74.7 E18-E14 S06-N05
      9 1092  752 (S10E01 L=107)   267     63        1 (1063/ 718-1123/ 779)  18.4 E03-W00 S11-S08
     10 1727  722 (S28W43 L=152)   422     35        1 (1702/ 698-1746/ 749) -90.0 W40-W45 S29-S26
     11 1010  621 (S17E09 L= 99)   213     67        0 (1001/ 583-1022/ 654)  68.2 E10-E07 S19-S14
     12  476  383 (S24E58 L= 50)  2092    259        3 ( 383/ 295- 568/ 438) -59.7 E65-E52 S31-S19
  

    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