http://www.k-wave.org/forum/topic/with-respect-to-forward-and-inverse-model Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:22:13 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/with-respect-to-forward-and-inverse-model#post-7869 Mon, 05 Oct 2020 16:08:50 +0000 spg2716 7869@http://www.k-wave.org/forum/ <p>Thank you very much Sir. </p> http://www.k-wave.org/forum/topic/with-respect-to-forward-and-inverse-model#post-7834 Sat, 19 Sep 2020 10:00:17 +0000 bencox 7834@http://www.k-wave.org/forum/ <p>Hi Soma,</p> <p>The 'inverse crime' refers to the situation when exactly the same set-up and model are used for both generating the data and performing the inversion. (See<br /> <a href="https://arxiv.org/abs/math-ph/0401050">https://arxiv.org/abs/math-ph/0401050</a>.) In any practical/experimental situation there will be noise in the data and uncertainty in other parameters such as the detector positions, and how sensitive the inversion is to these is crucial to how well is will perform in an application. The worry is that these aspects are not captured by a solution that commits the 'inverse crime'. On the other hand, performing an inversion that commits the inverse crime can be a useful check: if the inversion doesn't even work in this case there may be something wrong with the code, or perhaps the problem under investigation is very ill-posed. Common ways to ameliorate the inverse crime are to add noise to the simulated data, and use a different model for simulating the data and in the inversion (or at least using a different computational grid for the inversion).</p> <p>On your other question regarding the system matrix. You could construct it using k-Wave or some other model, or an analytical formula (Green's function), but it can get very large very easily. Each column will consist of a stack of time series at all the detector points generated by a point source at one grid point. So you could build it column-by-column by running lots of simulations, each with one point source at one grid point. (This will result in a full matrix, and it may be better to find a way to exploit some sparsity to reduce the memory required.) We don't usually construct the matrix explicitly, but if you do then it allows you to use all the usual tools of linear algebra to solve it.</p> <p>Hope that helps.</p> <p>Best wishes,<br /> Ben </p> http://www.k-wave.org/forum/topic/with-respect-to-forward-and-inverse-model#post-7820 Sat, 12 Sep 2020 19:29:21 +0000 spg2716 7820@http://www.k-wave.org/forum/ <p>Dear Bradley Sir,</p> <p>I have some conceptual doubt with respect to,<br /> 1. What is meant by inverse crime?<br /> 2. Just by changing the grid size how can we avoid the inverse crime ?<br /> 3. Please suggest reference or textbook to help me understand the matrix A. </p> <p> The forward model of PAT is represented as<br /> Ax = b, A = system matrix, x = initial pressure and b = measurement.</p> <p>Based on understanding gained - The system matrix is obtained by solving the governing photoacoustic wave propagation equation (either using Green's function approach or pseudospectral methods or finite difference). In kWave, the command- "kspaceFirstOrder2D" helps us create the A matrix and x to obtain b.<br /> (This becomes our Forward model)</p> <p>My intention is compute the system matrix and apply the SVE/SVD method to understand the behavior of it.</p> <p>In the following code , while doing reconstruction (Inverse model/process),<br /> you change the kgrid size in order to avoid the inverse crime but for the rest for the Time reversal example the same kgrid is used (does that mean inverse crime is committed and ignored). </p> <p>Reference Code:- (i) Attenuation Compensation Using Time Reversal Example<br /> (ii) 2D Time Reversal Reconstruction For A Line Sensor Example</p> <p>Please do correct me in case I am wrong </p> <p>Thank you for your time and consideration.</p> <p>Regards,<br /> Soma </p>