Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1617 1664 (N32W61 L= 30)   561    201        0 (1590/1627-1635/1700)  14.0 W56-W65 N30-N34
      2 1402 1537 (N30W37 L=  6)  1118    168        3 (1317/1485-1494/1583)  -3.6 W30-W46 N27-N32
      3  479 1481 (N44E36 L=293)   394     93        0 ( 440/1461- 519/1494)   7.1 E39-E31 N42-N44
      4  807 1200 (N20E11 L=318)   253    118        0 ( 751/1171- 865/1231) -35.5 E14-E07 N17-N22
      5 1543  917 (S09W31 L=359)  1296    278        5 (1485/ 804-1580/1065)  57.3 W25-W35 S14-N00
      6 1801  964 (S12W53 L= 21)  1332    292        3 (1722/ 827-1854/1065) -67.6 W47-W58 S21-S04
      7  220  922 (N10E60 L=269)  1094    217        1 ( 169/ 875- 273/ 991)  23.2 E67-E53 N07-N13
      8 1673  922 (S12W41 L=  9)   527    108        0 (1652/ 859-1698/ 975)  60.3 W38-W43 S15-S08
      9 1286  423 (S37W07 L=336)  1465    332        2 (1164/ 322-1357/ 554)  43.2 E03-W14 S43-S28
  

    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