Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1611 1460 (N28W49 L=301)  1373    240        4 (1547/1391-1687/1557) -67.4 W44-W54 N22-N34
      2  356 1419 (N32E53 L=199)   696    120        1 ( 295/1392- 405/1446)  -7.1 E62-E46 N29-N33
      3  256 1362 (N28E64 L=189)   442     50        1 ( 240/1334- 281/1385)  45.0 E66-E60 N25-N29
      4  318 1325 (N25E54 L=198)   801    103        0 ( 267/1308- 379/1348)   4.8 E60-E47 N24-N26
      5  672 1236 (N19E22 L=230)   226    105        0 ( 632/1207- 705/1284)  51.3 E25-E20 N17-N22
      6 1667  879 (S09W43 L=296)   714    180        2 (1636/ 808-1691/ 953)  60.9 W40-W45 S13-S04
      7 1748  856 (S11W51 L=303)   482     60        1 (1726/ 822-1775/ 883) -76.0 W48-W53 S13-S09
      8  591  545 (S24E34 L=219)   701    185        0 ( 566/ 455- 628/ 626)  81.3 E35-E31 S31-S18
      9  392  550 (S24E52 L=201)   434     74        1 ( 365/ 526- 423/ 580)  26.6 E55-E47 S25-S21
  

    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