Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  879 1864 (N43E47 L=184)   503     38        1 ( 859/1839- 895/1887) -51.3 E50-E43 N40-N45
      2 1241 1738 (N45E11 L=220)  1613    205        2 (1176/1668-1317/1787)  57.5 E14-E05 N37-N51
      3 1576 1398 (N33W24 L=255)   730    115        2 (1552/1340-1598/1460) -60.3 W21-W26 N29-N36
      4 1756 1344 (N37W42 L=273)  1548    292        1 (1687/1199-1814/1437) -38.2 W34-W49 N29-N41
      5 1336 1414 (N25W07 L=238)   222     68        0 (1309/1396-1374/1434)   0.0 W04-W09 N24-N26
      6  798  769 (S26E06 L=225)   400    135        0 ( 741/ 735- 860/ 801)  45.0 E08-E02 S29-S22
      7 1240  682 (S18W23 L=254)   639    266        1 (1149/ 648-1327/ 752)  14.0 W15-W27 S20-S13
      8 1004  465 (S38W18 L=250)   593    112        1 ( 970/ 418-1051/ 509) -83.7 W16-W20 S42-S33
      9 1146  339 (S41W36 L=267)   331     55        0 (1116/ 323-1173/ 361)  90.0 W34-W37 S42-S38
  

    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