Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1417 1705 (N35W47 L=347)   588     98        0 (1369/1675-1462/1730)   4.8 W41-W53 N33-N36
      2 1597 1661 (N28W61 L=  0)  1300    102        0 (1533/1641-1660/1680)  -7.1 W54-W70 N26-N29
      3 1560 1506 (N19W46 L=346)   286     45        0 (1536/1489-1588/1524) -45.0 W44-W48 N17-N20
      4  716 1331 (N25E14 L=286)  2336    594        3 ( 453/1258- 986/1393)  -7.7 E33-W04 N21-N29
      5 1795 1120 (S08W54 L=353)   495     42        1 (1773/1089-1821/1143) -63.4 W51-W56 S09-S05
      6 1442  823 (S18W22 L=321)   472    203        0 (1395/ 749-1470/ 917)  52.1 W17-W24 S21-S12
      7 1637  748 (S27W36 L=335)   262    130        0 (1610/ 694-1656/ 814)  66.8 W33-W37 S29-S22
      8 1277  717 (S21W10 L=309)   253     61        0 (1255/ 687-1295/ 757)  53.1 W07-W11 S22-S18
      9  521  389 (S28E52 L=247)  1057    193        0 ( 459/ 344- 582/ 462)  18.4 E60-E42 S30-S23
  

    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