Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1438 1538 (N26W42 L= 61)   641    156        1 (1369/1495-1517/1570) -41.2 W38-W46 N21-N29
      2 1239 1421 (N25W24 L= 42)   613    183        2 (1200/1343-1275/1505)  38.7 W18-W28 N21-N29
      3  863 1380 (N32E02 L= 17)   588    149        0 ( 802/1317- 913/1434) -68.2 E05-W00 N26-N36
      4 1827 1298 (N02W65 L= 83)  1333    188        3 (1776/1218-1876/1371) -68.2 W61-W69 S04-N06
      5  931 1208 (N20E01 L= 17)   253     53        1 ( 905/1185- 959/1232) -56.3 E02-W00 N18-N21
      6  771 1130 (N19E13 L=  5)   385    131        0 ( 705/1090- 807/1159) -56.3 E16-E10 N15-N21
      7  417 1051 (N20E40 L=339)   388    148        1 ( 361/1024- 472/1085)  11.3 E44-E34 N19-N21
      8  212  959 (N17E61 L=317)   903    189        0 ( 197/ 867- 233/1041)  84.8 E63-E57 N11-N22
      9 1660  752 (S25W36 L= 55)   446    103        0 (1641/ 707-1679/ 804) -90.0 W34-W37 S28-S21
     10  822  728 (S06E18 L=  0)   312    118        1 ( 789/ 676- 864/ 767)  33.7 E21-E15 S08-S04
     11  773  494 (S20E28 L=350)   783    294        0 ( 642/ 461- 863/ 557)  -3.4 E38-E21 S21-S16
     12 1049  484 (S26E11 L=  7)   683    219        2 ( 965/ 434-1148/ 542)   7.6 E18-E03 S28-S24
  

    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