Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  763 1769 (N31E45 L=256)   944    199        0 ( 716/1681- 817/1840) -61.9 E49-E40 N23-N38
      2  984 1772 (N38E30 L=270)   362     43        1 ( 956/1754-1012/1794)  45.0 E32-E28 N36-N40
      3 1006 1694 (N33E24 L=277)   563    124        1 ( 985/1625-1025/1749) -45.0 E27-E19 N28-N36
      4  528 1662 (N18E55 L=246)   492     83        0 ( 500/1624- 559/1701) -71.6 E56-E52 N14-N21
      5 1065 1352 (N13E07 L=293)  1332    319        3 (1037/1210-1090/1519) -57.1 E13-E02 N05-N22
      6  327 1124 (S18E48 L=253)   597    189        0 ( 275/1083- 377/1193)  54.5 E52-E45 S20-S13
      7  180  907 (S34E63 L=237)   720    119        1 ( 141/ 865- 226/ 954)  30.3 E71-E56 S37-S30
      8 1557  524 (S16W51 L=351)   368     55        1 (1533/ 500-1579/ 550) -21.8 W47-W52 S18-S15
      9 1614  426 (S19W63 L=  4)  2947    273       10 (1528/ 356-1725/ 493)  55.3 W56-W71 S25-S12
     10 1369  444 (S27W43 L=343)  1552    326        2 (1222/ 418-1513/ 482)  32.9 W32-W50 S32-S20
     11  746  242 (S67W17 L=317)  1510     67        1 ( 706/ 191- 791/ 294) -38.7 W09-W27 S72-S60
     12  921  254 (S56W29 L=330)   714     85        1 ( 897/ 212- 946/ 292) -32.5 W24-W35 S60-S53
  

    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