From arctica1963 at gmail.com Mon Nov 9 17:07:37 2020 From: arctica1963 at gmail.com (arctica1963) Date: Mon, 9 Nov 2020 09:07:37 -0700 (MST) Subject: [Scilab-users] Multiple regression on semi-log plot Message-ID: <1604938057702-0.post@n3.nabble.com> Hello, I am looking to determine multiple regression lines from a single power spectral dataset (log power vs radial wavenumber), and was wondering if it is feasible in Scilab to compute something similar to the attached plot? I did locate a Matlab code for finding a turning point in a plot and perhaps this may be a route to find the change in curvature? clear clf() //to find the "knee" dt = 0.01; t = 0:dt:1; y = exp(-10*t); //Compute first and second derivatives by finite differencing (centred) //yp = nan(size(y)); //ypp = nan(size(y)); yp=%nan*y; // similar behaviour to Matlab nan(size(y)) ypp=%nan*y; yp(2:$-1) = (y(3:$) - y(1:$-2))/(2*dt); // first derivative ypp(2:$-1) = (y(3:$) + y(1:$-2) - 2*y(2:$-1)) / (dt^2); // second derivative //Compute the curvature k = abs(ypp) ./ (1 + yp.^2).^1.5 //Find the maximum curvature point and plot it on the curve [kmax, idx] = max(k); plot(t, y, 'b', t(idx), y(idx), 'ro') //plot the curvature figure() plot(t, k) Any pointers or ideas would be welcome. Thanks Lester -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From stephane.mottelet at utc.fr Mon Nov 9 18:37:45 2020 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Mon, 9 Nov 2020 18:37:45 +0100 Subject: [Scilab-users] Multiple regression on semi-log plot In-Reply-To: <1604938057702-0.post@n3.nabble.com> References: <1604938057702-0.post@n3.nabble.com> Message-ID: Hi, This is an easy task that can be done by fitting a piecewise-affine function like this: y = a*(x-theta)+phi, for x < theta y = b*(x-theta)+phi, for x >= theta Here is an example : function y=fun(x, param) a = param(1); b = param(2); theta = param(3); phi = param(4); c = (x=theta)*b; y = c.*(x-theta)+phi; endfunction function r=residual(param, x, y, f) r = sum((y-f(x,param)).^2); endfunction function [r, dr, ind]=costf(p, ind, x, y, f) r = residual(p,x,y,f); dr = numderivative(list(residual,x,y,fun),p); endfunction x=linspace(0,4,50); theta = 1.5; phi = 1; a = grand(1,1,"unf",-2,2) b = grand(1,1,"unf",-2,2) y = fun(x,[a,b,theta,phi]) + rand(x,"normal")/5; [ropt,popt] = optim(list(costf,x,y,fun),[0,0,mean(x),mean(y)]); clf plot(x,y,'o',x,fun(x,popt),popt(3),popt(4),'xr') S. Le 09/11/2020 ? 17:07, arctica1963 a ?crit?: > Hello, > > I am looking to determine multiple regression lines from a single power > spectral dataset (log power vs radial wavenumber), and was wondering if it > is feasible in Scilab to compute something similar to the attached plot? > > I did locate a Matlab code for finding a turning point in a plot and perhaps > this may be a route to find the change in curvature? > > clear > clf() > //to find the "knee" > dt = 0.01; > t = 0:dt:1; > y = exp(-10*t); > //Compute first and second derivatives by finite differencing (centred) > //yp = nan(size(y)); > //ypp = nan(size(y)); > yp=%nan*y; // similar behaviour to Matlab nan(size(y)) > ypp=%nan*y; > yp(2:$-1) = (y(3:$) - y(1:$-2))/(2*dt); // first derivative > ypp(2:$-1) = (y(3:$) + y(1:$-2) - 2*y(2:$-1)) / (dt^2); // second derivative > //Compute the curvature > k = abs(ypp) ./ (1 + yp.^2).^1.5 > //Find the maximum curvature point and plot it on the curve > [kmax, idx] = max(k); > plot(t, y, 'b', t(idx), y(idx), 'ro') > //plot the curvature > figure() > plot(t, k) > > Any pointers or ideas would be welcome. > > Thanks > Lester > > > > > > -- > Sent from: https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Wed Nov 11 13:20:18 2020 From: p.muehlmann at gmail.com (P M) Date: Wed, 11 Nov 2020 13:20:18 +0100 Subject: [Scilab-users] IPCV imread Message-ID: Dear Experts, please improve my understanding of the IPCV function imread for following context: I have an image which is told to be a 3-channel image. At least with IrfanView I get a bit depth of 24 ... so I expect it to be an RBG image. However: The image contains - on purpose - only 2 colors. Now: img = imread(imPath); // reads image as boolean img = double(img); // convert to double. Question: Why does the imread-function read the image as a boolean image and not as a 3-channel image? Note: I am totally fine with the result of imread..just being curious. Thank you, Philipp -------------- next part -------------- An HTML attachment was scrubbed... URL: From chinluh.tan at bytecode-asia.com Thu Nov 12 12:32:10 2020 From: chinluh.tan at bytecode-asia.com (Chin Luh Tan) Date: Thu, 12 Nov 2020 19:32:10 +0800 Subject: [Scilab-users] IPCV imread In-Reply-To: References: Message-ID: <175bc3a2c19.c4a62d00253238.5452952409873861616@bytecode-asia.com> Hi,? These was a part in the code to check if the image read into Scilab only consist of 2 values, min and max, and it will be automatically converted to binary image.? If you're using Windows, consider to download version 4.1.2.3 from?https://github.com/tanchinluh/IPCV/tags?and use the extra input to force reading in 3 layers.? img = imread(imPath,?IMREAD_COLOR = 1); Hope this helps. Regards, Chin Luh ---- On Wed, 11 Nov 2020 20:20:18 +0800 P M wrote ---- Dear Experts, please improve my understanding of the IPCV function imread for following context: I have an image which is told to be a 3-channel image. At least with IrfanView I get a bit depth of 24 ... so I expect it to be an RBG image. However: The image contains - on purpose - only 2 colors. Now: img = imread(imPath);?????? // reads image as boolean img = double(img);??????????? // convert to double. Question:??????? Why does the imread-function read the image as a boolean image and not as a 3-channel image? Note:????????????? I am totally fine with the result of imread..just being curious. Thank you, Philipp _______________________________________________ users mailing list mailto:users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Thu Nov 12 16:48:27 2020 From: p.muehlmann at gmail.com (P M) Date: Thu, 12 Nov 2020 16:48:27 +0100 Subject: [Scilab-users] IPCV imread In-Reply-To: <175bc3a2c19.c4a62d00253238.5452952409873861616@bytecode-asia.com> References: <175bc3a2c19.c4a62d00253238.5452952409873861616@bytecode-asia.com> Message-ID: Quote: [...] These was a part in the code to check if the image read into Scilab only consist of 2 values, min and max, and it will be automatically converted to binary image. Quote End I thought about this, but did not expect it to be true. :-) Thanks for clarification. Philipp Am Do., 12. Nov. 2020 um 12:32 Uhr schrieb Chin Luh Tan < chinluh.tan at bytecode-asia.com>: > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > *Hi, These was a part in the code to check if the image read into Scilab > only consist of 2 values, min and max, and it will be automatically > converted to binary image. If you're using Windows, consider to download > version 4.1.2.3 from https://github.com/tanchinluh/IPCV/tags > and use the extra input to force > reading in 3 layers. img = imread(imPath, IMREAD_COLOR = 1);Hope this > helps.Regards,Chin Luh---- On Wed, 11 Nov 2020 20:20:18 +0800 P M > > wrote ----Dear > Experts,please improve my understanding of the IPCV function imread for > following context:I have an image which is told to be a 3-channel image.At > least with IrfanView I get a bit depth of 24 ... so I expect it to be an > RBG image.However: The image contains - on purpose - only 2 colors.Now:img > = imread(imPath); // reads image as booleanimg = > double(img); // convert to double.Question: Why does the > imread-function read the image as a boolean image and not as a 3-channel > image?Note: I am totally fine with the result of imread..just > being curious.Thank > you,Philipp_______________________________________________users mailing > list users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > * -------------- next part -------------- An HTML attachment was scrubbed... URL: From arctica1963 at gmail.com Sat Nov 14 13:40:53 2020 From: arctica1963 at gmail.com (arctica1963) Date: Sat, 14 Nov 2020 05:40:53 -0700 (MST) Subject: [Scilab-users] Multiple regression on semi-log plot In-Reply-To: References: <1604938057702-0.post@n3.nabble.com> Message-ID: <1605357653021-0.post@n3.nabble.com> Hello, Thanks for the idea and suggestions. Not too sure how to apply it, if you could give some pointers on the attached data and code. The ultimate idea is to get the slopes of the straight line segments. Many thanks, Lester clear clf() // Read data - wavelength (in km)), power, 1 standard deviation // Unknown data length; 3 columns -default space delimited // PSD_wavelength.dat from GMT grdfft radially averaged power spectra data = read('PSD_wavelength.dat',-1,3); wavelength = data(:,1); power = data(:,2); std_dev1 = data(:,3); ln_power = log(power); wavenumber = 1./wavelength; f=gcf(); //plot(wavenumber, ln_power) scatter(wavenumber, ln_power,'marker','.') a=gce().children; a.mark_mode = "on" a.mark_style = 0 a.mark_size_unit = "point" a.mark_size=3 xlabel ('wavenumber k (km-1)') ylabel ('Log (Power)') PSD_wavelength.dat -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From stephane.mottelet at utc.fr Sun Nov 15 23:04:51 2020 From: stephane.mottelet at utc.fr (stephane.mottelet at utc.fr) Date: Sun, 15 Nov 2020 23:04:51 +0100 Subject: [Scilab-users] Multiple regression on semi-log plot In-Reply-To: <1605357653021-0.post@n3.nabble.com> References: <1604938057702-0.post@n3.nabble.com> <1605357653021-0.post@n3.nabble.com> Message-ID: <20201115230451.Horde.gAA0kUnVJh1tckvESs5011w@webmail.utc.fr> Hello, You just have to replace "x" by "wavelength" and "y" by "ln_power". The slopes are the first two component of the optimal vector "popt" : clearclf()// Read data - wavelength (in km)), power, 1 standard deviation// Unknown data length; 3 columns -default space delimited // PSD_wavelength.dat from GMT grdfft radially averaged power spectra data = read('PSD_wavelength.dat',-1,3); wavelength = data(:,1); power = data(:,2); std_dev1 = data(:,3); ln_power = log(power); wavenumber = 1./wavelength; f=gcf(); //plot(wavenumber, ln_power) scatter(wavenumber, ln_power,'marker','.') a=gce().children; a.mark_mode = "on"a.mark_style = 0a.mark_size_unit = "point"a.mark_size=3 xlabel ('wavenumber k (km-1)')ylabel ('Log (Power)') function y=fun(x, param) a = param(1); b = param(2); theta = param(3); phi = param(4); c = (x=theta)*b; y = c.*(x-theta)+phi; endfunction function r=residual(param, x, y, f) r = sum((y-f(x,param)).^2); endfunction function [r, dr, ind]=costf(p, ind, x, y, f) r = residual(p,x,y,f); dr = numderivative(list(residual,x,y,fun),p); endfunction [ropt,popt] = optim(list(costf,wavenumber,ln_power,fun),[0,0,mean(wavenumber),mean(ln_power)]); plot(wavenumber,fun(wavenumber,popt),'r',popt(3),popt(4),'xr') arctica1963 a ?crit?: > Hello, > > Thanks for the idea and suggestions. Not too sure how to apply it, if you > could give some pointers on the attached data and code. The ultimate idea is > to get the slopes of the straight line segments. Many thanks, Lester > > clear > clf() > // Read data - wavelength (in km)), power, 1 standard deviation > // Unknown data length; 3 columns -default space delimited > > // PSD_wavelength.dat from GMT grdfft radially averaged power spectra > > data = read('PSD_wavelength.dat',-1,3); > > wavelength = data(:,1); > power = data(:,2); > std_dev1 = data(:,3); > > ln_power = log(power); > > wavenumber = 1./wavelength; > f=gcf(); > > //plot(wavenumber, ln_power) > > scatter(wavenumber, ln_power,'marker','.') > > a=gce().children; > a.mark_mode = "on" > a.mark_style = 0 > a.mark_size_unit = "point" > a.mark_size=3 > > xlabel ('wavenumber k (km-1)') > ylabel ('Log (Power)') > > PSD_wavelength.dat > > > -- > Sent from: > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html > _______________________________________________ > users mailing list > users at lists.scilab.orghttps://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From arctica1963 at gmail.com Mon Nov 16 09:49:21 2020 From: arctica1963 at gmail.com (arctica1963) Date: Mon, 16 Nov 2020 01:49:21 -0700 (MST) Subject: [Scilab-users] Multiple regression on semi-log plot In-Reply-To: <20201115230451.Horde.gAA0kUnVJh1tckvESs5011w@webmail.utc.fr> References: <1604938057702-0.post@n3.nabble.com> <1605357653021-0.post@n3.nabble.com> <20201115230451.Horde.gAA0kUnVJh1tckvESs5011w@webmail.utc.fr> Message-ID: <1605516561566-0.post@n3.nabble.com> Hello, Thanks for showing how it works, always good to get a proper understanding. Just wondering if the code can be adapted to find multiple straight line segments? Looking at my original test data, it appears to have "kinks" at wavenumber (0.05, 0.12 and 0.16), the last one is where the x plots. I am guessing that it needs some threshold or tolerance control as the angles at the kink points may be small or larger. Many thanks for the guidance. Lester -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From stephane.mottelet at utc.fr Mon Nov 16 09:55:50 2020 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Mon, 16 Nov 2020 09:55:50 +0100 Subject: [Scilab-users] Multiple regression on semi-log plot In-Reply-To: <1605516561566-0.post@n3.nabble.com> References: <1604938057702-0.post@n3.nabble.com> <1605357653021-0.post@n3.nabble.com> <20201115230451.Horde.gAA0kUnVJh1tckvESs5011w@webmail.utc.fr> <1605516561566-0.post@n3.nabble.com> Message-ID: <25a67dbb-44c5-22f1-04c3-32cf267a85f0@utc.fr> Hi, It is possible to fit a more general piecewise-affine model, but you will have to know in advance the number of "kinks" S. Le 16/11/2020 ? 09:49, arctica1963 a ?crit?: > Hello, > > Thanks for showing how it works, always good to get a proper understanding. > > Just wondering if the code can be adapted to find multiple straight line > segments? Looking at my original test data, it appears to have "kinks" at > wavenumber (0.05, 0.12 and 0.16), the last one is where the x plots. I am > guessing that it needs some threshold or tolerance control as the angles at > the kink points may be small or larger. > > Many thanks for the guidance. > Lester > > > > -- > Sent from: https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet From Christophe.Dang at sidel.com Mon Nov 16 10:26:07 2020 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Mon, 16 Nov 2020 09:26:07 +0000 Subject: [Scilab-users] {EXT} Re: Multiple regression on semi-log plot In-Reply-To: <25a67dbb-44c5-22f1-04c3-32cf267a85f0@utc.fr> References: <1604938057702-0.post@n3.nabble.com> <1605357653021-0.post@n3.nabble.com> <20201115230451.Horde.gAA0kUnVJh1tckvESs5011w@webmail.utc.fr> <1605516561566-0.post@n3.nabble.com> <25a67dbb-44c5-22f1-04c3-32cf267a85f0@utc.fr> Message-ID: Hello, > De : users De la part de St?phane Mottelet > Envoy? : lundi 16 novembre 2020 09:56 > > It is possible to fit a more general piecewise-affine model, > but you will have to know in advance the number of "kinks" Well, I think it is possible to automatically adjust the number of segments: 1. Perform the regression with only n = 1 line segment. 2. Increase the number of segments by 1, n = n + 1. Perform a new regression. You can detect the furthest point to know where to separate the segments. 3. Compare the quadratic errors between the current and the preceding regression. If the difference is small enough then end else go back to 2. Regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer General This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. From stephane.mottelet at utc.fr Mon Nov 16 10:45:52 2020 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Mon, 16 Nov 2020 10:45:52 +0100 Subject: [Scilab-users] {EXT} Re: Multiple regression on semi-log plot In-Reply-To: References: <1604938057702-0.post@n3.nabble.com> <1605357653021-0.post@n3.nabble.com> <20201115230451.Horde.gAA0kUnVJh1tckvESs5011w@webmail.utc.fr> <1605516561566-0.post@n3.nabble.com> <25a67dbb-44c5-22f1-04c3-32cf267a85f0@utc.fr> Message-ID: Hi, Le 16/11/2020 ? 10:26, Dang Ngoc Chan, Christophe a ?crit?: > Hello, > >> De : users De la part de St?phane Mottelet >> Envoy? : lundi 16 novembre 2020 09:56 >> >> It is possible to fit a more general piecewise-affine model, >> but you will have to know in advance the number of "kinks" > Well, I think it is possible to automatically adjust the number of segments: > > 1. Perform the regression with only n = 1 line segment. > > 2. Increase the number of segments by 1, n = n + 1. > Perform a new regression. > You can detect the furthest point to know where to separate the segments. > > 3. Compare the quadratic errors between the current and the preceding regression. > If the difference is small enough then end That's the more difficult part. Considering the provided data, which will never fit well this kind of model, the quadratic error will decrease quite regularly. If the real concern is about finding the points with maximum curvature, a non parametric model should be more adequate (Local regression, Kernel regression, splines, or wathever). S. > else go back to 2. > > Regards > > -- > Christophe Dang Ngoc Chan > Mechanical calculation engineer > > General > This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet From arctica1963 at gmail.com Fri Nov 20 08:40:39 2020 From: arctica1963 at gmail.com (arctica1963) Date: Fri, 20 Nov 2020 00:40:39 -0700 (MST) Subject: [Scilab-users] {EXT} Re: Multiple regression on semi-log plot In-Reply-To: References: <1604938057702-0.post@n3.nabble.com> <1605357653021-0.post@n3.nabble.com> <20201115230451.Horde.gAA0kUnVJh1tckvESs5011w@webmail.utc.fr> <1605516561566-0.post@n3.nabble.com> <25a67dbb-44c5-22f1-04c3-32cf267a85f0@utc.fr> Message-ID: <1605858039330-0.post@n3.nabble.com> Hello, Just a quick thought. But could a moving window over the wavenumber provide a simpler route to locating the kink points for computing the slope segments? Perhaps a window size of say 0.05 and shifted 0.025 with the existing methodology you outlined Stephane. Lester -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From p.muehlmann at gmail.com Fri Nov 20 10:03:39 2020 From: p.muehlmann at gmail.com (P M) Date: Fri, 20 Nov 2020 10:03:39 +0100 Subject: [Scilab-users] mesh2d - finding identical points Message-ID: Dear all, using the mesh2d - function I get an error: mesh2di: some points are identical. How to identify these points and remove the duplicates? There are roughly 500'000 coordinates.... Thank you, Philipp -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul18fr at gmail.com Fri Nov 20 10:22:35 2020 From: paul18fr at gmail.com (paul francedixhuit) Date: Fri, 20 Nov 2020 10:22:35 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: Hi Have a look to "unique" to identify (and remove) duplicates Paul Le ven. 20 nov. 2020 ? 10:04, P M a ?crit : > Dear all, > > using the mesh2d - function I get an error: mesh2di: some points are > identical. > > > How to identify these points and remove the duplicates? > > There are roughly 500'000 coordinates.... > > Thank you, > Philipp > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Fri Nov 20 10:24:20 2020 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Fri, 20 Nov 2020 09:24:20 +0000 Subject: [Scilab-users] {EXT} mesh2d - finding identical points In-Reply-To: References: Message-ID: Hello, > De : users De la part de P M Envoy? : > vendredi 20 novembre 2020 10:04 > > How to identify these points and remove the duplicates? You might try to use unique(). Maybe the following discussion can help: http://mailinglists.scilab.org/Scilab-users-Remove-duplicate-rows-and-sustain-original-order-of-rows-td4038110.html regards -- Christophe Dang Ngoc Chan Mechanical Calculation Engineer General This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. From stephane.mottelet at utc.fr Fri Nov 20 10:57:19 2020 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Fri, 20 Nov 2020 10:57:19 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: Hi Philipp, Glad to hear that some people use mesh2d. I know we have to improve the error message by getting the numbers at the Fortran level... There is also other errors which are not easy to fix without knowing the number of the problematic point (point too close to the boundary, for example). Can you file a bug on BZ ? S. Le 20/11/2020 ? 10:03, P M a ?crit?: > Dear all, > > using the mesh2d - function I get an error:???? mesh2di: some points > are identical. > > > How to identify these points and remove the duplicates? > > There are roughly 500'000 coordinates.... > > Thank you, > Philipp > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Fri Nov 20 11:30:58 2020 From: p.muehlmann at gmail.com (P M) Date: Fri, 20 Nov 2020 11:30:58 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: OK, some more background: I actually want to perform a delaunay triangulation on a set of X-Y-Z-coordinates. I know of CGLAB, but it seems only available to Scilab 6.0.x, while I am on Scilab 6.1.0 So the nearest thing to a delaunay triangulation I found is mesh2d... Because of the delaunay triangulation, I am not sure if the order of points may change. About unique: The result I get is actual of same size as the input matrix...I therefore do assume that each point is unique in the first place (which they actually should). What does repeat are individual x, y, z coordinates, but in different combinations. So the actual point is different. e.g.: x may repeat, but than the y-z-coordinates are different from point to point Thanks, Philipp Am Fr., 20. Nov. 2020 um 10:57 Uhr schrieb St?phane Mottelet < stephane.mottelet at utc.fr>: > Hi Philipp, > > Glad to hear that some people use mesh2d. I know we have to improve the > error message by getting the numbers at the Fortran level... There is also > other errors which are not easy to fix without knowing the number of the > problematic point (point too close to the boundary, for example). Can you > file a bug on BZ ? > > S. > Le 20/11/2020 ? 10:03, P M a ?crit : > > Dear all, > > using the mesh2d - function I get an error: mesh2di: some points are > identical. > > > How to identify these points and remove the duplicates? > > There are roughly 500'000 coordinates.... > > Thank you, > Philipp > > _______________________________________________ > users mailing listusers at lists.scilab.orghttps://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > -- > St?phane Mottelet > Ing?nieur de recherche > EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable > D?partement G?nie des Proc?d?s Industriels > Sorbonne Universit?s - Universit? de Technologie de Compi?gne > CS 60319, 60203 Compi?gne cedex > Tel : +33(0)344234688http://www.utc.fr/~mottelet > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at utc.fr Fri Nov 20 11:46:50 2020 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Fri, 20 Nov 2020 11:46:50 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: Hi again, Le 20/11/2020 ? 11:30, P M a ?crit?: > OK, > some more background: > > I actually want to perform a delaunay triangulation on a set of > X-Y-Z-coordinates. I suppose you mean a triangulation of scattered X-Y-Z but you call mesh2d with X-Y only ? > > I know of CGLAB, but it seems only available to Scilab 6.0.x, while I > am on Scilab 6.1.0 CGLAB is based on a very old version of CGAL and thus is very hard to maintain as is. There is a lot of work... > > So the nearest thing to a delaunay triangulation I found is mesh2d... Yeah. Suboptimal choice indeed... > Because of the delaunay triangulation, I am not sure if the order of > points may change. What do you mean ? > > About unique: > > The result I get is actual of same size as the input matrix...I > therefore do assume that each point is unique in the first place > (which they actually should). > > What does repeat are individual x, y, z coordinates, but in different > combinations. > So the actual point is different. > > > e.g.:??? x may repeat, but than the y-z-coordinates are different from > point to point Could you share the data (you can send me a private message) ? S. > > Thanks, > Philipp > > > > Am Fr., 20. Nov. 2020 um 10:57?Uhr schrieb St?phane Mottelet > >: > > Hi Philipp, > > Glad to hear that some people use mesh2d. I know we have to > improve the error message by getting the numbers at the Fortran > level... There is also other errors which are not easy to fix > without knowing the number of the problematic point (point too > close to the boundary, for example). Can you file a bug on BZ ? > > S. > > Le 20/11/2020 ? 10:03, P M a ?crit?: >> Dear all, >> >> using the mesh2d - function I get an error: mesh2di: some points >> are identical. >> >> >> How to identify these points and remove the duplicates? >> >> There are roughly 500'000 coordinates.... >> >> Thank you, >> Philipp >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > -- > St?phane Mottelet > Ing?nieur de recherche > EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable > D?partement G?nie des Proc?d?s Industriels > Sorbonne Universit?s - Universit? de Technologie de Compi?gne > CS 60319, 60203 Compi?gne cedex > Tel : +33(0)344234688 > http://www.utc.fr/~mottelet > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Fri Nov 20 12:05:41 2020 From: p.muehlmann at gmail.com (P M) Date: Fri, 20 Nov 2020 12:05:41 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: Am Fr., 20. Nov. 2020 um 11:47 Uhr schrieb St?phane Mottelet < stephane.mottelet at utc.fr>: > Hi again, > Le 20/11/2020 ? 11:30, P M a ?crit : > > OK, > some more background: > > I actually want to perform a delaunay triangulation on a set of > X-Y-Z-coordinates. > > I suppose you mean a triangulation of scattered X-Y-Z but you call mesh2d > with X-Y only ? > Yes, indeed. And there is the bug in my data. As x-y coordinates repeat, they do have a different z in the original data. But since I call mesh2d with x-y only there are indeed identical points. > > I know of CGLAB, but it seems only available to Scilab 6.0.x, while I am > on Scilab 6.1.0 > > CGLAB is based on a very old version of CGAL and thus is very hard to > maintain as is. There is a lot of work... > I am currently installing Scilab 6.0 to try CGLAB... > > So the nearest thing to a delaunay triangulation I found is mesh2d... > > Yeah. Suboptimal choice indeed... > > Because of the delaunay triangulation, I am not sure if the order of > points may change. > > What do you mean ? > The scattered points represent a 3D surface geometry. If the points will be resorted, I could imagine that the triangulation does not fit to the surface geometry anymore. > > About unique: > > The result I get is actual of same size as the input matrix...I therefore > do assume that each point is unique in the first place (which they actually > should). > > What does repeat are individual x, y, z coordinates, but in different > combinations. > So the actual point is different. > > > e.g.: x may repeat, but than the y-z-coordinates are different from > point to point > > Could you share the data (you can send me a private message) ? > > S. > > > Thanks, > Philipp > > > > Am Fr., 20. Nov. 2020 um 10:57 Uhr schrieb St?phane Mottelet < > stephane.mottelet at utc.fr>: > >> Hi Philipp, >> >> Glad to hear that some people use mesh2d. I know we have to improve the >> error message by getting the numbers at the Fortran level... There is also >> other errors which are not easy to fix without knowing the number of the >> problematic point (point too close to the boundary, for example). Can you >> file a bug on BZ ? >> >> S. >> Le 20/11/2020 ? 10:03, P M a ?crit : >> >> Dear all, >> >> using the mesh2d - function I get an error: mesh2di: some points are >> identical. >> >> >> How to identify these points and remove the duplicates? >> >> There are roughly 500'000 coordinates.... >> >> Thank you, >> Philipp >> >> _______________________________________________ >> users mailing listusers at lists.scilab.orghttps://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users >> >> -- >> St?phane Mottelet >> Ing?nieur de recherche >> EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable >> D?partement G?nie des Proc?d?s Industriels >> Sorbonne Universit?s - Universit? de Technologie de Compi?gne >> CS 60319, 60203 Compi?gne cedex >> Tel : +33(0)344234688http://www.utc.fr/~mottelet >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> >> > > _______________________________________________ > users mailing listusers at lists.scilab.orghttps://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > -- > St?phane Mottelet > Ing?nieur de recherche > EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable > D?partement G?nie des Proc?d?s Industriels > Sorbonne Universit?s - Universit? de Technologie de Compi?gne > CS 60319, 60203 Compi?gne cedex > Tel : +33(0)344234688http://www.utc.fr/~mottelet > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at utc.fr Fri Nov 20 12:15:33 2020 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Fri, 20 Nov 2020 12:15:33 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: In order to use mesh2d with your data, there must exist a plane (u,v) on which you could project your (x,y,z) data. The mesh obtained on projected data (u,v) would give you the triangulation of the (x,z,y) data. This won't be possible with mesh2d if you take for example scattered points on a sphere. However, any non-linear transformation "unfolding" the data could also help to use pure 2d tools like mesh2d. S. Le 20/11/2020 ? 12:05, P M a ?crit?: > > > Am Fr., 20. Nov. 2020 um 11:47?Uhr schrieb St?phane Mottelet > >: > > Hi again, > > Le 20/11/2020 ? 11:30, P M a ?crit?: >> OK, >> some more background: >> >> I actually want to perform a delaunay triangulation on a set of >> X-Y-Z-coordinates. > I suppose you mean a triangulation of scattered X-Y-Z but you call > mesh2d with X-Y only ? > > > Yes, indeed. And there is the bug in my data. > As x-y coordinates repeat, they do have a different z in the original > data. > But since I call mesh2d with x-y only there are indeed identical points. > >> >> I know of CGLAB, but it seems only available to Scilab 6.0.x, >> while I am on Scilab 6.1.0 > CGLAB is based on a very old version of CGAL and thus is very > hard? to maintain as is. There is a lot of work... > > > I am currently installing Scilab 6.0 to try CGLAB... > >> >> So the nearest thing to a delaunay triangulation I found is mesh2d... > Yeah. Suboptimal choice indeed... >> Because of the delaunay triangulation, I am not sure if the order >> of points may change. > What do you mean ? > > > The scattered points represent a 3D surface geometry. > If the points will be resorted, I could imagine that the triangulation > does not fit to the surface geometry anymore. > >> >> About unique: >> >> The result I get is actual of same size as the input matrix...I >> therefore do assume that each point is unique in the first place >> (which they actually should). >> >> What does repeat are individual x, y, z coordinates, but in >> different combinations. >> So the actual point is different. >> >> >> e.g.:??? x may repeat, but than the y-z-coordinates are different >> from point to point > > Could you share the data (you can send me a private message) ? > > S. > >> >> Thanks, >> Philipp >> >> >> >> Am Fr., 20. Nov. 2020 um 10:57?Uhr schrieb St?phane Mottelet >> >: >> >> Hi Philipp, >> >> Glad to hear that some people use mesh2d. I know we have to >> improve the error message by getting the numbers at the >> Fortran level... There is also other errors which are not >> easy to fix without knowing the number of the problematic >> point (point too close to the boundary, for example). Can you >> file a bug on BZ ? >> >> S. >> >> Le 20/11/2020 ? 10:03, P M a ?crit?: >>> Dear all, >>> >>> using the mesh2d - function I get an error:???? mesh2di: >>> some points are identical. >>> >>> >>> How to identify these points and remove the duplicates? >>> >>> There are roughly 500'000 coordinates.... >>> >>> Thank you, >>> Philipp >>> >>> _______________________________________________ >>> users mailing list >>> users at lists.scilab.org >>> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users >> >> -- >> St?phane Mottelet >> Ing?nieur de recherche >> EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable >> D?partement G?nie des Proc?d?s Industriels >> Sorbonne Universit?s - Universit? de Technologie de Compi?gne >> CS 60319, 60203 Compi?gne cedex >> Tel : +33(0)344234688 >> http://www.utc.fr/~mottelet >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> >> >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > -- > St?phane Mottelet > Ing?nieur de recherche > EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable > D?partement G?nie des Proc?d?s Industriels > Sorbonne Universit?s - Universit? de Technologie de Compi?gne > CS 60319, 60203 Compi?gne cedex > Tel : +33(0)344234688 > http://www.utc.fr/~mottelet > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Fri Nov 20 13:05:36 2020 From: p.muehlmann at gmail.com (P M) Date: Fri, 20 Nov 2020 13:05:36 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: I guess I do have such an u-v-plane. The background of the background : The task is to map an image onto a 3D geometry. The image includes reference points, which correspond to known x-y-z coordinates. So the image plane is (my guess) the u-v-plane. BR Philipp Am Fr., 20. Nov. 2020 um 12:16 Uhr schrieb St?phane Mottelet < stephane.mottelet at utc.fr>: > In order to use mesh2d with your data, there must exist a plane (u,v) on > which you could project your (x,y,z) data. The mesh obtained on projected > data (u,v) would give you the triangulation of the (x,z,y) data. This won't > be possible with mesh2d if you take for example scattered points on a > sphere. However, any non-linear transformation "unfolding" the data could > also help to use pure 2d tools like mesh2d. > > S. > Le 20/11/2020 ? 12:05, P M a ?crit : > > > > Am Fr., 20. Nov. 2020 um 11:47 Uhr schrieb St?phane Mottelet < > stephane.mottelet at utc.fr>: > >> Hi again, >> Le 20/11/2020 ? 11:30, P M a ?crit : >> >> OK, >> some more background: >> >> I actually want to perform a delaunay triangulation on a set of >> X-Y-Z-coordinates. >> >> I suppose you mean a triangulation of scattered X-Y-Z but you call mesh2d >> with X-Y only ? >> > > Yes, indeed. And there is the bug in my data. > As x-y coordinates repeat, they do have a different z in the original data. > But since I call mesh2d with x-y only there are indeed identical points. > > > >> >> I know of CGLAB, but it seems only available to Scilab 6.0.x, while I am >> on Scilab 6.1.0 >> >> CGLAB is based on a very old version of CGAL and thus is very hard to >> maintain as is. There is a lot of work... >> > > I am currently installing Scilab 6.0 to try CGLAB... > >> >> So the nearest thing to a delaunay triangulation I found is mesh2d... >> >> Yeah. Suboptimal choice indeed... >> >> Because of the delaunay triangulation, I am not sure if the order of >> points may change. >> >> What do you mean ? >> > > The scattered points represent a 3D surface geometry. > If the points will be resorted, I could imagine that the triangulation > does not fit to the surface geometry anymore. > >> >> About unique: >> >> The result I get is actual of same size as the input matrix...I therefore >> do assume that each point is unique in the first place (which they actually >> should). >> >> What does repeat are individual x, y, z coordinates, but in different >> combinations. >> So the actual point is different. >> >> >> e.g.: x may repeat, but than the y-z-coordinates are different from >> point to point >> >> Could you share the data (you can send me a private message) ? >> >> S. >> >> >> Thanks, >> Philipp >> >> >> >> Am Fr., 20. Nov. 2020 um 10:57 Uhr schrieb St?phane Mottelet < >> stephane.mottelet at utc.fr>: >> >>> Hi Philipp, >>> >>> Glad to hear that some people use mesh2d. I know we have to improve the >>> error message by getting the numbers at the Fortran level... There is also >>> other errors which are not easy to fix without knowing the number of the >>> problematic point (point too close to the boundary, for example). Can you >>> file a bug on BZ ? >>> >>> S. >>> Le 20/11/2020 ? 10:03, P M a ?crit : >>> >>> Dear all, >>> >>> using the mesh2d - function I get an error: mesh2di: some points are >>> identical. >>> >>> >>> How to identify these points and remove the duplicates? >>> >>> There are roughly 500'000 coordinates.... >>> >>> Thank you, >>> Philipp >>> >>> _______________________________________________ >>> users mailing listusers at lists.scilab.orghttps://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users >>> >>> -- >>> St?phane Mottelet >>> Ing?nieur de recherche >>> EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable >>> D?partement G?nie des Proc?d?s Industriels >>> Sorbonne Universit?s - Universit? de Technologie de Compi?gne >>> CS 60319, 60203 Compi?gne cedex >>> Tel : +33(0)344234688http://www.utc.fr/~mottelet >>> >>> _______________________________________________ >>> users mailing list >>> users at lists.scilab.org >>> http://lists.scilab.org/mailman/listinfo/users >>> >>> >> >> _______________________________________________ >> users mailing listusers at lists.scilab.orghttps://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users >> >> -- >> St?phane Mottelet >> Ing?nieur de recherche >> EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable >> D?partement G?nie des Proc?d?s Industriels >> Sorbonne Universit?s - Universit? de Technologie de Compi?gne >> CS 60319, 60203 Compi?gne cedex >> Tel : +33(0)344234688http://www.utc.fr/~mottelet >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> >> > > _______________________________________________ > users mailing listusers at lists.scilab.orghttps://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > -- > St?phane Mottelet > Ing?nieur de recherche > EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable > D?partement G?nie des Proc?d?s Industriels > Sorbonne Universit?s - Universit? de Technologie de Compi?gne > CS 60319, 60203 Compi?gne cedex > Tel : +33(0)344234688http://www.utc.fr/~mottelet > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at utc.fr Fri Nov 20 13:22:32 2020 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Fri, 20 Nov 2020 13:22:32 +0100 Subject: [Scilab-users] mesh2d - finding identical points In-Reply-To: References: Message-ID: <63e5e74c-39b1-9f38-3336-cc3c225c8cfb@utc.fr> Yes, just mesh the set of (u,v) coordinates corresponding to your (x,y,z) reference points. S. Le 20/11/2020 ? 13:05, P M a ?crit?: > I guess I do have such an u-v-plane. > > The background of the background : > > The task is to map an image onto a 3D geometry. > > The image includes reference points, which correspond to known x-y-z > coordinates. > > So the image plane is (my guess) the u-v-plane. > > BR > Philipp > > > > > Am Fr., 20. Nov. 2020 um 12:16?Uhr schrieb St?phane Mottelet > >: > > In order to use mesh2d with your data, there must exist a plane > (u,v) on which you could project your (x,y,z) data. The mesh > obtained on projected data (u,v) would give you the triangulation > of the (x,z,y) data. This won't be possible with mesh2d if you > take for example scattered points on a sphere. However, any > non-linear transformation "unfolding" the data could also help to > use pure 2d tools like mesh2d. > > S. > > Le 20/11/2020 ? 12:05, P M a ?crit?: >> >> >> Am Fr., 20. Nov. 2020 um 11:47?Uhr schrieb St?phane Mottelet >> >: >> >> Hi again, >> >> Le 20/11/2020 ? 11:30, P M a ?crit?: >>> OK, >>> some more background: >>> >>> I actually want to perform a delaunay triangulation on a set >>> of X-Y-Z-coordinates. >> I suppose you mean a triangulation of scattered X-Y-Z but you >> call mesh2d with X-Y only ? >> >> >> Yes, indeed. And there is the bug in my data. >> As x-y coordinates repeat, they do have a different z in the >> original data. >> But since I call mesh2d with x-y only there are indeed identical >> points. >> >>> >>> I know of CGLAB, but it seems only available to Scilab >>> 6.0.x, while I am on Scilab 6.1.0 >> CGLAB is based on a very old version of CGAL and thus is very >> hard? to maintain as is. There is a lot of work... >> >> >> I am currently installing Scilab 6.0 to try CGLAB... >> >>> >>> So the nearest thing to a delaunay triangulation I found is >>> mesh2d... >> Yeah. Suboptimal choice indeed... >>> Because of the delaunay triangulation, I am not sure if the >>> order of points may change. >> What do you mean ? >> >> >> The scattered points represent a 3D surface geometry. >> If the points will be resorted, I could imagine that the >> triangulation does not fit to the surface geometry anymore. >> >>> >>> About unique: >>> >>> The result I get is actual of same size as the input >>> matrix...I therefore do assume that each point is unique in >>> the first place (which they actually should). >>> >>> What does repeat are individual x, y, z coordinates, but in >>> different combinations. >>> So the actual point is different. >>> >>> >>> e.g.:??? x may repeat, but than the y-z-coordinates are >>> different from point to point >> >> Could you share the data (you can send me a private message) ? >> >> S. >> >>> >>> Thanks, >>> Philipp >>> >>> >>> >>> Am Fr., 20. Nov. 2020 um 10:57?Uhr schrieb St?phane Mottelet >>> >: >>> >>> Hi Philipp, >>> >>> Glad to hear that some people use mesh2d. I know we have >>> to improve the error message by getting the numbers at >>> the Fortran level... There is also other errors which >>> are not easy to fix without knowing the number of the >>> problematic point (point too close to the boundary, for >>> example). Can you file a bug on BZ ? >>> >>> S. >>> >>> Le 20/11/2020 ? 10:03, P M a ?crit?: >>>> Dear all, >>>> >>>> using the mesh2d - function I get an error:???? >>>> mesh2di: some points are identical. >>>> >>>> >>>> How to identify these points and remove the duplicates? >>>> >>>> There are roughly 500'000 coordinates.... >>>> >>>> Thank you, >>>> Philipp >>>> >>>> _______________________________________________ >>>> users mailing list >>>> users at lists.scilab.org >>>> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users >>> >>> -- >>> St?phane Mottelet >>> Ing?nieur de recherche >>> EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable >>> D?partement G?nie des Proc?d?s Industriels >>> Sorbonne Universit?s - Universit? de Technologie de Compi?gne >>> CS 60319, 60203 Compi?gne cedex >>> Tel : +33(0)344234688 >>> http://www.utc.fr/~mottelet >>> >>> _______________________________________________ >>> users mailing list >>> users at lists.scilab.org >>> http://lists.scilab.org/mailman/listinfo/users >>> >>> >>> >>> _______________________________________________ >>> users mailing list >>> users at lists.scilab.org >>> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users >> >> -- >> St?phane Mottelet >> Ing?nieur de recherche >> EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable >> D?partement G?nie des Proc?d?s Industriels >> Sorbonne Universit?s - Universit? de Technologie de Compi?gne >> CS 60319, 60203 Compi?gne cedex >> Tel : +33(0)344234688 >> http://www.utc.fr/~mottelet >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> >> >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > -- > St?phane Mottelet > Ing?nieur de recherche > EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable > D?partement G?nie des Proc?d?s Industriels > Sorbonne Universit?s - Universit? de Technologie de Compi?gne > CS 60319, 60203 Compi?gne cedex > Tel : +33(0)344234688 > http://www.utc.fr/~mottelet > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul18fr at gmail.com Fri Nov 20 14:48:46 2020 From: paul18fr at gmail.com (paul francedixhuit) Date: Fri, 20 Nov 2020 14:48:46 +0100 Subject: [Scilab-users] String and integer Message-ID: Hi all In ordre to change files names, I need to convert integers into strings but with zeros prior to it Example if i=1 i need 0001 and not 1 I don't know how to use something like "%4d" with "string" Thanks Paul -------------- next part -------------- An HTML attachment was scrubbed... URL: From Jean-Yves.Baudais at insa-rennes.fr Fri Nov 20 15:08:48 2020 From: Jean-Yves.Baudais at insa-rennes.fr (Jean-Yves Baudais) Date: Fri, 20 Nov 2020 15:08:48 +0100 (CET) Subject: [Scilab-users] String and integer In-Reply-To: References: Message-ID: <1989424254.375799.1605881328060.JavaMail.zimbra@insa-rennes.fr> Hi, A solution, of course not the best, --> if i<10, s=msprintf("000%d",i) end with iteration for high number, or with 3-floor(log10(i)) added zeros? --Jean-Yves ----- Original Message ----- > From: "paul francedixhuit" > To: "Users mailing list for Scilab" > Sent: Friday, 20 November, 2020 14:48:46 > Subject: [Scilab-users] String and integer > Hi all > > In ordre to change files names, I need to convert integers into strings but with > zeros prior to it > > Example if i=1 i need 0001 and not 1 > > I don't know how to use something like "%4d" with "string" > > Thanks > > Paul > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- -- Jean-Yves Baudais From denis.crete at thalesgroup.com Fri Nov 20 15:29:36 2020 From: denis.crete at thalesgroup.com (CRETE Denis) Date: Fri, 20 Nov 2020 14:29:36 +0000 Subject: [Scilab-users] String and integer In-Reply-To: References: Message-ID: <2d5d3fc1a6e546bc9c27863c595e986b@thalesgroup.com> Hello msprintf(?%04d?,1) seems to do the job HTH Denis De : users De la part de paul francedixhuit Envoy? : vendredi 20 novembre 2020 14:49 ? : Users mailing list for Scilab Objet : [Scilab-users] String and integer Hi all In ordre to change files names, I need to convert integers into strings but with zeros prior to it Example if i=1 i need 0001 and not 1 I don't know how to use something like "%4d" with "string" Thanks Paul -------------- next part -------------- An HTML attachment was scrubbed... URL: From david.cheze at cea.fr Fri Nov 20 15:01:25 2020 From: david.cheze at cea.fr (CHEZE David 227480) Date: Fri, 20 Nov 2020 14:01:25 +0000 Subject: [Scilab-users] String and integer In-Reply-To: References: Message-ID: Hi Paul, --> msprintf("%04d",3) ans = "0003" David De : users De la part de paul francedixhuit Envoy? : vendredi 20 novembre 2020 14:49 ? : Users mailing list for Scilab Objet : [Scilab-users] String and integer Hi all In ordre to change files names, I need to convert integers into strings but with zeros prior to it Example if i=1 i need 0001 and not 1 I don't know how to use something like "%4d" with "string" Thanks Paul -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul18fr at gmail.com Fri Nov 20 17:06:17 2020 From: paul18fr at gmail.com (paul francedixhuit) Date: Fri, 20 Nov 2020 17:06:17 +0100 Subject: [Scilab-users] String and integer In-Reply-To: References: Message-ID: Based on your advices, I've fix my issue. Thanks Paul ############################################### i = 1FileName = "MyFile"// PreviouslyNewFileName1 = FileName + string(i)// NowNewFileName2 = FileName + msprintf("%.4d",i)i = 666NewFileName3 = FileName + msprintf("%.4d",i) Le ven. 20 nov. 2020 ? 15:43, CHEZE David 227480 a ?crit : > Hi Paul, > > > > --> msprintf("%04d",3) > > ans = > > > > "0003" > > > > > > > > David > > > > *De :* users *De la part de* paul > francedixhuit > *Envoy? :* vendredi 20 novembre 2020 14:49 > *? :* Users mailing list for Scilab > *Objet :* [Scilab-users] String and integer > > > > Hi all > > > > In ordre to change files names, I need to convert integers into strings > but with zeros prior to it > > > > Example if i=1 i need 0001 and not 1 > > > > I don't know how to use something like "%4d" with "string" > > > > Thanks > > > > Paul > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Sat Nov 21 22:55:23 2020 From: p.muehlmann at gmail.com (P M) Date: Sat, 21 Nov 2020 22:55:23 +0100 Subject: [Scilab-users] converting octave script imremap to scilab Message-ID: Dear All, I am trying to map an image onto a 3D surface. The image contains reference points, which with known 2D pixel coordinates. The 3D equivalent of these reference points is as well known. The 3D shape is represented by a set of X-Y-Z coordinates. For mapping the image onto a delaunay-triangulated surface I found following function: "imremap"...which - as far as I understand - should do what I want. now: "imremap" uses "interp2"...and unfortunately both functions are in octave. imremap: https://searchcode.com/codesearch/view/9585363/ interp2: https://searchcode.com/codesearch/view/20327397/ Converting the functions to Scilab seems to become an elephant, especially for interp2. If one could avoid "interp2" , then "imremap" could be possible relatively easy converted to Scilab. Any ideas which Scilab function would work like interp2? Thank you, Philipp -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at utc.fr Sun Nov 22 08:55:07 2020 From: stephane.mottelet at utc.fr (=?utf-8?Q?St=C3=A9phane_Mottelet?=) Date: Sun, 22 Nov 2020 08:55:07 +0100 Subject: [Scilab-users] converting octave script imremap to scilab In-Reply-To: References: Message-ID: Hi, Does interp2 works by using the triangulation ? S. > Le 21 nov. 2020 ? 22:56, P M a ?crit : > > ? > Dear All, > > I am trying to map an image onto a 3D surface. > The image contains reference points, which with known 2D pixel coordinates. > The 3D equivalent of these reference points is as well known. > The 3D shape is represented by a set of X-Y-Z coordinates. > > For mapping the image onto a delaunay-triangulated surface I found following function: > > "imremap"...which - as far as I understand - should do what I want. > > now: "imremap" uses "interp2"...and unfortunately both functions are in octave. > > > imremap: https://searchcode.com/codesearch/view/9585363/ > > interp2: https://searchcode.com/codesearch/view/20327397/ > > > Converting the functions to Scilab seems to become an elephant, especially for interp2. > > If one could avoid "interp2" , then "imremap" could be possible relatively easy converted to Scilab. > > Any ideas which Scilab function would work like interp2? > > Thank you, > Philipp > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Sun Nov 22 22:48:29 2020 From: p.muehlmann at gmail.com (P M) Date: Sun, 22 Nov 2020 22:48:29 +0100 Subject: [Scilab-users] converting octave script imremap to scilab In-Reply-To: References: Message-ID: if you refer to interp2 from CGLAB toolbox....then no, I have no success with that. CGLAB call: vp = interp2(xp,yp,x,y,v)...were xp and yp are vectors while at octave: warped = grayinterp(im, XI, YI, interp, NA); // im seems to be a mxn matrix with:function [warped, valid] = grayinterp(im, XI, YI, interp, extrapval) if (strcmp(interp, "cubic")) warped = graybicubic(double(im), XI, YI, NA); else *warped = interp2(double(im), XI, YI, interp, NA);* endif valid = !isna(warped); warped(!valid) = extrapval;endfunction BR Philipp Am So., 22. Nov. 2020 um 08:55 Uhr schrieb St?phane Mottelet < stephane.mottelet at utc.fr>: > Hi, > > Does interp2 works by using the triangulation ? > > S. > > Le 21 nov. 2020 ? 22:56, P M a ?crit : > > ? > Dear All, > > I am trying to map an image onto a 3D surface. > The image contains reference points, which with known 2D pixel coordinates. > The 3D equivalent of these reference points is as well known. > The 3D shape is represented by a set of X-Y-Z coordinates. > > For mapping the image onto a delaunay-triangulated surface I found > following function: > > "imremap"...which - as far as I understand - should do what I want. > > now: "imremap" uses "interp2"...and unfortunately both functions are in > octave. > > > imremap: https://searchcode.com/codesearch/view/9585363/ > > > interp2: https://searchcode.com/codesearch/view/20327397/ > > > > Converting the functions to Scilab seems to become an elephant, especially > for interp2. > > If one could avoid "interp2" , then "imremap" could be possible > relatively easy converted to Scilab. > > Any ideas which Scilab function would work like interp2? > > Thank you, > Philipp > _______________________________________________ > users mailing list > users at lists.scilab.org > > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Mon Nov 23 09:08:49 2020 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Mon, 23 Nov 2020 08:08:49 +0000 Subject: [Scilab-users] {EXT} converting octave script imremap to scilab In-Reply-To: References: Message-ID: Hello Philipp, > De : users De la part de P M Envoy? : > samedi 21 novembre 2020 22:55 > > I am trying to map an image onto a 3D surface. If you are not "bound" to Scilab, then I suggest to have a look at Kitware Paraview https://www.paraview.org/ which can, as far as I remember, can import and export data before/after processing. Regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer General This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. From d.daloisio at enginsoft.com Mon Nov 23 17:12:30 2020 From: d.daloisio at enginsoft.com (Davide Daloisio) Date: Mon, 23 Nov 2020 17:12:30 +0100 Subject: [Scilab-users] CSV Building Message-ID: Good evening, I'm building a script to be able to create the matrix of two vectors that have different timesteps at the start. So I am writing a new matrix with two vectors that have the same timestep. In this specific case I have a time vector with a timestep (step) of about 2E-4 and then I want to create another vector with a constant speed (vel) equal to 2. To then couple them into a new matrix. When I run the script it gives me an error on line 7. Could someone give me a tip? [image: image.png] Davide Daloisio *Application Engineer * *SPDM - Simulation Process and Data Management* [image: EnginSoft] *LinkedIn*: https://www.linkedin.com/in/davidedaloisio/ *Office*: +39 0461 979338 *Mobile*: +39 347 8769382 *e-mail* : d.daloisio at enginsoft.com www.enginsoft.com www.caeconference.com ------------------------------------------ Please, consider your environmental responsibility! Before printing this e-mail ask yourself: "Do I need a hard copy?" -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image.png Type: image/png Size: 33795 bytes Desc: not available URL: From Christophe.Dang at sidel.com Tue Nov 24 09:35:01 2020 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Tue, 24 Nov 2020 08:35:01 +0000 Subject: [Scilab-users] CSV Building Message-ID: Hello Davide, > De : users De la part de Davide > Daloisio Envoy? : lundi 23 novembre 2020 17:13 > > So I am writing a new matrix with two vectors that have the same timestep. > In this specific case I have a time vector with a timestep (step) of > about 2E- When you're writing: A = B(C) Where B is a matrix, then C must be a matrix of Booleans (the cells you keep or reject) or of integers (the indices). So the code t_vel(1:step:l_vel) is not valid as step is not an integer. I'm not sure what you're trying to do but I guess you should try something like: new_t_vel = t_vel(1) + step*(0:l_vel); Hope this helps -- Christophe Dang Ngoc Chan Mechanical calculation engineer General This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. From p.muehlmann at gmail.com Tue Nov 24 13:55:36 2020 From: p.muehlmann at gmail.com (P M) Date: Tue, 24 Nov 2020 13:55:36 +0100 Subject: [Scilab-users] image insertion using scicv Message-ID: Dear Experts, I would like to implement a small image into a part of a big image. I would know how to do that with IPCV, but the small image is loaded using scicv "imread". The help to "Mat" states: - insertion: syntax is the same a the usual Scilab insertion, and accepts a Scilab matrix in input (not yet available). - concatenation: syntax is the same a the usual Scilab concatenation (not yet available). Is there any progress to make this available? Thank you, Philipp -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Tue Nov 24 14:15:44 2020 From: p.muehlmann at gmail.com (P M) Date: Tue, 24 Nov 2020 14:15:44 +0100 Subject: [Scilab-users] ipcv vs scicv Message-ID: Dear developers, as I have many scripts with "imread" I wonder if some time in the future duplicate functions names from IPCV and SCICV are going to be removed/changed. It's kind of disturbing to have to use atomsInstall/atomsRemove when using different scripts. I think this was discussed already in the past, though an update would be nice. Thank you, Philipp -------------- next part -------------- An HTML attachment was scrubbed... URL: From d.daloisio at enginsoft.com Tue Nov 24 14:51:23 2020 From: d.daloisio at enginsoft.com (Davide Daloisio) Date: Tue, 24 Nov 2020 14:51:23 +0100 Subject: [Scilab-users] CSV Building In-Reply-To: References: Message-ID: Thanks Christophe for your suggestion. I managed to solve. But later this error message appeared: exec('D:\users\A_DLO_AreaTecnica\A_DLO_R&D\R&D\REM\ROMCRF\to_DLO\Input\first_opt\script_interpolazione.sci', -1) at line 8 of executed file D:\users\A_DLO_AreaTecnica\A_DLO_R&D\R&D\REM\ROMCRF\to_DLO\Input\first_opt\script_interpolazione.sci inconsistent row/column dimensions Davide Daloisio *Application Engineer * *SPDM - Simulation Process and Data Management* [image: EnginSoft] *LinkedIn*: https://www.linkedin.com/in/davidedaloisio/ *Office*: +39 0461 979338 *Mobile*: +39 347 8769382 *e-mail* : d.daloisio at enginsoft.com www.enginsoft.com www.caeconference.com ------------------------------------------ Please, consider your environmental responsibility! Before printing this e-mail ask yourself: "Do I need a hard copy?" Il giorno mar 24 nov 2020 alle ore 09:35 Dang Ngoc Chan, Christophe < Christophe.Dang at sidel.com> ha scritto: > Hello Davide, > > > De : users De la part de Davide > > Daloisio Envoy? : lundi 23 novembre 2020 17:13 > > > > So I am writing a new matrix with two vectors that have the same > timestep. > > In this specific case I have a time vector with a timestep (step) of > > about 2E- > > When you're writing: > > A = B(C) > > Where B is a matrix, then C must be a matrix of Booleans (the cells you > keep or reject) or of integers (the indices). > > So the code > > t_vel(1:step:l_vel) > > is not valid as step is not an integer. > > I'm not sure what you're trying to do but I guess you should try something > like: > > new_t_vel = t_vel(1) + step*(0:l_vel); > > Hope this helps > > -- > Christophe Dang Ngoc Chan > Mechanical calculation engineer > > General > This e-mail may contain confidential and/or privileged information. If you > are not the intended recipient (or have received this e-mail in error), > please notify the sender immediately and destroy this e-mail. Any > unauthorized copying, disclosure or distribution of the material in this > e-mail is strictly forbidden. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Tue Nov 24 15:27:21 2020 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Tue, 24 Nov 2020 14:27:21 +0000 Subject: [Scilab-users] {EXT} Re: CSV Building In-Reply-To: References: Message-ID: Hello, > De : users De la part de Davide > Daloisio Envoy? : mardi 24 novembre 2020 14:51 > > inconsistent row/column dimensions Difficult to say without the code. You usually have this when you multiply matrices (the number of rows of the first one must match the number of columns of the second) or the like. -- Christophe Dang Ngoc Chan Mechanical calculation engineer General This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. From cfuttrup at gmail.com Tue Nov 24 17:22:50 2020 From: cfuttrup at gmail.com (Claus Futtrup) Date: Tue, 24 Nov 2020 17:22:50 +0100 Subject: [Scilab-users] ipcv vs scicv In-Reply-To: References: Message-ID: Hi all I imagine naming conflicts can occur across different ATOMS. How about a change where one has to specify which ATOMS module you're using when calling a function, such that both can coexist side-by-side ? I'm thinking Pythonic here, like if you import for example numpy, you can write import numpy as np ... from hereon, you write np.pi for the definition of pi (=3.14159 ... ), and so on. Maybe something for Scilab 7 (?) Cheers, Claus On 24-11-2020 14:15, P M wrote: > Dear developers, > > as I have many scripts with "imread" I wonder if some time in the > future duplicate functions names from IPCV and SCICV are going to be > removed/changed. > > It's kind of disturbing to have to use atomsInstall/atomsRemove when > using different scripts. > > I think this was discussed already in the past, though an update would > be nice. > > Thank you, > Philipp > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Tue Nov 24 19:46:37 2020 From: p.muehlmann at gmail.com (P M) Date: Tue, 24 Nov 2020 19:46:37 +0100 Subject: [Scilab-users] ipcv vs scicv In-Reply-To: References: Message-ID: Mh. Nice Idea. Maybe something like this is already possible? Install toolboxed but remove from autoloader. So scilab starts without toolboxed from atoms. The script includes than the command, which module shall be loaded. Thanks for the idea...i will try this. Philipp Am Dienstag, 24. November 2020 schrieb Claus Futtrup : > Hi all > > I imagine naming conflicts can occur across different ATOMS. > > How about a change where one has to specify which ATOMS module you're > using when calling a function, such that both can coexist side-by-side ? > > I'm thinking Pythonic here, like if you import for example numpy, you can > write import numpy as np ... from hereon, you write np.pi for the > definition of pi (=3.14159 ... ), and so on. > > Maybe something for Scilab 7 (?) > > Cheers, > Claus > > On 24-11-2020 14:15, P M wrote: > > Dear developers, > > as I have many scripts with "imread" I wonder if some time in the future > duplicate functions names from IPCV and SCICV are going to be > removed/changed. > > It's kind of disturbing to have to use atomsInstall/atomsRemove when using > different scripts. > > I think this was discussed already in the past, though an update would be > nice. > > Thank you, > Philipp > > > _______________________________________________ > users mailing listusers at lists.scilab.orghttp://lists.scilab.org/mailman/listinfo/users > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From p.muehlmann at gmail.com Wed Nov 25 04:10:17 2020 From: p.muehlmann at gmail.com (P M) Date: Wed, 25 Nov 2020 04:10:17 +0100 Subject: [Scilab-users] scicv - camera calibration Message-ID: Dear scicv developers, is camera calibration possible with scicv? Help states: The following feature domains are (or will be) covered: [...] - Camera calibration (calib3d) [...] I also do find a "opencv_calib3d2413.dll" somewhere in the scicv installation folders. However I do not find any scilab function pointing to the topic of camera calibration. Thank you, Philipp -------------- next part -------------- An HTML attachment was scrubbed... URL: From liviu9badea at gmail.com Thu Nov 26 14:22:26 2020 From: liviu9badea at gmail.com (Liviu Badea) Date: Thu, 26 Nov 2020 15:22:26 +0200 Subject: [Scilab-users] thermal barrier Message-ID: Hello, Please help me with this problem. I have a thermal barrier that has five layers I want to do in Scilab the following. Heat transfer; the stationary linear case to graphically represent the conductivity variation with the thickness of the layers and with the temperature. ReplyForward Hello, Please help me with this problem. I have a thermal barrier that has five layers I want to do in Scilab the following. Heat transfer; the stationary linear case to graphically represent the conductivity variation with the thickness of the layers and with the temperature. ReplyForward Compose: New Message [image: Minimize][image: Pop-out][image: Close] Compose: New Message [image: Minimize][image: Pop-out][image: Close] -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Thu Nov 26 15:11:27 2020 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Thu, 26 Nov 2020 14:11:27 +0000 Subject: [Scilab-users] {EXT} thermal barrier In-Reply-To: References: Message-ID: Hello, > De : users De la part de Liviu Badea > Envoy? : jeudi 26 novembre 2020 14:22 > > Heat transfer; the stationary linear case to graphically represent the > conductivity variation with the thickness of the layers and with the temperature. We would be happy to help you if you tell us exactly which problems you have. For the coding of the equation, you have to define a function: https://help.scilab.org/docs/6.1.0/en_US/function.html you'll have to define vectors of values for the thickness and temperature, see: https://help.scilab.org/docs/6.1.0/en_US/colon.html https://help.scilab.org/docs/6.1.0/en_US/linspace.html For the plot: https://help.scilab.org/docs/6.1.0/en_US/plot2d.html https://help.scilab.org/docs/6.1.0/en_US/plot3d.html Without further information, we can't help more. Please send us, if possible, the code you wrote and the errors that arise. Regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer General This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. From cfuttrup at gmail.com Thu Nov 26 19:19:35 2020 From: cfuttrup at gmail.com (Claus Futtrup) Date: Thu, 26 Nov 2020 19:19:35 +0100 Subject: [Scilab-users] fft Message-ID: Hi Scilabers I'm trying to understand how to use the FFT function in Scilab. I'm using v. 6.1.0 on a Windows 10 machine. This crashes Scilab: y = linspace(1,256,256) x = fft(y) BTW, the example in the online help (https://help.scilab.org/docs/6.1.0/en_US/fft.html) also crashes Scilab: sample_rate=1000; t = 0:1/sample_rate:0.6; N=size (t,'*'); //number of samples s=sin (2*%pi*50*t)+sin (2*%pi*70*t+%pi/4)+grand (1,N,'nor',0,1); y=fft (s); Cheers, Claus -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Thu Nov 26 19:24:07 2020 From: cfuttrup at gmail.com (Claus Futtrup) Date: Thu, 26 Nov 2020 19:24:07 +0100 Subject: [Scilab-users] fft In-Reply-To: References: Message-ID: <1ec27ad7-32db-c9a6-766e-31a4b6653de0@gmail.com> Hi there The 2D example from the help pages also crashes Scilab. Maybe my Scilab installation is broken - the FFT part? I can't remember, but I think I installed the Intel Math Kernel with everything. My computer is a Lenovo T480 with Intel Core i7 processor. I have the following installed: C:\Program Files\scilab-6.1.0>dir /S /B fft*.* C:\Program Files\scilab-6.1.0\bin\fftw C:\Program Files\scilab-6.1.0\bin\fftw.dll C:\Program Files\scilab-6.1.0\bin\fftw\fftw-readme.txt C:\Program Files\scilab-6.1.0\modules\fftw C:\Program Files\scilab-6.1.0\modules\fftw\etc\fftw.quit C:\Program Files\scilab-6.1.0\modules\fftw\etc\fftw.start Cheers, Claus On 26-11-2020 19:19, Claus Futtrup wrote: > > Hi Scilabers > > I'm trying to understand how to use the FFT function in Scilab. I'm > using v. 6.1.0 on a Windows 10 machine. > > This crashes Scilab: > y = linspace(1,256,256) x = fft(y) > > BTW, the example in the online help (https://help.scilab.org/docs/6.1.0/en_US/fft.html) also crashes Scilab: > > sample_rate=1000; t = 0:1/sample_rate:0.6; N=size > (t,'*'); //number > of samples s=sin > (2*%pi*50*t)+sin > (2*%pi*70*t+%pi/4)+grand > (1,N,'nor',0,1); > y=fft (s); > Cheers, > Claus -------------- next part -------------- An HTML attachment was scrubbed... URL: From fmiyara at fceia.unr.edu.ar Thu Nov 26 23:17:08 2020 From: fmiyara at fceia.unr.edu.ar (Federico Miyara) Date: Thu, 26 Nov 2020 19:17:08 -0300 Subject: [Scilab-users] fft In-Reply-To: <1ec27ad7-32db-c9a6-766e-31a4b6653de0@gmail.com> References: <1ec27ad7-32db-c9a6-766e-31a4b6653de0@gmail.com> Message-ID: <2cf9e8cc-e39f-8527-be5d-4fa5e9181f4c@fceia.unr.edu.ar> Claus, I recall that I had had a similar problem in certain occasion. The reason was that when installing Scilab I had selected an option to use a version of the FFT which, since it was much larger, I thought would be more complete, or faster... but it didn't work. Try to re-install Scilab and choose the default version of the FFT engine. It is good enough... and works! Regards, Federico Miyara On 26/11/2020 15:24, Claus Futtrup wrote: > Hi there > > The 2D example from the help pages also crashes Scilab. > > Maybe my Scilab installation is broken - the FFT part? I can't > remember, but I think I installed the Intel Math Kernel with > everything. My computer is a Lenovo T480 with Intel Core i7 processor. > I have the following installed: > > C:\Program Files\scilab-6.1.0>dir /S /B fft*.* > C:\Program Files\scilab-6.1.0\bin\fftw > C:\Program Files\scilab-6.1.0\bin\fftw.dll > C:\Program Files\scilab-6.1.0\bin\fftw\fftw-readme.txt > C:\Program Files\scilab-6.1.0\modules\fftw > C:\Program Files\scilab-6.1.0\modules\fftw\etc\fftw.quit > C:\Program Files\scilab-6.1.0\modules\fftw\etc\fftw.start > > Cheers, > Claus > > On 26-11-2020 19:19, Claus Futtrup wrote: >> >> Hi Scilabers >> >> I'm trying to understand how to use the FFT function in Scilab. I'm >> using v. 6.1.0 on a Windows 10 machine. >> >> This crashes Scilab: >> y = linspace(1,256,256) x = fft(y) >> >> BTW, the example in the online help (https://help.scilab.org/docs/6.1.0/en_US/fft.html) also crashes Scilab: >> >> sample_rate=1000; t = 0:1/sample_rate:0.6; N=size >> (t,'*'); //number >> of samples s=sin >> (2*%pi*50*t)+sin >> (2*%pi*70*t+%pi/4)+grand >> (1,N,'nor',0,1); >> y=fft (s); >> Cheers, >> Claus > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- El software de antivirus Avast ha analizado este correo electr?nico en busca de virus. https://www.avast.com/antivirus -------------- next part -------------- An HTML attachment was scrubbed... URL: From antoine.monmayrant at laas.fr Fri Nov 27 09:42:59 2020 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Fri, 27 Nov 2020 09:42:59 +0100 Subject: [Scilab-users] pathsep() returns ':' under linux ? Message-ID: <3c3af885-7259-a357-1eba-04c58dd2ede5@laas.fr> Hello all, I just found this today, on a linux machine (scilab 6.1.0 ubuntu 18.04 LTS): --> pathsep() ?ans? = ? ":" it definitely sounds like a bug to me, no? As anyone else experience some issues with pathsep ? Cheers, Antoine -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at utc.fr Fri Nov 27 09:58:19 2020 From: stephane.mottelet at utc.fr (=?utf-8?Q?St=C3=A9phane_Mottelet?=) Date: Fri, 27 Nov 2020 09:58:19 +0100 Subject: [Scilab-users] pathsep() returns ':' under linux ? In-Reply-To: <3c3af885-7259-a357-1eba-04c58dd2ede5@laas.fr> References: <3c3af885-7259-a357-1eba-04c58dd2ede5@laas.fr> Message-ID: <22CF94B4-BBCD-4E8D-AC48-3EC91B121F84@utc.fr> It?s the char to separate paths in the PATH environment variable. Different from filesep(). S. > Le 27 nov. 2020 ? 09:55, Antoine Monmayrant a ?crit : > > ? > Hello all, > > I just found this today, on a linux machine (scilab 6.1.0 ubuntu 18.04 LTS): > > --> pathsep() > ans = > > ":" > > it definitely sounds like a bug to me, no? > > As anyone else experience some issues with pathsep ? > > Cheers, > > Antoine > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From antoine.monmayrant at laas.fr Fri Nov 27 11:00:34 2020 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Fri, 27 Nov 2020 11:00:34 +0100 Subject: [Scilab-users] pathsep() returns ':' under linux ? In-Reply-To: <22CF94B4-BBCD-4E8D-AC48-3EC91B121F84@utc.fr> References: <3c3af885-7259-a357-1eba-04c58dd2ede5@laas.fr> <22CF94B4-BBCD-4E8D-AC48-3EC91B121F84@utc.fr> Message-ID: <005a481d-be70-2cbd-5dab-e971017649b8@laas.fr> Thank you all, I mixed up filesep() & pathsep(). Antoine On 27/11/2020 09:58, St?phane Mottelet wrote: > It?s the char to separate paths in the PATH environment variable. > Different from filesep(). > > S. > >> Le 27 nov. 2020 ? 09:55, Antoine Monmayrant >> a ?crit?: >> >> ? >> >> Hello all, >> >> I just found this today, on a linux machine (scilab 6.1.0 ubuntu >> 18.04 LTS): >> >> --> pathsep() >> ?ans? = >> >> ? ":" >> >> it definitely sounds like a bug to me, no? >> >> As anyone else experience some issues with pathsep ? >> >> Cheers, >> >> Antoine >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Fri Nov 27 16:35:11 2020 From: cfuttrup at gmail.com (Claus Futtrup) Date: Fri, 27 Nov 2020 16:35:11 +0100 Subject: [Scilab-users] fft In-Reply-To: <2cf9e8cc-e39f-8527-be5d-4fa5e9181f4c@fceia.unr.edu.ar> References: <1ec27ad7-32db-c9a6-766e-31a4b6653de0@gmail.com> <2cf9e8cc-e39f-8527-be5d-4fa5e9181f4c@fceia.unr.edu.ar> Message-ID: Hi Federico Thank you for your advice. I guess reinstalling Scilab is what I need to do. Best regards, Claus On 26-11-2020 23:17, Federico Miyara wrote: > > Claus, > > I recall that I had had a similar problem in certain occasion. The > reason was that when installing Scilab I had selected an option to use > a version of the FFT which, since it was much larger, I thought would > be more complete, or faster... but it didn't work. Try to re-install > Scilab and choose the default version of the FFT engine. It is good > enough... and works! > > Regards, > > Federico Miyara > > On 26/11/2020 15:24, Claus Futtrup wrote: >> Hi there >> >> The 2D example from the help pages also crashes Scilab. >> >> Maybe my Scilab installation is broken - the FFT part? I can't >> remember, but I think I installed the Intel Math Kernel with >> everything. My computer is a Lenovo T480 with Intel Core i7 >> processor. I have the following installed: >> >> C:\Program Files\scilab-6.1.0>dir /S /B fft*.* >> C:\Program Files\scilab-6.1.0\bin\fftw >> C:\Program Files\scilab-6.1.0\bin\fftw.dll >> C:\Program Files\scilab-6.1.0\bin\fftw\fftw-readme.txt >> C:\Program Files\scilab-6.1.0\modules\fftw >> C:\Program Files\scilab-6.1.0\modules\fftw\etc\fftw.quit >> C:\Program Files\scilab-6.1.0\modules\fftw\etc\fftw.start >> >> Cheers, >> Claus >> >> On 26-11-2020 19:19, Claus Futtrup wrote: >>> >>> Hi Scilabers >>> >>> I'm trying to understand how to use the FFT function in Scilab. I'm >>> using v. 6.1.0 on a Windows 10 machine. >>> >>> This crashes Scilab: >>> y = linspace(1,256,256) x = fft(y) >>> >>> BTW, the example in the online help (https://help.scilab.org/docs/6.1.0/en_US/fft.html) also crashes Scilab: >>> >>> sample_rate=1000; t = 0:1/sample_rate:0.6; N=size >>> (t,'*'); >>> //number of samples s=sin >>> (2*%pi*50*t)+sin >>> (2*%pi*70*t+%pi/4)+grand >>> (1,N,'nor',0,1); >>> y=fft (s); >>> Cheers, >>> Claus >> >> >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > > > > Libre de virus. www.avast.com > > > > <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2> > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: