Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1249 1582 (N41W14 L=337)   447     71        1 (1204/1560-1287/1600)  15.9 W10-W17 N38-N42
      2  514 1409 (N23E39 L=283)   310     53        1 ( 483/1392- 552/1429)  -8.1 E42-E35 N22-N24
      3 1154 1282 (N19W07 L=330)   559    241        0 (1058/1214-1241/1329) -26.6 W00-W12 N15-N21
      4  712 1232 (N13E21 L=301)   556    150        0 ( 652/1203- 771/1282) -14.0 E25-E17 N11-N16
      5 1651 1041 (N06W42 L=  5)   439     60        1 (1626/1009-1670/1071)  63.4 W40-W43 N04-N08
      6  683  722 (S19E21 L=302)   479    159        0 ( 647/ 659- 717/ 787)  58.0 E23-E18 S22-S14
      7 1613  704 (S15W43 L=  6)   354     59        1 (1589/ 673-1637/ 738)  76.0 W41-W44 S16-S12
      8 1709  649 (S18W54 L= 17)   352     32        0 (1692/ 628-1727/ 669) -26.6 W51-W56 S19-S16
      9  880  640 (S23E07 L=315)   344    141        0 ( 818/ 625- 937/ 661)   0.0 E11-E03 S23-S21
     10  341  524 (S36E57 L=265)  1072    133        2 ( 286/ 502- 404/ 553)   0.0 E66-E48 S37-S33
  

    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