Hi John,
No, i use directly Shapefile as initial_vgp_feature, the easiest way as we use mainly ESRI solution for our GIS data.
INPUT PARAMETERS
shp_with_plate_id = r"C:\Local\00_InputShp\Input_test.shp" # Input Shapefile with columns [PLATEID1, FROMAGE, TOAGE]
shp_rotated = r"C:\Local\01_OutputShp\Output_test.shp"
rot_file = r"C:\Local\T_ROT_GEM_jan21.rot"
age_r = 100
1 RECONSTRUCT THE SHAPEFILE
Load FC to rotate (if shapefile format, the pygplate feature collection can contain only 1 file)
vgp_features = pygplates.FeatureCollection(shp_with_plate_id)
Load rotation file as pygplate feature collection
rotation_features = pygplates.FeatureCollection(rot_file)
pygplates.reconstruct(vgp_features, rotation_features, shp_rotated, age_r)
del vgp_features, rotation_features2 REVERSE RECONSTRUCT THE RECONSTRUCTED SHAPEFILE
shp_to_derotate = shp_rotated # The shapefile to reverse reconstruct, is the shapefile reconstructed in first step
Load rotation file as pygplate feature collection
rotation_features = pygplates.FeatureCollection(rot_file)
Load FC to rotate (if shapefile format, the pygplate feature collection can contain only 1 file)
vgp_features = pygplates.FeatureCollection(shp_to_derotate)
pygplates.reverse_reconstruct(vgp_features, rotation_features, age_r)
del vgp_features, rotation_features