Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1689 (N41W42 L=292)   524    156        1 (1253/1648-1365/1708)   0.0 W35-W49 N38-N42
      2 1465 1671 (N35W54 L=304)   938    149        2 (1400/1635-1533/1704) -31.0 W49-W59 N31-N37
      3  542 1638 (N58E32 L=218)   531    128        2 ( 490/1595- 583/1671)  21.8 E38-E26 N55-N60
      4 1248 1463 (N28W25 L=275)   241     92        1 (1214/1433-1285/1482)   0.0 W21-W27 N26-N29
      5  340 1344 (N39E48 L=201)   816    142        2 ( 291/1277- 376/1377) -38.7 E55-E43 N33-N41
      6  625  901 (N07E27 L=223)   271     76        1 ( 608/ 863- 650/ 934)  56.3 E28-E25 N05-N09
      7  489  874 (N07E37 L=213)   316     77        2 ( 450/ 849- 522/ 902)   9.5 E39-E33 N06-N08
      8 1715  644 (S33W43 L=293)   476     77        1 (1693/ 617-1741/ 676) -45.0 W40-W46 S36-S30
      9 1557  545 (S36W27 L=277)   343     94        0 (1523/ 512-1585/ 584)  38.7 W24-W30 S37-S33
     10  433  556 (S11E51 L=199)   618     95        0 ( 388/ 533- 482/ 579)   0.0 E55-E45 S12-S10
     11  547  435 (S21E47 L=202)   934    151        2 ( 491/ 382- 626/ 474) -10.3 E51-E39 S24-S17
  

    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