# Plot pseudosection

**URL:** <https://simpeg.discourse.group/t/plot-pseudosection/114>\
**Category:** Electromagnetics\
**Created:** [June 5, 2020, 6:14am UTC](https://simpeg.discourse.group/t/plot-pseudosection/114 "2020-06-05T06:14:12Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![TravisJ](https://avatars.discourse-cdn.com/v4/letter/t/779978/32.png) [@TravisJ](https://simpeg.discourse.group/u/TravisJ)\
**Post date:** [June 5, 2020, 6:14am UTC](https://simpeg.discourse.group/t/plot-pseudosection/114/1 "2020-06-05T06:14:12Z")

</div>

Hi everyone,

I was wondering if it is possible to plot pseudosection with Simpeg? I have the apparent resistivity data, pseudo-depth and distance. I don’t need to calculate them again, I just basically need to plot pseudosection. However, I couldn’t do it. If it is possible or if you have any example, can you please help me? Thank you for your help.

---

<div class="post-metadata">

**Author:** ![lheagy](https://yyz2.discourse-cdn.com/free1/user_avatar/simpeg.discourse.group/lheagy/32/6_2.png) [@lheagy](https://simpeg.discourse.group/u/lheagy)\
**Post date:** [June 30, 2020, 3:07pm UTC](https://simpeg.discourse.group/t/plot-pseudosection/114/2 "2020-06-30T15:07:11Z")

</div>

Hi @TravisJ, Welcome! Apologies for our late response.

There are some utilities for plotting pseudosections in the resistivity module. We recently ran a tutorial at the Transform 2020 event and assembled some notebooks here: [https://github.com/simpeg/transform-2020-tutorial/](https://github.com/simpeg/transform-2020-tutorial/). In the first notebook, you will find an example that loads up some field data and plots pseudosections (both apparent resistivity and potentials).

If you have created a SimPEG survey and data object, the basic idea is

```python
from SimPEG.electromagnetics import resistivity as dc

# create a survey
survey = dc.Survey(source_list=source_list)

# create a data object
dc_data = Data(
    survey=survey,   
    dobs=potentials

# plot a pseudosection 
dc.utils.plot_pseudosection(
    dc_data, data_type="apparent resistivity"
)

```

Hope this helps!
