旋光效应可视检测液体的浓度（附MATLAB代码）

最新推荐文章于 2023-10-24 22:38:00 发布

liuz_notes

最新推荐文章于 2023-10-24 22:38:00 发布

阅读量3.8k

点赞数 5

分类专栏：杂项文章标签：信号处理算法边缘检测 matlab

本文链接：https://blog.csdn.net/liuz_notes/article/details/103477550

版权

杂项专栏收录该内容

11 篇文章 9 订阅

订阅专栏

目的

根据CCD相机采集到的旋光楔形图像，标定当前基于旋光可视检测的液体浓度光电检测系统实验系统参数；
然后由测得的比旋光度测定未知浓度葡萄糖溶液的浓度。

原理

一、旋光效应
当一束单一的平面偏振光通过某些晶体或物质的溶液时，其振动方向发生改变，此时光的振动面将以光的传播方向为轴线会旋转一定的角度，这种现象称为旋光现象，旋转的角度称为旋光度。物质的这种使偏振光的振动面旋转的性质叫做旋光性，具有旋光性的物质叫做旋光性物质或旋光物质。

二、角向偏振片显示原理
角向偏振显示器又称角向偏振片，其截面图如图1（a）所示。角向偏振片可以看作是一个线栅偏振，即由多个线型偏振器件构成的具有圆形同心传播轴的偏振器件。角向偏振片的每一个小区域都是一个线型偏振器，其传输轴的方向垂直于径向。当线偏振光通过角向偏振显示器之后，观察成像结果发现，在线偏振光的偏振方向上会出现两个相对的楔形暗区如图1(b)，对该图像进行处理，可以获得线偏振光的振动方向。

（a）

（b）

图1 (a)角向偏振显示器截面; (b)线偏振光经过角片器成像仿真图

三、系统原理
物质的旋光度实验光路如图4.2所示，由单色光源（钠光灯或激光器）发出的光，通过起偏镜（尼可尔棱镜）后，转变为平面偏振光（简称偏振光）。当偏振光通过样品管中的旋光物质时，由于旋光性，其振动平面会旋转一定角度，此时的偏振光在入射到角向偏振片后，成像到CCD相机上。拍摄到的图像的楔形暗区方向即为经过旋光物质作用后光束的偏振方向，与未经过旋光物质作用的光束的偏振方向的夹角即为样品的实测旋光度 $\theta$ 。
在这里插入图片描述

图2 实验光路

旋光度的大小除了取决于被测分子的立体结构外，还受到待测溶液的浓度、偏振光通过溶液的厚度（即样品管的长度）以及温度、所用光源的波长、所用溶剂等因素的影响，这些因素在测定结果中都要表示出来。常用比旋光度来表示物质的旋光性，比旋光度和旋光度的关系如下：
$\theta=\alpha\cdot c \cdot L$

上式中, $\theta$ 为所测得的旋光度；c为液层厚度为L (cm)的液体的浓度(g/ml)； $\alpha$ 为旋光物质的比旋光度，即单位浓度和单位长度下的旋光度，是旋光性物质的物理常数。通过测定旋光度，可以鉴定物质的纯度、测定溶液的浓度、密度和鉴别光学异构体。
实验系统中，在确定容器和旋光物质的情况下， $\alpha$ 与L 均为常数，上式可写为：
$\theta=K\cdot c$
其中K是当前实验系统待标定参数。

MATLAB处理程序

算法流程图如下：
在这里插入图片描述

未处理旋光图像

在这里插入图片描述

二值化图像

程序输出角度为3.0248°，根据公式θ=K⋅c及标定得到的K值计算得未知溶液浓度：
在这里插入图片描述

close all;
clear all;
[filename1, pathname1] = uigetfile({'*.tif';'*.jpg'; '*.bmp'});
if filename1==0, return, end 
oldpath=cd;  
cd(pathname1);
[filename2, pathname2] = uigetfile({'*.tif';'*.jpg'; '*.bmp'});
cd(oldpath);
if filename2==0, return, end 
FILENAME1 = [pathname1 filename1];
FILENAME2 = [pathname2 filename2];
SorImg1=imread(FILENAME1);
SorImg2=imread(FILENAME2);
subplot(1,2,1),imshow(SorImg1);
subplot(1,2,2),imshow(SorImg2);
threshold1=0.35;
threshold2=0.35 ;
figure(1);   %楔形暗区图像
subplot(1,2,1),[Orientation1,Centroid1]=RotAngFun(SorImg1,filename1,threshold1,threshold2);
subplot(1,2,2),[Orientation2,Centroid2]=RotAngFun(SorImg2,filename2,threshold1,threshold2);
Delta=Orientation1-Orientation2;  %单位是度
display(['the angle is ' num2str(Delta)]);
figure(2);  %外接椭圆长轴与x轴
subplot(1,2,1),hold on,PlotFun(Orientation1,Centroid1,size(SorImg1));
subplot(1,2,2),hold on,PlotFun(Orientation2,Centroid2,size(SorImg2));


function [Orientation,Centroid]=RotAngFun(SorImg,filename,threshold1,threshold2)
flag=strfind(filename,'-');
% GrayImg=rgb2gray(SorImg);
GrayImg = SorImg;
if(isempty(flag))
    BwImg=im2bw(GrayImg,threshold1);
else
    BwImg=im2bw(GrayImg,threshold2);
end
BwImg=bwareaopen(BwImg,1000);
BwImg=~bwareaopen(~BwImg,1000);
se=strel('disk',5);
BwImg=imopen(BwImg,se);
BwImg=imfill(BwImg,'holes');
imshow(BwImg);
S=regionprops(BwImg,'Orientation','Centroid');
Orientation=S.Orientation;
Centroid=S.Centroid;
end

function  PlotFun(Orientation,Centroid,ImgSize)
k=tan(deg2rad(Orientation));
b=-Centroid(2)-k*Centroid(1);
x=1:ImgSize(2);
y=k*x+b;
plot(x,-y,'r');hold on;
y1=Centroid(2)*ones(1,ImgSize(2));
plot(x,y1,'b');
title(num2str(Orientation));
end