滨海地区饮水型氟中毒和海水入侵十分普遍。 海水入侵引起地下水性质的改变对岩（土）氟迁移具有潜在影响, 但尚未有直接模拟实验证实, 其氟迁移机理也不明确。 用淡水和海（卤）水混合及配制溶液模拟海水入侵过程, 对含水砂层沉积物进行静态水-岩（土）氟迁移模拟实验, 对比分析沉积物FTIR特征, 探讨海水入侵对岩（土）氟迁移影响规律及机理。 结果表明, 含水砂层沉积物氟迁移量大小顺序为: 海水＞1∶1淡海＞淡水; 卤水＞1∶1淡卤＞淡水, 随海（卤）水混入程度增加, 含水砂层沉积物岩（土）氟迁移能力增强。 氟迁移能力随NaCl和NaHCO3浓度增加而增加, 随CaCl2浓度增加而降低。 对含水砂层沉积物傅里叶红外吸收光谱（FTIR）分析表明, 随NaCl和NaHCO3浓度增加, NaCl溶液中沉积物Si—O—Si键伸缩振动峰强度明显增强, 氟磷灰石弯曲振动峰减弱, O—H键特征峰变化不明显, 而NaHCO3溶液中沉积物羟基磷灰石弯曲振动峰和Si—O—Si键伸缩振动峰强度变化不明显, O—H键特征峰强度变化明显, 表明NaHCO3溶液主要通过OH-与F-离子交换, 而NaCl溶液是通过Si—O—Si键中Si—O置换来影响岩（土）氟迁移。 随CaCl2浓度增加, Si—O—Si键伸缩振动峰强度明显减弱, 氟磷灰石弯曲振动峰强度增强, 表明Ca2+强烈地抑制岩（土）氟析出。 同时, 随NaCl和NaHCO3浓度增加和CaCl2浓度减少, Si—F特征峰强度减弱, Si—O—Si键弯曲振动峰向低波数偏移。 由于大气CO2混入, 1 mol·L-1 CaCl2作用沉积物出现1 460和1 420 cm-1碳酸根特征峰。 由于碳酸根沉积作用, 1 mol·L-1的NaHCO3作用的沉积物在1 460 cm-1处峰强度增强, 且在875 cm-1出现新峰, 但是在淡海卤水和NaCl溶液作用的沉积物中无碳酸根特征峰, 表明含水砂层中无萤石（CaF2）溶解作用。 海水入侵引起的偏碱性、 高钠低钙的环境促进岩（土）氟迁移释放, 是滨海地区高氟地下水重要动力。

实现半色调印刷品原稿的光谱复制技术其前提要确定原稿所用油墨数目及油墨成份, 但目前应用于印刷品原稿原色油墨光谱预测的算法还有待研究, 且已有的基色色料光谱预测方法存在诸多弊端。 针对这一问题在非负矩阵分解算法基础上结合印刷品原稿光学特性, 创新的提出了一种基于约束条件非负矩阵分解的油墨光谱预测算法ISPNMF, 和对黑色油墨光谱预测结果优化的算法。 ISPNMF算法克服了基本非负矩阵分解有多重最优解和局部极小值的缺陷, 实现了预测算法唯一的全局最优解。 黑色油墨预测光谱优化的算法克服了彩色油墨对光线混合吸收给黑色油墨预测带来的干扰, 能优化得到逼近于实际黑色油墨光谱的预测值。 使用Konica Minolta C1085和HP indigo5600两台四色数码印刷机及其自身配备的墨粉和墨膏来摸拟不同品牌的油墨, 在230 g白卡纸上打印 IT8.7/3色标, 并使用X-rite i1 Pro2获取两样张的光谱反射率作为实验数据样本, 来探究并验证算法的准确性和实用性。 实验结果表明, 在印刷品原稿线性经验空间中能准确预测原稿所用原色油墨数目和油墨光谱, 且彩色油墨预测光谱与实际使用的油墨光谱相比其拟合度均高达99.9%, 光谱角距离均小于0.045, 黑色油墨的预测光谱经优化后与实际油墨光谱拟合度也高达99.9%。 这说明该算法不仅能实现对印刷品原稿原色油墨的准确预测, 而且可以精确匹配实际使用的原色油墨, 对实现印刷品原稿的光谱复制技术有重要意义。

在半色调印刷品原色油墨光谱预测技术中， 在原光谱反射率空间进行主成分分析得到的代表性基向量数目会大于实际所使用的原色油墨数目， 即光谱反射率空间并不适用光谱预测， 且主成分分析得到的基向量会出现负值， 没有物理意义。 针对上述问题， 建立了一个减色线性经验空间模型及其空间转换模型， 并探究了影响经验空间线性程度n值的因素， 通过实验及优化算法创新性找到确定最佳n值的方法， 并在该减色线性经验空间进行半色调原稿原色油墨的预测实验。 实验结果表明， 在不同n值下， 选取原始光谱反射率Rm和重构光谱反射率Rrecon的f范数平方值的最小值对应的n值做为建立线性经验空间确定最佳n值的方法是有效的； 为了将纸张和油墨类型对空间转换因数n的影响程度减到最小， 最终确定n值为3.5； 在减色线性经验空间进行印刷品原色油墨数目预测， 得到的代表性基向量数目恰好等于实际印刷使用的原色油墨数目4， 进行原色油墨光谱预测， 预测的4色油墨除K色外， 其他CMY色与实际原色油墨光谱相比拟合度GFC均大于99.9%。 即所提出的优化n值的新方法建立的减色线性经验空间是一个可作为半色调印刷品原色油墨数目预测和光谱预测的有效线性空间。

针对直接在光谱反射率空间, 对原稿颜色样本光谱的主成分分析会导致特征向量的数目超过真实物理维度(原稿所用基色色料)的数量, 以及特征向量和对应系数存在负值, 不能直接表示原稿基色色料的光谱特性和对应浓度等情况。 创新性的提出需根据原稿色料的光学特性建立一个完全线性的光谱空间, 并在该空间中使用带约束条件的非负矩阵分解实现对原稿基色数量和光谱形状进行预测的方法。 对此, 首先设计了一个实现对原稿基色色料光谱预测方法的总体研究方案和实现步骤, 再以透明色料原稿为例, 研究如何选择和构建一个符合其光学特性的光谱线性空间, 然后再在基本非负矩阵分解(BNMF)基础上提出针对基色色料光谱预测的有约束非负矩阵分解算法(SCNMF)。 针对BNMF算法会出现多重最优解, 为了提高预测精度, 使矩阵分解结果有明确的物理意义, 所提出的SCNMF算法需要满足4个约束条件: 非负性约束; 全加性约束; 平滑性约束; 稀疏性约束。 建立了满足约束条件的目标函数和迭代算法。 预测结果表明本文提出的新方法能有效的实现对原稿基色物理维度和基色色料光谱的准确预测。

针对光谱反射率空间并不是一个有效的基色色料线性混合空间的缺点，提出了建立基色色料线性混合空间的三项原则。根据不同原稿色料光学特性的不同，将原稿分为透明色料连续调原稿、不透明色料连续调原稿、半色调原稿(半透明色料)，并针对这三类原稿分别构建能真实反应原稿色料真实物理维度的光谱线性空间，以不透明色料原稿为例进行实验验证。实验结果表明，根据原稿色料光学特性建立的色料线性混合空间是一个有效的线性混合空间，在其上的光谱预测能真实反映原稿基色色料的数量和近似得到基色色料的光谱形状。

针对传统光谱降维方法其降维重构后的光谱数据仅是对原始光谱的数学逼近, 会出现光谱误差较小但颜色色差较大的缺点, 创新性的提出三种将人眼视觉特性与光谱降维相结合的方法。 其中, VPCA法直接将光谱光视效率函数加权到原始光谱上再进行降维, LMSPCA方法用LMS视稚响应构建加权矩阵对原始光谱加权后再进行降维, 在LMSPCA法中加权矩阵的构建有两种方式, 其主要区别在于视稚响应偏差的求取方式不同。 方式一中, L, M, S视稚响应偏差是各对应波长上的偏差取绝对值, 而方式二中, 其偏差是各对应波长上的偏差平方。 LMSPCAs法在LMSPCA法基础上再采用PCA(主成分分析)方法对损失的光谱进行降维。 实验结果表明VPCA法降维效果较差, LMSPCA法的两种加权矩阵降维效果接近, 皆可显著提高降维模型的色度精度, 但会降低模型的光谱精度, LMSPCAs法由于针对LMSPCA法因光谱加权引起的光谱损失再进行光谱补偿, 其在光谱精度、 色度精度以及变光照条件下的色差稳定性这三个方面都能较好地表征原始高维光谱反射率, 满足光谱颜色复制的要求。

针对传统光谱降维方法的缺点，提出一种基于矩阵理论的线性光谱降维方法。该方法在Cohen定义的矩阵R理论基础上提出修正矩阵，可将物体光谱反射率分解为基本光谱和同色异谱黑光谱，再分别对两者进行线性降维。基本光谱的三个基向量由颜色匹配函数与光源的乘积再进行正交化得到，其三个基向量的累积贡献率可达100%。同色异谱黑光谱的基向量由主成分分析推导，取前三个特征向量作为同色异谱黑光谱的基向量。实验结果表明，使用该方法构建的六维线性模型能在光谱和色度两个方面较好地表征原始高维光谱反射率，满足光谱颜色复制的要求。