Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  507 1608 (N18E51 L=162)   424     97        1 ( 489/1570- 523/1648) -59.0 E53-E49 N15-N20
      2  233 1404 (N00E69 L=144)   604    106        3 ( 203/1386- 257/1429)   0.0 E73-E65 N00-N02
      3  937  926 (S11E03 L=210)   632    245        0 ( 852/ 845-1049/ 981)  51.3 E05-W03 S18-S07
      4 1730  926 (N09W47 L=260)   408     90        2 (1708/ 888-1748/ 973) -45.0 W45-W49 N07-N11
      5  604  789 (S28E20 L=193)   353     59        1 ( 580/ 756- 634/ 824) -36.9 E22-E17 S29-S25
      6  468  680 (S38E29 L=183)   894    302        1 ( 362/ 660- 578/ 703)  25.3 E39-E20 S42-S33
      7  828  593 (S34E00 L=213)   784    213        0 ( 719/ 550- 920/ 661)  -7.1 E09-W07 S35-S32
      8 1289  525 (S25W31 L=244)  1633    459        1 (1063/ 455-1459/ 582)   9.5 W15-W45 S29-S22
      9 1686  462 (S17W65 L=278)  1456    383        3 (1618/ 433-1750/ 493)   6.0 W56-W75 S18-S14
  

    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