如何在matlab中绘制极坐标立体图中的x，y，z数据？

Question

我有来自 SAR 卫星仪器的数据，这些数据已经在极地立体投影中。数据带有 x,y 坐标，而不是 lat,lon。谁能提供有关如何在北半球方位角图上绘制 z 数据（归一化雷达截面）的见解？ x 和 y 变量的值范围均为

[-5498550, 5498450]

使用

imagesc

，我能够确认数据处于正确的投影位置，因为我可以看到格陵兰岛的轮廓。但是，我不确定如何在 matlab 中创建正确的地图轴或如何正确投影 z 数据。

Answer 1

你有绘图工具箱吗？如果是的话，试试这个：

xLimits = [-5498550, 5498450];
yLimits = xLimits

% grid
[xGrid, yGrid] = meshgrid(linspace(xLimits(1), xLimits(2), size(z, 2)), ...
                          linspace(yLimits(1), yLimits(2), size(z, 1)));

% Set up the map axes with an azimuthal projection centered on the North Pole
figure;
axesm('MapProjection', 'stereo', 'Origin', [90 0], 'FLatLimit', [60 90]);
axis off
setm(gca, 'MapLatLimit', [-90 90], 'MapLonLimit', [-180 180]);

% Display the data on the map
geoshow(xGrid, yGrid, z, 'DisplayType', 'surface');

mlabel on;
plabel on;
gridm on;