Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  601 1746 (N30E53 L=356)  1263    298        2 ( 559/1672- 647/1799) -23.6 E60-E44 N25-N33
      2  502 1614 (N19E52 L=357)  1745    398        4 ( 411/1486- 549/1669) -24.0 E62-E44 N11-N22
      3 1046 1537 (N27E13 L= 37)  1124    302        3 ( 933/1434-1180/1641) -12.8 E23-E01 N23-N31
      4  356 1551 (N12E62 L=347)  1397    289        5 ( 317/1481- 399/1637) -79.7 E66-E58 N07-N18
      5 1463 1502 (N37W16 L= 66)   560     92        1 (1428/1463-1507/1540) -23.2 W13-W20 N35-N39
      6  381 1456 (N07E55 L=355)  1183    291        5 ( 282/1396- 466/1521)  45.0 E62-E49 N01-N12
      7 1615 1357 (N32W30 L= 80)   350     92        0 (1570/1332-1651/1372)   0.0 W26-W33 N31-N33
      8  604 1252 (S01E30 L= 19)   704    262        2 ( 503/1220- 730/1279)  14.9 E37-E22 S03-N03
      9  700  574 (S39E09 L= 40)   342     69        1 ( 664/ 550- 742/ 598) -18.4 E12-E05 S39-S37
  

    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