Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1399 1694 (N38W50 L=299)   944    214        3 (1292/1665-1484/1715) -25.0 W41-W58 N34-N41
      2  542 1636 (N58E31 L=218)   592    107        1 ( 501/1591- 586/1671)  24.0 E36-E26 N54-N60
      3 1254 1464 (N28W26 L=275)   252     73        2 (1221/1434-1288/1484)   9.5 W22-W28 N26-N29
      4  346 1339 (N38E48 L=201)   865    160        0 ( 295/1278- 385/1378) -30.3 E54-E42 N33-N41
      5  637  903 (N07E26 L=223)   273     71        2 ( 617/ 869- 660/ 936)  56.3 E27-E24 N05-N08
      6  500  878 (N08E36 L=213)   295     64        0 ( 467/ 851- 530/ 899)  21.8 E38-E33 N06-N08
      7 1725  798 (S23W43 L=292)   828    164        3 (1677/ 737-1770/ 864)  49.4 W38-W48 S25-S18
      8 1712  648 (S33W42 L=291)   457     45        2 (1685/ 627-1733/ 676) -51.3 W39-W45 S35-S30
      9  445  561 (S11E49 L=200)   722    111        1 ( 396/ 533- 496/ 595)   0.0 E54-E44 S12-S08
     10  543  448 (S20E47 L=202)   964    186        2 ( 475/ 417- 635/ 498) -14.0 E51-E38 S22-S15
  

    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