Importance sampling

From GeopsyWiki
Jump to navigation Jump to search

This is an example based on the surface wave plugin. It can be run with any other inversion plugin.


  • Generate .report files as usual (command line or graphical user interface)
  • Run dinver with importance sampling mode (only command line), list several reports if needed.
 dinver -i DispersionCurve -importance-sampling inversion.report -ns 10000 -o 3_layers.report  -param 3_layers.param

You can use -resume to start this command several times appendingmodels to the same .report.


  • Plot the profile densities:
 gpdcreport 3_layers.report | gpprofile -vs -resample -d 4000 -n 40 | grep -v "#"  | awk '{print $2,$1}' |
 gphistogram -y-min 150 -y-max 3500 -y-sampling log -x-min 0 -x-max 4000 -x-count 41

If you have some artifacts on the plot (some depths appear with a much stronger density), adjust the -x-count to one more unit compared to the number of sub-depths (40 and 41 in this case). Vs is plotted down to 4 km with a step of 100 m in this case.