Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  877 1873 (N56E25 L=287)   942    102        2 ( 836/1833- 932/1897)  18.4 E30-E20 N51-N59
      2  910 1706 (N39E14 L=298)  1055    201        1 ( 848/1612- 953/1768) -47.5 E21-E09 N32-N44
      3 1109 1580 (N31W02 L=315)   331    118        1 (1090/1515-1130/1639) -77.5 W00-W03 N26-N35
      4  621 1566 (N27E32 L=280)   517    142        0 ( 557/1521- 663/1592) -24.4 E38-E27 N24-N29
      5 1450 1447 (N25W26 L=338)   280     53        0 (1432/1417-1471/1471)  53.1 W24-W27 N22-N26
      6 1204 1377 (N18W09 L=322)   245     59        1 (1175/1356-1244/1391)  26.6 W07-W11 N16-N18
      7  668  823 (S18E21 L=291)   310     87        1 ( 648/ 784- 692/ 866)  68.2 E22-E19 S20-S15
      8  460  763 (S23E37 L=275)   959    204        2 ( 389/ 691- 556/ 816)  25.0 E44-E29 S27-S19
      9  837  734 (S23E10 L=302)   228     63        0 ( 798/ 723- 872/ 748)  11.3 E12-E07 S23-S22
     10  521  525 (S38E37 L=276)   307     71        0 ( 479/ 499- 550/ 545)  26.6 E41-E33 S40-S36
     11  800  495 (S39E13 L=299)   705    101        3 ( 753/ 464- 848/ 542)  20.6 E17-E09 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