Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  761 1561 (N36E13 L= 61)   371     72        0 ( 727/1523- 785/1582) -51.3 E15-E11 N32-N37
      2  260 1356 (N27E56 L= 18)  2963    365        4 ( 224/1203- 321/1540) -43.8 E77-E47 N17-N41
      3 1366 1452 (N22W27 L=101)   817    206        0 (1265/1426-1474/1475) -19.7 W20-W34 N19-N24
      4 1865 1273 (N06W67 L=141)  2944    215        5 (1789/1212-1936/1356) -33.7 W58-W78 N02-N12
      5 1740  910 (S14W48 L=122)   320     23        0 (1717/ 888-1760/ 927) -14.0 W46-W50 S15-S13
      6  263  621 (S18E64 L= 10)   498     44        0 ( 247/ 594- 283/ 644)  15.9 E67-E59 S19-S16
      7  932  618 (S25E10 L= 64)   285     49        0 ( 903/ 601- 965/ 641) -33.7 E11-E08 S26-S23
      8 1325  447 (S41W17 L= 91)   689    100        2 (1275/ 422-1380/ 476)   0.0 W12-W21 S43-S38
      9 1219  414 (S43W08 L= 82)   439     60        0 (1193/ 392-1251/ 446) -59.0 W05-W10 S44-S39
     10 1131  382 (S44W00 L= 74)   419     58        1 (1104/ 363-1164/ 412)   0.0 E01-W03 S45-S41
  

    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