remove discs from lat/lon coverage

This is my code:

[xIntrinsic, yIntrinsic] = geographicToIntrinsic(r, discLat, discLon);
  [xc, yc] = polysplit(xIntrinsic, yIntrinsic);
  for ii = 1:(numel(xc)-1)
      xii = [xc{ii}; xc{ii+1}];
      yii = [yc{ii}; yc{ii+1}];
      k = convhull(xii,yii);
      xc{ii} = xii(k);
      yc{ii} = yii(k);
  end
  [xq, yq] = polyjoin(xc(1:end-1), yc(1:end-1));
  [latQ, lonQ] = intrinsicToGeographic(r, xq, yq);

When I plot, I see the discLat/discLon points (as circles) along with the bounding lines:

Is there an easy way to remove the discs so that all I see is the bounding outside lines? I'm using the areaint function on the latQ/lonQ points to get an area, and it is coming back way high, I believe because it includes those discs.

Thanks.

Answer 1

[Yuebin Zhou]

It seems the code is not complete according to the description. Can you put the entire of it online?

[Paul Huter]

Here is the script I am working from provided by Mathworks.