花見川流域を歩く HANAMIGAWA RYUIKI wo ARUKU

2025年7月11日金曜日

土坑SK082（千葉市有吉北貝塚）3Dモデル

Pit SK082 (Ariyoshikita Shell Mound, Chiba City) 3D Model

I created a 3D model of a pit with two holes dug at the bottom that overhang the wall. It seems worthwhile to investigate why the holes overhang the wall.

底面に、壁面に対してオーバーハングしている穴が2つ掘られた土坑の3Dモデルを作成しました。穴が壁面に対してオーバーハングしている理由の検討は価値があるように感じます。

1　土坑SK082（千葉市有吉北貝塚）3Dモデル

土坑SK082（千葉市有吉北貝塚）3Dモデル by arakiminoru on Sketchfab

土坑SK082（千葉市有吉北貝塚）3Dモデル

発掘調査報告書資料

3DF Zephyr v8.013でアップロード

3Dモデルの動画

3Dモデルの画像

2　メモ

自作Blenderアドオン「Bezier Z Interpolator」を使うことにより、精度の高い3Dモデリングが可能となり、3Dモデリング作業に対する納得感が少しずつ出てきました。

3Dモデルは底面の2つの穴のうち1つを除いたモデルを作成し、除いた穴のモデルを別につくり、ブーリアン（アドオンbool tool利用）で合成しました。

参考　資料の分析

土坑底面に、壁面に対してオーバーハングしている穴が2つ掘られています。2つの穴が底面に対してオーバーハングしている理由の検討は価値があるように直観します。方位に関連して風除けとか日射量などに関連するのでしょうか？そもそもこの土坑の利用の仕方自体が自分には全くわかっていません。当面は、考古学的考察ではなく、3Dモデリング技術習得に焦点を当てた活動を展開することにします。

縄文遺構実測図の縁取線（地形変換線）の高度線形補間Blenderアドオン

Blender add-on for linear interpolation of the border lines (terrain transformation lines) of Jomon ruins survey maps

The border lines (terrain transformation lines) of Jomon ruins survey maps are not contour lines, but change in altitude. I created a Blender add-on that linearly interpolates the altitude distribution of the entire border line from several known point altitudes and displays it in the Blender3D viewport. An essential basic tool for 3D modeling of Jomon ruins has been completed.

縄文遺構実測図の縁取線（地形変換線）は等高線ではなく、高度が変化します。判明している幾つかのポイント高度から縁取線全体の高度分布を線形補間して、Blender3Dビューポートに表示するBlenderアドオンを作成しました。縄文遺構3Dモデリングの必須基礎ツールが完成しました。

1　縁取線高度線形補間Blenderアドオンの機能

このアドオンでは、3Dビューポートの真上からのビュー（テンキー7）で描いた縁取線（べジェ曲線）の頂点インデックスに2つ以上の既知の高度数値を手入力します。アドオンは、入力された既知高度数値に基づき、それ以外の頂点インデックスに線形補間（内挿、外挿）で高度数値を生成し、縁取線頂点高度が3Dビューポートで移動します。

2　縁取線高度線形補間Blenderアドオンの使い方

2-1　縁取線の作成・表示

Blender3Dビューポートに縄文遺構実測図を真上からのビュー（テンキー7）でスケール通りに配置します。この実測図を下敷きにして目的の縁取線をBezier Toolkitを使ってクリックによりべジェ曲線で生成します。

生成したべジェ曲線を編集モードで選択しておきます。

なお、閉じたべジェ曲線の場合は「オブジェクトデータプロパティ→アクティブスプライン→ループのUチェックを外す」操作で開いたべジェ曲線にしておきます。

2-2　アドオンでの高度線形補間

2-2-1　3Dビューポートにおけるインデックス（番号）の表示

アドオンの「インデックス表示切替」をクリックして、3Dビューポートの縁取線（べジェ曲線）頂点にインデックス（番号）を表示します。

2-2-2　初期化

「初期化」をクリックして、ID欄（頂点番号欄）とZ値欄（標高欄）をアドオン画面に表示します。

ID欄は全て表示されていますが、Z値欄は空欄となっています。

2-2-3　既知高度の手入力

実測図から判明している既知高度を関連する頂点のZ値欄に手入力します。補間のために、2つ以上の値を入力する必要があります。

既知高度を入力した様子（54頂点に対して既知高度を4頂点に手入力）

既知高度を入力した様子（高度は全て0となっている）

2-2-4　補間とその反映

「補間して反映」をクリックすると、Z値欄の空欄部分に線形補間された数値（高度）が生成されるとともに、3Dビューポートでべジェ曲線頂点が補間高度に移動します。

線形補間した様子（X座標、Y座標は不動）

線形補間した様子（頂点のZ座標が補間高度に移動している）

参考　土坑実測図（土坑肩線の高度が読みとれるのは4点）

2-2-5　やり直し

不都合があれば、「初期化」をクリックして、Z値欄に既知高度を入力し直し、「補間して反映」をクリックします。

2-2-6　確定

不都合がなければ、「確定」で結果を確定します。

元が閉じたべジェ曲線であった場合は、「オブジェクトデータプロパティ→アクティブスプライン→ループのUをチェックする」で開いたべジェ曲線から閉じたべジェ曲線に戻します。

3　メモ

高度線形補間を行った縁取線を使うことによって、縄文遺構3Dモデリングの精度を上げることができます。

このアドオンはBlenderバージョン4.4.3を対象にしています。

4　感想

このアドオンは「べジェ曲線の頂点インデックス（番号）表示」BlenderPythonスクリプト及び、「内挿と外挿を含む線形補間」Pythonスクリプトを最初に作成し、その2つを主な要素として新たに構成したPythonスクリプトです。それぞれの段階でのChatGPTとのやりとりはいつも以上に回数が多くなりました。恐らく、このアドオンの機能を最初からChatGPTに求めても、完成に至らなかったと想像します。

5　Pythonスクリプト

bl_info = {
    "name": "Bezier Z Interpolator",
    "blender": (4, 4, 3),
    "category": "Curve",
}

import bpy
import blf
from bpy_extras import view3d_utils
import numpy as np

handler = None

def draw_callback_px(self, context):
    obj = context.active_object
    if not obj or obj.type != 'CURVE' or obj.mode != 'EDIT':
        return
    region = context.region
    rv3d = context.space_data.region_3d
    curve = obj.data

    font_id = 0
    blf.size(font_id, 14)
    for spline in curve.splines:
        if spline.type != 'BEZIER':
            continue
        for i, bez in enumerate(spline.bezier_points):
            co_world = obj.matrix_world @ bez.co
            screen_coord = view3d_utils.location_3d_to_region_2d(region, rv3d, co_world)
            if screen_coord:
                blf.position(font_id, screen_coord.x, screen_coord.y, 0)
                blf.draw(font_id, str(i))

def register_draw_handler():
    global handler
    if handler is None:
        handler = bpy.types.SpaceView3D.draw_handler_add(
            draw_callback_px, (None, bpy.context), 'WINDOW', 'POST_PIXEL'
        )

def unregister_draw_handler():
    global handler
    if handler is not None:
        bpy.types.SpaceView3D.draw_handler_remove(handler, 'WINDOW')
        handler = None

def interpolate_z(values):
    ids = np.array([i for i, v in enumerate(values) if v != ""])
    zs = np.array([float(v) for v in values if v != ""])
    if len(ids) < 2:
        return values
    f = np.interp(range(len(values)), ids, zs)
    return [str(round(z, 4)) for z in f]

class CurveZEntry(bpy.types.PropertyGroup):
    z: bpy.props.StringProperty(name="Z", default="")

class CURVE_PT_z_interpolator_panel(bpy.types.Panel):
    bl_label = "Z値補間ツール"
    bl_idname = "CURVE_PT_z_interpolator_panel"
    bl_space_type = 'VIEW_3D'
    bl_region_type = 'UI'
    bl_category = "Bezier Z"

    def draw(self, context):
        layout = self.layout
        scene = context.scene
        layout.operator("curve.toggle_index_display", text="インデックス表示切替")
        layout.operator("curve.z_initialize", text="初期化")

        layout.label(text="Z値入力:")
        for i, entry in enumerate(scene.bezier_z_entries):
            row = layout.row()
            row.label(text=f"ID {i}")
            row.prop(entry, "z", text="")

        layout.operator("curve.z_interpolate", text="補間して反映")
        layout.operator("curve.z_confirm", text="確定")

class CURVE_OT_z_initialize(bpy.types.Operator):
    bl_idname = "curve.z_initialize"
    bl_label = "Zエントリ初期化"

    def execute(self, context):
        obj = context.active_object
        if not obj or obj.type != 'CURVE':
            self.report({'WARNING'}, "カーブオブジェクトを選択してください")
            return {'CANCELLED'}
        context.scene.bezier_z_entries.clear()
        for spline in obj.data.splines:
            if spline.type == 'BEZIER':
                for _ in spline.bezier_points:
                    context.scene.bezier_z_entries.add()
        self.report({'INFO'}, "Zエントリを初期化しました")
        return {'FINISHED'}

class CURVE_OT_z_interpolate(bpy.types.Operator):
    bl_idname = "curve.z_interpolate"
    bl_label = "Z値補間"

    def execute(self, context):
        entries = context.scene.bezier_z_entries
        z_list = [e.z for e in entries]
        z_interp = interpolate_z(z_list)

        obj = context.active_object
        if obj and obj.type == 'CURVE':
            index = 0
            mw = obj.matrix_world
            imw = mw.inverted()

            for spline in obj.data.splines:
                if spline.type != 'BEZIER':
                    continue

                bez_points = spline.bezier_points
                n = len(bez_points)

                for i, bez in enumerate(bez_points):
                    if index >= len(z_interp):
                        break

                    # グローバル座標でZ補正
                    world_co = mw @ bez.co
                    world_co.z = float(z_interp[index])
                    bez.co = imw @ world_co
                    entries[index].z = z_interp[index]
                    index += 1

                # 始点・終点のハンドルをたたむ
                if n >= 2:
                    p_start = spline.bezier_points[0]
                    p_end = spline.bezier_points[-1]

                    p_start.handle_left = p_start.co.copy()
                    p_start.handle_right = p_start.co.copy()
                    p_end.handle_left = p_end.co.copy()
                    p_end.handle_right = p_end.co.copy()

        return {'FINISHED'}

class CURVE_OT_z_confirm(bpy.types.Operator):
    bl_idname = "curve.z_confirm"
    bl_label = "Z値確定"

    def execute(self, context):
        self.report({'INFO'}, "Z値を確定しました")
        return {'FINISHED'}

class CURVE_OT_toggle_index_display(bpy.types.Operator):
    bl_idname = "curve.toggle_index_display"
    bl_label = "インデックス表示切替"

    def execute(self, context):
        if handler is None:
            register_draw_handler()
            self.report({'INFO'}, "インデックス表示を開始")
        else:
            unregister_draw_handler()
            self.report({'INFO'}, "インデックス表示を停止")
        return {'FINISHED'}

classes = (
    CurveZEntry,
    CURVE_PT_z_interpolator_panel,
    CURVE_OT_z_initialize,
    CURVE_OT_z_interpolate,
    CURVE_OT_z_confirm,
    CURVE_OT_toggle_index_display,
)

def register():
    for cls in classes:
        bpy.utils.register_class(cls)
    bpy.types.Scene.bezier_z_entries = bpy.props.CollectionProperty(type=CurveZEntry)
    # register() では active_object にアクセスしない！

def unregister():
    for cls in reversed(classes):
        bpy.utils.unregister_class(cls)
    unregister_draw_handler()
    del bpy.types.Scene.bezier_z_entries

if __name__ == "__main__":
    register()

2025年7月9日水曜日

Blenderべジェ曲線の頂点インデックス（番号）表示

Blender Bezier curve vertex index (number) display

Blender Bezier curve vertex index (number) display is now possible. This allows efficient and accurate manipulation of Bezier curves. This BlenderPython has been completed for the first time in a month.

Blenderべジェ曲線の頂点インデックス（番号）表示ができるようになりました。これでべジェ曲線の操作が効率的かつ正確にできるようになります。1カ月ぶりにこのBlenderPythonが完成しました。

1　Blenderべジェ曲線の頂点インデックス（番号）表示

Blenderべジェ曲線の頂点インデックス（番号）表示

この例は閉じたべジェ曲線で、0～77のインデックスが表示されています。

2　Blenderバージョン4.4.3用の頂点インデックス表示BlenderPythonスクリプト

import bpy
import blf
from bpy_extras import view3d_utils

handler = None

def draw_callback_px(self, context):
    obj = context.active_object
    if not obj or obj.type != 'CURVE' or obj.mode != 'EDIT':
        return

    region = context.region
    rv3d = context.space_data.region_3d
    curve = obj.data

    font_id = 0
    blf.size(font_id, 14)

    for spline in curve.splines:
        if spline.type != 'BEZIER':
            continue

        for i, bez in enumerate(spline.bezier_points):
            if not bez.select_control_point:
                continue

            # 3D位置 → スクリーン座標
            co_world = obj.matrix_world @ bez.co
            screen_coord = view3d_utils.location_3d_to_region_2d(region, rv3d, co_world)

            if screen_coord:
                blf.position(font_id, screen_coord.x, screen_coord.y, 0)
                blf.draw(font_id, str(i))

def register_draw_handler():
    global handler
    if handler is None:
        handler = bpy.types.SpaceView3D.draw_handler_add(
            draw_callback_px, (None, bpy.context), 'WINDOW', 'POST_PIXEL'
        )
        print("インデックス表示を開始しました。")

def unregister_draw_handler():
    global handler
    if handler is not None:
        bpy.types.SpaceView3D.draw_handler_remove(handler, 'WINDOW')
        handler = None
        print("インデックス表示を停止しました。")

class CURVE_OT_toggle_index_display(bpy.types.Operator):
    bl_idname = "curve.toggle_index_display"
    bl_label = "ベジェインデックス表示切り替え"

    def execute(self, context):
        if handler is None:
            register_draw_handler()
            self.report({'INFO'}, "表示開始")
        else:
            unregister_draw_handler()
            self.report({'INFO'}, "表示停止")
        return {'FINISHED'}

def register():
    bpy.utils.register_class(CURVE_OT_toggle_index_display)

def unregister():
    bpy.utils.unregister_class(CURVE_OT_toggle_index_display)

if __name__ == "__main__":
    register()
    bpy.ops.curve.toggle_index_display()

3　頂点インデックス（番号）表示の意義

Blenderで、XY座標平面に縄文遺構の幾つかの縁取線をトレースし、それらの縁取線を該当する標高に移動し、その縁取線からアドオンCurves To Meshで3Dモデルを作成します。しかし、縁取線は等高線ではなく、標高が変化する線分です。その標高変化を表現するには、標高が判明している少数の頂点の間をトポロジー的に補間する必要があります。その補間作業のために、標高判明頂点を把握する必要があり、どうしても、全体の頂点インデックスが必要になります。

このBlenderPythonスクリプトにより、その頂点インデックスが3Dビューポートで一見して判るようになりました。

縄文遺構3Dモデル構築作業にとって重要な基礎ツールの一つができました。

4　メモ

このBlenderPythonスクリプトは1カ月程前にChatGPTにより作成チャレンジして、半分失敗しました。当時作成したスクリプトはBlenderバージョン3では使えましたが、バージョン4では使えませんでした。失敗した理由は、ChatGPTにBlenderバージョンを明示しなかったからであると推察します。

今回はBlenderバージョン4.4.3で使うことを前提に指示しました。ChatGPTとのやりとりの中で、Blenderバージョン4.4の特性に関する扱いが幾つか出てきました。このような経緯から、Blenderのバージョンを明示しなければ、今回も失敗していた可能性があります。

べジェ曲線カラー化でやる気が出る

Coloring Bezier curves motivates me

I'm working hard to learn the basics of Blender in order to proceed with the project to create a 3D model of Jomon ruins. While coloring Bezier curves in the Blender 3D viewport, I found that they were easier to operate than I expected, and I felt motivated. I've decided to master Bezier curve manipulation techniques.

縄文遺構3Dモデル化プロジェクトを進めるために、Blender基礎技術習得に励んでいます。その中でBlender3Dビューポートでべジェ曲線に色を付けてたところ、思いのほか操作しやすくなり、「やる気」が湧き出てきました。べジェ曲線操作技術を極めることにします。

1　Blender3Dビューポートにおけるべジェ曲線の表現

オブジェクトモード　デフォルトのべジェ曲線（非選択）表現（黒線でほとんど見分けられない）

オブジェクトモード　新たに塗色したべジェ曲線（非選択）表現

オブジェクトモード　デフォルトのべジェ曲線（選択）表現

編集モード　新たに塗色したべジェ曲線（非選択）表現

編集モード　デフォルトのべジェ曲線（選択）表現

オブジェクトモード、編集モード共に、べジェ曲線（カーブ）は非選択デフォルトでは黒色となり、背景と識別しにくく、とても扱いづらいものです。選択すると判りやすくなるので、これまでは非選択での見にくさは無意識に我慢してきました。しかし、べジェ曲線の扱いについて意識的に学習を深める中で、思い切って色を変えてみました。オブジェクトモードでは黄緑に、編集モードでは水色にしました。色を変えた途端に選択していなくてもその存在が明瞭に判るようになりました。

同時に、べジェ曲線の扱いをとことん極めてみたいという欲望が生まれました。「やる気」が突然湧き出して来たということです。自分の心理機構がかなり単純であることに驚きました。その心理分析ははいつか深めることにして、今は何はともあれべジェ曲線の特性を極めて、縄文遺構3Dモデル作成技術の基礎を構築することにします。

Blenderでのべジェ曲線色変更画面

プリファレンス→テーマ→3Dビューポート→ワイヤーフレーム（及びワイヤ編集）

2　べジェ曲線の操作技術メモ

最近確認したべジェ曲線操作技術をメモします。

2-1　べジェ曲線を閉じる（ループをつくる）

2-1-1　Bezier Toolkitによる操作

自分はべジェ曲線はアドオンBezier Toolkitを使って描いています。このアドオンの操作では閉じる時にはALTキーを押しながらクリックすると最短にある頂点にスナップして、べジェ曲線を閉じることができます。

2-1-2　Fキーによる辺連結

既にある開いたべジェ曲線を閉じるときには両端の頂点を選択してFキーを押せば、閉じることができます。

2-2　べジェ曲線を開く

2-2-1　辺の削除

閉じたべジェ曲線の隣接する2つの頂点を選択して、X→セグメントで2つの頂点の間の辺が削除されます。

2-2-2　アクティブスプラインの操作

べジェ曲線が選択された状態で、「オブジェクトデータプロパティ→アクティブスプライン→ループのUチェックを外す」操作をすると、閉じたべジェ曲線の一部の辺が表示されなくなり、開いたべジェ曲線になります。Uチェックを入れると再び閉じたべジェ曲線になります。

この操作の説明はBlenderマニュアルにあるのですが、自分にはその意味がさっぱり判りません。Blender初心者の醍醐味です。なぜそのような機能があるのかという意味が判らないのは初心者である今だけだと思います。

2-3　べジェ曲線を分割する

2-3-1　コピー、削除の繰り返し

べジェ曲線を複数コピーして、それぞれ別の場所の複数頂点を選択してX→セグメントで削除すれば、結果としてべジェ曲線を分割できます。

2-3-2　Bezier Toolkitによる操作

閉じたべジェ曲線が1つの平面に存在している場合、その同じ平面に別の閉じたべジェ曲線を描くいてその二つに交点があると、2つのべジェ曲線は交点により切断されてバラバラの開いたべジェ曲線になります。不用なべジェ曲線を削除すれば、当該べジェ曲線を分割したことになります。Illustratorのパスファインダーのような機能がBezier Toolkitにあるということです。

2-4　ハンドルタイプの変更

縄文遺構を実測図からトレースする時は操作の単純化、正確化のためにべジェ曲線ハンドルタイプをベクトル（カクカクの線分）にしています。モデリングのためのメッシュ化ではハンドルタイプを整列にして滑らかな線分で遺構を表現することにします。

2-5　Z座標値の補間

自分は現在、数本の閉じたべジェ曲線から、アドオンCurves To Meshでメッシュ生成することを3Dモデル作成の基本としています。この操作になかで、閉じたべジェ曲線のZ座標値の補間が重要です。その補間方法について、BlenderPythonスクリプトを作成したましたので、次の記事で紹介します。

2025年7月5日土曜日

縄文陥穴に設置された欺瞞装置の想定

Assumption of a deception device installed in the Jomon Pitfall

The deception device installed in the Jomon Pitfall at Ariyoshikita Shell Mound was assumed from the actual measurement drawing and added to the 3D model.

有吉北貝塚の縄文陥穴に設置された欺瞞装置を実測図から想定して、3Dモデルに追記しました。

1　有吉北貝塚の縄文陥穴に設置された欺瞞装置想定

有吉北貝塚の縄文陥穴に設置された欺瞞装置想定 by arakiminoru on Sketchfab

有吉北貝塚の縄文陥穴に設置された欺瞞装置想定

元3Dモデルは「陥穴TP005（千葉市有吉北貝塚）3Dモデル」

想定した欺瞞装置

想定の根拠とした実測図に表現された四隅掘り込み跡

3DF Zephyr v8.013でアップロード

3Dモデルの動画

3Dモデルの画像

2　メモ

実測図に表現される四隅掘り込み跡を陥穴欺瞞装置と解釈しました。具体的には掘り込み跡として表現される事象を土中に打ち込まれた4本の杭跡と解釈しました。4本の杭上面が、小枝や草や土から構成される欺瞞地表を支える2本の棒の支えであると解釈しました。

2025年7月3日木曜日

ブログ過去記事リストの更新（4293記事）と過去記事リスト改善策検討

Updating the blog's past article list (4293 articles) and considering ways to improve the past article list

It's been a year since I updated the past article list for the blog "Walking the Hanami River Basin" (4293 articles). There's a limit to not managing a blog. I apologize to everyone. As part of a consideration of ways to improve the past article list page, I created a Python script to set up category-based classification tags.

ブログ「花見川流域を歩く」の過去記事リストの更新（4293記事）を1年ぶりにしました。ブログを管理しないテイタラクにも程があるというものです。皆様に謝ります。過去記事リストページの改善策検討の一環として、カテゴリー別大分類タグ設置のPythonスクリプトを作成してみました。

1　ブログ「花見川流域を歩く」の過去記事リストの更新

ブログ「花見川流域を歩く」の過去記事リストを更新しました。過去記事総数は4293記事となりました。

過去記事リスト

ブログを毎月管理しないというテイタラクにも程があるというものです。皆様に謝ります。

過去記事リストページ

2　過去記事リストページの改善策検討

2-1　カテゴリ別タグ導入

過去記事リストページの改善策を検討しました。その一つとしてカテゴリ別タグの導入（例：#考古学 #Blender #古代交通）について検討してみました。全記事を10〜15の大分類でタグ付けし、タグで記事にスピーディーに到達できるようにする工夫です。4293と記事数が多くても読者が関連トピックをすぐに絞り込めるようにしてみたいと思います。

2-2　カテゴリ別大分類タグの試案（ChatGPT作成）

4293記事表題リストからカテゴリ別大分類タグ試案をChatGPTにつくってもらいました。大雑把なイメージでいえば、このような感じの大分類タグをつくり、過去記事リストを検索できるようにしたのです。

このカテゴリ別大分類タグの試案はChatGPTが作成したもので、自分の頭脳で思考したものではありません。あくまでも見本です。

【1】花見川・地形・自然観察

地形学的視点、河川争奪、地名、地形3Dモデル、地質、景観変遷

【2】考古遺跡研究（加曽利・有吉・大膳野など）

千葉市内の遺跡における考古学的観察と記録

【3】縄文土器・石器観察と3Dモデル

土器・石器・貝製品・装飾品等の3D観察記録、型式分類

【4】土偶・埴輪・古墳研究

縄文〜古墳時代の人物・動物表現物、埴輪や古墳の構造・象徴性分析

【5】3Dモデルと視覚化技術（Blender・点群・密度図など）

3D分布モデル、密集度ヒートマップ、Blender活用、視覚化工夫

【6】遺物データベース・分析・統計処理

遺物台帳電子化、QGIS・PostgreSQL・Pythonを用いた定量分析

【7】儀礼考古学・象徴研究・仮説構築

コンテクスト分析、儀礼・祭祀・象徴行動の推定

【8】海外地形・歴史・文化との比較

中央アジア・イラン等の地形や古代構造との比較研究

【9】考古学切手・文化財に関する趣味活動

考古モチーフ切手、文化財モデリング、展示会観覧

【10】技術メモ・Python/Blenderスクリプト

作業効率化、技術的工夫、スクリプトの実用例

【11】学習・講座・読書・理論深化

書籍・雑誌・講演会などから得た考古学的知見の紹介・思索

【12】博物館・特別展・現地観覧記

展示会・特別展・発掘現場等の訪問と所感

【13】月間／年間ふりかえり・メモ

活動総括、習得のふりかえり、ブログ運営記念

【14】AI・デジタル技術との対話・試行

ChatGPT、OCR、AI画像生成など新技術との実験的対話

【15】その他・雑記・技術遊び

技術遊びや生活小話など、分類が難しいユニーク記事

2-3　カテゴリ別大分類タグ（試案）に基づく4293記事の自動分類Pythonスクリプト

よき友であるChatGPTにカテゴリ別大分類タグ（試案）に基づく4293記事の自動分類Pythonスクリプトを作成してもらいました。1つの記事に複数のタグを付ける多重タグ制のスクリプトです。

import csv

# 分類カテゴリとキーワード辞書（簡易版）
category_keywords = {
    "花見川・地形・自然観察": ["花見川", "地形", "河川", "堀割", "印旛", "干渉色", "古地図"],
    "考古遺跡研究": ["貝塚", "遺跡", "遺構", "有吉", "加曽利", "大膳野", "北斜面"],
    "縄文土器・石器観察と3Dモデル": ["土器", "石器", "観察記録3Dモデル", "破片", "深鉢", "小型", "両耳壺"],
    "土偶・埴輪・古墳研究": ["土偶", "埴輪", "古墳", "遮光器", "人形塚", "武人", "祖霊"],
    "3Dモデルと視覚化技術": ["3D", "Blender", "ヒートマップ", "点群", "分布図", "可視化", "geometry nodes"],
    "遺物データベース・分析・統計処理": ["台帳", "QGIS", "PostgreSQL", "統計", "データベース", "頻度分布", "標高分布"],
    "儀礼考古学・象徴研究・仮説構築": ["儀礼", "祭祀", "骨", "装飾品", "散乱", "ミニチュア", "コンテクスト"],
    "海外地形・歴史・文化との比較": ["イラン", "ルート砂漠", "岩塩", "中央アジア", "サルチャハン", "アラル海"],
    "考古学切手・文化財に関する趣味活動": ["切手", "考古学切手", "ツタンカーメン", "スキタイ", "文化財"],
    "技術メモ・Python/Blenderスクリプト": ["スクリプト", "Python", "メモ", "テクスチャ", "非破壊", "干渉色"],
    "学習・講座・読書・理論深化": ["講座", "読書", "学習", "理論", "方法論", "物質性"],
    "博物館・特別展・現地観覧記": ["博物館", "観覧", "展示", "展", "発掘調査", "展示状況"],
    "月間／年間ふりかえり・メモ": ["ふりかえり", "活動", "周年", "明けまして", "通過"],
    "AI・デジタル技術との対話・試行": ["ChatGPT", "OCR", "生成AI", "AI", "感動"],
    "その他・雑記・技術遊び": ["遊ぶ", "楽しむ", "雑記", "故障", "トピック", "Photoshop", "Illustrator"]
}

# タイトルを読み込む（1記事1行）
with open("articles.txt", "r", encoding="utf-8") as f:
    lines = [line.strip() for line in f if line.strip()]

# 出力準備
with open("classified_articles.csv", "w", newline="", encoding="utf-8") as csvfile:
    writer = csv.writer(csvfile)
    writer.writerow(["記事タイトル", "分類タグ"])

    for title in lines:
        matched_categories = []
        for category, keywords in category_keywords.items():
            if any(keyword in title for keyword in keywords):
                matched_categories.append(category)
        writer.writerow([title, ", ".join(matched_categories) if matched_categories else "未分類"])

2-4　4293記事の自動分類結果

4293記事の自動分類はcsvファイルで書き出されます。

4293記事の自動分類結果（部分）

この結果をブログページに実装するにはまだ幾つかステップが必要ですが、それがChatGPTの指南でできることが確信できたので、自分にとってとても素晴らしい活動になりました。

2-5　カテゴリ別大分類タグ作成Pythonスクリプトについて

今回の作業ではカテゴリ別大分類タグの試案を直接ChatGPTに作成してもらいましたが、それをPythonスクリプトで行い、そのPythonスクリプトを自分好みに改善するプロセスが大切であると考えます。追って、それをChatGPTに作成してもらうことにします。

2-6　カテゴリ別大分類タグ導入以外の改善策

次のような改善策が考えられますが、追って詳しく検討することにします。

・ページ専用のキーワード検索機能の実装

・テーマ別年表・マップ表示

・シリーズ記事の目次ページ作成

・画像付きインデックス作成

3　参考　ブログテーマの変遷

ChatGPTに4293記事表題からブログテーマの変遷をまとめてもらいました。あくまでもChatGPTが作成したもので、自分の頭脳の思考は一切関わっていない資料です。

全てのブログ記事は自分が書いたものといっても、それを全部咀嚼しなおして変遷を整理するには一定の時間が必要となります。その時間を捻出することは、現在自分の興味優先順位からして事実上不可能です。従って、この資料は自分にとっても意義があるものといわざるをえません。

●ブログ「花見川流域を歩く」テーマの時系列的まとめ（ChatGPT作成）

1　 2011～2013年：地域自然・地形・歴史への関心

初期は花見川流域の自然観察（桜・鳥など）や「地形・歴史（地峡、河川争奪、古代交通）」への探究が中心。

「花見川地峡」や「印旛沼筋河川争奪」の自然地史の考察が大きなテーマに。

読者参加型ではなく、筆者の興味に基づく地道な観察・考察を継続。

2　 2014～2017年：地形・遺跡とデジタル技術の融合

地理情報（GIS、3D地形）、地名、古代遺跡の関係性を深く追究。

特に「花見川」「加曽利貝塚」「大膳野南貝塚」など地域遺跡に焦点を当てる。

歴史研究とデジタルマッピング（3Dモデル、等高線解析等）が結びつく。

3　 2018～2020年：考古学的視点と実証分析の深化

竪穴住居・土器・人骨など遺構や遺物に関する検討を重ねる。

観察・仮説・検証のサイクルが形成され、研究的関心が強まる。

単なる趣味ブログから、市民考古学者的な探究活動に発展。

4　 2021～2023年：Blender・Pythonなど技術の本格活用

3D空間での遺物分布分析や可視化が主軸に。

「有吉北貝塚北斜面貝層」を継続的に分析し、点群・メッシュ・ヒートマップなどの高度表現を試みる。

BlenderやPythonによる独自ツール開発・習熟が目立つ。

5　 2024～2025年：考古学・視覚化・異文化地形への展開

イランや中央アジアの地形分析・3D可視化への興味拡大。

「儀礼考古学」「物質性」「類推法」などの概念学習から理論的深化も図る。

考古学切手、古墳地割線、埴輪の3D観察など文化財全般へ視野拡大。

興味の深め方は、単なる観察→モデル化→定量分析→抽象的思考と高度化している。

2025年7月1日火曜日

BlenderのprimitiveCUBEのデフォルトの大きさが2m×2m×2mである理由

Why is the default size of Blender's primitiveCUBE 2m x 2m x 2m?

I investigated why the default size of Blender's primitiveCUBE is 2m x 2m x 2m, not 1m x 1m x 1m. I understood that this is related to visibility and operability (visual sensation on a computer and operability with hand muscles) when constructing and generating a CUBE as a solid object at the origin of 3D space.

BlenderのprimitiveCUBEのデフォルトの大きさが1m×1m×1mではなく、2m×2m×2mである理由を調べました。3D空間原点で立体オブジェクトとしてのCUBEを構築生成する場合の視認・操作性（パソコンでの視認感覚・手筋肉操作性）が関わっていることを理解しました。

1　CUBE大きさに関するケアレスミス

有吉北貝塚北斜面貝層の3D座標取得遺物をBlender3DビューポートにプロットするBlenderPythonスクリプトを状況に応じて幾つか作成しています。遺物すべてを5㎝×5㎝×5㎝のCUBEで表現しています。

最近、プロット高速化改善策を施したBlenderPythonスクリプトを新たにつくり走らせたところ、遺物の大きさ（CUBEの大きさ）が異なってしまうミスを犯しました。ミスの原因は、BlenderにおけるCUBEのデフォルトの大きさは2m×2m×2mであると理解していましたが、primitiveCUBEと指定すれば、それはprimitiveCUBE(size=1)を意味すると勘違いしていたからでした。

このケアレスミスをキッカケに、BlenderではなぜCUBEのデフォルトの大きさ、つまりprimitiveCUBEの大きさが2m×2m×2mであり、区切りのよい単位と自分が考えた1m×1m×1mのCUBE（つまりprimitiveCUBE(size=1)）ではないのか、その理由をwebやChatGPTで調べてみました。

2　Blenderにおいて bpy.ops.mesh.primitive_cube_add() で生成されるCUBEのデフォルトサイズが size=2（= 2m×2m×2m）である理由

理由は3つほどあるようです。

2-1　Blenderの内部基準（Unit Scale）が「1 Blender Unit = 1m」であるため

Blenderではデフォルトの単位系は「1 Blender Unit = 1m」です。

その上で、size=2 という引数は「中心から各方向に1m（つまり全体で2m）」という意味になります。

Blenderでは空間原点（あるいは中心）との対応でオブジェクトを認識し操作しますから、原点（あるいは中心）とCUBEの大きさを考える時、次のような関係になります。

size=2 → -1m 〜 +1m の立方体（合計 2m）

size=1 → -0.5m 〜 +0.5m の立方体（合計 1m）

つまり、原点（あるいは中心）を意識し、それとの関係でCUBEの大きさを捉える時、size=2 という引数が「1 Blender Unit = 1m」に対応することになり、最も基礎的大きさであると考えることができます。

2-2　分かりやすい視認性を優先した設計

Blenderでは、画面立上げと同時に立方体（CUBE）が追加され、表示されます。この際、1m立方体だと小さすぎて3Dビューポートで目立たないことがあります。初心者ユーザーでも「見てすぐわかる形」が重要であり、そのため、デフォルトでやや大きめの 2mサイズを選んでいると考えられます。

2-3　Blender初期バージョンとの互換性

Blenderの初期のバージョンでは、実世界単位（mやcm）ではなく「Blender Unit」が独立した単位として使われていました。その当時、size=2（-1〜+1）で作られたCUBEは、原点を中心に対称な形で作るための定番でした。後から単位系（Metric）に対応した際に、古いスクリプトとの互換性を保つため、「デフォルトsize=2」が維持された可能性があります。

3　感想

BlenderのprimitiveCUBEのデフォルトサイズがsize=2である理由を知ることができました。その根本は原点（中心）から「1 Blender Unit = 1m」のサイズのCUBEという点にあります。それは、3D空間のなかでCUBEを構築生成する行為（パソコン画面で3D空間を認識して、マウスやキーボードを手筋肉で操作してCUBEをつくる活動）があってはじめて生まれた観念であると理解しました。腑に落ちました。

Blender3DビューポートにおけるprimitiveCUBE(size=1)とprimitiveCUBE(size=2)

2025年6月ブログ活動のふりかえり

Looking back on blog activities in June 2025

I looked back on the activities of the blog "Walking the Hanami River Basin" in June 2025.

Almost every day in June, I worked hard on inputting data for the Ariyoshi Kita Shell Mound in the morning, and in the afternoon, I enjoyed creating 3D models of the ruins and 3D models of the Iranian topography.

ブログ「花見川流域を歩く」の2025年6月活動をふりかえりました。

6月はほぼ毎日、午前中は有吉北貝塚入力作業に励み、午後は遺構3Dモデル作成やイラン地形3Dモデル作成を楽しみました。

1　ブログ「花見川流域を歩く」

・2025年6月の記事数は27です。

・ほぼ毎日午前中は有吉北貝塚北斜面貝層の遺物平面図から遺物の平面座標読取デジタル化作業に励みました。作業日誌記事を3編書きました。

・午後は有吉北貝塚の遺構（陥穴など）3Dモデル作成（記事6編）及びイラン地形3Dモデル作成（記事13編）を楽しみました。

・考古学切手記事1編を書きました。

2　ブログ「花見川流域を歩く　自然・風景編」

・早朝散歩記事を1編書きました。

3　2025年6月活動の特徴

3-1　有吉北貝塚北斜面貝層遺物平面座標読取デジタル化作業の本格スタート

有吉北貝塚北斜面貝層遺物平面座標読取デジタル化作業の本格スタートを切ることができました。毎日午前中は有無を言わせずに入力作業を行っています。地道な継続作業なくして6万4000件の遺物3Dデータベースは作成できません。（手入力作業による数万件以上のデータベース作成作業は千葉県遺跡分布データベース、千葉県小字分布データベースなどで既体験しているので、地道な作業の大切さは身に染みて理解しています。）

6月は作業の効率化に注力しました。これから本格化する作業をいかに効率化するか、作業詳細面でのPythonスクリプト活用の工夫改善を多数行いました。

3-2　遺構3Dモデル作成プロジェクトスタート

考古における3Dモデルという観点で自分の技術をふりかえると、遺物をフォトグラメトリで3Dモデルにする技術だけでは跛行感がぬぐい切れません。そこで、発掘調査報告書掲載実測図（平面図、断面図）から遺構（竪穴住居跡、土坑、陥穴、炉穴、…）の3Dモデル作成技術を身に着けるプロジェクトをスタートさせました。

6月は自分のBlenderによる3Dモデル作成技術がいかに稚拙であるか実感する日々となりました。その稚拙感は物事のスタートでいつも感じる、その分野専門技術にはじめて接触した時の実感です。しかし、よくよく考えてみると、現時点より技術が下がることはあり得ないので、後は専門技術が少しづつ向上するだけですから、未来（と言っていもごく短い未来しかありませんが）が楽しみです。

このような感想をくどくどと記事に書くことははばかりました。

3-3　イラン地形3Dモデル作成完結

自分最初の会社（応用地理を標榜する日本最初の建設コンサルタント会社）で一緒だった中家さんから雑誌地理4月号（イラン特集号）を3月に謹呈していただきました。2カ月遅れとなりましたが、その雑誌論説の感想文装いで、イラン関連記事を書きました。しかし、記事内実は自分のイラン地形3Dモデル作成を楽しむ活動となってしまいました。中家さんや著者の皆さんにはお詫びします。6月に戦争が勃発し、イランに特別の興味を持たざるを得ない状況下での地形3Dモデル作成作業は、いつまでも印象に残ると思います。

3-4　考古学切手記事

考古学切手記事1編を書きました。この記事を書く中で、考古学切手ストックをあらためて眺めてみると、自分の考えていた以上に豊なコレクションになっていることに気が付きました。今後考古学切手蒐集趣味活動を発展させることにします。まずはいままで蒐集したコレクションの本格整理に取り掛かることにします。

4　2025年7月活動の展望

7月は6月活動を継続します。午前中は有吉北貝塚入力作業に励み、午後は縄文遺構3Dモデル作成を楽しむことにします。午前中活動では作業効率化のためのPythonスクリプト作成技術向上、午後活動ではBlender技術向上を目指します。

参考

ブログ「花見川流域を歩く」2025年6月記事

ブログ「花見川流域を歩く　自然・風景編」2025年6月記事

晴れ・川霧

2025年6月　Sketchfabに投稿した3Dモデル

2025年6月　YouTubeに投稿した動画

2025年6月　ブログ「花見川流域を歩く」投稿記事に掲載した画像

2025年6月30日月曜日

イランのサルチャハン付近背斜構造地形（ノーマルマップバージョン）

Anticlinal structure near Sarchahan, Iran (normal map version)

After reading the article "Can you drink alcohol in Iran?" in the April issue of Geography magazine, I learned about the alcohol-free social situation in Iran and the taste of non-alcoholic beer.

Although not related to this article, I created a normal map version of the 3D model of the anticlinal structure and enjoyed the mysterious appearance of it being illuminated in the spotlight.

雑誌地理4月号記事「イランで酒は飲めるのか？」を読んで、イラン禁酒社会状況とノンアルビールの味について学習しました。

この記事とは関係しませんが、背斜構造地形3Dモデルのノーマルマップバージョンを作成して、スポットライトが当たるような不思議な様子を楽しみました。

1　イランのサルチャハン付近背斜構造地形（ノーマルマップバージョン）

イランのサルチャハン付近背斜構造地形（ノーマルマップバージョン） by arakiminoru on Sketchfab

イランのサルチャハン付近背斜構造地形（ノーマルマップバージョン）

Anticlinal structure topography near Sarchahan, Iran (Normal Map Version)

Vertical ratio: ×1

Attributions

Generator: DEM Net Elevation API - https://elevationapi.com

Mesh reduction: geometry3Sharp (GradientSpace, Ryan Schmidt) BSL 1.0 License - https://github.com/gradientspace/geometry3Sharp/blob/master/LICENSE

Digital Elevation Model: AW3D30 OpenTopography - https://opentopography.org/

Imagery: Esri World Imagery - https://services.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer

テクスチャはheightmap、normalmap、texture画像の合成

3Dモデルの動画

3Dモデルの画像

2　メモ

雑誌地理4月号記事今村遼平著「イランで酒は飲めるのか？」を読んで、イラン禁酒社会の状況とイランのノンアルコールビールの味がとてもまずいことを学習しました。

テクスチャ画像の作成

イランのサルチャハン付近背斜構造地形（干渉色バージョン）

Anticlinal structure near Sarchahan, Iran (interference color version)

I read the article "Iran, the Earthquake Country" in the April issue of Geography magazine and learned about the current state of earthquakes in Iran and disaster prevention measures.

Although not related to this article, I enjoyed creating an interference color version of a 3D model of an anticlinal structure using an interference color conversion tool (Map Art Institute).

雑誌地理4月号記事「地震国イラン」を読んで、イラン地震事情と防災対策の様子について学習しました。

この記事とは関係しませんが、背斜構造地形3Dモデルの干渉色バージョンを、干渉色変換ツール（地図アート研究所）を利用して作成し、楽しみました。

1　イランのサルチャハン付近背斜構造地形（干渉色バージョン）

イランのサルチャハン付近背斜構造地形（干渉色バージョン） by arakiminoru on Sketchfab

イランのサルチャハン付近背斜構造地形（干渉色バージョン）

Anticlinal structure topography near Sarchahan, Iran (Interference color version)

Vertical ratio: ×1

Attributions

Generator: DEM Net Elevation API - https://elevationapi.com

Mesh reduction: geometry3Sharp (GradientSpace, Ryan Schmidt) BSL 1.0 License - https://github.com/gradientspace/geometry3Sharp/blob/master/LICENSE

Digital Elevation Model: AW3D30 OpenTopography - https://opentopography.org/

Imagery: Esri World Imagery - https://services.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer

テクスチャ画像の干渉色化 - https://ymd5022002.github.io/map-art-jp/tools/index.html

3Dモデルの動画

3Dモデルの画像

2　メモ

雑誌地理4月号尾上篤生著「地震国イラン」を読んで、イラン地震事情と防災対策の様子について学習しました。木のないイランでは土で家をつくることが多く、地震被害は日本とくらべ特段に多いとのことです。

テクスチャ画像の干渉色化は次のステップで行いました。

テクスチャ画像の干渉色化

干渉色変換ツールを無償で提供している地図アート研究所（代表　やまだこーじ　さん）に感謝します。

登録: 投稿 (Atom)

私の散歩論

2025年7月11日金曜日

2025年7月9日水曜日

2025年7月5日土曜日

2025年7月3日木曜日

2025年7月1日火曜日

2025年6月30日月曜日