1

pykrige を使用して 3D で土壌ボリュームをクリグしようとしています。通常のクリギングを使用していますが、データのクリギングが完了すると、それを ASCII ファイルに書き込む方法がわかりません。コードは次のとおりです。

from pykrige.ok3d import OrdinaryKriging3D
from pykrige.uk3d import UniversalKriging3D
import numpy as np
from matplotlib import pyplot as plt
import pandas as pd
import pykrige.kriging_tools as kt

data_in = pd.read_excel("Krig3D_data.xlsx")
data = data_in.to_numpy()

gridx = np.arange(392050, 392250, 5)
gridy = np.arange(5138395, 5138595,5)
gridz = np.arange(0.0, 5, 0.25)
                     
ok3d = OrdinaryKriging3D(
    data[:, 0], data[:, 1], data[:, 2], data[:, 3], variogram_model="linear"
)

k3d1, ss3d = ok3d.execute("grid", gridx, gridy, gridz)

クリグ データを ascii ファイルにエクスポートして、別の場所でスライスして視覚化できるようにする必要があります。

ありがとう！

4

1 に答える 1