Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1685 1404 (N08W53 L=314)  1017    207        4 (1647/1328-1711/1476)  53.1 W49-W57 N03-N11
      2  510  914 (N09E33 L=229)   389    110        1 ( 469/ 866- 543/ 948) -71.6 E35-E31 N05-N11
      3  653  815 (N00E26 L=235)   476    120        1 ( 585/ 793- 707/ 838) -14.0 E30-E22 S01-N02
      4  278  772 (N05E56 L=206)   994    153        2 ( 231/ 742- 331/ 797) -26.6 E60-E50 N03-N07
      5 1262  751 (S17W07 L=269)   200     55        1 (1229/ 736-1292/ 769)   0.0 W04-W09 S18-S16
      6 1888  743 (S37W58 L=319)   846    190        1 (1863/ 727-1924/ 769) -29.1 W54-W63 S39-S34
      7 1319  723 (S21W10 L=272)   198     50        0 (1294/ 707-1348/ 743)   0.0 W08-W12 S21-S19
      8  402  578 (S08E52 L=210)   425     71        0 ( 377/ 542- 424/ 603) -63.4 E54-E50 S10-S06
      9  450  437 (S18E55 L=207)  4091    692       16 ( 341/ 353- 597/ 525) -51.8 E65-E43 S25-S11
     10  716  313 (S32E42 L=220)   338     65        1 ( 698/ 284- 735/ 342)  33.7 E44-E38 S33-S29
  

    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