地質調査所研究資料集、 no. 139

岩石・鉱物の反射スペクトルデータとその表示ソフトウエア (PCSPEC)
Reflectance data of rocks and minerals, and its display software

浦井 稔 (地質調査所)
Minoru URAI (Geological Survey of Japan)

内容紹介

可視~短波長赤外域における岩石・鉱物の反射スペクトルはリモートセンシングによる資源探査の基礎資料として不可欠である。地質調査所では、野外電磁特性の研究として岩石・鉱物の反射スペクトル測定ならびに、X線回折法による鉱物の同定や化学分析を行ない、これをデータベースに蓄積してきた。その成果の一つとして反射スペクトル・カタログを出版した(浦井ほか、1989)。その後、この反射スペクトル・カタログのデータをデジタル・データで使用したいとの要望が強まったため、反射スペクトル・データをその表示用ソフトウエアPCSPECとともに、フロッピーディスクとして公開することとした。これらのデータおよびソフトウエアはNEC製のパーソナルコンピュータで使用することを前提として作成されたが、条件さえ整えば、他のパーソナルコンピュータでも使用可能である。

頁数:16

受理日:1989年12月20日

目 次

1 はじめに
2 表示ソフトウエア PCSPEC
2.1 ハードウエア・ソフトウエア環境
2.2 PCSPEC使用の手引き
2.3 画面のハードコピー
2.4 他のパーソナル・コンピュータへの移植
3 岩石・鉱物の反射スペクトルデータの利用法
3.1 試料の記述、化学分析値およびX線回折結果の利用
3.2 反射スペクトル・データの利用
Appendix 1 各ファイルのフォーマット
Appendix 2 試料の記述、化学分析値およびX線回折結果の利用例
Appendix 3 反射スペクトル・データの利用例

1 はじめに

可視~短波長赤外域における岩石・鉱物の反射スペクトルはリモートセンシングによる資源探査の基 礎資料として不可欠である。地質調査所では、野外電磁特性の研究として岩石・鉱物の反射スペクトル測 定ならびに、X線回折法による鉱物の同定や化学分析を行ない、これをデータベースに蓄積してきた。そ の成果の一つとして反射スペクトル・カタログを出版した(浦井ほか、1989)。その後、この反射スペクトル・カタログのデータをデジタル・データで使用したいとの要望が強 まったため、反射スペクトル・データをその表示用ソフトウエアPCSPECとともに、フロッピーディスクと して公開することとした。
この反射スペクトル・データには反射スペクトル・カタログとして出版したデータのうち、X 線チャートを除く全てのデータが含まれている。すなわち、約80個の粘土鉱物と約30個の炭酸塩岩試料 における試料の記述、0.5um-2.5umにおける反射スペクトル、化学分析値およびX線回折結果である。反射 スペクトル・データの利用法およびそのプログラム例を示した。データについての詳細は反射スペクト ル・カタログ(浦井ほか、1989)を参照されたい。
表示用ソフトウエアPCSPECは、反射スペクトル・カタログのデータのうち、X線チャートを除く全ての データを表示することができる。
これらのデータおよびソフトウエアはNEC製のパーソナルコンピュータで使用することを前提として作成されたが、条件さえ整えば、他のパーソナルコンピュータでも使用可能である。

2 表示ソフトウエア PCSPEC

表示用ソフトウエアPCSPECは、反射スペクトル・カタログのデータのうち、X線チャートを除く全てのデータを表示することができる。

2.1 ハードウエア・ソフトウエア環境

PCSPECは以下のハードウエア・ソフトウエア環境で作成された。従って、PCSPECは以下のハードウエア ・ソフトウエア環境があれば十分動作するが、これが満たされない場合でも動作する可能性はある。

ハードウエア環境
CPU PC-9801VM2
ディスプレー PC-KD852
プリンタ PC-PR201HC
メモリ 640KB
ソフトウエア環境
OS MS-DOS v.3.1

2.2 PCSPEC使用の手引き

反射スペクトル・データおよびPCSPECは2枚の5インチ2HDのフロピーディスクで提供される。まず、PC -9801VM2をMS-DOS v.3.1で立ちあげる。次に、コントロールディスクをドライブB: に挿入して、
A> B:PCSPEC
と入力する。すると、Fig.1の画面が表示される。しばらくすると、
Remove Control Disk and insert Data Disk to drive B:
と出るので、データディスクをドライブB:に挿入して、ENTERキーを押す。すると、Fig.2の画面が表示さ れる。目的の試料に対応する番号を入力する。次のページを表示したい場合はENTERキーを、プログラム を終了したい場合は-1をそれぞれ入力する。目的の試料に対応する番号を入力した場合は、Fig.3に示す 試料の記述、化学分析値およびX線回折結果が表示される。ここで、 ENTERキーを押せば、Fig.4に示す残りの試料の記述、化学分析値およびX線回折結果が表示される。次に、 ENTERキーを押せば、Fig.5に示す0.5um-2.5umにおける反射スペクトルが表示される。 さらに、ENTERキーを押せば、試料番号を選択するメニュー(Fig.2)に戻る。
2.3 画面のハードコピー

文字画面のハードコピーは、COPYキーを使用してとることができる。グラフィック画面のハードコピ ーは、Fig.5に示す画面でPと入力することによってとれる。
2.4 他のパーソナル・コンピュータへの移植

NEC PC9801以外のパーソナル・コンピュータで、反射スペクトル・データおよびPCSPECを使用するため には、データの変換およびプログラムの変換を行なう必要がある。PCSPECは、マイクロソフト社のクイ ックベーシックVersion 4.2の文法で記述されている。当該パーソナル・コンピュータがマイクロソフト社のクイックベーシックをもつ場合、プログラム変換作業は比較的簡単であろう。

1) データの変換

当該パーソナル・コンピュータのOSがMS-DOSでない場合は、次に示すファイルのデータ変換を行なわ なければならない。
コントロールディスク
PCSPEC.BAS PCSPECのソースコード
PCSPEC.CNT 試料の記述、化学分析値およびX線回折結果のデータ
データディスク
*.WVL スペクトルデータの波長部
*.REF スペクトルデータの反射強度部
これらは、市販のユーティリティ・プログラムで変換できる。

2)プログラムの変換

PCSPEC.BASはマイクロソフト社のクイックベーシックで記述されているため、当該パーソナル・コン ピュータでマイクロソフト社のクイックベーシックを使用して、PCSPEC.BASをコンパイルしなければな らない。当該パーソナル・コンピュータがマイクロソフト社のクイックベーシックを持たない場合、他の ベーシック言語を使用することも可能であるが、グラピィック部分のプログラム変換は困難であろう。

3 岩石・鉱物の反射スペクトルデータの利用法

反射スペクトル・データおよびPCSPECは2まいの5インチ2HDのフロピーディスクで提供される。これには以下のファイルが含まれる。


コントロールディスク
PCSPEC.EXE PCSPECの実行形式プログラム
PCSPEC.BAS PCSPECのソースコード
PCSPEC.CNT 試料の記述、化学分析値およびX線回折結果のデータ


データディスク
*.WVL スペクトルデータの波長部
*.REF スペクトルデータの反射強度部


PCSPEC.CNTは各試料のデータが1レコードとして全各試料のデータが格納されている(正確には、各試料 のデータが1レコードとしてX線回折で得られた鉱物の個数存在する)。一方、*.WVLや *.REFは各試料のデータが1ファイルとして格納されている。各ファイルのフォーマットを Appendix 1に示す。

3.1 試料の記述、化学分析値およびX線回折結果の利用

試料の記述、化学分析値およびX線回折結果を読むためには、Appendix 1およびPCSPEC.BASの RDATRサブルーチンを使用する。具体的には、Appendix 2を参照のこと。

3.2 反射スペクトル・データの利用

反射スペクトル・データを読むためには、Appendix 1およびPCSPEC.BASのGETSPTサブルーチンを使用する。具体的には、Appendix 3を参照のこと。

参考文献

浦井 稔・佐藤 功・二宮芳樹・古宇田亮一・宮崎芳徳・山口 靖(1989) 可視から短波長赤外域における 岩石・鉱物の反射スペクトル・カタログ,地質調査所, pp. 367.

Fig.1 PCSPEC initial menu

Fig.2 Sample selection menu

Fig.3 Sample discription, chemical analysis and X-ray diffraction results

Fig.4 Remeined sample discription, chemical analysis and X-ray diffraction results

Fig.5 Example of spectrum display
Appendix 1 File format
PCSPEC.CNT
Feild no.  Item               Name*
Sample record
1       Sample ID             sid$
2       Sample name           name$
3       Color                 color$
4       Sampling date         sdate$
5       Nation                nation$
6       Prefecture            pref$
7       Locality              site$
8       Latitude              latd%
9       Latitude              latm
10      Longitude             lond%
11      Longitude             lonm
12      Sample shape          shape$
13      Alteration            alt$
Lab. measurement record
14      Measurement ID        lid$
15      Sample ID             lsid$
16      Date                  ldate$
17      Equipment             lequ$
18      Reference             lref$
19      Parameter             lprm$
20      Treatment             ltre$
21      Lamp                  lamp$
Chemical composition record
22      Measurement ID        cid$
23      Sample ID             csid$
24      Date                  cdate$
25      Equipment             cequ$
26      Parameter             cprm$
27      SiO2                  sio2
27      TiO2                  tio2
28      Al2O3                 Al2O3
29      Fe2O3                 Fe2O3
31      FeO                   feo
32      MnO                   mno
33      MgO                   mgo
34      CaO                   cao
35      Na2O                  na2o
36      K2O                   K2O
37      P2O5                  P2O5
38      S                     s
39      C                     c
40      ZrO2                  ZrO2
41      H2O+                  h2op
42      H2O-                  h2om
43      CO2                   co2
44      SO3                   so3
45      LOI                   loi
46      BaO                   bao
Mineral composition record
47      Sample ID             nsid$
48      Chemical Rec.         ncid$
49      Quartz                q
50      Corumdum              co
51      Andalusite            ad
52      Microcline            mi
53      Albite                ab
54      Anorthite             an
55      Kaolinite             ka
56      Sericite              se
57      Chlorite              ch
58      Gibbsite              gi
59      Montmorillonite       mo
60      Pyrophyllite          pp
61      Diaspore              dp
62      Alunite               al
63      Zircon                z
64      Halite                hl
65      Thernardite           th
66      Calcite               cc
67      Magnesite             mg
68      Siderite              sd
69      Limonite              li
70      Magnetite             mt
71      Hematite              he
72      Chromite              cm
73      Pyrite                pr
74      Ilmenite              il
75      Rutile                ru
76      Apatite               ap
77      Fluorite              fr
78      Enstatite             en
79      Others                ot
80      Condition             cond$
81      Comment               com$
Xray measurement record
82      Sample ID             xsid$
83      Measurement ID        xid$
84      Date                  xdate$
85      Equipment             equ$
86      Parameter             xprm$
87      Treatment             xtre$
Detected minerals
88      Measurement ID        xxid$
89      Detected minreal      xmin$
*: Names are corresponed to the variable names of the BASIC program.
Appendix 2 Example program to read control data
COMMON SHARED cnt$, dath$, datt$, viewx, viewy, windx, windy, windx1, windx2
no% = 120
mno% = 10
cnt$ = "b:pcspec.cnt"
dath$ = "b:"
datt$ = ".ref"
viewx = 639
viewy = 399
windx1 = 3
windy1 = 1
windx = 30
windy = 20
REM $DYNAMIC
DIM SHARED sid(no%) AS STRING * 20, name$(no%), color$(no%), sdate(no%) AS STRING * 10, nation$(no%)
DIM SHARED pref$(no%), site$(no%), latd%(no%), latm(no%), lond%(no%)
DIM SHARED lonm(no%), shape$(no%), alt$(no%)
DIM SHARED lid$(no%), lsid(no%) AS STRING * 12, ldate(no%) AS STRING * 10, lequ(no%) AS STRING * 12, lref$(no%)
DIM SHARED lprm$(no%), ltre$(no%), lamp$(no%)
DIM SHARED cid$(no%), csid(no%) AS STRING * 12, cdate(no%) AS STRING * 10, cequ(no%) AS STRING * 12, cprm$(no%)
DIM SHARED sio2(no%) AS SINGLE, tio2(no%), Al2O3(no%), Fe2O3(no%), feo(no%)
DIM SHARED mno(no%), mgo(no%), cao(no%), na2o(no%), K2O(no%), P2O5(no%)
DIM SHARED s(no%), c(no%), ZrO2(no%), h2op(no%), h2om(no%), co2(no%)
DIM SHARED so3(no%), loi(no%), bao(no%)
DIM SHARED nsid(no%) AS STRING * 12, ncid$(no%), q(no%), co(no%), ad(no%), mi(no%)
DIM SHARED ab(no%), an(no%), ka(no%), se(no%), ch(no%), gi(no%)
DIM SHARED mo(no%), pp(no%), dp(no%), al(no%), z(no%), hl(no%), th(no%)
DIM SHARED cc(no%), mg(no%), sd(no%), li(no%), mt(no%), he(no%)
DIM SHARED cm(no%), pr(no%), il(no%), ru(no%), ap(no%), fr(no%)
DIM SHARED en(no%), ot(no%), cond$(no%), com$(no%)
DIM SHARED xid$(no%), xsid(no%) AS STRING * 12, xdate(no%) AS STRING * 10, xequ(no%) AS STRING * 12, xprm$(no%)
DIM SHARED xtre$(no%)
DIM SHARED xxid$(no%), xmin(no%, mno%) AS STRING * 20, xmno%(no%)
DIM SHARED wvl(1151), igrt%(9)
DIM x(1151), y(1151), h(800), quo(800), idep%(400)
REM $STATIC
OPEN cnt$ FOR INPUT AS #1
i% = 1
j% = 1
CALL rdatr(i%)
xmno%(i%) = 1
CLS
CALL title
DO WHILE NOT EOF(1) AND i% < no%
LOCATE 23, 1
PRINT " Reading control file ---"; i%
i% = i% + 1
CALL rdatr(i%)
IF sid$(i%) = sid$(i% - 1) THEN
xmno%(i% - 1) = xmno%(i% - 1) + 1
xmin$(i% - 1, xmno%(i% - 1)) = xmin$(i%, 1)
i% = i% - 1
ELSE
xmno%(i%) = 1
END IF
LOOP
CLOSE #1
rno% = i%

SUB rdatr (i%)
INPUT #1, sid$(i%), name$(i%), color$(i%), sdate$(i%), nation$(i%)
INPUT #1, pref$(i%), site$(i%), latd%(i%), latm(i%), lond%(i%)
INPUT #1, lonm(i%), shape$(i%), alt$(i%)
INPUT #1, lid$(i%), lsid$(i%), ldate$(i%), lequ$(i%), lref$(i%)
INPUT #1, lprm$(i%), ltre$(i%), lamp$(i%)
INPUT #1, csid$(i%), cid$(i%), cdate$(i%), cequ$(i%), cprm$(i%)
INPUT #1, sio2(i%), tio2(i%), Al2O3(i%), Fe2O3(i%), feo(i%)
INPUT #1, mno(i%), mgo(i%), cao(i%), na2o(i%), K2O(i%), P2O5(i%)
INPUT #1, s(i%), c(i%), ZrO2(i%), h2op(i%), h2om(i%), co2(i%)
INPUT #1, so3(i%), loi(i%), bao(i%)
INPUT #1, nsid$(i%), ncid$(i%), q(i%), co(i%), ad(i%), mi(i%)
INPUT #1, ab(i%), an(i%), ka(i%), se(i%), ch(i%), gi(i%)
INPUT #1, mo(i%), pp(i%), dp(i%), al(i%), z(i%), hl(i%), th(i%)
INPUT #1, cc(i%), mg(i%), sd(i%), li(i%), mt(i%), he(i%)
INPUT #1, cm(i%), pr(i%), il(i%), ru(i%), ap(i%), fr(i%)
INPUT #1, en(i%), ot(i%), cond$(i%), com$(i%)
INPUT #1, xsid$(i%), xid$(i%), xdate$(i%), xequ$(i%), xprm$(i%)
INPUT #1, xtre$(i%)
INPUT #1, xxid$(i%), xmin$(i%, 1)
END SUB
Appendix 3 Example program to read reflectance data
SUB getspt (sname$, x(), y(), n%, wmin, wmax, sno$, comm$)
STATIC wname$
OPEN dath$ + sname$ + datt$ FOR INPUT AS #1
LINE INPUT #1, a$
INPUT #1, b$
FOR i% = 1 TO 1151: INPUT #1, y(i%): NEXT i%
LINE INPUT #1, a$
sno$ = MID$(a$, 1, 18)
comm$ = MID$(a$, 19, 60)
CLOSE #1
IF wname$ <> b$ THEN
wname$ = b$
OPEN dath$ + wname$ FOR INPUT AS #1
LINE INPUT #1, a$
FOR i% = 1 TO 1151: INPUT #1, wvl(i%): wvl(i%) = wvl(i%) / 10000: NEXT i%
FOR i% = 1 TO 8: INPUT #1, igrt%(i%): NEXT i%
CLOSE #1
igrt%(9) = 800
END IF
n1s = 1
n1e = igrt%(2)
n2s = igrt%(4)
n2e = igrt%(6)
n3s = igrt%(8)
n3e = igrt%(9)
j% = 1
FOR i% = n1s TO n1e
IF wvl(i%) > wmin AND wvl(i%) < wmax THEN
x(j%) = wvl(i%)
y(j%) = y(i%)
j% = j% + 1
END IF
NEXT i%
FOR i% = n2s TO n2e
IF wvl(i%) > wmin AND wvl(i%) < wmax THEN
x(j%) = wvl(i%)
y(j%) = y(i%)
j% = j% + 1
END IF
NEXT i%
FOR i% = n3s TO n3e
IF wvl(i%) > wmin AND wvl(i%) < wmax THEN
x(j%) = wvl(i%)
y(j%) = y(i%)
j% = j% + 1
END IF
NEXT i%
n% = j% - 1
FOR i% = 1 TO n%
IF y(i%) > 150 THEN y(i%) = 150
IF y(i%) < 0 THEN y(i%) = 0
NEXT i%
END SUB
最新情報

この研究資料集のフロッピーディスクに収録したファイルは、以下の場所よりダウンロードして使用することができます。

自己解凍形式 LZH圧縮形式 ZIP圧縮形式

文およびファイルの引用例

浦井 稔(1989) 岩石・鉱物の反射スペクトルデータとその表示ソフトウェア(PCSPEC),地質調査所研究資料集, no. 139,地質調査所.