DefineProjection_managementでPRJが作成されませんでした

Question

私たちはAutoCAD DWGファイルをSDEフィーチャクラスにインポートするpythonスクリプトを持っています。これらのDWGファイルは様々な第三者によって作成され、インポートのために送信されます。インポートプロセスの一つのステップはDWGのPRJファイルを作成することです。ほとんどの場合、これは問題なく動作しますが、約5％のDWGでPRJファイルの作成時に問題が発生します。PRJの作成にはarcpy.DefineProjection_managementを使用しており、常に実行されますが、必ずしもPRJを作成しません。作成されない場合でもエラーは発生しません。コードは単にPRJを作成せずに静かに進行します。何が問題か知っている方や同様の問題を経験した方はいらっしゃいますか？

Pythonバージョン: 2.7 と 3.13.7

Arcpyバージョン: 10.9.1 と 3.6.4

SDEバージョン: 10.9.1

def project(strFullWorkingFilePath, dblScaleFactor):

    try:

# スケールファクター: 描画を移動するために偽の北方向と東方向に掛けられます。

        # 空間参照計算にはfloatではなくDecimalデータ型を使用します。

        sfi = Decimal(dblScaleFactor)

        sf = Decimal(1.0)\/Decimal(sfi)

        fe_4204 = Decimal(1968500.0) * Decimal(sf)

        fn_4204 = Decimal(13123333.33333333) * Decimal(sf)

# CADファイルのPRJファイルに格納される座標系文字列を構築します。

        Coordinate_System = "PROJCS['NAD_1983_StatePlane_Texas_South_Central_FIPS_4204_Feet',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Lambert_Conformal_Conic'],PARAMETER['False_Easting'," + str(fe_4204) + "],PARAMETER['False_Northing'," + str(fn_4204) + "],PARAMETER['Central_Meridian',-99.0],PARAMETER['Scale_Factor'," + str(sf) + "],PARAMETER['Standard_Parallel_1',28.38333333333333],PARAMETER['Standard_Parallel_2',30.28333333333333],PARAMETER['Latitude_Of_Origin',27.83333333333333],UNIT['Foot_US',0.3048006096012192]]"

# ツール: Define ProjectionはCADファイルと並んで正しいPRJファイルを作成します。

        arcpy.DefineProjection_management(strFullWorkingFilePath, Coordinate_System)

Tonyalmeida · Answer

まず、"D:\Working\Drawing.dwg\Polyline" を指していることを確認してください

prj が次のようなもので作成されていることを確認してください

import os
import arcpy

prj_path = os.path.splitext(strFullWorkingFilePath)[0] + ".prj"

arcpy.AddMessage("========================================")
arcpy.AddMessage("DWG: {}".format(strFullWorkingFilePath))
arcpy.AddMessage("DWG 存在: {}".format(os.path.exists(strFullWorkingFilePath)))
arcpy.AddMessage("PRJ: {}".format(prj_path))
arcpy.AddMessage("PRJ 存在 BEFORE: {}".format(os.path.exists(prj_path)))

try:
    arcpy.DefineProjection_management(
        strFullWorkingFilePath,
        Coordinate_System
    )

arcpy.AddMessage("DefineProjection は正常に返されました。")

except Exception as ex:
    arcpy.AddError("DefineProjection に失敗しました: {}".format(ex))
    arcpy.AddError(arcpy.GetMessages(2))
    raise

arcpy.AddMessage("PRJ 存在 AFTER: {}".format(os.path.exists(prj_path)))

if os.path.exists(prj_path):
    arcpy.AddMessage(
        "PRJ サイズ: {} バイト".format(os.path.getsize(prj_path))
    )
else:
    arcpy.AddError(
        "DefineProjection は正常に返されましたが、PRJ が存在しません！"
    )

次のように投影法を定義します。

arcpy.DefineProjection_management(
    strFullWorkingFilePath,
    Coordinate_System
)

prj_path = os.path.splitext(strFullWorkingFilePath)[0] + ".prj"

if not os.path.exists(prj_path):
    arcpy.AddWarning(
        "DefineProjection は完了しましたが、PRJ は作成されませんでした: {}".format(prj_path)
    )
    raise RuntimeError(
        "PRJ ファイルが作成されませんでした: {}".format(prj_path)
    )

arcpy.AddMessage("PRJ が正常に作成されました: {}".format(prj_path))

NAD 1983 StatePlane Texas South Central FIPS 4204 Feet の一般的な Esri WKID は 2278 です。

import arcpy
from decimal import Decimal

sfi = Decimal(str(dblScaleFactor))
sf = Decimal("1") / sfi

fe_4204 = Decimal("1968500.0") * sf
fn_4204 = Decimal("13123333.33333333") * sf

sr = arcpy.SpatialReference(2278)

print(sr.name)
print(sr.factoryCode)