function [ PC1V ] = PGSTA( IPC1 )
%calculating CV vector, 8 ways
%IPC1は29(y)*15(x)の行列
%a=y, b=xであることに留意

enda = size(IPC1, 1);
endb = size(IPC1, 2);

for a=2:enda-1;
for b=2:endb-1;

if IPC1(a,b)<IPC1(a+1,b), u1(a,b)=1i;
else u1(a,b)=0.0;
end

if IPC1(a,b)<IPC1(a-1,b), u2(a,b)=-1i;
else u2(a,b)=0.0;
end

if IPC1(a,b)<IPC1(a,b+1), u3(a,b)=1;
else u3(a,b)=0.0;
end

if IPC1(a,b)<IPC1(a,b-1), u4(a,b)=-1;
else u4(a,b)=0.0;
end

if IPC1(a,b)<IPC1(a+1,b+1), u5(a,b)=(0.5)^(0.5)+(0.5)^(0.5)*1i;
else u5(a,b)=0.0;
end

if IPC1(a,b)<IPC1(a+1,b-1), u6(a,b)=-(0.5)^(0.5)+(0.5)^(0.5)*1i;
else u6(a,b)=0.0;
end

if IPC1(a,b)<IPC1(a-1,b-1), u7(a,b)=-(0.5)^(0.5)-(0.5)^(0.5)*1i;
else u7(a,b)=0.0;
end

if IPC1(a,b)<IPC1(a-1,b+1), u8(a,b)=(0.5)^(0.5)-(0.5)^(0.5)*1i;
else u8(a,b)=0.0;
end


end

end

PC1V=u1+u2+u3+u4+u5+u6+u7+u8;

%PC1Xの最終行、列に0行、0列を加える
% PC1V = [PC1V zeros(enda-1, 1)];
% PC1V = (PC1V zeros(1, endb));
% 
% Rp=real(PC1V);
% 
% Ip=imag(PC1V);
% 
% quiver(XI,YI,Rp,Ip);

%出力されるベクトルのエッジはダミーデータなので削除すること

end






