Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1105 1849 (N55E01 L=285)   413     68        1 (1077/1825-1126/1873)  66.8 E03-W00 N51-N58
      2  450 1620 (N30E50 L=236)   356     72        2 ( 421/1587- 479/1655)  90.0 E52-E49 N27-N33
      3  945 1566 (N29E10 L=276)   238     76        0 ( 900/1552- 985/1577)   0.0 E12-E06 N27-N29
      4 1525 1446 (N26W31 L=317)   262     82        0 (1483/1425-1563/1467) -18.4 W28-W34 N24-N26
      5  820 1139 (N01E13 L=273)   181     61        0 ( 799/1112- 844/1165)  63.4 E14-E11 N00-N02
      6  748  705 (S26E16 L=270)   223     67        0 ( 720/ 678- 783/ 726)  36.9 E17-E13 S27-S24
      7  881  501 (S39E06 L=280)   283     56        0 ( 846/ 488- 921/ 514)   0.0 E09-E03 S39-S37
      8 1131  466 (S39W13 L=299)  1020    173        3 (1061/ 426-1187/ 512) -16.7 W07-W17 S41-S35
      9 1601  355 (S41W64 L=350)  1004    143        2 (1565/ 330-1662/ 392)  21.8 W56-W70 S42-S38
  

    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