# pyGPlates: restore points to paleomap

**URL:** <https://discourse.gplates.org/t/pygplates-restore-points-to-paleomap/339>\
**Category:** Q&A\
**Created:** [9 July 2021 22:55 UTC](https://discourse.gplates.org/t/pygplates-restore-points-to-paleomap/339 "2021-07-09T22:55:44Z")\
**Posts on this page:** 1\
**Showing post:** 5

<div class="post-metadata">

**Author:** ![john.cannon](https://discourse.gplates.org/letter_avatar_proxy/v4/letter/j/3e96dc/32.png) [@john.cannon](https://discourse.gplates.org/u/john.cannon)\
**Post date:** [12 July 2021 05:19 UTC](https://discourse.gplates.org/t/pygplates-restore-points-to-paleomap/339/5 "2021-07-12T05:19:06Z")

</div>

> [@caowenrong](#):
>
> Following the example [import geometries](https://www.gplates.org/docs/pygplates/sample-code/pygplates_import_geometries_and_assign_plate_ids.html), I am able to restore the present-day points to their paleo-positions and generate associated gpmlz and shape files.

Great! That was quick.

> [@caowenrong](#):
>
> Another native question is how to append properties associated with points and restore them back together with the lat and long to a geological time? For example, I have SiO2 wt% for some plutons having lat and long and I want to see how the SiO2 was spatially distributed at 100 Ma?

You can add SiO2 as an extra column in your text file (after lon and lat) and then set that as a property on each point feature that you’re creating.

Have a look at [this thread](https://discourse.gplates.org/t/issue-with-add-a-new-property-to-a-feature-in-pygplates/156/12) for a discussion on how to do this. Since it’s not a [standard GPlates property](https://www.gplates.org/docs/gpgim/) there are some issues that that thread talks about (like being able to load these imported features into GPlates). But essentially you could do something like:

```auto
feature.add(
    pygplates.PropertyName.create_gpml('SiO2'),
    pygplates.XsDouble(SiO2_value),
    pygplates.VerifyInformationModel.no)

```

…and that would be fine if you stay in pyGPlates land where you could later get back the value from the reconstructed feature (after having reconstructed with [reconstruct](https://www.gplates.org/docs/pygplates/generated/pygplates.reconstruct.html) using the `group_with_feature` option) with something like:

```auto
SiO2_value = reconstructed_feature.get_value(
    pygplates.PropertyName.create_gpml('SiO2')).get_double()

```

However, as mentioned in the thread, loading into GPlates would be a problem, in which case the thread shows how to use a Shapefile attribute (eg, `feature.set_shapefile_attribute('SiO2', SiO2_value)`).

---

_[View the full topic](https://discourse.gplates.org/t/pygplates-restore-points-to-paleomap/339)._