Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1385 1635 (N46W24 L=150)   387     86        1 (1356/1608-1416/1655)  26.6 W20-W26 N43-N48
      2  538 1488 (N25E40 L= 87)  2308    531        6 ( 407/1354- 627/1579) -23.7 E54-E29 N17-N30
      3 1222 1511 (N35W09 L=136)   407    112        2 (1187/1467-1255/1561) -56.3 W06-W11 N31-N38
      4 1669 1337 (N27W44 L=170)   297     77        1 (1646/1319-1705/1358) -21.8 W41-W47 N25-N28
      5 1716 1195 (N18W47 L=174)   867    131        3 (1673/1141-1745/1248) -60.3 W43-W49 N14-N21
      6  934 1129 (N07E06 L=120)   200     67        1 ( 904/1105- 956/1165)  63.4 E08-E05 N05-N09
      7  430 1095 (S00E40 L= 87)  1879    378        3 ( 327/1024- 587/1157)  11.3 E47-E27 S05-N03
      8  575  725 (S22E27 L= 99)   227     64        0 ( 544/ 712- 607/ 741)  11.3 E29-E24 S22-S20
      9  388  674 (S27E43 L= 84)   337    120        0 ( 358/ 644- 420/ 715)  39.8 E46-E39 S29-S24
     10  720  557 (S32E17 L=110)   304     57        2 ( 693/ 538- 750/ 584)  26.6 E18-E14 S33-S29
  

    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