1

I am very new to QGIS development, but not GIS. My goal is to create a copy of an existing layer (created within the QGIS GUI) from one SpatiaLite file to another, then add new attributes and populate records to the copy.

I have successfully used the QgsVectorFileWriter.writeAsVectorFormatV2() to create the copy, but checking the validity of the copy returns false and when adding to the map has a broken link in the TOC. The original layer is empty and in 4326 crs. I have done quite a bit of searching but can't track down what the issue is.

I have tried adding a transform in the map properties, I have digitized to a polygon in the copied layer and to see if creating an extent might help. I have checked the encoding. Just not sure where to go from here.

Here is the code:


uri = QgsDataSourceUri()
db_path = os.path.join(
 "C:\\",
 "Users",
 getpass.getuser(),
 "Documents",
 "temp_gis_data",
 "temp_gis_data.sqlite",
)
db_path2 = os.path.join(
 "C:\\",
 "Users",
 getpass.getuser(),
 "Documents",
 "temp_gis_data",
 "temp_gis_data2.sqlite",
)
uri.setDatabase(db_path)
schema = ""
table = "test_poly"
geom_column = "geom"
uri.setDataSource(schema, table, geom_column)
driverName = "SpatiaLite"
vlayer = QgsVectorLayer(uri.uri(), table, driverName)
print(vlayer.isValid()) --> True
QgsProject.instance().addMapLayer(vlayer)
options = QgsVectorFileWriter.SaveVectorOptions()
options.driverName = driverName 
options.layerName = "test_poly"
if not os.path.exists(db_path2):
 print("exists false")
 options.actionOnExistingFile = QgsVectorFileWriter.CreateOrOverwriteFile
 options.EditionCapability = QgsVectorFileWriter.CanAddNewLayer
else:
 options.actionOnExistingFile = QgsVectorFileWriter.CreateOrOverwriteLayer
transform_context = QgsProject.instance().transformContext()
err = QgsVectorFileWriter.writeAsVectorFormatV2(
 vlayer, db_path2,transform_context, options
)
print(err)
#update layer with data from db
uri = QgsDataSourceUri()
uri.setDatabase(db_path2)
print(uri.uri())
schema = ""
table = "test_poly"
geom_column = "geom"
uri.setDataSource(schema, table, geom_column)
vlayer2 = QgsVectorLayer(uri.uri(), table, driverName)
QgsProject.instance().addMapLayer(vlayer2) --> broken link to source, but can confirm a layer was created
print(vlayer2.isValid()) -->False
PolyGeo
65.5k29 gold badges115 silver badges350 bronze badges
asked Dec 24, 2020 at 18:23

1 Answer 1

0

The geometry column name was changed to "geometry" in the copy from "geom" which I set in the original layer. I updated the geometry column name parameter and it fixed my problem.

PolyGeo
65.5k29 gold badges115 silver badges350 bronze badges
answered Dec 24, 2020 at 18:55

Your Answer

Draft saved
Draft discarded

Sign up or log in

Sign up using Google
Sign up using Email and Password

Post as a guest

Required, but never shown

Post as a guest

Required, but never shown

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

Start asking to get answers

Find the answer to your question by asking.

Ask question

Explore related questions

See similar questions with these tags.