関数read_latlonDEM_to_3d マニュアル

(The documentation of function read_latlonDEM_to_3d)

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 is read from a single file. The user specifies the file name.

◆形式(Format)

##include <3d_data/readwrite.h>
inline struct _3d_data read_latlonDEM_to_3d(const char ∗inputfile)

◆引数(Arguments)

inputfile

入力ファイル名。 GSIdem2latlonDataコマンドから出力されたファイル、またはそれと同じ書式を持つテキストファイルとする。
The name of an input file, which must be a text file that has the same format as the output files from GSIdem2latlonData command.

◆戻り値(Return value)

引数inputfileで指定したファイルから読み込んだ数値標高モデル。 2次元データ(\(N_z=1\)の3次元データ) であり以下のように与える。
A digital elevation model read from the file specified by argument inputfile. 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)。なお標高が未定義値(undefined)の場合には-99999999.99を代入する。
The value of 3-D data at each location: ground surface elevation (m). If the elevation is undefined, -99999999.99 is inserted.

◆使用例(Example)

struct _3d_data topography=read_latlonDEM_to_3d ("/dat/DEM_GSI/5337/FG-GML-5337-27-dem10b-20250101.dat");

作成する3次元データ(変数topography)において、例えば北緯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\) となる。
The coordinates are converted to \(x\) and \(y\) in the 3-D data (variable topography) 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\).