private string MAxValue(string[] fileNames, int bandNo, float zoom)
{
foreach (string f in fileNames)
{
if (!File.Exists(f))
{
PrintInfo("所选择的数据:\"" + f + "\"不存在。");
return(null);
}
}
Dictionary <string, FilePrdMap> filePrdMap = new Dictionary <string, FilePrdMap>();
for (int i = 0; i < fileNames.Length; i++)
{
FilePrdMap map = new FilePrdMap(fileNames[i], zoom, new VaildPra(float.MinValue, float.MaxValue), new int[] { bandNo });
if (map.BandCount < 1)
{
PrintInfo("请选择正确的数据进行最大值合成。");
return(null);
}
filePrdMap.Add("ndviFile" + i.ToString(), map);
}
ITryCreateVirtualPrd tryVPrd = new TryCreateVirtualPrdByMultiFile();
IVirtualRasterDataProvider prd = tryVPrd.CreateVirtualRasterPRD(ref filePrdMap);
if (prd == null)
{
throw new Exception("数据间无相交部分,无法创建虚拟数据提供者!");
}
IRasterPixelsVisitor <float> visitor = new RasterPixelsVisitor <float>(new ArgumentProvider(prd, null));
List <int> bands = new List <int>();
for (int i = 0; i < fileNames.Length; i++)
{
bands.Add(filePrdMap["ndviFile" + i.ToString()].StartBand);
}
PrintInfo("正在进行最大值合成,请稍后...!");
IPixelFeatureMapper <float> result = new MemPixelFeatureMapper <float>("0MAX", prd.Width * prd.Height, new Size(prd.Width, prd.Height), prd.CoordEnvelope, prd.SpatialRef);
visitor.VisitPixel(bands.ToArray(),
(idx, values) =>
{
result.Put(idx, values.Max() * zoom);
}
);
RasterIdentify rid = new RasterIdentify(fileNames);
IInterestedRaster <float> iir = new InterestedRaster <float>(rid, new Size(prd.Width, prd.Height), prd.CoordEnvelope, prd.SpatialRef);
iir.Put(result);
iir.Dispose();
return(iir.FileName);
}
public override IExtractResult Make(Action <int, string> progressTracker, IContextMessage contextMessage)
{
_contextMessage = contextMessage;
if (_argumentProvider == null)
{
return(null);
}
if (_argumentProvider.GetArg("AlgorithmName") == null)
{
PrintInfo("参数\"AlgorithmName\"为空。");
return(null);
}
string algorith = _argumentProvider.GetArg("AlgorithmName").ToString();
if (algorith != "0MAX")
{
PrintInfo("指定的算法\"" + algorith + "\"没有实现。");
return(null);
}
string[] fileNames = GetStringArray("SelectedPrimaryFiles");
if (fileNames == null || fileNames.Count() == 0)
{
PrintInfo("请选择参与最大值合成的数据！");
return(null);
}
foreach (string f in fileNames)
{
if (!File.Exists(f))
{
PrintInfo("所选择的数据:\"" + f + "\"不存在。");
return(null);
}
}
Dictionary <string, FilePrdMap> filePrdMap = new Dictionary <string, FilePrdMap>();
for (int i = 0; i < fileNames.Length; i++)
{
FilePrdMap map = new FilePrdMap(fileNames[i], 1, new VaildPra(Int16.MinValue, Int16.MaxValue), new int[] { 1 });
if (map.BandCount < 1)
{
PrintInfo("请选择正确的数据进行最大值合成。");
return(null);
}
filePrdMap.Add("ndviFile" + i.ToString(), map);
}
ITryCreateVirtualPrd tryVPrd = new TryCreateVirtualPrdByMultiFile();
IVirtualRasterDataProvider prd = tryVPrd.CreateVirtualRasterPRD(ref filePrdMap);
if (prd == null)
{
throw new Exception("数据间无相交部分,无法创建虚拟数据提供者!");
}
IRasterPixelsVisitor <float> visitor = new RasterPixelsVisitor <float>(new ArgumentProvider(prd, null));
List <int> bands = new List <int>();
for (int i = 0; i < fileNames.Length; i++)
{
bands.Add(filePrdMap["ndviFile" + i.ToString()].StartBand);
}
Size size = new Size(prd.Width, prd.Height);
IPixelFeatureMapper <Int16> result = new MemPixelFeatureMapper <Int16>("0MAX", prd.Width * prd.Height, size, prd.CoordEnvelope, prd.SpatialRef);
visitor.VisitPixel(bands.ToArray(),
(idx, values) =>
{
result.Put(idx, (Int16)values.Max());
}
);
result.SetDispaly(false);
return(result);
}
