Skip to content

Project to Cross Section tool issue where OrientationPoints not outputting features with apparent dips calculated to the feature class of the Original GeMS database #137

@PhilMillerNMT

Description

@PhilMillerNMT

The tool completes successfully but is not outputting the features to the CrossSection"Output name token" feature class after creating the scratch feature class "CSA_ProfilePartsOrientationPoints_located". It creates the empty feature class—CS"Output name token"_OrientationPoints—with unknown SR, in the original gdb. However, that is the last step in the script that I can verify completes successfully. The tool states it completes successfully. The tool simply is not outputting the MZ profile apparent dip orientation point features into the original gdb in CrossSection"Output name token"\CS"Output name token"_OrientationPoints. The warning message output is:

"9could not create feature from objectid {in_row[oid_i]} in {loc_points}"
SIDEBAR—Could we also fix this "9" typo? Yes I can on my local copy of the script, but it would be nice to get this corrected for everyone.

I have attached the output tool history log.
ProjectMapDataToCrossSection.txt

So I think that means it is failing in this step of the code:

convert from cartesian to geographic angle

        csAzi = cartesianToGeographic(in_row[n("LOC_ANGLE")])
        vals[n("LocalCSAzimuth")] = csAzi
        vals[n("DistanceFromSection")] = in_row[n("Distance")]
        if isOrientationData:
            # vals[n("MapAzimuth")] = in_row[n("Azimuth")]
            if isAxial(in_row[n("Type")]):
                appInc, oblique = apparentPlunge(
                    in_row[n("Azimuth")], in_row[n("Inclination")], csAzi
                )
                dip_direction = in_row[n("Azimuth")]
            else:
                appInc, appIncVE, oblique = apparentDip(
                    in_row[n("Azimuth")], in_row[n("Inclination")], csAzi
                )
                dip_direction = in_row[n("Azimuth")] + 90
                # inclination_direction = bearing_sum(in_row[n("Azimuth")], 90)
                if dip_direction > 360:
                    dip_direction = dip_direction - 360
            # plotAzi = plotAzimuth(inclinationDirection, csAzi, appInc)
            plotAzi = symbol_rotation(dip_direction, csAzi, appIncVE)

            rowObliquity = round(oblique, 2)
            vals[n("DipDirection")] = round(dip_direction, 2)
            vals[n("ApparentInclination")] = round(appInc, 2)
            vals[n("ApparentInclinationVE")] = round(appIncVE, 2)
            vals[n("PlotRotation")] = round(plotAzi, 2)
            vals[n("Obliquity")] = rowObliquity

        # {out_name}_ID is last item
        vals[-1] = f"{id_pref}{str(i)}"
        try:
            out_rows.insertRow(vals)
        except Exception as e:
            guf.addMsgAndPrint(
                f"9could not create feature from objectid {in_row[oid_i]} in {loc_points}",
                1,
            )
            guf.addMsgAndPrint(e)

    # change _ID field to ID field
    alter_id(fc_name, flds, out_path)

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions