Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  531 1742 (N26E59 L=329)  1109    108        3 ( 491/1689- 558/1790) -48.8 E65-E54 N21-N29
      2  417 1685 (N20E66 L=322)  1126    105        1 ( 366/1651- 452/1713)  45.0 E71-E62 N17-N22
      3 1136 1639 (N33E10 L= 19)   445    117        0 (1118/1585-1151/1702) -58.0 E12-E07 N28-N37
      4 1552 1649 (N49W22 L= 51)   365     49        0 (1535/1616-1570/1677) -90.0 W20-W24 N45-N51
      5 1608 1423 (N32W28 L= 57)  1629    299        1 (1519/1303-1692/1527) -32.0 W20-W36 N26-N38
      6 1700 1218 (N22W38 L= 67)   347     53        0 (1667/1194-1728/1241)   0.0 W35-W40 N21-N23
      7  476  954 (S23E33 L=356)  1009    215        4 ( 395/ 900- 580/1001)  29.7 E39-E25 S27-S18
      8  890  714 (S28E00 L= 29)   303    170        0 ( 818/ 686- 973/ 735)  32.0 E04-W04 S31-S24
      9  565  676 (S39E22 L=  6)   461    133        1 ( 526/ 624- 611/ 723)  53.1 E25-E19 S42-S34
     10  386  444 (S61E50 L=338)   696    112        0 ( 362/ 406- 420/ 494)  74.1 E57-E44 S63-S56
  

    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