Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  764 1645 (N26E36 L=281)  1408    379        6 ( 669/1495- 858/1787) -28.1 E53-E23 N18-N34
      2 1050 1739 (N41E23 L=294)   352     43        1 (1021/1721-1074/1760)  26.6 E24-E20 N39-N43
      3 1093 1653 (N36E15 L=302)   868    136        0 (1056/1580-1134/1715) -38.7 E20-E09 N31-N40
      4 1570 1669 (N57W27 L=343)  1206    137        1 (1517/1634-1635/1693)  29.5 W19-W42 N50-N63
      5  267  950 (S26E52 L=265)  1422    294        1 ( 196/ 832- 370/1022)  34.4 E61-E41 S35-S20
      6  665  693 (S32E14 L=303)   707    232        1 ( 567/ 675- 782/ 712)  21.8 E21-E06 S35-S28
      7 1184  527 (S28W25 L=342)   819    152        1 (1140/ 458-1216/ 598) -66.0 W22-W26 S32-S23
      8  517  385 (S61E26 L=290)  2310    313        4 ( 399/ 310- 622/ 461) -12.3 E44-E10 S65-S56
      9 1323  412 (S31W40 L=357)   808    281        1 (1203/ 374-1428/ 450)  11.9 W31-W50 S33-S27
  

    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