Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1303 1765 (N49E03 L= 51)   492     89        0 (1266/1726-1330/1788) -26.6 E06-W00 N45-N51
      2  815 1706 (N29E33 L= 20)   353     68        0 ( 795/1674- 837/1740) -38.7 E35-E30 N26-N31
      3 1618 1648 (N52W31 L= 85)   422     18        0 (1596/1629-1640/1666)  29.7 W27-W34 N49-N54
      4 1305 1559 (N32W04 L= 57)   228     37        0 (1287/1538-1326/1580) -33.7 W02-W05 N30-N33
      5  958 1402 (N13E13 L= 41)   323    144        2 ( 943/1335- 973/1466) -58.0 E15-E10 N08-N16
      6 1687 1174 (N19W38 L= 91)   288    132        0 (1665/1139-1717/1238) -51.3 W35-W40 N17-N22
      7  202 1025 (S23E60 L=353)  2366    255        3 ( 129/ 924- 289/1090)  24.4 E72-E50 S29-S18
      8 1368  387 (S34W45 L= 99)   525    102        1 (1333/ 367-1418/ 422)   0.0 W40-W49 S35-S32
      9 1163  135 (S57W66 L=119)   776     84        0 (1148/  91-1178/ 174)  -9.2 W56-W87 S59-S53
  

    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