Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  874 1652 (N29E28 L=233)  1500    343        5 ( 766/1518- 949/1726) -23.2 E39-E18 N22-N33
      2  676 1624 (N22E40 L=221)  1317    234        3 ( 595/1539- 747/1708)  85.9 E42-E37 N15-N29
      3 1367 1392 (N26W11 L=272)  1871    437        5 (1257/1187-1424/1554) -49.6 W00-W18 N14-N34
      4  321 1040 (S19E46 L=215)   438    116        0 ( 271/1007- 373/1065)  36.9 E50-E41 S22-S16
      5 1669 1010 (N12W40 L=301)   263     74        0 (1657/ 973-1679/1051) -78.7 W39-W40 N09-N14
      6  319  901 (S28E45 L=216)   418     90        0 ( 274/ 890- 374/ 916)  11.3 E49-E39 S28-S26
      7 1178  855 (S09W13 L=274)   347    147        0 (1111/ 815-1234/ 884)   0.0 W08-W17 S10-S08
      8  647  435 (S51E10 L=251)  1653    301        1 ( 527/ 360- 764/ 525) -11.3 E23-W02 S53-S47
  

    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