function [ Area, output ] = reading_ums( inputfile, logornot, Sx, Sy )
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here
%   开始读取文件
file_id=fopen(inputfile,'rb');
fseek(file_id,10,'bof');
data=fread(file_id,'uint');
fclose(file_id);
%   产生横坐标
x=1:1:length(data);
%   绘图&面积计算初始化
Area=0;
if nargin==1
    logornot=0;
    Sx=1;
    Sy=length(data);
elseif nargin==2
    Sx=1;
    Sy=length(data);
end
%   绘图
if logornot==0
    plot(x,data);
    output=data;
elseif logornot==1
    index=find(data~=0);
    logdata=log(data);
    plot(x(index),logdata(index));
    output=logdata;
end
%   计算面积
for count=Sx:Sy
    Area=Area+data(count);
end
end