Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1383 1507 (N22W37 L=288)   528    102        3 (1338/1469-1425/1546)   0.0 W33-W42 N20-N23
      2 1176 1507 (N29W25 L=275)   676    143        1 (1104/1485-1251/1530) -36.9 W21-W29 N25-N32
      3  917 1491 (N36W08 L=258)   262     76        0 ( 898/1456- 944/1521)  21.8 W05-W10 N33-N37
      4  125 1056 (N28E73 L=177)   916     66        0 (  86/1040- 164/1079)  -6.0 E85-E65 N27-N29
      5 1723  884 (S22W40 L=291)   287     40        1 (1700/ 859-1741/ 907)  45.0 W37-W42 S23-S20
      6  475  692 (N01E42 L=208)   812    192        0 ( 443/ 636- 534/ 760) -74.1 E44-E38 S02-N05
      7 1210  366 (S39E09 L=241)   545    139        2 (1154/ 332-1272/ 410)   5.7 E13-E03 S40-S37
      8  613  247 (S30E60 L=190)   981     99        3 ( 579/ 205- 652/ 278)  71.6 E65-E55 S34-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