title: How to output a raster layer calculated from a probability function using 2 other raster layers
body: <p>We have two raster layers with the same resolution and want to an output layer of the same resolution. However, the equation is quite complex. The layers contain information about wind shape and scale (weibull k and A) within a cell.  First, we must find the wind speed probability using the probability function:</p>
  <p>P(S| k, A) = (k/A) * ((S/A)^(k-1)) *e^(-(S/A)^k)</p>
  <p>S = wind speed</p>
  <p>k = shape</p>
  <p>A = scale</p>
  <p>This distribution of wind speeds is thus the wind speed on the x axis and the probability of the y axis. This is then multiplied by the power curve (per hour) for a given turbine model. The power curve tells how much power is generated (y axis) at a given wind speed (x axis). By multiplying the probability distribution of wind speeds by the power curve (* 8760 hour per year), the result should be a probability (y axis) distribution of expected power output (x axis). By summing the expected power outputs, the result is the estimated annual energy output per cell. This output needs to be a raster file for GIS. We are using QGIS but information pertaining the ArcGIS also helpful.</p>
  <p>The question is: What is the best was to perform this series of equations? In Python or R, what packages would be useful here? We have no idea where to start.</p>