You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
84 lines
2.1 KiB
Matlab
84 lines
2.1 KiB
Matlab
function V1_cross_orient_tuning
|
|
grating_wavel=26;
|
|
|
|
pxsize=65;
|
|
iterations=6;
|
|
|
|
grating_angles=[0:22.5:179];
|
|
diff_angles=[22.5:22.5:90];
|
|
node=5;
|
|
|
|
[wFFon,wFFoff]=dim_conv_V1_filter_definitions;
|
|
|
|
da=0;
|
|
for diff_angle=diff_angles
|
|
da=da+1;
|
|
figure(da),clf;
|
|
for test=[1:3]
|
|
ga=0;
|
|
for grating_angle=grating_angles
|
|
ga=ga+1;
|
|
if test==1
|
|
I=image_orientation_square(pxsize,grating_wavel,2,grating_angle);
|
|
elseif test==2
|
|
I=image_orientation_square(pxsize,grating_wavel,2,grating_angle+diff_angle);
|
|
else
|
|
I=image_cross_orientation_square(pxsize,grating_wavel,2,grating_angle,diff_angle,1);
|
|
end
|
|
I=I*0.5+0.5;
|
|
[Ion,Ioff]=preprocess_image(I);
|
|
|
|
[a,b]=size(I);
|
|
%plot original image
|
|
maxsubplot(4,length(grating_angles),ga+(test-1)*length(grating_angles)),
|
|
imagesc(I,[0,1]);
|
|
axis('equal','tight'), set(gca,'XTick',[],'YTick',[],'FontSize',11);
|
|
drawnow;
|
|
|
|
%perform competition
|
|
[y,ron,roff,eon,eoff,Y]=dim_conv_on_and_off(wFFon,wFFoff,Ion,Ioff,iterations);
|
|
y=mean(Y,4);
|
|
resp(test,ga)=y(ceil(a/2),ceil(b/2),node);
|
|
|
|
if ga==node
|
|
respt(:,test,da)=Y(ceil(a/2),ceil(b/2),node,:);
|
|
end
|
|
end
|
|
end
|
|
figure(da),clf;
|
|
subplot(2,2,1);
|
|
plot(grating_angles,resp(1,:),'k--o','LineWidth',3,'MarkerSize',9,'MarkerFaceColor','w')
|
|
hold on
|
|
plot(grating_angles,resp(2,:),'b--o','LineWidth',3,'MarkerSize',9,'MarkerFaceColor','w')
|
|
plot(grating_angles,resp(3,:),'r-s','LineWidth',3,'MarkerSize',9,'MarkerFaceColor','w')
|
|
axis([0,180,0,0.01]);
|
|
set(gca,'XTick',[0:45:180],'YTick',[0:0.01:0.01],'FontSize',20,'Box','off');
|
|
xlabel('Orientation (degrees)'),ylabel('Response')
|
|
|
|
maxsubplot(11,11,104)
|
|
I=image_cross_orientation_square(pxsize,grating_wavel,2,90,diff_angle,1);
|
|
I=I*0.5+0.5;
|
|
imagesc(I,[0,1]);
|
|
axis('equal','tight')
|
|
set(gca,'XTick',[],'YTick',[]);
|
|
cmap=colormap('gray');
|
|
drawnow;
|
|
end
|
|
figure(5),clf
|
|
|
|
da=find(diff_angles==67.5);
|
|
subplot(2,2,1);
|
|
plot(respt(:,:,da))
|
|
|
|
da=find(diff_angles==45);
|
|
subplot(2,2,2);
|
|
plot(respt(:,:,da))
|
|
|
|
da=find(diff_angles==22.5);
|
|
subplot(2,2,3);
|
|
plot(respt(:,:,da))
|
|
|
|
da=find(diff_angles==90);
|
|
subplot(2,2,4);
|
|
plot(respt(:,:,da))
|