%处理架构的脚本%曲线的首次读取扫描图像a=imread（'scan_6523211656_1.jpg'）；b=a（230:530158:558,1）；【I J】=发现（b<150）；jmin=最小值（J）；jmax=最大值（J）；newI=零（1，jmax-jmin+1）；newJ=零（1，jmax jmin+1）；计数=1；对于ind=jmin：jmaxvals=查找（J==ind）；如果numel（vals）>0newI（计数）=舍入（总和（I（vals））/numel（val））；newJ（计数）=ind；计数=计数+1；结束结束newJ=新J（1:count-1）；newI=新I（1:count-1）；新I=300-newI；图（newJ，newI）；fid=fopen（'silh.dat'，'w'）；fprintf（fid，“%3d”，newJ）；fprintf（fid，'\n'）；fprintf（fid，'%3d'，newI）；fclose（fid）；