Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  838 1370 (N28E10 L=152)   509    242        2 ( 718/1349- 940/1387)   3.6 E18-E02 N26-N29
      2 1516 1296 (N17W36 L=198)   488    199        1 (1428/1274-1590/1361) -32.5 W30-W41 N15-N22
      3  700 1038 (N08E20 L=142)   401    107        1 ( 662/1015- 750/1072) -18.4 E22-E16 N06-N10
      4 1900 1016 (S05W69 L=232)  1080     73        1 (1859/ 991-1944/1039) -10.0 W63-W80 S07-S03
      5 1804  827 (S15W57 L=219)   580    116        1 (1764/ 790-1837/ 851) -18.4 W52-W61 S18-S13
      6  694  701 (S13E24 L=139)   427    163        0 ( 631/ 678- 746/ 745)   7.1 E27-E19 S14-S10
      7  273  704 (S12E59 L=103)   390     56        0 ( 259/ 676- 288/ 734)  76.0 E60-E57 S13-S09
      8 1337  604 (S24W18 L=180)  1833    488        2 (1180/ 438-1499/ 732)  36.5 W06-W29 S35-S16
      9 1744  634 (S28W54 L=216)   516     87        0 (1713/ 598-1777/ 663)  26.6 W50-W57 S29-S25
     10  589  333 (S39E46 L=117)   535     72        0 ( 571/ 296- 603/ 377)  35.0 E50-E40 S42-S35
     11 1060  267 (S49E05 L=157)   519     81        0 (1030/ 243-1094/ 299) -31.0 E07-E01 S51-S45
  

    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