Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1351 1716 (N40W48 L=  6)  1242    171        3 (1287/1668-1418/1752) -54.5 W43-W52 N35-N43
      2 1644 1556 (N22W61 L= 18)   864    106        0 (1625/1512-1671/1606)  50.2 W57-W67 N18-N24
      3  643 1373 (N36E18 L=299)   963    257        3 ( 523/1342- 759/1398) -17.5 E28-E09 N32-N39
      4  961 1317 (N25W03 L=320)   589    207        2 ( 876/1259-1038/1350) -41.2 E02-W06 N20-N28
      5  913 1146 (N16E04 L=313)   369    176        0 ( 873/1099- 969/1176) -31.0 E06-E01 N13-N18
      6 1901 1094 (S11W67 L= 25)   719     78        3 (1870/1055-1919/1128) -90.0 W62-W70 S13-S09
      7 1844  996 (S15W57 L= 14)  1339    198        3 (1813/ 922-1882/1085) -70.0 W53-W61 S20-S09
      8  310  921 (N13E51 L=267)  1459    244        5 ( 263/ 819- 353/1025) -72.9 E54-E45 N06-N19
      9  145  914 (N14E72 L=246)  1017     89        5 ( 128/ 874- 169/ 954)  59.0 E75-E67 N11-N16
     10 1129  836 (S06W02 L=320)   270    126        0 (1076/ 802-1177/ 856) -36.9 W00-W04 S08-S04
     11 1463  629 (S27W19 L=337)   264     41        0 (1438/ 611-1482/ 651)  33.7 W17-W20 S28-S25
     12  805  536 (S18E24 L=293)   431    125        0 ( 745/ 524- 864/ 550) -14.0 E28-E20 S19-S15
  

    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