Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1688 1478 (N12W58 L=340)  1090    198        3 (1615/1406-1759/1527) -57.5 W53-W61 N05-N16
      2  490 1426 (N45E28 L=254)   318     68        0 ( 463/1397- 512/1460)  36.9 E30-E25 N42-N46
      3 1677 1426 (N10W54 L=336)   436     97        2 (1644/1381-1709/1455) -78.7 W51-W56 N06-N11
      4 1693 1185 (S03W46 L=328)   462    129        1 (1642/1162-1738/1203) -29.7 W41-W48 S05-S01
      5  915 1091 (N13E04 L=277)  1197    441        3 ( 747/1049-1117/1116) -15.3 E15-W07 N08-N17
      6  344  860 (N13E48 L=234)  2839    553        7 ( 183/ 784- 535/ 940) -23.2 E62-E34 N06-N19
      7 1536  865 (S16W26 L=308)   475     71        1 (1507/ 835-1575/ 897) -59.0 W23-W28 S18-S13
      8 1398  666 (S24W13 L=294)   737    195        1 (1338/ 581-1433/ 758)  48.4 W06-W16 S27-S18
      9  796  644 (S09E23 L=259)   202     76        0 ( 773/ 610- 827/ 675)  31.0 E25-E20 S10-S07
     10  780  561 (S13E27 L=255)   267     56        1 ( 756/ 542- 818/ 585) -63.4 E28-E24 S15-S11
     11  431  353 (S20E64 L=218)  2940    293        1 ( 355/ 242- 537/ 451) -87.0 E67-E61 S30-S11
  

    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