Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  784 1649 (N40E28 L=123)   794    129        1 ( 728/1606- 838/1684)  -9.5 E34-E22 N36-N42
      2 1422 1580 (N42W26 L=178)   923    139        1 (1374/1531-1480/1620) -32.0 W21-W31 N39-N45
      3  366 1517 (N25E60 L= 91)  1795    171        5 ( 335/1432- 393/1597) -32.9 E71-E53 N19-N30
      4 1639 1497 (N38W47 L=198)   472     85        0 (1602/1462-1672/1527) -38.7 W43-W49 N35-N40
      5  837 1402 (N23E17 L=135)   596    203        0 ( 809/1316- 884/1496)  85.2 E17-E14 N17-N29
      6  461  835 (S16E36 L=115)  2948    551        6 ( 283/ 627- 604/1010)  44.0 E58-E26 S32-S03
      7 1107  778 (S13W08 L=159)   267     78        1 (1091/ 740-1124/ 822) -68.2 W06-W08 S15-S10
      8 1786  700 (S12W62 L=213)  1249    224        4 (1708/ 641-1838/ 780)  36.0 W54-W68 S16-S06
      9 1479  711 (S14W33 L=185)   488    159        0 (1413/ 688-1556/ 729)   5.2 W28-W39 S14-S12
     10 1419  526 (S26W34 L=185)   541    120        2 (1373/ 487-1470/ 558)  56.3 W30-W37 S29-S23
     11  574  404 (S46E33 L=118)   571     79        1 ( 541/ 378- 613/ 428)  26.6 E37-E27 S48-S43
     12 1434  385 (S37W42 L=194)  1199    114        4 (1396/ 335-1472/ 428) -66.8 W38-W47 S40-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