Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  786 1708 (N34E30 L=254)   940    210        1 ( 736/1601- 819/1781) -42.9 E38-E23 N26-N39
      2  478 1635 (N25E51 L=232)   643     66        3 ( 440/1616- 520/1660)   0.0 E55-E46 N24-N26
      3 1334 1636 (N38W12 L=295)   262     71        0 (1304/1611-1370/1659) -31.0 W09-W14 N36-N39
      4 1254 1595 (N33W06 L=290)   308     94        0 (1227/1547-1277/1640) -50.2 W03-W08 N30-N36
      5  378 1536 (N17E55 L=229)   968    140        3 ( 345/1476- 402/1593) -63.4 E59-E51 N13-N21
      6  173  993 (S17E63 L=220)  3970    376        4 ( 120/ 854- 268/1164)  55.6 E73-E53 S25-S06
      7  592  846 (S22E25 L=258)  1628    489        4 ( 410/ 755- 794/ 989)  34.7 E39-E13 S29-S11
      8  888  829 (S19E05 L=278)   657    299        0 ( 829/ 760- 981/ 872)  12.5 E09-W00 S23-S16
      9 1451  768 (S14W31 L=314)   302     60        1 (1435/ 738-1468/ 805) -45.0 W29-W32 S15-S12
     10 1057  709 (S24W07 L=291)   205     40        1 (1037/ 687-1084/ 732)   0.0 W05-W08 S25-S22
     11  273  710 (S34E56 L=228)   337     72        0 ( 236/ 694- 309/ 726)   0.0 E60-E51 S34-S32
     12 1541  584 (S23W43 L=327)  1650    411        1 (1387/ 477-1672/ 726)  59.0 W36-W48 S32-S12
     13  750  282 (S61E11 L=272)  1098    143        4 ( 667/ 250- 813/ 310)   0.0 E22-E02 S63-S58
  

    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