26 views

Skip to first unread message

Apr 16, 2023, 2:35:42 AMApr 16

to

Let us say we want the gradient of a Time series. For finite time series we

will be losing one data point at each end using most most methods. We

subtract each data point from its adjacent data point the one to either so

I'd let's say it's going to be the right side consistency consistently.

This gives us the time series which is the gradient of that original signal.

Let us say we want the gradient of a two-dimensional spatial series

otherwise known as an image. Buying analogous process I'm subtracting and

getting something like the gradient.

Symbolically the gradient is the partial derivative in each direction for

any number of dimensions.

Well I have been playing with the two dimensional spatial series data or

image format in mathcad and found a function called cfft which computes the

complex fast fourier transform of data in 2 dimensions. The icfft function

reverses that result to restore an image. Lowpass and high-pass filtering

are easy but not spectacularly easy in the cfft domain. They are

spectacular easy in the one-dimensional case.

Pixel data in two dimensions in space in grayscale our conventionally a

number of bits say 8 bits for 256 levels of Gray. And are analogous in

color spaces. Edge detection produces a result similar to an artist's

pencil drawing.

Well I collect pencil drawings and also pencil drawing software.

And I use mathcad.

I usually use the variable name capital m for matrix containing two

dimensional grayscale image data in mathcad. I compute the gradient several

different ways in method.

I have a new way to compute the gradient in mathcad and symbolically.

icfft (i * cfft (M))= grad(M).

Discussion?

Douglas

will be losing one data point at each end using most most methods. We

subtract each data point from its adjacent data point the one to either so

I'd let's say it's going to be the right side consistency consistently.

This gives us the time series which is the gradient of that original signal.

Let us say we want the gradient of a two-dimensional spatial series

otherwise known as an image. Buying analogous process I'm subtracting and

getting something like the gradient.

Symbolically the gradient is the partial derivative in each direction for

any number of dimensions.

Well I have been playing with the two dimensional spatial series data or

image format in mathcad and found a function called cfft which computes the

complex fast fourier transform of data in 2 dimensions. The icfft function

reverses that result to restore an image. Lowpass and high-pass filtering

are easy but not spectacularly easy in the cfft domain. They are

spectacular easy in the one-dimensional case.

Pixel data in two dimensions in space in grayscale our conventionally a

number of bits say 8 bits for 256 levels of Gray. And are analogous in

color spaces. Edge detection produces a result similar to an artist's

pencil drawing.

Well I collect pencil drawings and also pencil drawing software.

And I use mathcad.

I usually use the variable name capital m for matrix containing two

dimensional grayscale image data in mathcad. I compute the gradient several

different ways in method.

I have a new way to compute the gradient in mathcad and symbolically.

icfft (i * cfft (M))= grad(M).

Discussion?

Douglas

Apr 30, 2023, 3:46:53 PMApr 30

to

On Sunday, April 16, 2023 at 1:35:42=E2=80=AFAM UTC-5, Douglas Goncz

Physics and other fields.

The usual way to approach this is:

(1) model the point cloud as a function (that's interpolation),

(2) take the relevant derivatives of that function.

In ALGLIB, the "interpolation" module has the

"spline1d", "spline2d" and "spline3d" packages.

The spline 2D includes routines to initialize splines from point clouds

(both point grids, and with more recent revisions, unordered point clouds).

spline2dbuild*()

and to configure the splines

spline2dset*()

to calculuate interpolated values off the spline (singly or in batch)

spline2dcalc*()

to do affine transformations on both the input and output spaces

(including resampling)

spline2dlintrans*(), spline2dresample*()

and (tada!) to differentiate

spline2ddiff*().

The RBF module is closer to AI-level modelling and the original

ALGLIB developer has been harping on it, from day one, as a better

method for interpolation; but I can't say too much about it.

A curated version of the C++ free edition of ALGLIB (and a work in

progress) can be found on GitHub, currently under

https://github.com/LydiaMarieWilliamson/ALGLIB_cpp

which may be later integrated with, and into, curated versions Mesa

(translated into C++), OpenCV (and other image 2D to 3D and 3D modelling

software) as well as serving, for AI application, as an engine for

revised versions of LLaMA or GPT.

[[Mod. note -- Other good sources of free interpolation code include

https://www.netlib.org

https://www.gnu.org/software/gsl/doc/html/interp.html

-- jt]]

A.A.S. M.E.T. 1990 wrote:

> Let us say we want the gradient of a Time series.

[...]
> Let us say we want the gradient of a Time series.

> Let us say we want the gradient of a two-dimensional spatial series

> otherwise known as an image.

Well-known problems for scientific and numeric programming for
> otherwise known as an image.

Physics and other fields.

The usual way to approach this is:

(1) model the point cloud as a function (that's interpolation),

(2) take the relevant derivatives of that function.

In ALGLIB, the "interpolation" module has the

"spline1d", "spline2d" and "spline3d" packages.

The spline 2D includes routines to initialize splines from point clouds

(both point grids, and with more recent revisions, unordered point clouds).

spline2dbuild*()

and to configure the splines

spline2dset*()

to calculuate interpolated values off the spline (singly or in batch)

spline2dcalc*()

to do affine transformations on both the input and output spaces

(including resampling)

spline2dlintrans*(), spline2dresample*()

and (tada!) to differentiate

spline2ddiff*().

The RBF module is closer to AI-level modelling and the original

ALGLIB developer has been harping on it, from day one, as a better

method for interpolation; but I can't say too much about it.

A curated version of the C++ free edition of ALGLIB (and a work in

progress) can be found on GitHub, currently under

https://github.com/LydiaMarieWilliamson/ALGLIB_cpp

which may be later integrated with, and into, curated versions Mesa

(translated into C++), OpenCV (and other image 2D to 3D and 3D modelling

software) as well as serving, for AI application, as an engine for

revised versions of LLaMA or GPT.

[[Mod. note -- Other good sources of free interpolation code include

https://www.netlib.org

https://www.gnu.org/software/gsl/doc/html/interp.html

-- jt]]

Reply all

Reply to author

Forward

0 new messages

Search

Clear search

Close search

Google apps

Main menu