/// <summary>
/// 光谱基线校正
/// </summary>
/// <param name="xyDataList">第一项是X，后面是Y</param>
/// <param name="method">基线校正方法.</param>
/// <returns></returns>
public static List <double[]> BaselineCorrect(List <double[]> xyDataList, BaseLineMethod method)
{
if (xyDataList.Count < 2)
{
ErrorString = "Invalid parameters";
return(null);
}
var xyDatas = CommonMethod.CombineSpectrumDatas(xyDataList);
if (xyDatas == null)
{
ErrorString = CommonMethod.ErrorString;
return(null);
}
if (BaselineCorrect(xyDatas, xyDataList.Count, xyDataList[0].Length, method) == false)
{
return(null);
}
return(CommonMethod.SplitSpectrumDatas(xyDatas, xyDataList.Count));
}
/// <summary>
/// AdvancedIdentify定性分析(光谱的数据点数量必须相同, X轴也必须相同)
/// </summary>
/// <param name="apiList">API原始光谱插值后，一列X轴，一列Y轴</param>
/// <param name="inferList">干扰组分原光谱插值后，一列X轴，多列Y轴</param>
/// <param name="sampleList">待分析物原始光谱插值后，一列X轴，一列Y轴</param>
/// <param name="coefMethod">相关方法，0为标准，1为增强</param>
/// <param name="baseline">基线校正方法</param>
/// <param name="derivatives">导数平滑方法(有基线校正，忽略导数平滑)</param>
/// <param name="smoothPoints">平滑窗口大小(5-99)</param>
/// <param name="polyDegree">导数的阶数,2,3,4</param>
/// <param name="firstX">起始计算X轴，-1表示使用API光谱的firstX</param>
/// <param name="lastX">结束计算X轴，-1表示使用API光谱的lastX</param>
/// <returns></returns>
public static double Identify(List <double[]> apiList, List <double[]> inferList, List <double[]> sampleList, CoefficientMethod coefMethod,
BaseLineMethod baseline = BaseLineMethod.None, DerivativeMethod derivatives = DerivativeMethod.None,
int smoothPoints = 7, int polyDegree = 2, double firstX = double.MinValue, double lastX = double.MaxValue)
{
if (apiList.Count != 2 || inferList.Count < 2 || sampleList.Count != 2)
{
ErrorString = "Invalid parameters";
return(0);
}
if (smoothPoints > 99 || smoothPoints < 2 || polyDegree < 2 || polyDegree > 4)
{
ErrorString = "Invalid parameters";
return(0);
}
//需要截取一部分光谱区间
if (firstX != double.MinValue && lastX != double.MaxValue)
{
apiList = CommonMethod.GetRangeData(apiList, firstX, lastX);
inferList = CommonMethod.GetRangeData(inferList, firstX, lastX);
sampleList = CommonMethod.GetRangeData(sampleList, firstX, lastX);
}
firstX = apiList[0][0];
lastX = apiList[0][apiList[0].Length - 1];
int dataCount = apiList[0].Length;
//格式不相同，全部兼容到api光谱
if (!CommonMethod.IsSameXDatas(inferList[0], firstX, lastX, dataCount))
{
inferList = Ai.Hong.Algorithm.PreProcessor.SplineCubicInterpolation(inferList, firstX, lastX, dataCount);
}
if (!CommonMethod.IsSameXDatas(sampleList[0], firstX, lastX, dataCount))
{
sampleList = Ai.Hong.Algorithm.PreProcessor.SplineCubicInterpolation(sampleList, firstX, lastX, dataCount);
}
//限定2种计算系数的方法，1=无权重, 2=强度
coefMethod = coefMethod != CoefficientMethod.Normal ? CoefficientMethod.Weight : coefMethod;
//拷贝为数组，好在Ai.Hong.Algorithm中处理
var srcApiDatas = Ai.Hong.Algorithm.CommonMethod.CombineSpectrumDatas(apiList);
var apiDatas = Ai.Hong.Algorithm.CommonMethod.CombineSpectrumDatas(apiList);
var inferDatas = Ai.Hong.Algorithm.CommonMethod.CombineSpectrumDatas(inferList);
var sampleDatas = Ai.Hong.Algorithm.CommonMethod.CombineSpectrumDatas(sampleList);
if (apiDatas == null || inferDatas == null || sampleDatas == null)
{
ErrorString = CommonMethod.ErrorString;
return(0);
}
//需要基线校正
if (baseline != BaseLineMethod.None)
{
BaselineCorrect(apiDatas, apiList.Count, dataCount, baseline);
BaselineCorrect(inferDatas, inferList.Count, dataCount, baseline);
BaselineCorrect(sampleDatas, sampleList.Count, dataCount, baseline);
}
else if (derivatives != DerivativeMethod.None) //导数平滑
{
coefMethod = coefMethod == CoefficientMethod.Normal ? coefMethod : CoefficientMethod.SquareWeight; //纪南的程序里面，导数后，增强方法计算系数的方式不一样（使用强度平方）
PreProcessor.SGSmoothAndDerivative(apiDatas, apiList.Count, dataCount, smoothPoints, polyDegree, (int)derivatives);
PreProcessor.SGSmoothAndDerivative(inferDatas, inferList.Count, dataCount, smoothPoints, polyDegree, (int)derivatives);
PreProcessor.SGSmoothAndDerivative(sampleDatas, sampleList.Count, dataCount, smoothPoints, polyDegree, (int)derivatives);
//纪南的导数没有除以X轴步长
double stepx = apiDatas[1] - apiDatas[0];
for (int i = dataCount; i < apiDatas.Length; i++)
{
apiDatas[i] *= stepx;
}
stepx = inferDatas[1] - inferDatas[0];
for (int i = dataCount; i < inferDatas.Length; i++)
{
inferDatas[i] *= stepx;
}
stepx = sampleDatas[1] - sampleDatas[0];
for (int i = dataCount; i < sampleDatas.Length; i++)
{
sampleDatas[i] *= stepx;
}
}
return(Identify(srcApiDatas, apiDatas, inferDatas, inferList.Count, sampleDatas, coefMethod));
}
/// <summary>
/// 光谱基线校正
/// </summary>
/// <param name="spectrumDatas">第一行为X数据，后面几行为Y数据</param>
/// <param name="specRows">数据的行数</param>
/// <param name="specCols">数据的列数</param>
/// <param name="method">基线校正方法</param>
/// <returns></returns>
public static bool BaselineCorrect(double[] spectrumDatas, int specRows, int specCols, BaseLineMethod method)
{
if (spectrumDatas == null || spectrumDatas.Length != specRows * specCols)
{
ErrorString = "Invalid parameters";
return(false);
}
if (method == BaseLineMethod.None)
{
return(true);
}
try
{
if (CommonMethod.Is64BitVersion())
{
return(ADBaselineCorrect64(spectrumDatas, specRows, specCols, (int)method));
}
else
{
return(ADBaselineCorrect32(spectrumDatas, specRows, specCols, (int)method));
}
}
catch (Exception ex)
{
ErrorString = ex.Message;
return(false);
}
}
