Exporting shapes in a different coordinate system

[Alexey Pasumansky]

Hello Thibaud,

You need to loop through all the vertices of every shapes and convert the coordinates from EPSG::4326 to EPSG::3948. Then change chunk.shapes.crs.

[Alexey Pasumansky]

I think it should be something like the following:

Code: [Select]

chunk = PhotoScan.app.document.chunk
shapes = chunk.shapes
source = shapes.crs
out = PhotoScan.CoordinateSystem("EPSG::3948")
for shape in chunk.shapes:
      shape.vertices = [PhotoScan.CoordinateSystem.transform(v, source, out) for v in shape.vertices]
shapes.crs = out


[Thibaud Capra]

Hello Alexey,

Thanks for your precious help, as usual!

I'm also trying to export my shapes, using chunk.exportShapes as you recommended. My question is, if I have multiple shape layers, will it create one shape file per layer or one merged shape file?

I'd like to keep my shapes separated according to the shape layers, using maybe the shape layer's name to identify them. Is it possible?

[Alexey Pasumansky]

Hello Thibaud,

In chunk.exportShapes() function there's optional "groups" argument that can be used, if you wish to export only certain layers. Otherwise, all shape layers will be saved to the same file.

[Thibaud Capra]

Hello Alexey,

I don't get the shape group creation part.
Are groups the corresponding layers by default? Or do I have to assign them first to each different ShapeGroup key? If need be, how do I do that?
For example, I made a small test project with two layers I want to export (see below) as two different shape files.

For now, I have the following code:

Code: [Select]

import os, PhotoScan

doc = PhotoScan.app.document
chunk = doc.chunk
shapes = chunk.shapes
source = shapes.crs
out = PhotoScan.CoordinateSystem("EPSG::3948")
type_chunk = ['Green']

for chunk in doc.chunks:
    if not any(type_chunk[0] in chunk.label for type_chunk in type_chunk):
        for shape in chunk.shapes:
            shape.vertices = [PhotoScan.CoordinateSystem.transform(v, source, out) for v in shape.vertices]
        shapes.crs = out
    print("""Les SHP de """ + chunk.label + """ ont été transformées avec succès !""")

for chunk in doc.chunks:   
    for shape in chunk.shapes:
        if not any(type_chunk[0] in chunk.label for type_chunk in type_chunk):
            chunk.exportShapes(path = "C:\\PFE_CAPRA\\Essais_Shapes\\Green.shp", groups = ShapeGroup(key = 1))
            chunk.exportShapes(path = "C:\\PFE_CAPRA\\Essais_Shapes\\Approche.shp", groups = ShapeGroup(key = 2))

Thanks

[Alexey Pasumansky]

Hello Thibaud,

I would rather use something like the following:

Code: [Select]

for layer in chunk.shapes.groups:
      chunk.exportShapes(path, groups = layer)


[Thibaud Capra]

Hello Alexey,

Code: [Select]

for layer in chunk.shapes.layers:
      chunk.exportShapes(path, groups = layer)


I think this layer export method works in 1.3.x, I am currently using PhotoScan 1.2.6 Build 2834... I have the following error message:

Traceback (most recent call last):
File "C:\PFE_CAPRA\Essais_Shapes\Script_Shapes.py", line 18, in <module>
for layer in chunk.shapes.layers:
AttributeError: 'PhotoScan.Shapes' object has no attribute 'layers'

[Alexey Pasumansky]

Hello Thibaud,

It was my mistake, should be chunk.shapes.groups.

[Thibaud Capra]

Hello Alexey,

I think this layer grouping is too recent, as I have the following error now:

     chunk.exportShapes(path, groups = layer)
 TypeError: 'groups' is an invalid keyword argument for this function

Or is it some error in my code?

Code: [Select]

import os, PhotoScan

doc = PhotoScan.app.document
chunk = doc.chunk
shapes = chunk.shapes

path = "C:\\PFE_CAPRA\\Essais_Shapes\\Shapes"
for chunk in doc.chunks:   
    for layer in chunk.shapes.groups:
        chunk.exportShapes(path, groups = layer)

[Alexey Pasumansky]

Hello Thibaud,

In the version 1.2.6 exportShapes() has only path and items arguments, so it seems that the code can be applied to version 1.3 only.

[Thibaud Capra]

Hello Alexey,

Too bad, I'll see if I change the rest of the scripts I wrote to adapt it to 1.3.x to allow me to export automatically all my shapes, or stick to what I have for now and export my shapes the old-fashioned way, aka manually!

As usual, thanks for the great support!

[Alexey Pasumansky]

Hello Thibaud,

Alternative way in 1.2.6 could be removing all shape layers, except the required, export, then re-open of the project. Or duplicating the original chunk for each layer export and then removing the duplicate.

[Thibaud Capra]

Hello Alexey,

The problem is, I have shapes in multiple chunks, would this work in a for loop doing roughly this:

Code: [Select]

for chunk in doc.chunks:
    Delete all layers but one based on layer labels
    Export layer
    Reopen project
    Export next layer the same way
    Loop through chunks
Additional difficulty (but I found a way around for this one), it musn't loop through all the chunks but only those with some specific keyword in the chunk label.
Wouldn't closing the project reset the progress in the various chunks?

I'm aware of the fact that this is a very specific question, but my aim is to automatize as much of the processing as possible!

[Alexey Pasumansky]

Hello Thibaud,

Then you can consider using chunk duplicates:

for chunk in list(doc.chunks):
     loop though the index of the shape layers
           temp_chunk = chunk.copy()
           remove all layers but one
           export layer
           remove temp_chunk

[Thibaud Capra]

Hello Alexey,

Can I remove layers based on their labels? If so, how?
I'm using keywords in the labels to determine where to export the file. (The same way I use keywords in chunk labels to process different LODs in one batch)