Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1583 1579 (N48W32 L=147)   359     54        1 (1555/1565-1615/1598)   8.1 W28-W35 N46-N49
      2  712 1466 (N17E30 L= 85)   776    155        4 ( 644/1418- 772/1517)   0.0 E35-E25 N14-N20
      3  648 1406 (N12E32 L= 83)   300     86        1 ( 608/1369- 675/1436)  90.0 E35-E30 N09-N14
      4  998 1373 (N18E09 L=106)   727    185        3 ( 926/1332-1088/1414)  29.1 E13-E04 N14-N20
      5  700  748 (S26E14 L=101)   894    293        1 ( 597/ 693- 800/ 827)  47.5 E20-E08 S31-S18
      6 1525  451 (S24W52 L=167)   685     87        2 (1495/ 420-1558/ 489)  78.7 W49-W56 S26-S21
      7  598  446 (S50E19 L= 96)   336     43        1 ( 572/ 432- 631/ 467)   9.5 E21-E15 S51-S47
      8  830  383 (S48W04 L=119)  1093    217        1 ( 759/ 358- 922/ 434)  17.1 E02-W11 S51-S43
      9 1383  324 (S36W50 L=165)   433     64        0 (1354/ 301-1417/ 347)  78.7 W48-W52 S38-S33
  

    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