Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  513 1580 (N15E51 L= 17)   521    104        0 ( 490/1540- 549/1621) -24.0 E54-E45 N13-N17
      2  317 1490 (N06E64 L=  3)   489     84        1 ( 291/1451- 333/1535) -63.4 E66-E62 N03-N09
      3  372 1506 (N08E59 L=  8)   437     55        0 ( 350/1480- 400/1524)  26.6 E61-E56 N06-N09
      4 1473 1428 (N32W18 L= 85)   348     89        1 (1439/1391-1500/1457) -18.4 W14-W20 N30-N34
      5  142 1077 (S21E71 L=356)   621     76        1 ( 128/1043- 162/1113)  45.0 E73-E67 S23-S18
      6  408  759 (S36E36 L= 31)   687    209        2 ( 329/ 717- 494/ 790)  20.6 E44-E28 S39-S32
      7  680  628 (S37E11 L= 56)  1814    356        6 ( 563/ 539- 808/ 747)  56.3 E19-E04 S46-S25
      8 1701  584 (S09W60 L=127)  2054    271        6 (1589/ 502-1801/ 653)  52.4 W53-W72 S17-S03
      9 1484  558 (S17W43 L=110)   303     94        1 (1451/ 524-1522/ 597)  78.7 W41-W45 S19-S14
     10 1389  419 (S28W45 L=112)  1007    208        1 (1302/ 366-1479/ 463)  61.4 W41-W49 S33-S22
     11 1243  361 (S36W39 L=106)   461     72        0 (1218/ 333-1274/ 395) -23.2 W35-W42 S37-S34
  

    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