Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1589 1638 (N26W59 L=184)  1111    185        3 (1535/1564-1632/1698)  78.7 W56-W62 N20-N30
      2 1723 1379 (N08W57 L=181)   333     84        2 (1701/1353-1739/1408)   9.5 W53-W59 N06-N09
      3 1786 1146 (S07W53 L=178)  1134    272        2 (1705/1025-1825/1229) -74.1 W48-W56 S14-N00
      4 1639  651 (S32W32 L=157)  1371    251        1 (1527/ 553-1732/ 736)  18.4 W21-W40 S35-S28
      5 1122  619 (S19E04 L=121)   237    115        0 (1064/ 607-1180/ 631) -23.2 E07-W00 S20-S17
      6  811  534 (S17E25 L=100)   435    128        2 ( 776/ 493- 851/ 570)  18.4 E27-E21 S19-S14
      7  708  542 (S14E31 L= 94)   228     35        0 ( 685/ 524- 733/ 556)  18.4 E32-E29 S15-S13
      8  928  474 (S23E20 L=105)   554    131        0 ( 872/ 442- 996/ 525) -58.0 E22-E16 S26-S18
      9  681  438 (S20E37 L= 87)   755    187        2 ( 614/ 364- 723/ 517) -90.0 E39-E35 S26-S13
     10  433  435 (S16E57 L= 67)   439     75        1 ( 404/ 409- 454/ 458) -45.0 E59-E55 S18-S14
     11  522  398 (S20E52 L= 73)  1608    341        1 ( 474/ 329- 585/ 447)  36.0 E58-E45 S25-S16
  

    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