Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  857 1825 (N62W25 L=317)   751    113        2 ( 804/1802- 922/1844) -18.4 W18-W32 N58-N65
      2  707 1767 (N64W04 L=296)   449     62        0 ( 688/1732- 729/1797)  29.1 W00-W09 N60-N66
      3 1175 1737 (N42W38 L=331)   879    199        3 (1095/1683-1266/1779)  -6.3 W29-W47 N39-N43
      4  386 1603 (N61E45 L=247)  1218    190        1 ( 311/1550- 463/1651)   7.1 E59-E31 N57-N63
      5 1513 1596 (N23W51 L=343)   570    203        0 (1455/1536-1578/1635) -69.4 W45-W53 N17-N25
      6  525 1177 (N27E27 L=265)   839    302        2 ( 394/1122- 649/1226) -31.4 E37-E18 N20-N32
      7 1808  943 (S21W49 L=341)  2987    454        6 (1694/ 780-1899/1115)  25.7 W36-W63 S27-S14
      8  283  906 (N16E52 L=241)   473    151        1 ( 243/ 875- 349/ 960) -41.2 E54-E46 N13-N20
      9 1371  874 (S11W15 L=308)   350     82        0 (1357/ 830-1386/ 918)  78.7 W13-W16 S13-S08
     10 1426  567 (S32W12 L=304)   423     82        1 (1392/ 545-1477/ 587) -18.4 W09-W15 S33-S29
     11 1196  404 (S35E08 L=284)  1358    274        1 (1082/ 348-1320/ 473)   3.0 E17-W02 S37-S32
     12  819  170 (S42E51 L=242)  1872    169        6 ( 743/ 126- 894/ 236)  11.3 E63-E38 S46-S37
  

    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