Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  501 1735 (N30E64 L=344)  1076    105        1 ( 485/1682- 515/1790) -20.9 E79-E58 N26-N34
      2 1498 1600 (N46W22 L= 71)   302     65        0 (1470/1571-1524/1625) -26.6 W19-W25 N44-N47
      3 1648 1484 (N41W37 L= 86)   394     51        0 (1623/1449-1670/1512) -45.0 W34-W39 N38-N43
      4  525 1448 (N12E42 L=  7)   398    103        1 ( 504/1409- 541/1505) -39.8 E45-E39 N10-N15
      5 1207 1233 (N14W06 L= 55)   253     32        0 (1190/1207-1230/1258)  90.0 W05-W07 N11-N15
      6  851 1135 (N01E12 L= 37)  1037    291        1 ( 772/1024- 933/1257)  75.1 E14-E09 S06-N09
      7 1565 1190 (N19W30 L= 79)   388     98        0 (1515/1159-1609/1214)   0.0 W26-W34 N17-N19
      8 1653  789 (S02W45 L= 94)   853    141        2 (1606/ 748-1704/ 849)  60.3 W41-W48 S05-N02
      9  608  572 (S39E19 L= 30)   236     74        0 ( 578/ 544- 634/ 605)  56.3 E21-E16 S41-S35
  

    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