Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1111 1573 (N38W01 L=298)   312     72        0 (1084/1551-1155/1593)  26.6 W00-W04 N36-N39
      2 1255 1526 (N36W13 L=309)   253     57        0 (1229/1504-1285/1540)  11.3 W10-W15 N34-N37
      3 1715 1484 (N37W53 L=349)  1090    199        2 (1641/1438-1782/1519) -22.6 W45-W58 N33-N38
      4  546 1445 (N24E38 L=259)   346     39        0 ( 523/1422- 572/1470)  56.3 E39-E36 N22-N25
      5  651 1438 (N25E29 L=267)   453     98        0 ( 604/1414- 694/1461) -15.9 E33-E26 N23-N25
      6 1616 1317 (N24W40 L=336)  1640    283        4 (1511/1229-1716/1386) -29.4 W31-W47 N19-N28
      7  245  986 (S07E56 L=241)  1496    212        2 ( 216/ 874- 269/1091) -69.0 E60-E53 S14-N00
      8  761  733 (S18E15 L=282)   292     45        1 ( 733/ 710- 790/ 757)  26.6 E16-E12 S20-S16
      9 1804  657 (S16W64 L=  1)  1086    212        4 (1763/ 627-1841/ 702)  23.2 W59-W70 S17-S12
  

    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