Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  456 1666 (N34E50 L=132)  1158    257        1 ( 413/1630- 514/1717) -26.6 E57-E43 N31-N38
      2 1589 1421 (N25W36 L=218)   753    156        1 (1521/1393-1658/1447)   0.0 W30-W41 N22-N26
      3  606 1371 (N14E29 L=153)   210     60        1 ( 578/1338- 629/1394) -45.0 E30-E27 N11-N15
      4  373 1268 (N07E45 L=138)   673    107        2 ( 350/1219- 402/1327) -81.9 E46-E42 N03-N10
      5 1807 1295 (N19W54 L=236)   368     35        1 (1785/1267-1826/1321) -71.6 W51-W55 N17-N20
      6 1909 1118 (N10W65 L=248)  1451    223        2 (1864/1029-1943/1188) -68.2 W59-W71 N04-N14
      7 1274 1003 (S04W15 L=197)   973    329        1 (1215/ 853-1310/1153) -80.0 W10-W17 S13-N04
      8 1644  948 (S04W41 L=223)   399    136        0 (1588/ 904-1694/ 970)  32.0 W36-W44 S07-S02
      9 1645  804 (S12W43 L=225)   347    122        1 (1611/ 784-1684/ 818)  21.8 W40-W45 S13-S11
     10  403  772 (S23E41 L=141)   573    183        0 ( 360/ 708- 468/ 817)  37.9 E45-E35 S27-S20
     11 1793  640 (S20W61 L=243)  3752    481        3 (1658/ 527-1921/ 721)  31.0 W51-W76 S28-S13
     12 1678  646 (S21W49 L=232)   302     64        1 (1661/ 620-1701/ 672) -36.9 W47-W51 S22-S19
     13 1444  418 (S39W39 L=221)   716    206        2 (1345/ 396-1536/ 450)  26.6 W31-W45 S41-S35
  

    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