private void Erode(int width, int height, float seaLevel, Overworld.MapData[,] heightMap, int numRains, int rainLength, int numRainSamples, float[,] buffer)
{
float remaining = 1.0f - Progress.Value - 0.2f;
float orig = Progress.Value;
for (int x = 0; x < width; x++)
{
for (int y = 0; y < height; y++)
{
buffer[x, y] = heightMap[x, y].Height;
}
}
for (int i = 0; i < numRains; i++)
{
LoadingMessage = "Erosion " + i + "/" + numRains;
Progress.Value = orig + remaining * ((float)i / (float)numRains);
Vector2 currentPos = new Vector2(0, 0);
Vector2 bestPos = currentPos;
float bestHeight = 0.0f;
for (int k = 0; k < numRainSamples; k++)
{
int randX = PlayState.Random.Next(1, width - 1);
int randY = PlayState.Random.Next(1, height - 1);
currentPos = new Vector2(randX, randY);
float h = Overworld.GetHeight(buffer, currentPos);
if (h > bestHeight)
{
bestHeight = h;
bestPos = currentPos;
}
}
currentPos = bestPos;
const float erosionRate = 0.99f;
Vector2 velocity = Vector2.Zero;
for (int j = 0; j < rainLength; j++)
{
Vector2 g = Overworld.GetMinNeighbor(buffer, currentPos);
float h = Overworld.GetHeight(buffer, currentPos);
if (h < seaLevel || g.LengthSquared() < 1e-12)
{
break;
}
Overworld.MinBlend(Overworld.Map, currentPos, erosionRate * Overworld.GetValue(Overworld.Map, currentPos, Overworld.ScalarFieldType.Erosion), Overworld.ScalarFieldType.Erosion);
velocity = 0.1f * g + 0.7f * velocity + 0.2f * MathFunctions.RandVector2Circle();
currentPos += velocity;
}
}
}
private void Weather(int width, int height, float T, Vector2[] neighbs, float[,] buffer)
{
for (int x = 0; x < width; x++)
{
for (int y = 0; y < height; y++)
{
buffer[x, y] = Overworld.Map[x, y].Height * Overworld.Map[x, y].Faults;
}
}
int weatheringIters = 10;
for (int iter = 0; iter < weatheringIters; iter++)
{
for (int x = 0; x < width; x++)
{
for (int y = 0; y < height; y++)
{
Vector2 p = new Vector2(x, y);
Vector2 maxDiffNeigh = Vector2.Zero;
float maxDiff = 0;
float totalDiff = 0;
float h = Overworld.GetHeight(buffer, p);
float lowestNeighbor = 0.0f;
for (int i = 0; i < 4; i++)
{
float nh = Overworld.GetHeight(buffer, p + neighbs[i]);
float diff = h - nh;
totalDiff += diff;
if (diff > maxDiff)
{
maxDiffNeigh = neighbs[i];
maxDiff = diff;
lowestNeighbor = nh;
}
}
if (maxDiff > T)
{
Overworld.AddValue(Overworld.Map, p + maxDiffNeigh, Overworld.ScalarFieldType.Weathering, (float)(maxDiff * 0.4f));
Overworld.AddValue(Overworld.Map, p, Overworld.ScalarFieldType.Weathering, (float)(-maxDiff * 0.4f));
}
}
}
for (int x = 0; x < width; x++)
{
for (int y = 0; y < height; y++)
{
Vector2 p = new Vector2(x, y);
float w = Overworld.GetValue(Overworld.Map, p, Overworld.ScalarFieldType.Weathering);
Overworld.AddHeight(buffer, p, w);
Overworld.Map[x, y].Weathering = 0.0f;
}
}
}
for (int x = 0; x < width; x++)
{
for (int y = 0; y < height; y++)
{
Overworld.Map[x, y].Weathering = buffer[x, y] - Overworld.Map[x, y].Height * Overworld.Map[x, y].Faults;
}
}
}
