Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  519 1623 (N58E23 L= 97)   406     66        0 ( 492/1605- 557/1641)  -7.1 E27-E18 N56-N60
      2 1123 1598 (N35W26 L=146)  1151    228        1 (1016/1566-1232/1624) -36.9 W19-W32 N30-N39
      3  289 1313 (N42E49 L= 71)   495     68        1 ( 255/1283- 318/1342) -35.5 E55-E45 N39-N44
      4 1689  998 (S15W40 L=161)  1920    481        3 (1561/ 838-1774/1171)  33.0 W26-W48 S20-S05
      5 1350 1000 (S04W17 L=138)   400    123        0 (1322/ 944-1385/1053) -90.0 W16-W18 S08-N00
      6  219  791 (N11E63 L= 58)   915    139        1 ( 194/ 725- 246/ 846)  82.9 E65-E60 N06-N14
      7 1421  530 (S35W11 L=131)  1134    210        1 (1379/ 444-1452/ 633)  69.0 W05-W13 S40-S27
      8  551  326 (S24E57 L= 63)   410     45        0 ( 523/ 313- 581/ 341)  -9.5 E60-E54 S24-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