3D Plot on a Cylindrical Surface

CYL3D plots 3D data on a cylinder.

Description

CYL3D(Zin,theta_min,theta_max,Rho,H,meshscale) plots the 3D profile Zin as a mesh plot on a cylindrical surface of radius Rho, with mesh size determined by meshscale, for equally spaced
angles between theta_min and theta_max. H is the cylinder height.
 
CYL3D(Zin,...,meshscale,plotspec) plots the 3D profile Zin with a plot type specification. If plotspec = 'surf' a standard MATLAB
surface is plotted, whereas 'mesh', 'surfc' or 'meshc' will plot mesh, surface with countour, or mesh with contour, respectively.

CYL3D(Zin,...,meshscale,interpspec) plots the 3D profile Zin with the interpolation specification, interpspec, which can be one of 'spline', 'linear', 'nearest' or 'cubic'. The default interpolation is linear if not specified.

CYL3D(Zin,...,meshscale,Zscale) plots the 3D profile Zin with the data scaling factor, Zscale. This allows you to scale the peaks and troughs of the data on the surface if the radius is
relatively large. Zscale larger than 1 magnifies the data range correspondingly. Zscale defaults to 1 if not specified.

[Xout,Yout,Zout,Cmap] = CYL3D(Zin,...) returns output values corresponding to the Cartesian positions (Xout,Yout,Zout) with colour mapping, Cmap.

This function was developed under Ver. 7 R14SP1 and has been updated to be compatible with earlier Matlab versions, however there is no guarantee of compatibility.

If you do have any problems using the function please contact the author.

• MATLAB release: MATLAB 7.0.1 (R14SP1)

[John Hazel]

Absolutely useless. Producess errors even in examples:
>> P = peaks(20);
cyl3d(P,-pi/2,pi,20,20,.4,'spline',1.5);
??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> pol2cart at 12
RCST = r*cos(theta*rho);

Error in ==> cyl3d at 464
    [Xout(j,:) Yout(j,:) Zout(j,:)] = pol2cart(theta(j,:),Rho+Zscale*T(j,:),Z(j,:));

[J De Freitas]

John,
I'm disappointed it didn't work for you and I fully understand your frustration. It would be useful though if you could send me details of the Matlab Release you use and together we can do some detective work to sort out your problems. Cheers.

[Cristian Nastase]

 

[Cristian Nastase]

Very good. It has minor problems with Matlab 6.5. Took me 5 mins to correct it. All the examples work now. Most probably because v7.0 is more tolerant. Please find attached the PATCH. Question: How can I change it to do the same on a arbitrary inclined plane (or arbitrary surface in 3D).
--- diff
517c517
< h=colorbar; set(h,'YTickLabel',stepLabel);
---
> colorbar('YTickLabel',stepLabel);
571,572c571,572
< x(1:num,nLevel) = C(1,2:num+1)';
< y(1:num,nLevel) = C(2,2:num+1)';
---
> x(1:num,nLevel) = C(1,2:num+1);
> y(1:num,nLevel) = C(2,2:num+1);

[akbar sohoulipour]

thankyou

[Bouziane Yebka]