1
$\begingroup$

I recently stumbled into this question: Say, for example, we are doing mean-free-path sampling (distance sampling) in a scattering medium and to keep it simple, let's only consider the single scattering case, but I'm only interested in part of the region, as shown in the following figure. It won't feel like it is a problem if I just discard all the samples outside of region. enter image description here However, when I revisit the mathematical expression of the expectation, I found that: $$ \mathbb{E}(\hat{I}) = \int_0^{+\inf} \underbrace{\frac{f(x)}{p(x)}}_{\text{estimator}} \underbrace{p(x)}_{\text{PDF}} \text{d}x $$ Yeah, if the PDF term in the denominator of the estimator and the PDF of the sample $x$ are the same, these two terms will cancel each other out and we are left with the integral we want. Yet, discarding all the samples outside of the interested region is seemed to alter the PDF of the samples. I think that, for now, all the recorded samples have a conditional PDF: $$ p(x|x\text{ in region}) $$ If we are using the original estimator, for example, if we sample $x$ from an exponential distribution, we will normally have this estimator: $$ \frac{f(x)}{\sigma_t \exp(-\sigma_t x)} $$ Since we are expecting that when we sample $x$ according to $-\log(1 - x)/\sigma_t$,the sample $x$ actually follows the same distribution, but it seems that when we discard samples, these two PDF terms can not be cancelled out.

My question is:

  • Does this mean that we can should modify the PDF in the estimator according to the region we are interested? Am I right about this?
  • What should be the correct distribution of the samples? A exponential distribution "confined" in a small region?
  • What is the case for more bounces? If the one-bounce case is biased, does this mean that the multi-bounce case is also biased?
$\endgroup$
1
  • $\begingroup$ We can break down the result calculated when the interested region is infinity into non-overlapping segments. So, if every segment is biased, why they can be summed to an unbiased result? $\endgroup$ Aug 11, 2023 at 11:53

1 Answer 1

2
$\begingroup$

I think so. I did some numerical experiments using Monte-Carlo intergation.

discard no-discard within-range uniform GT
$x^2,x\in[0, 2)$ 14.705318 2.66772274 2.6668968 2.666666
$\sin x, x\in[0, \pi)$ 7.4154146 1.99990323 2.0004143 2
$\text{sigmoid},x\in[0, 1)$ 6.5161143 0.62023634 0.6201093 0.620115
$1 / (x+1)^2,x\in[0, 1)$ 5.2540569 0.49955521 0.4999410 0.5

discard means that we discard the samples that are outside of the range, no-discard means that those out-of-range samples will be kept track of. within-range uniform means that we sample in that range with uniform distribution so it does not have out-of-range samples in the first place. The original sampling approach is exponential which means that it samples in $[0, \infty)$ therefore will have out-of-range samples.

$\endgroup$
1
  • $\begingroup$ Thanks for these experiments! I think they are able to verify my thoughts. $\endgroup$ Aug 22, 2023 at 9:15

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.