関数read_latlonDEM_to_3d_multi マニュアル

(The documentation of function read_latlonDEM_to_3d_multi)

Last Update: 2025/12/22

◆機能・用途(Purpose)

緯度・経度に基づく格子点で標高が与えられた数値標高モデルを struct _3d_data型構造体に読み込む。緯度経度で範囲を指定する。数値標高モデルは複数ファイルにまたがっていても良い。
Read a digital elevation model, composed of the elevations at latitude- and longitude-based grid nodes, into a struct _3d_data-type structure. The data to be read are specified by latitude and longitude ranges, and can be separated in multiple files.

入力ファイル(数値標高モデルのデータ)は GSIdem2latlonDataコマンドから出力されたファイル、またはそれと同じ書式を持つテキストファイルとする。緯度・経度範囲に応じた複数のファイルに分かれていても良いが以下の3条件を満たさなければならない。
The file formats of input data (digital elevation model) must be identical to the file format of the output files from GSIdem2latlonData command. The data can be divided as multiple files, as long as they satisfy the following requirements.

個々のファイルがカバーする緯度経度範囲は日本の地域メッシュ区分 ( 関数JapanMeshCode2latlonRangeのマニュアル参照) における1次〜3次メッシュのいずれかと一致する。
The spatial range in the latitude and longitude of each file must be identical to that of one of the 1st, 2nd, or 3rd regional meshes of Japan; see the documentation of function JapanMeshCode2latlonRange for the regional meshes.
ファイル名はそのメッシュコードを含んでいる。
Its file name must contain the mesh code.
ファイル名のうちメッシュコード以外の部分は全てのファイルで共通である。
The names of all files must be identical except for the mesh code.

◆形式(Format)

##include <3d_data/readwrite.h>
inline struct _3d_data read_latlonDEM_to_3d_multi
(const char ∗inputfile_name_with_pattern,
const double Nmin,const double Nmax, const double Emin,const double Emax)

◆引数(Arguments)

inputfile_name_with_pattern	入力ファイル名。数値標高モデルのファイル名であり、以下をパターンとして含んでいて良い。 The name of the input files that contain the data of the digital elevation model. This argument may consist of the following patterns. %mesh1 1次メッシュコード(4桁) The 1st mesh code (4 digits). %mesh2 2次メッシュコード(2桁) The 2nd mesh code (2 digits). %mesh3 3次メッシュコード(2桁) The 3rd mesh code (2 digits).
Nmin	数値標高モデルを読み込みたい緯度範囲の下限(度)。 The lower limit latitude (degrees) of the spatial range to read the digital elevation model.
Nmax	数値標高モデルを読み込みたい緯度範囲の上限(度)。 The upper limit latitude (degrees) of the spatial range to read the digital elevation model.
Emin	数値標高モデルを読み込みたい経度範囲の下限(度)。 The lower limit longitude (degrees) of the spatial range to read the digital elevation model.
Emax	数値標高モデルを読み込みたい経度範囲の上限(度)。 The upper limit longitude (degrees) of the spatial range to read the digital elevation model.

◆戻り値(Return value)

引数Nmin, Nmax, Emin, Emaxで指定した範囲をカバーする数値標高モデル。 [Nmin, Nmax], [Emin, Emax]の範囲ではなくそれらを包含するファイルのデータ範囲全体を読み込む。 2次元データ(\(N_z=1\)の3次元データ) であり以下のように与える。
A digital elevation model that covers the spatial range specified by arguments Nmin, Nmax, Emin, and Emax. The output 3-D data is not truncated at Nmin, Nmax, Emin, and Emax but includes the entire data ranges of the corresponding files. This is a 2-D data (3-D data with \(N_z=1\)) and defined as follows.

\(x\)成分: 経度。数値標高モデルの座標が度・分・秒単位で与えられている場合であっても実数への換算に伴う丸め誤差を生じないように度・分部分を含めた秒単位の実数 (すなわち\(3600\times\)度\(+60\times\)分\(+\)秒) で与える。
The \(x\)-component represents the longitude, given as a real number in the second unit that includes the degreen and minute parts in the calculation (\(3600\times\) degree \(+60\times\) minute \(+\) second) to avoid rounding errors in the conversion of the coordinates from degree:minute:second to a real number.
\(y\)成分: 緯度。\(x\)成分と同じ与え方とする。
The \(y\)-component represents the longitude, given in the same manner as the \(x\)-component.
\(z\)成分: 使用しないダミー値であり0とする。標高ではない点に留意。 3次元データの\(z\)成分は\(x\), \(y\)と独立に決まる「独立変数としての座標」であるのに対し、数値標高モデルにおける標高は\((x,y)\)の関数(従属変数)である。
The \(z\)-component is dummy zero value that is not used. Note that the \(z\)-component is not an elevation; the \(z\)-component of a 3-D data is an independent variable for a coordinate, but a function of \((x,y)\) (a dependent variable) that is the case in a digital elevation model.
各位置での3次元データの値: 地表面の標高(m)。
The value of 3-D data at each location: ground surface elevation (m).

◆使用例(Example)

struct _3d_data topography=read_latlonDEM_to_3d_multi ("/dat/DEM_GSI/%mesh1/FG-GML-%mesh1-%mesh2-dem10b-20250101.dat",35.12345,36.78901,137.65432,138.90123);

この例では /dat/DEM_GSI/????/FG-GML-????-??-dem10b-20250101.dat (????は1次メッシュコード、 ??は2次メッシュコード) のパターンにマッチするファイルを北緯35.12345-36.78901度、東経137.65432-138.76543度の範囲で探索して数値標高モデルを読み込む。作成する3次元データにおいて、例えば北緯35度30分36秒の位置は\(y=35\times 3600+30\times 60+36=127836\)、東経137度54分7.2秒の位置は\(x=137\times 3600+54\times 60+7.2=496447.2\) となる。
In this example, the digital elevation model is read from files that match a pattern of /dat/DEM_GSI/????/FG-GML-????-??-dem10b-20250101.dat, where ???? and ?? are the 1st and 2nd mesh codes, respectively, in a spatial range of 35.12345-36.78901N, 137.65432-138.76543E. The coordinates are converted to \(x\) and \(y\) in the 3-D data as, for example, from 35°30’36”N to \(y=35\times 3600+30\times 60+36=127836\) and from 137°54’7.2” to \(x=137\times 3600+54\times 60+7.2=496447.2\).