Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1022 1544 (N35W17 L=147)   724    244        2 ( 931/1518-1151/1563) -31.6 W11-W24 N30-N38
      2 1632  943 (S16W34 L=164)   790    249        0 (1556/ 877-1697/1058)  19.7 W27-W41 S18-S11
      3 1208  928 (S04W07 L=138)   527    103        1 (1176/ 876-1246/ 983) -82.9 W06-W08 S08-N00
      4  259  892 (N16E55 L= 75)   449     63        2 ( 229/ 866- 287/ 926) -31.0 E58-E52 N14-N18
      5 1398  852 (S14W16 L=147)   251     76        0 (1370/ 819-1417/ 891)  53.1 W14-W18 S15-S11
      6 1620  770 (S26W30 L=160)   598    204        3 (1586/ 702-1655/ 852)  66.0 W27-W33 S30-S21
      7 1422  583 (S31W12 L=142)   257     44        0 (1403/ 562-1442/ 607)  33.7 W10-W13 S32-S29
      8 1299  491 (S33W01 L=131)  1234    262        3 (1245/ 397-1353/ 589)  77.9 E05-W05 S39-S25
  

    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