Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1600 1587 (N19W57 L= 53)  1839    320        3 (1496/1542-1726/1618) -29.1 W48-W67 N12-N22
      2 1725 1247 (S03W50 L= 46)   874    139        0 (1700/1191-1761/1322)  74.1 W47-W53 S06-N01
      3  715 1023 (N14E18 L=338)   222     80        0 ( 681/ 999- 760/1040)   0.0 E20-E14 N12-N14
      4  223  877 (N16E58 L=298)   390     45        0 ( 206/ 848- 241/ 898)  45.0 E59-E55 N14-N17
      5  261  762 (N09E58 L=298)  1632    397        3 ( 146/ 707- 366/ 810) -17.7 E72-E50 N03-N12
      6 1554  622 (S33W22 L= 18)   267     70        0 (1537/ 586-1567/ 658)  53.1 W20-W23 S35-S30
      7  938  518 (S21E20 L=336)   381    131        0 ( 908/ 466- 959/ 579)  45.0 E23-E16 S23-S17
      8  545  312 (S24E56 L=300)  2066    332        3 ( 507/ 220- 589/ 405)  59.0 E61-E50 S32-S17
      9 1054  250 (S43E26 L=330)   345     45        1 (1028/ 229-1074/ 272)   0.0 E28-E23 S45-S41
  

    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