Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1079 1752 (N39E19 L= 62)   527     90        3 (1033/1727-1117/1771)   9.5 E22-E16 N36-N41
      2 1249 1686 (N39E05 L= 76)   975    180        4 (1207/1641-1308/1724) -24.0 E09-W00 N34-N42
      3 1501 1642 (N46W16 L= 97)   449     57        1 (1472/1619-1531/1667)  -9.5 W12-W18 N43-N47
      4  877 1540 (N19E22 L= 58)   681    193        0 ( 822/1470- 941/1620) -26.6 E28-E16 N15-N22
      5 1337 1409 (N23W08 L= 89)   691    312        4 (1297/1253-1371/1522) -59.0 W03-W13 N14-N29
      6 1754 1228 (N25W42 L=122)  1410    260        2 (1658/1140-1834/1326) -15.5 W32-W50 N22-N28
      7  238  943 (S29E52 L= 29)   713    177        4 ( 187/ 905- 291/ 996)  16.7 E57-E47 S31-S25
      8  713  459 (S50E05 L= 76)   580    137        2 ( 656/ 428- 789/ 501)   4.8 E10-W02 S50-S46
      9  864  207 (S65W27 L=108)   600     97        2 ( 839/ 172- 891/ 233) -22.6 W21-W33 S68-S62
  

    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