Complex number matrix in OpenCV (Java)
Hello
I am working on image data and the end goal is to use fourier transform on a product of the data and some complex matrix.
CV_64FC1 1000x1000 matrix of grayscaled image data (I have this)
Secondary 1000x1000 matrix of complex numbers (I don't have this)
The first matrix will multiply with the second complex matrix and after that I will use a fourier transform on it
Is there a way of doing this in OpenCV?
edit:
// image matrix
Bitmap bitmap_image = ...
Mat matrix = new Mat(1000, 1000, CvType.CV_64FC1);
Utils.bitmapToMat(bitmap_image, matrix);
Imgproc.cvtColor(matrix, matrix, Imgproc.COLOR_RGBA2GRAY);
matrix.convertTo(matrix, CvType.CV_64FC1, 1.0/255.0); // type when initialised is ignored
// secondary matrix
Mat secondary_matrix = new Mat(1000, 1000, CvType.CV_64FC1);
for (int i=0; i<1000; i++) {
for (int j=0; j<1000; j++) {
// ... do some stuff to get value
secondary matrix.put(i,j,value);
}
}
// now need to make secondary matrix complex by multiplying by some function of i
edit2:
Context: I have a hologram image in OpenCV matrix form (single channel), I have another matrix filled with real numbers (single channel) - this matrix represents the sampling period in the spatial domain. This latter matrix isn't used directly, it passes through an e^(i*...) function to make the values complex. Finally, a discrete fourier transform is applied on the per-element product of the image data and this complex matrix I just described.
please show your coding attempt !
@berak, silly of me to forget to show it - i've updated it now. The secondary matrix will be made complex, then multiplied with the first matrix, and then I will use OpenCV DFT
"according to the Apache Commons complex class" -- forget that. you'll have to do with opencv means, not apache, maybe it's just CvType.CV_64FC2 ?
still it's unclear, why you're trying to multiply a 1 channel Mat with a 2 channel one. context, again ? (why are you trying to do this, and what outcome do you expect ?)
@berak, Where do you see the 2 channel matrix? I am converting code from MATLAB to Java, in MATLAB it's really easy to have a matrix for the image, a matrix of complex numbers, and then I multiply them together using dot product and then pass them to the 2D fourier transform function. My outcome is to take a hologram and apply fresnel reconstruction to it, I could add more context if you like but I think it's not necessary to understand the issue
" Where do you see the 2 channel matrix?" -- in your question. a complex Mat should be mapped to either CV_64FC2 or CV_32FC2 in opencv
"I think it's not necessary to understand the issue" -- it is ! basically, you're asking - how would anyone else write this code, given X constraints, so you need to be explicit about the latter.
@berak, So I need 2 channels for the imaginary and real part each. How would I do something like this then:
How would I multiply those two matrices then? Should the first one have its single channel duplicated so it has two channels similar to the second matrix?
you cannot multiply a 1 channel Mat with a 2 channel one. again, context ?
(also, careful, matlab has seperate planes for channels, while opencv has it interleaved. porting code from there is always daunting)
it probably helps, if you can give a more mathsy (mathjax?) receipe, and abstract away from terrible matlab code
also, is this a per-element or a matrix multiplication ? (be concise, and don't expect any silly matlab knowledge here, it's far beyond the scope of this)
@berak, I don't know what more context you want... I am trying to perform a hologram reconstruction using fresnel algorithms: Image (1) ->Real Matrix, predicted diffraction of light (2) ->Real Matrix, Fresnel algorithm (3): e^(-ipi/λz) * (2) -> complex(2) ... reconstruct -> fourier transform of (1)(3)
per-element, MATLAB might be silly but at least it works ;)
(2) is the interesting part. show, how you do this
after that, it'll be probably like
dft(src,dst,COMPLEX_INPUT);
. we'll also have to discuss, where the "center" of it is.(again, matlab != opencv)
(and again, the code you show , does not reflect any of it)