projector/poly_projector/poly_projector.py

Refer to the Projector Example

# $Id$
 
# Copyright (C) 2008-2014, Roman Lygin. All rights reserved.
# Copyright (C) 2014-2026, CADEX. All rights reserved.
 
# This file is part of the Manufacturing Toolkit software.
 
# You may use this file under the terms of the BSD license as follows:
 
# Redistribution and use in source and binary forms, with or without
# modification, are permitted provided that the following conditions are met:
# * Redistributions of source code must retain the above copyright notice,
#   this list of conditions and the following disclaimer.
# * Redistributions in binary form must reproduce the above copyright notice,
#   this list of conditions and the following disclaimer in the documentation
#   and/or other materials provided with the distribution.
 
# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
# AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
# ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
# LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
# CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
# SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
# CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
# ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
# POSSIBILITY OF SUCH DAMAGE.
 
import os
import sys
from pathlib import Path
 
import mtk.MTKCore as mtk
 
sys.path.append(os.path.abspath(os.path.dirname(Path(__file__).resolve()) + "/../../"))
sys.path.append(os.path.abspath(os.path.dirname(Path(__file__).resolve()) + "/../../helpers/"))
 
import LicenseHelper
 
import mtk_license as license
 
def main(theSource: str, theOutFolder: str) -> int:
    aKey = license.Value()
 
    if not LicenseHelper.SetupRuntimeKey() :
        return 1
        
    try:
        mtk.LicenseManager.Activate(aKey)
    except mtk.LicenseError as anException:
        mtk.LicenseManager.Deactivate()
        print("Failed to activate Manufacturing Toolkit license: " + anException.what())
        return 1
 
    aModel = mtk.ModelData_Model()
    aReader = mtk.ModelData_ModelReader()
 
    # Reading the file
    if not aReader.Read(mtk.UTF16String(theSource), aModel):
        mtk.LicenseManager.Deactivate()
        print("Failed to read the file " + theSource)
        return 1
 
    print("Model: " + str(aModel.Name()) + "\n")
 
    aProjectionDataCollection = [
        (mtk.UTF16String("X"), mtk.Geom_Direction.XDir()),
        (mtk.UTF16String("Y"), mtk.Geom_Direction.YDir()),
        (mtk.UTF16String("Z"), mtk.Geom_Direction.ZDir()),
    ]
 
    for aProjectionName, aProjectionDirection in aProjectionDataCollection:
        anOutFilePath = mtk.UTF16String(theOutFolder + "projection_" + str(aProjectionName) + ".stl")
 
        if ComputeProjection(aModel, aProjectionDirection, aProjectionName, anOutFilePath):
            print("The output projection is written: " + str(anOutFilePath) + "\n")
        else:
            print("Failed to save the result: " + str(anOutFilePath) + "\n")
 
    mtk.LicenseManager.Deactivate()
 
    return 0
 
if __name__ == "__main__":
    if len(sys.argv) != 3:
        print("Usage: " + sys.argv[0] + " <input_file> <output_folder>, where:")
        print("    <input_file> is a name of the file to be read")
        print("    <output_folder> is a folder to save projections to (must end with '/' or '\\')")
        sys.exit(1)
 
    aSource = os.path.abspath(sys.argv[1])
    aOutFolder = os.path.abspath(sys.argv[2])
 
    sys.exit(main(aSource, aOutFolder))
 
def PrintProjectionInfo(thePart: mtk.ModelData_Part,
                        theDirectionName: mtk.UTF16String,
                        theProjection: mtk.Projector_Projection) -> None:
    print(f"Part [{thePart.Name()}], projection {theDirectionName}:")
    if theProjection:
        print(f"    area            = {theProjection.Area()} mm")
        print(f"    outer perimeter = {theProjection.OuterPerimeter()} mm")
    else:
        print("    undefined")
    print("")
 
class ProjectionComputer(mtk.ModelData_ModelElementVoidVisitor):
    def __init__(self, theDirection: mtk.Geom_Direction, theDirectionName: mtk.UTF16String):
        super().__init__()
        self.myDirection = theDirection
        self.myDirectionName = theDirectionName
        self.myProjector = mtk.Projector_PolyProjector()
        self.myPartProjections = []
 
    def SaveProjection(self, theFileName: mtk.UTF16String) -> bool:
        if len(self.myPartProjections) == 0:
            return False
 
        aPart = mtk.ModelData_Part(mtk.UTF16String("Projections"))
 
        for aProjection in self.myPartProjections:
            aMeshBody = mtk.ModelData_MeshBody(aProjection.Mesh())
            aPart.AddBody(aMeshBody)
 
        anOutModel = mtk.ModelData_Model(mtk.UTF16String("Projector"))
        anOutModel.AddRoot(aPart)
 
        aWriter = mtk.ModelData_ModelWriter()
        return aWriter.Write(anOutModel, theFileName)
 
    def VisitPart(self, thePart: mtk.ModelData_Part):
        aProjection = self.myProjector.Perform(thePart, self.myDirection)
        PrintProjectionInfo(thePart, self.myDirectionName, aProjection)
        if aProjection:
            self.myPartProjections.append(aProjection)
 
def ComputeProjection(theModel: mtk.ModelData_Model,
                      theDirection: mtk.Geom_Direction,
                      theDirectionName: mtk.UTF16String,
                      theOutFileName: mtk.UTF16String) -> bool:
    aComputer = ProjectionComputer(theDirection, theDirectionName)
    theModel.Accept(aComputer)
    return aComputer.SaveProjection(theOutFileName)