Archivo:Rohn Ellipsograph a.gif
Contenido de la página no disponible en otros idiomas.
Apariencia
De Wikipedia, la enciclopedia libre
Rohn_Ellipsograph_a.gif (500 × 500 píxeles; tamaño de archivo: 2,81 MB; tipo MIME: image/gif, bucleado, 200 frames, 8,0s)
Este es un archivo de Wikimedia Commons, un depósito de contenido libre hospedado por la Fundación Wikimedia. Más abajo se reproduce su página de descripción con la información sobre su origen y licencia. |
Resumen
DescripciónRohn Ellipsograph a.gif |
English: Karl Rohn Ellipsograph (similar to the HAFF-Ellipsograph)
Deutsch: Ellipsenzeichner (Ellipsenzirkel) nach Karl Rohn (ähnlich dem HAFF-Ellipsograph).
Er funktioniert auf Basis der Papierstreifenkonstruktion 1. Art und verwendet hierfür die Idee der gleitenden Kreise.
|
Fecha | |
Fuente | Trabajo propio |
Autor | Jahobr |
Otras versiones |
|
GIF desarrollo InfoField | |
Código fuente InfoField | MATLAB codefunction Rohn_Ellipsograph()
% source code for drawing Rohn_Ellipsograph
% produces a GIF and a SVG
%
% 2017-05-28 Jahobr
p = 1;
nFrames = 200;
xySize = 500;
scaleReduction = 2; % reduction for nice antialiasing
angleSmall = linspace(0,2*pi,nFrames+1); % define gear position in frames
angleSmall = angleSmall(1:end-1); % remove last frame, it would be double
angleSmall = circshift(angleSmall',round(nFrames*0.82)); % shift a view frames to get good visibility of all lines in frame 1
[pathstr,fname] = fileparts(which(mfilename)); % save files under the same name and at file location
figHandle = figure(15677755); clf
set(figHandle, 'Units','pixel');
set(figHandle, 'position',[1 1 [xySize xySize]*scaleReduction]); % big start image for antialiasing later [x y width height]
%set(figHandle, 'GraphicsSmoothing','on') % requires at least version 2014b
axesHandle = axes;
hold(axesHandle,'on')
set(axesHandle,'position',[-0.05 -0.05 1.1 1.1]); % stretch axis bigger as figure, easy way to get rid of ticks [x y width height]
xyLim = p*[-2.73 2.73];
xlim(xyLim); ylim(xyLim);
axis equal; drawnow;
diamScale = 1.4;
col.E = round([1 0 0 ]*255)/255; % red
col.Hori = round([0 0.7 0 ]*255)/255; % green
col.Vert = round([0.2 0.2 1 ]*255)/255; % blue
col.HoriDark = round([0 0.6 0 ]*255)/255; % green
col.VertDark = round([0 0 1 ]*255)/255; % blue
reducedRGBimage = uint8(ones(xySize,xySize,3,nFrames)); % allocate
for iFrame = 1:nFrames
currentAngle = angleSmall(iFrame);
cla(axesHandle) % fresh frame
disc(p*cos(currentAngle),0,p*diamScale,currentAngle,col.Hori,col.Hori,4,'-') % green disc
plotBox(1.25*diamScale*[p,-p], diamScale*p*[1.02 1.25],col.HoriDark,col.HoriDark,6) % horizontal bar
plotBox(1.25*diamScale*[p,-p],-diamScale*p*[1.02 1.25],col.HoriDark,col.HoriDark,6) % horizontal bar
disc(0,p*sin(currentAngle),p*diamScale,currentAngle,col.Vert,col.Vert,4,'-') % blue disc outline
plotBox( diamScale*p*[1.02 1.25],1.25*diamScale*[p,-p],col.VertDark,col.VertDark,6) % vertical bar
plotBox(-diamScale*p*[1.02 1.25],1.25*diamScale*[p,-p],col.VertDark,col.VertDark,6) % vertical bar
plot([0 0],[p,-p],'-.','linewidth',3,'color',col.VertDark) % vertical line
plot([p,-p],[0 0],'-.','linewidth',3,'color',col.HoriDark) % horizontal line
plot(p*cos(currentAngle),0,'.','MarkerSize',40,'linewidth',4,'color',col.Hori*0.7) % small dot to indicate the precise position
plot(0,p*sin(currentAngle),'.','MarkerSize',40,'linewidth',4,'color',col.Vert*0.7) % small dot to indicate the precise position
ellTwoOffset = 0.2;
% ellipse(0,0,ellTwoOffset,p-ellTwoOffset,col.E,4);
% plot(ellTwoOffset*cos(currentAngle),(p-ellTwoOffset)*sin(currentAngle),'.','MarkerSize',40,'linewidth',3,'color',col.E);
ellipse(0,0,p-ellTwoOffset,ellTwoOffset,col.E,4);
plot((p-ellTwoOffset)*cos(currentAngle),ellTwoOffset*sin(currentAngle),'.','MarkerSize',40,'linewidth',3,'color',col.E);
%% save animation
xlim(xyLim); ylim(xyLim); drawnow % set axis limits
if iFrame == 1 % save SVG
if ~isempty(which('plot2svg'))
plot2svg(fullfile(pathstr, [fname '_a.svg']),figHandle) % by Juerg Schwizer, See http://www.zhinst.com/blogs/schwizer/
else
disp('plot2svg.m not available; see http://www.zhinst.com/blogs/schwizer/');
end
end
f = getframe(figHandle);
reducedRGBimage(:,:,:,iFrame) = imReduceSize(f.cdata,scaleReduction); % the size reduction: adds antialiasing
end
map = createImMap(reducedRGBimage,64,[1 1 1;struct2map(col)]); % colormap
im = uint8(ones(xySize,xySize,1,nFrames)); % allocate
for iFrame = 1:nFrames
im(:,:,1,iFrame) = rgb2ind(reducedRGBimage(:,:,:,iFrame),map,'nodither');
end
imwrite(im,map,fullfile(pathstr, [fname '_a.gif']),'DelayTime',1/25,'LoopCount',inf) % save gif1
disp([fname '_a.gif has ' num2str(numel(im)/10^6 ,4) ' Megapixels']) % Category:Animated GIF files exceeding the 50 MP limit
%%
function disc(x,y,r,rot,colFa,colEd,linWi,linSt)
% x coordinates of the center
% y coordinates of the center
% r is the radius of the circle
% rot rotation of the slit in "rad"
% colFa facecolor
% colEd edgecolor
% linWi linewidth
% linSt linestyle
angleOffPoints = linspace(0,2*pi,300);
xc = x + r*sin(angleOffPoints);
yc = y + r*cos(angleOffPoints);
xs = 0.9* [r r -r -r r r]; % x slit
ys = 0.05*[0 -r -r r r 0]; % y slit
rotM = [cos(-rot) -sin(-rot); sin(-rot) cos(-rot)];
vecTemp = rotM*[xs; ys]; % rotate slit
xs = vecTemp(1,:)+x;
ys = vecTemp(2,:)+y;
[boolX,boolY] = polybool('subtraction', xc,yc, xs,ys);
boolX = boolX(~isnan(boolX)); % I do not know why there is an NaN
boolY = boolY(~isnan(boolY)); % I do not know why there is an NaN
patch(boolX,boolY,'k','linewidth',4,'edgecolor','none','facecolor',colFa,'FaceAlpha',.3); % ,'FaceAlpha',0.2% transparencc does onl work in 'renderer','OpenGL' but that fails in combination with "getframe"
plot(xc,yc,'linewidth',linWi,'linestyle',linSt,'color',colEd);
plot(xs,ys,'linewidth',linWi,'linestyle',linSt,'color',colEd);
xl = [0 0];
yl = r*[0.1 0.9];
vecTemp = rotM*[xl; yl];
plot( vecTemp(1,:)+x, vecTemp(2,:)+y,'k','linewidth',linWi,'linestyle',linSt,'color',colEd);
plot(-vecTemp(1,:)+x,-vecTemp(2,:)+y,'k','linewidth',linWi,'linestyle',linSt,'color',colEd);
function h = ellipse(x,y,a,b,col,linw)
% x coordinates of the center
% y coordinates of the center
% a radius1
% b radius2
angleOffPoints = linspace(0,2*pi,300);
xe = x + a*cos(angleOffPoints);
ye = y + b*sin(angleOffPoints);
h = plot(xe,ye,'-','linewidth',linw,'color',col);
function plotBox(x,y,colFa,colEd,linw)
% x [left right]
% y [top bottom]
% colFa face color [r g b]
% colEd edge color [r g b]
% linw line width
xs = [x(1) x(1) x(2) x(2) x(1) x(1)]; % x
ys = [mean(y) y(1) y(1) y(2) y(2) mean(y)]; % y joint in the middle of an edge to get nice corners
patch(xs,ys,colFa,'EdgeColor',colEd,'LineWidth',linw,'FaceAlpha',.3); %
function map = struct2map(RGB)
% RGB: struct of depth 1 with [r g b] in each field
fNames = fieldnames(RGB);
nNames = numel(fNames);
map = NaN(nNames,3); % allocate
for iName = 1:nNames
map(iName,:) = RGB.(fNames{iName}); %
end
function im = imReduceSize(im,redSize)
% Input:
% im: image, [imRows x imColumns x nChannel x nStack] (unit8)
% imRows, imColumns: must be divisible by redSize
% nChannel: usually 3 (RGB) or 1 (grey)
% nStack: number of stacked images
% usually 1; >1 for animations
% redSize: 2 = half the size (quarter of pixels)
% 3 = third the size (ninth of pixels)
% ... and so on
% Output:
% imNew: unit8([imRows/redSize x imColumns/redSize x nChannel x nStack])
%
% an alternative is : imNew = imresize(im,1/reduceImage,'bilinear');
% BUT 'bicubic' & 'bilinear' produces fuzzy lines
% IMHO this function produces nicer results as "imresize"
[nRow,nCol,nChannel,nStack] = size(im);
if redSize==1; return; end % nothing to do
if redSize~=round(abs(redSize)); error('"redSize" must be a positive integer'); end
if rem(nRow,redSize)~=0; error('number of pixel-rows must be a multiple of "redSize"'); end
if rem(nCol,redSize)~=0; error('number of pixel-columns must be a multiple of "redSize"'); end
nRowNew = nRow/redSize;
nColNew = nCol/redSize;
im = double(im).^2; % brightness rescaling from "linear to the human eye" to the "physics domain"; see youtube: /watch?v=LKnqECcg6Gw
im = reshape(im, nRow, redSize, nColNew*nChannel*nStack); % packets of width redSize, as columns next to each other
im = sum(im,2); % sum in all rows. Size of result: [nRow, 1, nColNew*nChannel]
im = permute(im, [3,1,2,4]); % move singleton-dimension-2 to dimension-3; transpose image. Size of result: [nColNew*nChannel, nRow, 1]
im = reshape(im, nColNew*nChannel*nStack, redSize, nRowNew); % packets of width redSize, as columns next to each other
im = sum(im,2); % sum in all rows. Size of result: [nColNew*nChannel, 1, nRowNew]
im = permute(im, [3,1,2,4]); % move singleton-dimension-2 to dimension-3; transpose image back. Size of result: [nRowNew, nColNew*nChannel, 1]
im = reshape(im, nRowNew, nColNew, nChannel, nStack); % putting all channels (rgb) back behind each other in the third dimension
im = uint8(sqrt(im./redSize^2)); % mean; re-normalize brightness: "scale linear to the human eye"; back in uint8
function map = createImMap(imRGB,nCol,startMap)
% createImMap creates a color-map including predefined colors.
% "rgb2ind" creates a map but there is no option to predefine some colors,
% and it does not handle stacked images.
% Input:
% imRGB: image, [imRows x imColumns x 3(RGB) x nStack] (unit8)
% nCol: total number of colors the map should have, [integer]
% startMap: predefined colors; colormap format, [p x 3] (double)
imRGB = permute(imRGB,[1 2 4 3]); % step1; make unified column-image (handling possible nStack)
imRGBcolumn = reshape(imRGB,[],1,3,1); % step2; make unified column-image
fullMap = double(permute(imRGBcolumn,[1 3 2]))./255; % "column image" to color map
[fullMap,~,imMapColumn] = unique(fullMap,'rows'); % find all unique colores; create indexed colormap-image
% "cmunique" could be used but is buggy and inconvenient because the output changes between "uint8" and "double"
nColFul = size(fullMap,1);
nColStart = size(startMap,1);
disp(['Number of colors: ' num2str(nColFul) ' (including ' num2str(nColStart) ' self defined)']);
if nCol<=nColStart; error('Not enough colors'); end
if nCol>nColFul; warning('More colors than needed'); end
isPreDefCol = false(size(imMapColumn)); % init
for iCol = 1:nColStart
diff = sum(abs(fullMap-repmat(startMap(iCol,:),nColFul,1)),2); % difference between a predefined and all colores
[mDiff,index] = min(diff); % find matching (or most similar) color
if mDiff>0.05 % color handling is not precise
warning(['Predefined color ' num2str(iCol) ' does not appear in image'])
continue
end
isThisPreDefCol = imMapColumn==index; % find all pixel with predefined color
disp([num2str(sum(isThisPreDefCol(:))) ' pixel have predefined color ' num2str(iCol)]);
isPreDefCol = or(isPreDefCol,isThisPreDefCol); % combine with overall list
end
[~,mapAdditional] = rgb2ind(imRGBcolumn(~isPreDefCol,:,:),nCol-nColStart,'nodither'); % create map of remaining colors
map = [startMap;mapAdditional];
|
Licencia
Yo, el titular de los derechos de autor de esta obra, la publico en los términos de la siguiente licencia:
Este archivo está disponible bajo la licencia Creative Commons Dedicación de Dominio Público CC0 1.0 Universal. | |
La persona que ha asociado una obra a este documento lo dedica al dominio público mediante la cesión mundial de sus derechos bajo la ley de derechos de autor y todos los derechos legales adyacentes propios de dicha, en el ámbito permitido por ley. Puedes copiar, modificar, distribuir y reproducir el trabajo, incluso con objetivos comerciales, sin pedir aprobación del autor.
http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse |
Elementos representados en este archivo
representa a
Algún valor sin elemento de Wikidata
28 may 2017
image/gif
Historial del archivo
Haz clic sobre una fecha y hora para ver el archivo tal como apareció en ese momento.
Fecha y hora | Miniatura | Dimensiones | Usuario | Comentario | |
---|---|---|---|---|---|
actual | 15:54 24 oct 2017 | 500 × 500 (2,81 MB) | Jahobr | graph smoothing, thick border | |
14:32 28 may 2017 | 500 × 500 (1,98 MB) | Jahobr | User created page with UploadWizard |
Usos del archivo
La siguiente página usa este archivo:
Uso global del archivo
Las wikis siguientes utilizan este archivo:
- Uso en fr.wikipedia.org