Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  458 1661 (N58E53 L= 22)  1516    183        3 ( 364/1636- 525/1690)   2.8 E85-E42 N55-N60
      2 1704 1447 (N18W59 L=135)   477     99        0 (1658/1410-1740/1477) -50.2 W56-W62 N15-N21
      3  718 1432 (N37E15 L= 60)   523    153        4 ( 645/1414- 795/1448)  -9.5 E21-E09 N34-N38
      4  371 1061 (N17E44 L= 31)   564    132        1 ( 339/ 997- 402/1108) -66.8 E47-E41 N12-N19
      5 1459  688 (S22W23 L= 98)   768    242        1 (1387/ 605-1526/ 780)  39.3 W17-W28 S25-S16
      6 1172  603 (S22W03 L= 78)   675    288        0 (1080/ 534-1262/ 715)  31.6 E04-W09 S25-S16
      7 1715  651 (S30W46 L=121)   414     58        3 (1700/ 619-1737/ 681) -68.2 W43-W48 S32-S27
      8  685  581 (S16E29 L= 46)   461    136        1 ( 659/ 527- 729/ 622)  35.5 E32-E25 S19-S13
      9 1194  483 (S31W02 L= 77)   929    345        2 (1062/ 419-1337/ 560)  13.4 E08-W13 S33-S27
     10 1644  533 (S38W40 L=116)   375     56        0 (1614/ 519-1673/ 552)   9.5 W37-W43 S38-S36
     11 1422  452 (S38W20 L= 95)   581    188        1 (1356/ 402-1501/ 507)  21.0 W13-W26 S41-S35
     12  786  444 (S27E27 L= 48)   559     92        3 ( 765/ 399- 811/ 491)  81.9 E28-E24 S30-S23
  

    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