Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1338 1769 (N55W39 L=163)   529    103        0 (1287/1750-1384/1790)   7.1 W32-W48 N53-N57
      2 1301 1551 (N37W23 L=147)   373     88        1 (1270/1508-1331/1578) -68.2 W21-W25 N34-N39
      3  559 1424 (N30E34 L= 90)  1976    424        2 ( 410/1285- 692/1535) -27.3 E52-E21 N20-N37
      4 1050 1262 (N19W02 L=126)   499    221        0 (1004/1184-1132/1311)  45.0 W00-W07 N13-N21
      5  731  973 (N01E18 L=106)  1511    447        2 ( 645/ 816- 863/1103)  48.6 E24-E09 S08-N09
      6 1462  912 (S04W28 L=152)   267     68        0 (1443/ 886-1489/ 945) -45.0 W26-W29 S05-S02
      7  608  636 (S20E29 L= 95)  1421    312        3 ( 520/ 505- 686/ 731)  45.0 E39-E22 S29-S13
      8 1392  517 (S30W26 L=150)  3793    755        8 (1088/ 357-1596/ 730)  34.4 W03-W39 S42-S16
      9  928  669 (S18E07 L=117)   242     52        1 ( 909/ 648- 952/ 693)  45.0 E08-E05 S19-S16
     10 1703  642 (S23W51 L=175)   407     66        0 (1690/ 605-1715/ 678) -78.7 W49-W52 S25-S20
     11  376  632 (S21E49 L= 75)   361     63        2 ( 356/ 600- 399/ 660)  38.7 E52-E46 S23-S19
     12  981  487 (S31E04 L=120)  1421    479        2 ( 845/ 384-1196/ 637)  35.8 E15-W10 S39-S20
  

    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