Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1615 1649 (N32W58 L=158)   530     59        0 (1581/1623-1643/1676) -45.0 W55-W61 N30-N34
      2 1477 1368 (N15W33 L=133)   442    119        0 (1461/1315-1491/1436)  69.4 W31-W35 N11-N19
      3  952 1385 (N22E01 L= 99)   477    173        1 ( 867/1373-1045/1400) -16.7 E05-W05 N20-N23
      4  232 1093 (N12E55 L= 45)   940    151        4 ( 187/1031- 276/1155) -45.0 E60-E50 N07-N16
      5 1346  734 (S22W17 L=117)   369    103        0 (1311/ 698-1385/ 774)  38.7 W14-W20 S23-S19
      6  925  708 (S19E10 L= 90)   192     51        0 ( 906/ 681- 947/ 733)  56.3 E11-E07 S20-S17
      7  888  649 (S22E13 L= 87)   491    194        0 ( 806/ 617- 961/ 686)  -5.7 E17-E07 S24-S19
      8  581  545 (S25E37 L= 63)   572    114        1 ( 531/ 530- 643/ 566) -20.6 E40-E32 S26-S23
      9  993  386 (S41E11 L= 89)   973    165        0 ( 920/ 350-1063/ 422)   9.5 E17-E05 S43-S39
  

    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