Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  394 1531 (N31E49 L=160)   384     29        2 ( 371/1507- 417/1552)  90.0 E51-E46 N28-N32
      2 1528 1438 (N22W35 L=244)   302    102        0 (1480/1405-1567/1469) -33.7 W31-W37 N20-N24
      3 1692 1324 (N15W47 L=256)   345     59        1 (1667/1301-1719/1350) -26.6 W44-W48 N14-N17
      4  180 1201 (N10E63 L=146)  4020    283        9 ( 131/1059- 235/1319) -82.9 E71-E56 N01-N17
      5  486 1295 (N14E35 L=174)   657    213        0 ( 400/1261- 583/1315) -12.1 E41-E27 N11-N15
      6  681  831 (S14E22 L=187)   899    275        0 ( 634/ 713- 712/ 962)  71.6 E26-E19 S21-S06
      7 1268  864 (S13W15 L=224)   198     46        1 (1240/ 845-1292/ 891) -33.7 W13-W16 S13-S11
      8 1239  791 (S17W13 L=222)   288     99        1 (1193/ 780-1291/ 809)  -9.5 W10-W16 S17-S16
      9 1651  449 (S39W56 L=265)   815     99        2 (1593/ 431-1703/ 472)   4.1 W48-W64 S40-S37
     10  682  457 (S39E28 L=181)   302     80        0 ( 643/ 446- 732/ 468)   0.0 E31-E23 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