Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  643 1715 (N29E47 L=139)  1204    159        3 ( 586/1667- 700/1745)  29.7 E51-E43 N25-N31
      2  365 1577 (N14E63 L=124)  3432    399        4 ( 304/1504- 426/1659)  76.0 E70-E57 N08-N20
      3 1026  825 (S15W05 L=192)   310     55        3 ( 997/ 794-1060/ 851)  45.0 W03-W06 S16-S12
      4 1265  717 (S15W22 L=209)   485    136        0 (1199/ 649-1322/ 778)  74.7 W20-W23 S20-S09
      5 1620  729 (S05W45 L=232)   682    120        2 (1594/ 661-1643/ 778) -74.1 W42-W47 S08-S01
      6  735  561 (S39E06 L=180)  1109    322        2 ( 593/ 533- 903/ 586)  17.7 E17-W05 S42-S34
      7  320  531 (S53E48 L=138)   795    205        2 ( 267/ 511- 371/ 551)  14.9 E55-E40 S55-S50
      8 1460  523 (S21W42 L=229)   292     93        1 (1410/ 507-1496/ 549)  45.0 W39-W44 S22-S18
      9 1456  404 (S27W49 L=235)  2744    784        4 (1176/ 340-1656/ 447)  11.6 W27-W71 S35-S22
  

    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