アンテナパターンスカイマップ(antenna pattern skymap)

jbname = 'test';
detectors = {'H1','L1','V1','LCGT'};
phiV=[-180:4:180];
thetaV=[-90:4:90];

outdatasuffix = ['fpfcts_',jbname,'.mat'];
fpfctsfile = fullfile('./',['fpfcts',outdatasuffix]);
[Fp,Fc,Ts] = genfpfcts(detectors,phiV,thetaV,'');


load /opt/workspace/myfiles/utils/coast.mat

imagesc(phiV,thetaV,sqrt(Fp(:,:,1).^2)); axis xy;
hold on;plot(long,lat,'w');hold off;
xlabel('longitude'); ylabel('latitude');
set(gca, 'XTick', -180:45:180);
set(gca, 'XTickLabel', {'180 W', '135 W', '90 W', '45 W', '0', '45 E', ...
        '90 E', '135 E', '180 E'});
set(gca, 'YTick', -90:45:90);
set(gca, 'YTickLabel', {'90 S', '45 S', '0', '45 N', '90 N'});
colorbar('EastOutside')
%colorbar('horiz');
title('Fp');
%shg

savedataname = ['antennaSkymap',detectors{1},detectors{2},detectors{3},detectors{4}];