Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  798 1505 (N30E11 L=194)   666    123        1 ( 741/1475- 856/1532) -20.6 E15-E07 N27-N32
      2 1648 1178 (N03W42 L=247)   377     56        0 (1625/1149-1672/1218)  63.4 W39-W43 N01-N05
      3 1866 1110 (S02W62 L=267)  1054    178        2 (1819/1069-1909/1145) -15.3 W56-W67 S04-N00
      4  702 1105 (N06E19 L=186)   216     51        0 ( 678/1074- 722/1133)  45.0 E20-E17 N04-N07
      5 1873  963 (S11W61 L=267)   724    133        2 (1846/ 923-1899/1004)  39.8 W58-W65 S13-S08
      6 1745  927 (S12W48 L=253)   431     66        2 (1717/ 899-1769/ 952) -21.8 W45-W50 S14-S10
      7  416  846 (S07E41 L=164)  1164    267        1 ( 339/ 779- 500/ 906) -45.0 E48-E35 S11-S03
      8  859  867 (S10E11 L=194)   190     46        0 ( 847/ 838- 877/ 900) -76.0 E11-E10 S11-S07
      9 1698  758 (S22W45 L=250)  1379    241        1 (1627/ 654-1774/ 870)  52.4 W39-W50 S28-S15
     10 1598  786 (S20W36 L=241)   330     53        0 (1577/ 758-1624/ 811)  33.7 W34-W38 S21-S18
     11  823  581 (S27E17 L=188)   418     60        1 ( 791/ 547- 860/ 601)   0.0 E19-E14 S29-S25
     12 1546  408 (S46W39 L=245)   889    106        1 (1502/ 372-1585/ 440)   7.1 W34-W42 S49-S43
  

    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