chart.m

这是matlab在地球物理数据处理方面的源码

load shotinfo.mat

for i = 1:85
	plot(shotinfo(1,i),i,'r*')

	hold on
	plot([shotinfo(2,i):shotinfo(2,i)+47],i*ones(1,48),'go')
	hold on
end

xlabel('Station Number')
ylabel('Shot Number')


numshot = 85;

mincdp = (min(shotinfo(2,:)') + min(shotinfo(1,:)'))/2;
maxcdp = (max(shotinfo(2,:)'+47) + max(shotinfo(1,:)'))/2;

cdps =  mincdp:0.5:maxcdp;

cmps = zeros(length(cdps),numshot);
temp = 1:48;  % Always 48 shots in order here.
for i = 1:numshot
	if i == 19 | i == 44
	cmps(2*i-2:2*i-1+46,i) = temp';
	elseif i > 50 & i < 66
	cmps(2*i-1+10:2*i-1+47+10,i) = temp'; 
	else
	cmps(2*i-1:2*i-1+47,i) = temp';
	end
end

stack = zeros(length(cdps),2);
stack(:,1) = cdps';
for i = 1:length(cdps)
	for j = 1:85
	

	if cmps(i,j) ~= 0
		stack(i,2) = stack(i,2) + 1;
	plot(stack(i,1),j,'k.')
	end
	end
end		



plot(stack(:,1),stack(:,2))