contourm fails to identify closed contour lines in some cases

1 次查看（过去 30 天）
Alexander
Alexander 2015-8-1
编辑： Alexander 2015-8-3
What I have is a 2-D geolocated data grid F for hundreds of timesteps. Though it is a global dataset, I am interested in the coordinates of a certain contour line (value=2), which encircles the main parts of the north polar region.
To get these coordinates I have to make a contourm plot, there is no low-level function in the mapping toolbox (as contourc in the x/y-plane). I get a contour matrix C which contains the lat/lon coordinates of every closed contour line (or not closed if on the boundary).
The data of each contourline are separated in C, so they can be processed individually.
The code looks like that:
lat=90:0.75:-90;
lon=0:0.75:359.25;
figure
axesm('MapProjection','stereo','Origin',[90 0 0],'MapLatLimit',[27 90])
C = contourm(lat,lon,F,[2 2]);
In certain cases the result is as expected, shown here. Note that the contour encircles the northpole.
But in other cases (here 12 hours later) a normally closed contour line is splitted in two separated sequences in C, for illustration I have plotted these two parts.
As it seems this happens if the contour line do not encircle the pole. In these cases the first point of the green line starts with lon=0 and as it seems the algorithm to identify a single contour stops when a point reaches lon=0 again, not taking into account the points that follow. However it should be one contour line.
As I said I am mainly interested in the coordinates not in the plot. Since I am dealing with many fields I cannot effort to check in every case whether the contour is splitted and if so which sequences in C I need to join.
If a developer of the mapping toolbox read this: Is this a known bug?
Or is there a workaround for this problem?
I am grateful for any hints or suggestions.
Alex

请先登录，再回答此问题。

回答（1 个）

Chad Greene
Chad Greene 2015-8-1
I've had this problem with contourm on stereographic projections. After you've initialized the map, try this:
[X,Y] = mfwdtran(lat,lon);
[C,h] = contour(X,Y,F,[2 2]);
  3 个评论
显示 1更早的评论
Chad Greene
Chad Greene 2015-8-1
Oh interesting. I'm not sure what the problem is. Perhaps an issue that Matlab doesn't know 360 degrees longitude is the same thing as 0 degrees? Most of my polar datasets are equally spaced in polar stereographic coordinates, not equally spaced in geo coordinates, so that may be the issue.
A clunky fix, but perhaps you can convert C to x,y,z coordinates with C2xyz and color the contour lines manually by plotting with plot?
Alexander
Alexander 2015-8-1
编辑：Alexander 2015-8-3
Meanwhile I tried to find the code for the algorithm of finding the contour lines in contourm. I didn't manage but I think no matter if the data are projected or not everything is based on the contourc function. A short description of the Contouring Algorithm can be found in the Matlab Help.
So I think you're right. It seems that Matlab cannot deal with Longitudes in a right manner if one has a global data set.
The fact that it doesn't work neither with contourm nor contour (after re-mapping on x/y) seems to confirm this.
A clunky fix, but perhaps you can convert C to x,y,z coordinates with C2xyz and color the contour lines manually by plotting with plot?
The problem is I need the coordinates of the line rather than the plots.
But thanks anyway.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Mapping Toolbox 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by