Deep burn wound flap transplantation fitting method and system

Abstract

This invention relates to the field of burn surgery technology, specifically to a method and system for adapting skin flaps for deep burn wounds. The method includes the following steps: acquiring shearing information, extracting the rotation trajectory and contact area, determining trajectory continuity and identifying interference trends, analyzing spatial relationships to locate discontinuous segments and constructing mapping boundaries, dividing segments, adjusting the order, extracting deviations, aligning the final segment morphology, and generating burn wound skin flap adaptation results. In this invention, through shearing trajectory direction analysis and contact interference area identification, a tissue stress distribution map is generated and discontinuous path segments are extracted. Boundary mapping is established to achieve rotation path segmentation, forming a control sequence with directional consistency to enhance path stability. Combined with error offset information, trajectory execution deviations are dynamically corrected, achieving precise alignment of the final segment rotation direction with the fitting boundary, thus improving rotation control accuracy, path response capability, and tissue adaptation effect.

Description

Technical Field

[0001] This invention relates to the field of burn surgery technology, and in particular to a method and system for adapting skin flaps for deep burn wounds. Background Technology

[0002] Burn surgery encompasses the clinical management and surgical reconstruction techniques for burns to the skin and deep tissues caused by heat, electric current, chemicals, radiation, and other factors. Its core aspects include wound debridement, infection control, tissue protection, wound coverage, and reconstructive repair. This field typically integrates debridement, skin grafting, flap transplantation, tissue expansion, and skin substitutes to repair damaged tissues and restore function. Especially in cases of large-area, deep burns, blood-supply flap transplantation is necessary to repair refractory wounds and reduce scarring and functional impairment. The clinical technical requirements include precise matching of the transplant site and donor site, ensuring blood supply, and complex postoperative adaptation procedures. Traditional methods for adapting skin flaps for deep burns involve transplanting blood-supply skin flaps to effectively cover and repair the wound during deep burn treatment. This is achieved through preoperative empirical measurement of the wound area and manual trimming of the flap's shape and size to match the recipient wound in area, contour, and thickness. It often relies on the surgeon's clinical experience to determine the donor site selection and design path, and then gradually adjusts the flap's direction, rotation angle, and suture edges during surgery to complete the adaptation process. This method primarily achieves flap-wound alignment through physical measurements of wound length and width, estimation of the flap's physiological curvature and rotation angle, and layer-by-layer suturing during surgery.

[0003] Current technologies rely on empirical measurements and manual cutting for flap transplantation and fitting. They lack effective monitoring of tissue shearing direction and contact dynamics during rotation, making it difficult to identify directional abrupt changes and tissue interference areas during path execution. They also lack a precise trajectory segmentation mechanism to address discontinuities in the rotation path and a real-time quantification method for rotation errors, which limits the accuracy and efficiency of path adjustment. Furthermore, the lack of a clear boundary alignment strategy at the path termination stage leads to unstable fitting, accumulation of path matching errors, and decreased wound fitting accuracy, affecting flap tissue function recovery and postoperative fitting results. Summary of the Invention

[0004] To achieve the above objectives, the present invention employs the following technical solution: a method for adapting skin flaps for deep burn wounds, comprising the following steps: S1: Collect shear sensing information of the flap donor area edge and flap base area, extract the shear direction trajectory and shear contact distribution area during flap rotation, determine the direction continuity of the shear direction trajectory, identify the interference trend of the shear contact distribution area, and generate a stress interference tissue distribution map. S2: Analyze the spatial overlap between the force-affected tissue distribution map and the original rotation trajectory of the flap, judge whether the change in shear direction angle at adjacent times is continuous, locate the discontinuous segments in the rotation trajectory, construct the intervention area mapping boundary, and generate a path segmentation and positioning map; S3: Construct a rotation trajectory structure based on the path segmentation and positioning map, divide the rotation control segments with related directions, adjust the connection order and direction consistency, and generate a segmented rotation trajectory sequence; S4: Call the segmented rotation trajectory sequence, determine its consistency with the flap stress state, re-extract the shearing direction and shearing contact distribution area after simulation execution, compare the error offset, and generate a segment execution offset map; S5: Based on the fragment, perform an offset map, adjust the rotation direction of the flap's final segment and the termination boundary of the fitting path, align the tissue contact morphology of the final segment of the rotation trajectory with the edge of the burn wound, and generate the burn wound flap transplantation adaptation result.

[0005] As a further aspect of the present invention, the stress interference tissue distribution map includes shear direction distribution characteristics, shear contact area boundary, and shear stress concentration area; the path segmentation and positioning map includes trajectory segmentation nodes, interference area contour boundary, and spatial overlap determination area; the segmented rotation trajectory sequence includes rotation control sub-segments, direction consistency identifiers, and segment connection sequence information; the segment execution offset map includes shear direction offset value, contact area offset range, and rotation angle error distribution; and the burn wound flap transplantation adaptation result includes termination adhesion boundary, tissue contact adaptation morphology, and rotation path end segment docking status.

[0006] As a further aspect of the present invention, the interference trend of the shear contact distribution area refers to the change pattern of the direction and intensity of shear stress in the shear contact area during rotation, which involves spatial displacement or fracture.

[0007] As a further aspect of the present invention, the direction-associated rotation control segment refers to a flap rotation path segment that maintains the same shearing direction and is continuous with the direction of the preceding and following trajectory segments.

[0008] As a further aspect of the present invention, the specific steps of S1 are as follows: S101: Collect shear sensing information of the flap donor area edge and the base region, extract shear displacement and time series information, perform vectorization processing, and judge whether the shear direction angle change is continuous at adjacent time points to obtain the shear direction continuity judgment result. S102: Based on the shear direction continuity judgment result, collect the contact pressure distribution and time index information, cluster the contact areas at different times, construct a two-dimensional contact matrix, calculate the pressure gradient change and evolution trend, and obtain the shear contact interference change feature group. S103: Based on the shear direction continuity judgment result and the shear contact interference change feature group, match the abnormal shear direction and the interference enhancement region, map the corresponding unit to the donor area and the substrate area contour, and generate a stress interference tissue distribution map.

[0009] As a further aspect of the present invention, the specific steps of S2 are as follows: S201: Based on the force-interference tissue distribution map and the original rotation trajectory of the flap, retrieve the spatial overlapping area, calculate the spatial intersection of the trajectory point and the interference area, and obtain the trajectory interference overlapping area index set. S202: Call the index set of the trajectory interference overlapping area, determine the position of trajectory interruption based on the angle change and displacement change of the trajectory points, and extract the corresponding boundary points to obtain the set of non-continuous trajectory boundary segments; S203: Based on the set of non-continuous trajectory boundary segments and the set of indexes of overlapping areas of trajectory interference, construct a closed boundary contour, and construct an image matrix based on the image coordinate system corresponding to the force interference tissue distribution map, perform region labeling mapping and superimpose it into the image matrix to generate a path segmentation and positioning map.

[0010] As a further aspect of the present invention, the specific steps of S3 are as follows: S301: Based on the path segmentation and positioning map, obtain the spatial coordinate sequence and sequential index information of the segmented path segments, perform distance calculation and direction vector extraction on the endpoint coordinates of adjacent path segments, and perform association judgment according to the endpoint connectivity judgment threshold to generate a multi-segment rotation trajectory structure. S302: Based on the multi-segment rotation trajectory structure, calculate the direction vector sequence inside each trajectory segment, compare the direction angles between adjacent segments, and associate and group the trajectory segments according to a preset direction consistency judgment threshold to obtain the rotation control segment sequence; S303: Based on the rotation control segment sequence, adjust the connection order index and direction vector of the segments, perform vector reversal processing on segments with opposite directions, rearrange the segment order index, and combine and output the adjusted segments in a serialization manner to generate a segmented rotation trajectory sequence.

[0011] As a further aspect of the present invention, the specific steps of S4 are as follows: S401: Call the rotation control segment in the segmented rotation trajectory sequence, detect the rotation angle parameter and displacement parameter when the corresponding segment is executed, and simultaneously collect the flap shearing direction data and shearing contact distribution area data. Make a consistency judgment on the rotation parameters and shearing direction vector, obtain the matching state record corresponding to each segment, and generate a rotation force consistency judgment sequence. S402: Based on the rotational force consistency determination sequence, after each rotation is executed, the shearing direction data and shearing contact distribution area data are re-acquired, the difference in direction vector and the offset of contact area at the same spatial index position before and after execution are compared, the offset and offset direction are calculated and quantized, and the segment offset error feature set is obtained. S403: Based on the segment offset error feature set, map the spatial index position and offset of the rotation control segment to a unified image coordinate system, perform regional marking and matrix superposition processing on the offset direction and magnitude, perform image encoding and integration according to the segment order, and generate a segment execution offset map.

[0012] As a further aspect of the present invention, the specific steps of S5 are as follows: S501: Based on the spatial position and direction parameters of the last rotating segment in the segment execution offset map, extract the end direction vector corresponding to the last rotating segment, and calculate the angle between the last segment direction vector and the normal vector of the end boundary of the fitting path. If the angle exceeds the preset direction adjustment threshold, perform direction rotation matrix adjustment to obtain the last segment direction calibration parameter group. S502: Call the terminal direction calibration parameter group and the flap tissue edge contour coordinate data, perform position offset processing on the terminal path under the calibration direction, and retrieve the overlapping area index in the burn wound edge contour. The overlap between the two is scored and calculated according to the boundary point distribution matching to obtain the path termination alignment boundary index set. S503: Based on the path termination alignment boundary index set, perform vertex sequence pairing and boundary closure processing on the tissue contact area and wound edge contour of the flap rotation end segment, and jointly encode the aligned structure with the original trajectory coordinate set to establish the burn wound flap transplantation adaptation result.

[0013] A deep burn wound flap transplantation fitting system includes: The shear induction analysis module is used to achieve S1: collecting shear induction information of the flap donor area edge and flap base area, extracting the shear direction trajectory and shear contact distribution area during flap rotation, determining the direction continuity of the shear direction trajectory, identifying the interference trend of the shear contact distribution area, and generating a stress interference tissue distribution map; The trajectory intervention and positioning module is used to achieve S2: analyze the spatial overlap between the force-interference tissue distribution map and the original rotation trajectory of the flap, judge whether the shear direction angle change is continuous at adjacent time points, locate the discontinuous segment in the rotation trajectory, construct the intervention area mapping boundary, and generate a path segmentation and positioning map; The segmented trajectory reconstruction module is used to implement S3: construct a rotation trajectory structure based on the path segmentation and positioning map, divide the rotation control segments with associated directions, adjust the connection order and direction consistency, and generate a segmented rotation trajectory sequence; The offset analysis module is used to implement S4: call the segmented rotation trajectory sequence, determine the consistency with the flap force state, re-extract the shearing direction and shearing contact distribution area after simulation execution, compare the error offset, and generate a segment execution offset map; The terminal adaptation optimization module is used to implement S5: based on the fragment execution offset map, adjust the rotation direction of the flap end segment and the termination overlap boundary of the fitting path, align the tissue contact morphology of the end segment of the rotation trajectory with the edge of the burn wound, and generate the burn wound flap transplantation adaptation result.

[0014] Compared with the prior art, the advantages and positive effects of the present invention are as follows: In this invention, by analyzing the shear trajectory direction and identifying the contact interference area, a tissue stress distribution map is generated and discontinuous segments of the path are extracted. Boundary mapping is established to achieve rotation path segmentation, forming a control sequence with directional consistency to enhance path stability. Combined with error offset information, trajectory execution deviation is dynamically corrected to achieve precise alignment between the final rotation direction and the fitting boundary, thereby improving rotation control accuracy, path response capability, and tissue adaptation effect. Attached Figure Description

[0015] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0016] Figure 1 This is a schematic diagram of the steps of the present invention; Figure 2 This is a detailed schematic diagram of S1 of the present invention; Figure 3 This is a detailed schematic diagram of S2 of the present invention; Figure 4 This is a detailed schematic diagram of S3 of the present invention; Figure 5 This is a detailed schematic diagram of S4 of the present invention; Figure 6 This is a detailed schematic diagram of S5 of the present invention; Figure 7 This is a system module diagram of the present invention. Detailed Implementation

[0017] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0018] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0019] In the embodiments of this invention, the terms "image" and "picture" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, their intended meanings are consistent. Similarly, the terms "of," "corresponding (relevant)," and "corresponding" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, their intended meanings are consistent.

[0020] In this embodiment of the invention, sometimes a subscript such as W1 may be written in a non-subscript form such as W1. When the difference is not emphasized, the meaning they express is the same.

[0021] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0022] Please see Figure 1 This invention provides a method for adapting skin flaps for deep burn wounds, comprising the following steps: S1: By deploying a shear sensing device at the edge of the flap donor area and the flap base area, the shear direction trajectory and shear contact distribution area during the flap rotation process are extracted, the shear direction trajectory is judged to determine the direction continuity, the interference trend of the shear contact distribution area is identified, and the distribution map of the stress-affected tissue is output. S2: Based on the distribution map of the tissue under force interference, analyze the spatial overlap between the original rotation trajectory of the flap and the interference area, locate the discontinuous segments in the trajectory and construct the boundary mapping of the intervention area to generate a path segmentation and localization map; S3: Based on the path segmentation and positioning map, construct an operable multi-segment rotation trajectory structure, divide the path into rotation control segments with directional association, adjust the connection order and directional consistency of the control segments, and output the segmented rotation trajectory sequence; S4: Call the rotation control segment in the segmented rotation trajectory sequence, determine the consistency between each rotation operation and the force state of the flap, re-extract the shearing direction and shearing contact distribution area after each simulation execution, compare the error offset information generated during execution, and output the segment execution offset map. S5: Based on the fragment execution offset map, adjust the rotation direction of the flap terminal segment and the termination boundary of the fitting path, and align the tissue contact morphology of the flap terminal segment with the edge of the burn wound to generate the burn wound flap transplantation adaptation result.

[0023] In this embodiment, based on wound images, 3D scan data, and CTA / MRA vascular imaging data, the donor and recipient areas of the flap are spatially reconstructed to construct a unified 3D geometric model. After completing the initial model construction, the elastic parameters, thickness distribution, and vascular course information of the flap tissue are introduced into the model as constraints. By setting the rotation center, rotation angle range, and path boundaries, the movement process of the flap transferring from the donor area to the wound is parameterized. Furthermore, based on the 3D model, combined with tissue shear force sensing information and contact pressure distribution data, the stress state of the flap at each moment during rotation and adhesion is gradually simulated to obtain the corresponding shear direction change sequence and contact area evolution characteristics.

[0024] Specifically, by continuously sampling the direction vectors, displacement changes, and contact area boundaries of each path node during the simulation, locations of abrupt directional changes and areas of force interference are identified and mapped onto the 3D model to form path intervention marker areas. Based on these marker areas, the original design path is reconstructed in segments, adjusting the connection order and directional consistency of each segment to obtain a segmented rotation trajectory sequence that satisfies organizational mechanical constraints. After path reconstruction, each rotation segment is virtually executed, and the shear direction, contact area, and positional offset before and after execution are compared and analyzed to quantify the path execution error, and the rotation parameters are iteratively corrected accordingly.

[0025] After path optimization, the final rotation path is matched with the wound edge contour. By calculating the angle between the final direction vector and the wound boundary normal vector, the final direction is calibrated and adjusted, and the fitting boundary is aligned to ensure a stable fit between the flap tip and the wound edge. Finally, the simulation-optimized path planning results are applied to the actual surgical procedure, achieving consistency between preoperative virtual design and intraoperative operation, and improving the spatial matching accuracy and tissue adaptation effect of flap transplantation.

[0026] The force interference tissue distribution map includes shear direction distribution characteristics, shear contact area boundary, and shear stress concentration area; the path segmentation and positioning map includes trajectory segmentation nodes, interference area outline boundary, and spatial overlap determination area; the segmented rotation trajectory sequence includes rotation control sub-segments, direction consistency identifier, and segment connection sequence information; the segment execution offset map includes shear direction offset value, contact area offset range, and rotation angle error distribution; and the burn wound flap transplantation adaptation results include termination adhesion boundary, tissue contact adaptation morphology, and rotation path end segment docking status.

[0027] Please see Figure 2 The specific steps of S1 are as follows: S101: Collect shear sensing information of the flap donor area edge and the base region, extract shear displacement and time series information, perform vectorization processing, and judge whether the shear direction angle change is continuous at adjacent time points to obtain the shear direction continuity judgment result. First, shear force sensor arrays are deployed at the edges of the supply area and the base area. Each sensor records shear force changes along two main directions (usually horizontal and vertical) at 1-second intervals. By continuously recording the shear force changes at each node and combining this with the initial position coordinates of the sensors, the shear displacement of each node per unit time is calculated, forming a shear displacement time series. Further, vectors are constructed by combining the displacement directions of adjacent time nodes, forming a shear direction vector sequence on the time axis. At each time point, the corresponding shear direction angle is extracted, and the change of this angle over time is continuously recorded to determine if there are any abrupt changes. By calculating the difference in direction angles between two adjacent time points item by item, the angle change value is compared with the predicted value. A threshold for continuous shearing direction is set for comparison. This threshold is set based on the fluctuation of the directional angle change under the actual stable state of flap shearing, and is usually set to 15 degrees. If there are more than five sudden changes in direction within a continuous time period, that is, the number of times the directional angle change exceeds the threshold is five or more, it is considered that the continuity of shearing direction has been disrupted. In actual operation, if 9 sensor nodes are set up at the edge of the donor area, and within a 10-second sampling period, the 3rd, 4th, 5th, 6th and 7th nodes record continuous directional changes of 20 degrees, 18 degrees, 22 degrees, 25 degrees and 21 degrees respectively, all exceeding the set threshold more than five times, then it is determined that there is a shearing direction continuity interruption in the area corresponding to these nodes, thus obtaining the shearing direction continuity judgment result.

[0028] S102: Based on the shear direction continuity judgment result, collect contact pressure distribution and time index information, cluster the contact area at each time, construct a two-dimensional contact matrix, calculate the pressure gradient change and evolution trend, and obtain the shear contact interference change feature group. An array of pressure sensing units is deployed in the contact area between the donor area and the basal tissue. Each unit samples at a fixed frequency, recording the contact pressure value at each moment, forming a time-based pressure data sequence. For the pressure distribution at each moment, a two-dimensional contact pressure matrix containing all sensing points is constructed. The pressure data set of the contact area at the same time point is spatially clustered. The clustering process is based on the similarity of pressure values ​​between adjacent units. If the pressure difference between two units does not exceed 3 kPa, they are classified into the same contact area. This results in multiple spatially continuous contact blocks at each time point. By tracking the average pressure change of the same contact block across multiple consecutive time points, its temporal sequence is identified. The evolution trend in the column is analyzed, and then the pressure change gradient is analyzed. By calculating the pressure difference between each unit, the pressure gradient at the current time point is obtained and compared with the pressure gradient of the corresponding unit at the previous time point. The gradient change is calculated. If the gradient change exceeds the preset disturbance identification threshold, it is determined that the contact block has interference behavior. The disturbance identification threshold is based on the normal contact pressure change rate of the tissue, for example, it is set to 10 kPa per second. If a contact block has a pressure gradient change greater than this threshold more than twice within 3 consecutive seconds, it is determined that the contact block has shear contact interference. Finally, the contact blocks that meet the disturbance judgment conditions at each time point are classified and organized to form a shear contact interference change feature group.

[0029] S103: Based on the shear direction continuity judgment result and the shear contact interference change feature group, match the abnormal shear direction and the interference enhancement area, map the corresponding unit to the donor area and the base area contour, and generate a stress interference tissue distribution map; First, the coordinates of the identified shear direction interruption points are extracted from the shear direction judgment results to form a set of suspected abnormal region nodes. Simultaneously, the coordinates of contact region elements with drastic pressure changes and enhanced interference are extracted from the shear contact interference feature set to construct a set of interference-enhanced regions. Cross-matching is performed on the two sets to determine whether the shear direction interruption point falls inside or on the edge of any interference-enhanced contact block. If there is spatial overlap, the node is considered to have both directional abrupt changes and pressure disturbance characteristics, and is classified as a stress-interference-enhanced element. Through geometric mapping relationships, this type of element is further classified... The location is mapped onto the original donor area or basal area tissue structure outline map to generate a tissue stress interference map. For each stress-enhancing unit in the mapped area, a 3-row, 3-column neighboring unit area is extended from it as the center and uniformly identified as the stress interference range area. For example, if the unit with coordinates (23, 35) satisfies the interference enhancement condition, then all units in the range from (22, 34) to (24, 36) around it are marked as interference areas and marked with specific gray levels or colors in the tissue map. The spatial locations of multiple interference points are integrated to complete the generation of the final stress interference tissue distribution map.

[0030] Please see Figure 3 The specific steps of S2 are as follows: S201: Based on the distribution map of the tissue under force interference and the original rotation trajectory of the flap, retrieve the spatial overlapping area, calculate the spatial intersection of the trajectory points and the interference area, and obtain the index set of trajectory interference overlapping area; First, the two-dimensional trajectory coordinate sequence of the flap's original rotation path during surgery is standardized and unified with the distribution of the tissue affected by force interference. Figure 1 In the image coordinate system, the original trajectory is represented as a continuous set of points, each with position coordinates (x, y) and rotation order. Then, the marked interference enhancement units in the force interference distribution map are traversed point by point, recording the spatial boundary range of all force-affected areas. A spatial overlap detection mechanism based on coordinate matching is established to determine whether each trajectory point falls within the boundary of any force interference unit. When the trajectory point coordinates are within the boundary of the interference area, spatial overlap is determined. The sequential number of the trajectory points is combined to form a trajectory point index. Finally, all points related to the interference are output as a set of index numbers. Trajectory points with spatial intersection in the region constitute the trajectory interference overlap region index set. For example, if the original trajectory contains 100 trajectory points, and 5 interference units are identified in the interference distribution map, their boundaries correspond to rectangular regions such as (20, 20)-(25, 25) and (30, 30)-(34, 34) in the image coordinates, when the trajectory point numbered 15, 16, 48, 49, 72, 73, and 74 fall within these regions, they are recorded as interference overlap points. Then the output index set is {15, 16, 48, 49, 72, 73, 74}, thus obtaining the trajectory interference overlap region index set.

[0031] S202: Call the track interference overlapping area index set, determine the location of the track interruption based on the angle change and displacement change of the track points, and extract the corresponding boundary points to obtain the set of non-continuous track boundary segments; First, based on the sequence number of the original trajectory points, calculate the angle difference between each trajectory point and its adjacent points to determine if there is an angle abrupt change. Specifically, construct a direction vector from two adjacent trajectory points and iterate through all trajectory points, comparing the angle change between the previous and current direction vectors. If the angle change exceeds a set continuity threshold (e.g., 20 degrees), record the trajectory point as an angle abrupt change point. Simultaneously, monitor the displacement changes of the trajectory points, i.e., determine if there are abnormal fluctuations in the Euclidean distance between the current and previous points. If the displacement suddenly increases beyond a set threshold (e.g., 15 pixels), it is also recorded as a displacement abrupt change point. Combining the results of angle and displacement abrupt changes, determine the location where the trajectory continuity is interrupted. If a trajectory point simultaneously satisfies the angle... If the abrupt change in displacement or the presence of the same characteristics at multiple consecutive points are identified, the segment is determined to be a trajectory interruption interval. The starting and ending points of the interruption segment are then extracted as boundary points. In the above steps, the judgment range is limited by the aforementioned trajectory interference overlap area index set. That is, the above judgment process is only performed on trajectory points located within the interference area, thereby effectively reducing the processing burden on non-critical areas. For example, if the interference index set is {15, 16, 48, 49, 72, 73, 74}, and the angle change between points 48 and 49 is found to be 26 degrees and the displacement change is 18 pixels, both exceeding the set threshold, then points 48 and 49 are extracted as a boundary segment. Finally, multiple such boundary pairs are formed and sorted for output, resulting in a set of non-continuous trajectory boundary segments.

[0032] S203: Based on the set of non-continuous trajectory boundary segments and the index set of trajectory interference overlapping areas, construct a closed boundary contour, and construct an image matrix based on the image coordinate system corresponding to the force interference tissue distribution map. Perform region labeling mapping and superimpose it into the image matrix to generate a path segmentation and positioning map. First, the coordinates of the start and end points of all discontinuous trajectory segments are extracted, and then connected according to the order relationship between adjacent trajectory segments to construct closed bounding boxes around these segments. For each pair of boundary start and end points, the boundary is extended to the edge points of the surrounding interference area according to the actual direction of the trajectory. All relevant trajectory nodes are connected and arranged sequentially to form a closed boundary contour. Then, this contour is used as a labeled region for image space mapping processing. The image pixel coordinates covered by this contour are assigned specific label values, and a bitmap-level overlay operation is performed with the original organized image matrix to generate a segmented image with path interference regions. In the matrix, the label value of each pixel position records the region affiliation in integer form. Non-path regions are assigned a value of 0, and interference path regions are assigned a value of 1 or a higher level for hierarchical identification. Finally, a path segmentation and localization map with complete annotation, continuous coverage, and clear boundaries is formed. For example, if a non-continuous trajectory segment is composed of the start and end points (48, 49), and the closed contour constructed around the interference region in the image corresponds to the coordinate range (120, 130) - (140, 150), then all pixels in this range in the image matrix are assigned a value of 1, and other regions are assigned a value of 0. The superimposed image is the final path segmentation and localization map.

[0033] Please see Figure 4 The specific steps of S3 are as follows: S301: Based on the path segmentation and positioning map, obtain the spatial coordinate sequence and sequential index information of the segmented path segments, perform distance calculation and direction vector extraction on the endpoint coordinates of adjacent path segments, and perform association judgment according to the endpoint connectivity judgment threshold to generate a multi-segment rotation trajectory structure. First, the marked independent trajectory regions in the path segmentation map are extracted into several individual path segments. Each segment is bounded by its starting and ending points, forming a coordinate sequence of continuous coordinate points. These coordinates are then labeled with their sequential numbers within the overall path structure. The coordinates of the endpoints of all path segments are extracted. The distance between the ending point of an adjacent path segment and the starting point of the next path segment is compared segment by segment. The Euclidean distance between the two points is calculated and compared to a preset endpoint connectivity threshold. This threshold is set based on the minimum permissible breakpoint distance during flap rotation, for example, 20 pixels. If the distance between the endpoints of two path segments is less than or equal to 20 pixels, they are considered connectable path segments, and a connectivity marker operation is performed. Simultaneously, the vector formed by the two endpoints is used as the path segment direction vector, and this direction information is recorded for subsequent trajectory structure determination. All path segments that satisfy connectivity are sequentially combined to form continuous rotation segments, constructing a multi-segment rotation trajectory structure. For example, if four path segments are extracted from the segmentation image, numbered T1 to T4, with endpoint coordinates of (50, 60)-(70, 80), (71, 81)-(90, 100), (91, 100)-(110, 115), and (130, 130)-(145, 145), the endpoint distances between T1 and T4 are 1.41, 1.00, and 2.24 pixels, respectively, all less than the 20-pixel threshold, and can be determined as connected segments. The distance between T4 and T5 is approximately 28.28 pixels, exceeding the threshold, and does not satisfy the connectivity condition. Finally, T1, T2, T3, and T4 are combined into a multi-segment trajectory structure, while T5 is an independent segment, generating a multi-segment rotation trajectory structure.

[0034] S302: Based on the multi-segment rotation trajectory structure, calculate the direction vector sequence inside each trajectory segment, compare the direction angles between adjacent segments, and associate and group the trajectory segments according to the preset direction consistency judgment threshold to obtain the rotation control segment sequence; For each trajectory segment, extract its coordinate sequence composed of continuous points. Calculate the connecting direction vector for any pair of adjacent coordinate points, record the direction changes along the entire path, and form a direction vector sequence within that segment. For any two adjacent trajectory segments, calculate the angle between the average direction vectors of each segment. If the angle is less than a set direction consistency judgment threshold, the two trajectory segments are considered to have the same direction. This threshold is set based on the allowable range of actual organizational rotation direction offset, for example, 30 degrees. When the direction angle between two segments is between 0 and 30 degrees, they are considered to be segments in the same direction. For trajectories that meet the conditions... The segments are associated and grouped to form a trajectory control unit with a unified direction. All segment pairs are traversed step by step and their directional relationships are determined. All trajectory segments with similar directions are merged into the same group to form one or more rotation control segment sequences. For example, if the direction of segment T1 is 45 degrees, the direction of T2 is 50 degrees, the direction of T3 is 48 degrees, and the direction of T4 is 120 degrees, then the directional differences of T1, T2, and T3 are 5 degrees and 3 degrees, respectively, which are all less than the 30-degree threshold and are assigned to the same directional control group. However, the angle between T4 and T3 is 72 degrees, which exceeds the threshold and is assigned to a new group. Finally, the rotation control segment sequence is obtained.

[0035] S303: Based on the rotation control segment sequence, adjust the connection order index and direction vector of the segments, perform vector reversal processing on segments with opposite directions, rearrange the segment order index, and combine and output the adjusted segments in a serialization manner to generate a segmented rotation trajectory sequence. First, the trajectory segments within each control group are arranged according to the original trajectory index order. Then, the direction vectors of all adjacent trajectory segments are compared one by one to determine if the direction vector of the current segment is opposite to that of the previous segment. The criterion for opposite directions is that the included angle is close to 180 degrees, i.e., the dot product of the two direction vectors is less than -0.9. If the direction is determined to be opposite, the direction vector of that trajectory segment is reversed, i.e., its start and end positions are swapped to make its direction consistent with the previous segment. After performing the direction correction, all trajectory segments are re-indexed and adjusted according to the connectivity order between them, starting from the original sequence. The segments are reordered to ensure the continuity of the path in physical space. Finally, all the trajectory segments that have been aligned and ordered are spliced ​​together according to the new index to form a complete serialized rotation trajectory. For example, if the directions of segments T1, T2, and T3 are 45 degrees, 225 degrees, and 47 degrees respectively, and the angle between the direction of T2 and T1 is 180 degrees, which meets the reverse standard, a reversal operation needs to be performed. The start and end points of T2 are swapped and the direction is adjusted to 45 degrees, which is consistent with the direction of T1 and T3. Then, they are renumbered as T1, T2, and T3 to complete the serialization and generate a segmented rotation trajectory sequence.

[0036] Please see Figure 5 The specific steps of S4 are as follows: S401: Call the rotation control segment in the segmented rotation trajectory sequence, detect the rotation angle parameter and displacement parameter when the corresponding segment is executed, and simultaneously collect the flap shearing direction data and shearing contact distribution area data. Make a consistency judgment on the rotation parameters and shearing direction vector, obtain the matching state record corresponding to each segment, and generate a rotation force consistency judgment sequence. First, the starting and ending coordinates of each segment are read sequentially according to the order of the segments in the trajectory sequence. The rotation angle is calculated based on the difference between the starting and ending positions, and is expressed as the angle between the line connecting the two points and the horizontal axis in the coordinate system. Simultaneously, the total displacement of the trajectory points in that segment is obtained; this displacement is the sum of the cumulative Euclidean distances between all points on the path of that segment, and is used as the actual displacement parameter of that segment. Then, the shear sensor array is simultaneously activated. Within the actual execution time window of each rotation segment, shear direction vector data is acquired from the donor flap tissue, and the shear contact distribution area data on the contact surface is recorded. The shear direction vector is represented by the shear force direction per unit time at each sampling point, and is a set of two-dimensional unit vectors. The shear contact distribution area data consists of the pressure distribution matrix at the same moment that is above the shear response threshold. The cell set is then used to perform a vector-by-vector consistency check on the rotation angle parameter of each rotation segment and the shear direction vector collected within the corresponding time period. By analyzing the angle between the rotation direction and the shear direction, if the absolute angle is less than the consistency threshold of 30 degrees, it is considered that the execution direction of the segment is consistent with the shear direction; otherwise, it is inconsistent. The judgment result is recorded as "1" or "0" to indicate whether it is consistent. For example, for the third rotation trajectory, its rotation angle is 45 degrees, and the average value of the collected shear direction vector is 50 degrees. The angle between the two is 5 degrees, which meets the consistency judgment standard and is recorded as 1. However, the angle of the fourth segment is 60 degrees, which exceeds the threshold and is recorded as 0. After completing the direction consistency judgment for all trajectory segments, a matching status record for each segment is formed, generating a rotation force consistency judgment sequence.

[0037] S402: Based on the rotational force consistency judgment sequence, after each rotational execution, the shearing direction data and shearing contact distribution area data are re-acquired. The difference in direction vector and the offset of the contact area at the same spatial index position before and after execution are compared. The offset and offset direction are calculated and quantized to obtain the segment offset error feature set. After each rotation, shear direction data and shear contact distribution area data are re-acquired. During execution, for each control segment, a new round of shear vector and contact area data acquisition is performed immediately after its completion. The current shear direction vector is extracted and compared one-to-one with the shear direction vector recorded before the segment execution. The difference in direction vectors at the same spatial position is calculated, specifically the absolute value of the difference in direction angle. If the difference is greater than a set threshold of 20 degrees, it is marked as a shear direction abrupt change point. Then, the cell positions of the shear contact distribution area before and after execution are compared, and the changes in cell indexes that have shifted position are recorded. The offset is calculated as the magnitude of the difference in coordinates between the old and new indices in the image coordinate system, and the offset direction is extracted as the direction angle of the offset vector. Finally, the offset and offset direction of all points are normalized. The process involves uniform processing and quantization encoding. The encoding method divides the offset into 5-pixel levels and assigns numbers to the directions in 8 quadrants (with a step size of 45 degrees). A two-dimensional encoding value is generated for each offset point. For example, if the offset level is 2 and the direction is located in the 135° quadrant, the point is encoded as (2, 3). For instance, if the shearing direction of the fifth segment of the trajectory is 60 degrees before execution and becomes 85 degrees after execution, the direction difference is 25 degrees, which exceeds the threshold. This is then counted as a sudden change point. At the same time, the shearing contact area is detected. The original location was (45, 60), and the new location is (48, 64). The offset is 5 units, the direction is approximately 53 degrees, the corresponding level is 1, the direction quadrant is 2, and the encoding is (1, 1). After completing the offset encoding point by point, the encoding information of all offset points of the segment is summarized to form the segment offset error feature set.

[0038] S403: Based on the segment offset error feature set, map the spatial index position and offset of the rotation control segment to a unified image coordinate system, perform regional marking and matrix superposition processing on the offset direction and magnitude, perform image encoding and integration according to the segment order, and generate a segment execution offset map. First, a coordinate mapping table is established, binding the spatial starting position index of each trajectory segment to its segment number. For the offset error points collected before and after the segment's execution, the offset direction and offset magnitude information in its encoding are read, and the position coordinates before the offset are found in the unified image coordinate system. Centered on this position, the offset magnitude is extended along the specified direction on the image matrix according to the offset direction, and the pixel label value of the target coordinates is modified accordingly. The offset direction is assigned as quadrant encoding, and the offset magnitude value is multiplied by 10 as the pixel grayscale value, forming a region feature jointly labeled with direction and magnitude. Next, this region is written into the corresponding segment region layer in the target image matrix. If multiple segments have overlapping areas, they are superimposed in the order of the segments. The segment executed later can overwrite the encoding value of the previous segment. A new marker layer is added to the image matrix to record the index label of each segment. For example, if an offset point in the 6th segment is encoded as (3, 2), that is, the offset level is 3, the corresponding amplitude is 15 pixels, and the direction quadrant is the 2nd quadrant (the direction angle range is 45-90 degrees), then the original position is offset by 15 units in that direction, and the pixel value 30 is written to the image matrix at the target position, and the quadrant value 2 is written to the direction layer. The writing operation of the offset area of ​​all segments is completed point by point. After the segments are superimposed layer by layer in the order of the segments, the integrated image output is completed, and the segment execution offset map is generated.

[0039] Please see Figure 6 The specific steps of S5 are as follows: S501: Based on the spatial position and orientation parameters of the last rotating segment in the fragment execution offset map, extract the end orientation vector corresponding to the last rotating segment, and calculate the angle between the last orientation vector and the normal vector of the end boundary of the fitting path. If the angle exceeds the preset orientation adjustment threshold, perform orientation rotation matrix adjustment and obtain the last orientation calibration parameter set. First, locate the last segment of the rotation trajectory sequence in the segment execution offset map, i.e., the final segment. Construct a direction vector by using the termination coordinate point of this segment and the previous trajectory point. This vector reflects the final motion direction of the rotation final segment and is recorded as the end direction vector. Simultaneously, read the boundary contour coordinate point sequence that fits the path termination boundary. Calculate the local normal vector direction of two adjacent coordinate points along the outer edge of the path at the boundary termination point. This direction is perpendicular to the boundary tangent direction. Then, calculate the angle between the final segment direction vector and this normal vector. The angle formed by the two vectors is used as the criterion. If the calculated result is greater than the set direction adjustment threshold, it is determined that there is a direction deviation. This threshold is set as the maximum allowable rotation. The rotation error angle is typically set to 20 degrees. Depending on whether this threshold is exceeded, a direction rotation matrix adjustment operation is performed. This involves correcting the direction of the final segment's direction vector through reverse rotation compensation. The correction angle is the remaining amount obtained by subtracting the direction adjustment threshold from the included angle value. The direction is based on the direction of motion of the rotating segment as the positive reference, and is corrected clockwise or counterclockwise. The corresponding direction calibration parameter set is recorded based on the correction direction and angle. For example, if the final segment's direction vector is 60 degrees and the termination boundary normal vector is 100 degrees, then the included angle is 40 degrees, exceeding the threshold of 20 degrees. A rotation adjustment is required, with a correction angle of 20 degrees and a clockwise rotation direction. The direction calibration parameter set is (clockwise, 20). The final segment direction calibration parameter set is then obtained.

[0040] S502: Call the terminal direction calibration parameter group and the flap tissue edge contour coordinate data, perform position offset processing on the terminal path under the calibration direction, and retrieve the overlapping area index in the burn wound edge contour. The overlap between the two is scored and calculated according to the boundary point distribution matching to obtain the path termination alignment boundary index set. First, using the starting point of the final segment of the flap contour as a reference, the direction vector of the final segment is rotated. Following the rotation direction and angle recorded in the direction calibration parameter set, all coordinate points in the final segment path are rotated by the same angle, generating a coordinate sequence of the final segment path after direction adjustment. Then, this corrected path is overlaid with the flap tissue edge contour data. The overlapping area in the image coordinate system is extracted, and the coordinates of all boundary points of the final segment path are scanned to find points that coincide with all coordinate points in the burn wound edge contour within a 2-pixel range. Their index information is recorded, establishing a set of overlapping area indexes. Finally, the final segment of the path and the wound edge are matched... The matching process involves calculating a score, specifically evaluating the relative positions of each pair of points in the overlapping region index, and calculating the ratio of the length of the continuous distribution segment to the total number of boundary points. This ratio reflects the degree of overlap. If overlapping points appear in a concentrated distribution area in the path boundary and the number of matching points exceeds 30%, it is judged as a valid alignment region. The boundary point numbers corresponding to this region are included in the path termination alignment boundary index set. For example, if there are 60 points in the final path boundary and 22 matching points are continuously distributed, the overlap ratio is 36.7%, which exceeds the preset 30% threshold. Then, the set of indexes of these 22 points is the path termination alignment boundary index set.

[0041] S503: Based on the path termination alignment boundary index set, perform vertex sequence pairing and boundary closure processing on the tissue contact area and wound edge contour of the flap rotation end segment, and jointly encode the aligned structure with the original trajectory coordinate set to establish the burn wound flap transplantation adaptation result. First, extract the path boundary points and wound boundary points from the alignment boundary index set, and pair them according to their spatial position. Each pair of paired points constitutes a vertex mapping unit. Construct boundary connection segments for the connection paths between each mapping unit. If the distance between any two paired points exceeds a set boundary spacing threshold (e.g., 5 pixels), supplement the intermediate connection points using linear interpolation to form a continuous boundary closed chain. This chain is used to connect the final rotation path with the boundary of the wound edge, forming a complete closed structure. Then, expand the coordinates of all points within this structure according to the point numbering format of the original trajectory coordinate set. The encoding is performed in the format of (original trajectory number, supplementary marker bit). The original trajectory number is used to track the source path information, and the supplementary marker bit is used to identify that the point belongs to a closed pair structure. For example, if the original trajectory number is T_end.45, the corresponding closed supplementary point is T_end.45-B. By combining the coordinate points of all closed structures with the original trajectory points, a complete joint coordinate set of trajectory and structure is constructed. Finally, the joint coordinate set is mapped to the image coordinate system of the burn area, and a regional matching identifier with the original surface contour is established to form a structural adaptation mapping relationship of tissue contact and adhesion area, and to establish the adaptation result of burn wound flap transplantation.

[0042] Please see Figure 7 A deep burn wound flap transplantation fitting system, including: The shear induction analysis module is used to achieve S1: collecting shear induction information of the flap donor area edge and flap base area, extracting the shear direction trajectory and shear contact distribution area during flap rotation, determining the direction continuity of the shear direction trajectory, identifying the interference trend of the shear contact distribution area, and generating a stress interference tissue distribution map; The trajectory intervention and positioning module is used to achieve S2: analyze the spatial overlap between the force-affected tissue distribution map and the original rotation trajectory of the flap, judge whether the change of shear direction angle at adjacent time points is continuous, locate the discontinuous segments in the rotation trajectory, construct the intervention area mapping boundary, and generate a path segmentation and positioning map; The segmented trajectory reconstruction module is used to implement S3: construct a rotation trajectory structure based on the path segmentation and positioning map, divide the rotation control segments with related directions, adjust the connection order and direction consistency, and generate a segmented rotation trajectory sequence; The offset analysis module is used to implement S4: call the segmented rotation trajectory sequence, judge its consistency with the flap stress state, re-extract the shearing direction and shearing contact distribution area after simulation, compare the error offset, and generate the segment execution offset map; The terminal adaptation optimization module is used to implement S5: based on the fragment execution offset map, adjust the rotation direction of the flap end segment and the termination overlap boundary of the fitting path, align the tissue contact morphology of the end segment of the rotation trajectory with the edge of the burn wound, and generate the burn wound flap transplantation adaptation result.

[0043] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for deep burn wound flap grafting adaptation, characterized in that, Includes the following steps: S1: Collect shear sensing information of the flap donor area edge and flap base area, extract the shear direction trajectory and shear contact distribution area during flap rotation, determine the direction continuity of the shear direction trajectory, identify the interference trend of the shear contact distribution area, and generate a stress interference tissue distribution map. S2: Analyze the spatial overlap between the force-affected tissue distribution map and the original rotation trajectory of the flap, judge whether the change in shear direction angle at adjacent times is continuous, locate the discontinuous segments in the rotation trajectory, construct the intervention area mapping boundary, and generate a path segmentation and positioning map; S3: Based on the path segmentation and positioning map, construct a rotation trajectory structure, divide the rotation control segments with related directions, adjust the connection order and direction consistency, and generate a segmented rotation trajectory sequence; S4: Call the segmented rotation trajectory sequence, determine its consistency with the flap stress state, re-extract the shearing direction and shearing contact distribution area after simulation execution, compare the error offset, and generate a segment execution offset map; S5: Based on the fragment, perform an offset map, adjust the rotation direction of the flap's final segment and the termination boundary of the fitting path, align the tissue contact morphology of the final segment of the rotation trajectory with the edge of the burn wound, and generate the burn wound flap transplantation adaptation result.

2. The deep burn wound flap grafting adaptation method according to claim 1, characterized in that, The stress-interference tissue distribution map includes shear direction distribution characteristics, shear contact area boundaries, and shear stress concentration areas. The path segmentation and positioning map includes trajectory segmentation nodes, interference area contour boundaries, and spatial overlap determination areas. The segmented rotation trajectory sequence includes rotation control sub-segments, direction consistency identifiers, and segment connection sequence information. The segment execution offset map includes shear direction offset values, contact area offset ranges, and rotation angle error distribution. The burn wound flap transplantation adaptation results include termination adhesion boundaries, tissue contact adaptation morphology, and rotation path end segment docking status.

3. The deep burn wound flap grafting adaptation method according to claim 1, wherein, The interference trend of the shear contact distribution area refers to the change pattern of the direction and intensity of shear stress in the shear contact area during rotation, which may result in spatial displacement or fracture.

4. The deep burn wound flap grafting adaptation method according to claim 1, wherein, The direction-associated rotation control segment refers to a flap rotation path segment that maintains the same shearing direction and is continuous with the directions of the preceding and following trajectory segments.

5. The deep burn wound flap grafting adaptation method according to claim 1, wherein, The specific steps of S1 are as follows: S101: Collect shear sensing information of the flap donor area edge and the base region, extract shear displacement and time series information, perform vectorization processing, and judge whether the shear direction angle change is continuous at adjacent time points to obtain the shear direction continuity judgment result. S102: Based on the shear direction continuity judgment result, collect the contact pressure distribution and time index information, cluster the contact areas at different times, construct a two-dimensional contact matrix, calculate the pressure gradient change and evolution trend, and obtain the shear contact interference change feature group. S103: Based on the shear direction continuity judgment result and the shear contact interference change feature group, match the abnormal shear direction and the interference enhancement region, map the corresponding unit to the donor area and the substrate area contour, and generate a stress interference tissue distribution map.

6. The deep burn wound flap grafting adaptation method according to claim 1, wherein, The specific steps of S2 are as follows: S201: Based on the force-interference tissue distribution map and the original rotation trajectory of the flap, retrieve the spatial overlapping area, calculate the spatial intersection of the trajectory point and the interference area, and obtain the trajectory interference overlapping area index set. S202: Call the index set of the trajectory interference overlapping area, determine the position of trajectory interruption based on the angle change and displacement change of the trajectory points, and extract the corresponding boundary points to obtain the set of non-continuous trajectory boundary segments; S203: Based on the set of non-continuous trajectory boundary segments and the set of indexes of overlapping areas of trajectory interference, construct a closed boundary contour, and construct an image matrix based on the image coordinate system corresponding to the force interference tissue distribution map, perform region labeling mapping and superimpose it into the image matrix to generate a path segmentation and positioning map.

7. The deep burn wound flap grafting adaptation method according to claim 1, wherein, The specific steps for S3 are as follows: S301: Based on the path segmentation and positioning map, obtain the spatial coordinate sequence and sequential index information of the segmented path segments, perform distance calculation and direction vector extraction on the endpoint coordinates of adjacent path segments, and perform association judgment according to the endpoint connectivity judgment threshold to generate a multi-segment rotation trajectory structure. S302: Based on the multi-segment rotation trajectory structure, calculate the direction vector sequence inside each trajectory segment, compare the direction angles between adjacent segments, and associate and group the trajectory segments according to a preset direction consistency judgment threshold to obtain the rotation control segment sequence; S303: Based on the rotation control segment sequence, adjust the connection order index and direction vector of the segments, perform vector reversal processing on segments with opposite directions, rearrange the segment order index, and combine and output the adjusted segments in a serialization manner to generate a segmented rotation trajectory sequence.

8. The deep burn wound flap grafting adaptation method according to claim 1, wherein, The specific steps of S4 are as follows: S401: Call the rotation control segment in the segmented rotation trajectory sequence, detect the rotation angle parameter and displacement parameter when the corresponding segment is executed, and simultaneously collect the flap shearing direction data and shearing contact distribution area data. Make a consistency judgment on the rotation parameters and shearing direction vector, obtain the matching state record corresponding to each segment, and generate a rotation force consistency judgment sequence. S402: Based on the rotational force consistency determination sequence, after each rotation is executed, the shearing direction data and shearing contact distribution area data are re-acquired, the difference in direction vector and the offset of contact area at the same spatial index position before and after execution are compared, the offset and offset direction are calculated and quantized, and the segment offset error feature set is obtained. S403: Based on the segment offset error feature set, map the spatial index position and offset of the rotation control segment to a unified image coordinate system, perform regional marking and matrix superposition processing on the offset direction and magnitude, perform image encoding and integration according to the segment order, and generate a segment execution offset map.

9. The method for adapting skin flaps for deep burn wounds according to claim 1, characterized in that, The specific steps of S5 are as follows: S501: Based on the spatial position and direction parameters of the last rotating segment in the segment execution offset map, extract the end direction vector corresponding to the last rotating segment, and calculate the angle between the last segment direction vector and the normal vector of the end boundary of the fitting path. If the angle exceeds the preset direction adjustment threshold, perform direction rotation matrix adjustment to obtain the last segment direction calibration parameter group. S502: Call the terminal direction calibration parameter group and the flap tissue edge contour coordinate data, perform position offset processing on the terminal path under the calibration direction, and retrieve the overlapping area index in the burn wound edge contour. The overlap between the two is scored and calculated according to the boundary point distribution matching to obtain the path termination alignment boundary index set. S503: Based on the path termination alignment boundary index set, perform vertex sequence pairing and boundary closure processing on the tissue contact area and wound edge contour of the flap rotation end segment, and jointly encode the aligned structure with the original trajectory coordinate set to establish the burn wound flap transplantation adaptation result.

10. A flap transplantation and fitting system for deep burn wounds, characterized in that, The system is used to implement the deep burn wound flap transplantation and adaptation method according to any one of claims 1-9, the system comprising: The shear induction analysis module is used to achieve S1: collecting shear induction information of the flap donor area edge and flap base area, extracting the shear direction trajectory and shear contact distribution area during flap rotation, determining the direction continuity of the shear direction trajectory, identifying the interference trend of the shear contact distribution area, and generating a stress interference tissue distribution map; The trajectory intervention and positioning module is used to achieve S2: analyze the spatial overlap between the force-interference tissue distribution map and the original rotation trajectory of the flap, judge whether the shear direction angle change is continuous at adjacent time points, locate the discontinuous segment in the rotation trajectory, construct the intervention area mapping boundary, and generate a path segmentation and positioning map; The segmented trajectory reconstruction module is used to implement S3: construct a rotation trajectory structure based on the path segmentation and positioning map, divide the rotation control segments with associated directions, adjust the connection order and direction consistency, and generate a segmented rotation trajectory sequence; The offset analysis module is used to implement S4: call the segmented rotation trajectory sequence, determine the consistency with the flap force state, re-extract the shearing direction and shearing contact distribution area after simulation execution, compare the error offset, and generate a segment execution offset map; The terminal adaptation optimization module is used to implement S5: based on the fragment execution offset map, adjust the rotation direction of the flap end segment and the termination overlap boundary of the fitting path, align the tissue contact morphology of the end segment of the rotation trajectory with the edge of the burn wound, and generate the burn wound flap transplantation adaptation result.

Images