Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  701 1733 (N63E12 L=112)   344     49        0 ( 680/1709- 725/1753)  26.6 E14-E08 N60-N64
      2  477 1334 (N35E34 L= 90)   337     80        0 ( 438/1304- 513/1357) -29.7 E37-E30 N33-N37
      3 1762  952 (S14W48 L=172)  1549    367        2 (1722/ 864-1825/1060)  59.0 W44-W55 S19-S07
      4 1669 1012 (S08W41 L=165)   348     58        1 (1652/ 977-1690/1039)  56.3 W39-W42 S09-S06
      5 1894  914 (S21W65 L=189)   607     72        0 (1866/ 880-1920/ 953)  14.0 W60-W70 S22-S18
      6 1517  821 (S15W27 L=151)   229     59        0 (1502/ 789-1534/ 853)  63.4 W25-W28 S17-S13
      7 1082  654 (S16E04 L=120)   645    283        0 (1045/ 549-1175/ 714)  35.0 E08-W02 S22-S12
      8  430  542 (S13E51 L= 73)   408     75        0 ( 409/ 504- 452/ 578)  35.5 E54-E47 S15-S10
      9 1511  446 (S42W24 L=148)   337     43        0 (1491/ 424-1531/ 471)  56.3 W22-W26 S43-S39
     10 1409  371 (S45W14 L=138)   417     65        2 (1392/ 338-1435/ 407)  74.1 W12-W16 S48-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