public RasterOp(RasterOp op)
{
_raster = new float?[op.Width, op.Height];
for (int i = 0; i < Width; i++)
{
for (int j = 0; j < Height; j++)
{
_raster[i, j] = op.Read(i, j);
}
}
}
public RasterOp(RasterOp op)
{
_raster = new float?[op.Width, op.Height];
for (int i = 0; i < Width; i++)
{
for (int j = 0; j < Height; j++)
{
_raster[i, j] = op.Read(i, j);
}
}
}
public void Overlay(RasterOp raster, Func <float, float, float> func)
{
if (Width != raster.Width || Height != raster.Height)
{
throw new ArgumentException("栅格大小不一致，无法叠加操作");
}
for (int i = 0; i < Width; i++)
{
for (int j = 0; j < Height; j++)
{
if (Read(i, j).HasValue&&raster.Read(i, j).HasValue)
{
float basic = Read(i, j).GetValueOrDefault();
float overlay = raster.Read(i, j).GetValueOrDefault();
Write(i, j, func(basic, overlay));
}
}
}
}
/// <summary>
/// 计算一个高铁城市的在高铁未通车的时间成本
/// </summary>
/// <returns>时间成本操作</returns>
private RasterOp CalculationCity(string cityName)
{
RasterOp rasterOp=new RasterOp(_rasterReader);
City city = _allCities.First(item => item.Name == cityName);
Postion pos = _rasterReader.Coordinate(city.XCoord,city.YCoord);
return rasterOp.Calculator(pos);
}
/// <summary>
/// 计算一个城市对其中的影响
/// </summary>
/// <param name="basic"></param>
/// <param name="city"></param>
/// <param name="value"></param>
private void Overlay(RasterOp basic, RasterOp city,float value)
{
for (int i = 0; i < basic.Width; i++)
{
for (int j = 0; j < basic.Height; j++)
{
if (basic.Read(i,j).HasValue)
{
var timeCost = (float)city.Read(i, j);
basic.Write(i,j,(float)basic.Read(i,j)+timeCost*value);
}
}
}
}
//private void Wirte(RasterOp op, string foldePath, string name)
//{
// RasterWriter writer = new RasterWriter(foldePath, name,Info);
// op.WriteRaster(writer,"TIFF");
//}
/// <summary>
///
/// </summary>
/// <param name="basic"></param>
/// <param name="rasterFolderPath"></param>
/// <param name="fileName"></param>
/// <param name="cityName"></param>
private void Overlay(RasterOp basic, string rasterFolderPath,string fileName,string cityName)
{
RasterReader reader = new RasterReader(rasterFolderPath, fileName + ".tif");
RasterOp op=new RasterOp(reader);
Overlay(basic, op, CityValues[cityName]);
}
private void Overlay(RasterOp basic, RasterOp city,string cityName)
{
for (int i = 0; i < basic.Width; i++)
{
for (int j = 0; j < basic.Height; j++)
{
if (basic.Read(i, j).HasValue)
{
var timecost = (float) city.Read(i, j);
float shorestTimeCost = ShortestTimeCost(i, j, _cityPos[cityName]);
if (Math.Abs(shorestTimeCost) < 10e-5)
{
continue;
}
basic.Write(i, j, (float)basic.Read(i, j) + timecost/shorestTimeCost * CityValues[cityName]); }
}
}
}
private void Run(RasterReader reader, CalculatorCity city, string folderPath)
{
RasterOp rasterOp=new RasterOp(reader);
Postion pos = reader.Coordinate(city.XCoord,city.YCoord);
var result = rasterOp.Calculator(pos);
RasterWriter writer=new RasterWriter(folderPath,city.Name+"_高铁未通车",reader.RasterInfo);
result.WriteRaster(writer, "TIFF");
}
private void Run(object p)
{
Hashtable para=p as Hashtable;
ContinuousWait wait = para["wait"] as ContinuousWait;
var folderPath = para["folderPath"].ToString();
try
{
string preWorkSpace = Path.GetDirectoryName(_preRasterFilePath);
string preFileName = Path.GetFileNameWithoutExtension(_preRasterFilePath);
string aftWorkSpace = Path.GetDirectoryName(_aftRasterFilePath);
string aftFileName = Path.GetFileNameWithoutExtension(_aftRasterFilePath);
RasterReader preReader = new RasterReader(preWorkSpace, preFileName + ".tif");
RasterReader aftReader = new RasterReader(aftWorkSpace, aftFileName + ".tif");
RasterOp preOp = new RasterOp(preReader);
RasterOp aftOp = new RasterOp(aftReader);
RasterOp res = preOp.Clone();
res.Reset();
for (int i = 0; i < preOp.Width; i++)
{
for (int j = 0; j < preOp.Height; j++)
{
if (preOp.Read(i, j).HasValue)
{
float orgin = (float) preOp.Read(i, j);
float now = (float) aftOp.Read(i, j);
if (Math.Abs(orgin) > 10e-5)
{
float rate = (now - orgin)/orgin;
res.Write(i, j, rate);
}
}
}
}
RasterWriter writer=new RasterWriter(folderPath,RasterName,preReader.RasterInfo);
res.WriteRaster(writer,"TIFF");
para["ret"] = true;
}
catch (Exception e)
{
_log.Error(e.GetType()+e.Message+e.StackTrace);
para["ret"] = false;
}
finally
{
wait.CloseWait();
}
}
private void Overlay(RasterOp basic, RasterOp city, string cityName)
{
for (int i = 0; i < basic.Width; i++)
{
for (int j = 0; j < basic.Height; j++)
{
if (basic.Read(i, j).HasValue)
{
float timeCost = (float)city.Read(i, j);
if (Math.Abs(timeCost) < 10e-5)
{
if (_popDic.ContainsKey(cityName))
{
float tt = (float)(3*Math.Log10(_popDic[cityName]*10));
basic.Write(i,j,(float)basic.Read(i,j)+CityValues[cityName]/tt);
}
}
else
{
basic.Write(i, j, (float)basic.Read(i, j) + CityValues[cityName] / timeCost);
}
}
}
}
}
public void Overlay(RasterOp raster,Func<float,float,float> func)
{
if (Width!=raster.Width||Height!=raster.Height)
{
throw new ArgumentException("栅格大小不一致，无法叠加操作");
}
for (int i = 0; i < Width; i++)
{
for (int j = 0; j < Height; j++)
{
if (Read(i,j).HasValue&&raster.Read(i,j).HasValue)
{
float basic =Read(i, j).GetValueOrDefault();
float overlay =raster.Read(i, j).GetValueOrDefault();
Write(i,j,func(basic,overlay));
}
}
}
}
protected void Wirte(RasterOp op, string foldePath, string name)
{
RasterWriter writer = new RasterWriter(foldePath, name, Info);
op.WriteRaster(writer, "TIFF");
}
private void Init()
{
var firstCity = Cities.First(item => item.IsSelected);
RasterReader reader = new RasterReader(HighTrainNoFolderPath, firstCity.Name + "_高铁未通车" + ".tif");
Info = reader.RasterInfo;
RasterOp = new RasterOp(reader);
RasterOp.Reset();
}
private bool WriteTimeCostRaster(string folderPath)
{
try
{
//读取两个文件
var landReader = new RasterReader(folderPath, LandUse + ".tif");
var roadReader = new RasterReader(folderPath, TrafficRoad + ".tif");
RasterOp landOp = new RasterOp(landReader);
RasterOp roadOp = new RasterOp(roadReader);
//取最小值
landOp.Overlay(roadOp,RasterAlgorithm.MinPixel);
////写入文件
RasterWriter writer = new RasterWriter(folderPath, TimeCostName, landReader.RasterInfo);
landOp.WriteRaster(writer,"TIFF");
return true;
}
catch (Exception e)
{
_log.Error(e.Message + e.StackTrace);
}
return false;
}
