Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  471 1609 (N49E36 L=235)   351     62        0 ( 439/1595- 499/1629)   0.0 E39-E32 N48-N50
      2 1168 1491 (N26W18 L=289)   716    163        1 (1112/1438-1247/1531)  15.3 W13-W24 N23-N27
      3  456 1435 (N36E34 L=237)   692    175        2 ( 383/1404- 532/1461) -11.3 E41-E26 N33-N39
      4  734 1111 (N11E16 L=255)   195     55        1 ( 710/1094- 765/1129) -33.7 E16-E13 N09-N12
      5 1851  841 (S25W58 L=329)   780    155        0 (1814/ 770-1886/ 915)  45.0 W54-W62 S28-S20
      6 1615  594 (S36W35 L=306)   269     44        1 (1585/ 574-1634/ 613) -14.0 W32-W36 S37-S35
      7  353  492 (S19E61 L=209)  2755    287        6 ( 254/ 388- 454/ 582) -59.0 E68-E54 S26-S11
      8  891  522 (S26E18 L=252)   259     98        0 ( 853/ 487- 926/ 565)  26.6 E21-E15 S28-S24
  

    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