colour.m

细胞生长结构可视化工具箱-MATLAB Toolbox1999.zip

% Calculate a non-linear greyscale colour bar which maximises contrast

NoC = 64;	% No Colours to be used

freq = zeros(NoC,1);

scale = max(max(cmap1/(NoC-1)));

ii=0;
xi=minx:(maxx-minx)/NoPts:maxx;
yi=maxy:-(maxy-miny)/NoPts:miny;
for i=minx:(maxx-minx)/NoPts:maxx
	ii=ii+1;
	jj=NoPts+2;
	for j=maxy:-(maxy-miny)/NoPts:miny
		jj=jj-1;
		Ip=[i,j];
		% Find enclosing triangle
		tri = findtri(Ip,wvis,C);
		if (tri ~= [])	% if is within network
			index = floor(cmap1(jj,ii)./scale) + 1;
			freq(index) = freq(index)+1;
		end;
	end
end

cumfreq=[];
for i=1:NoC
	cumfreq(i) = sum(freq(1:i));
end

plot(freq);
hold on;
plot(cumfreq);