自動整理反射率的匯出資料,並計算反射比數值

2021-10-06 14:12:31 字數 2990 閱讀 1209

'''

匯入pandas模組,xlwt模組

將txt檔案中的資料按照;進行分割,寫入list列表中。

'''import pandas as pd

import matplotlib as mp

import matplotlib.pyplot as plt

filename =

"d:\python temp\分光光度計資料\pel_data\export\sample21.樣品.asc"

#將list列表的格式轉換為dateframe格式,直接寫入excel檔案中。

with

open

(filename,

'r+'

, encoding=

'ansi'

)as f:

list =

[i[:-1

].split

('\t'

)for i in f.

readlines()

] #列表生成器,將檔案每行資料按上述方法處理後放入列表

dataframe = pd.

dataframe

(list)

dataframe.

to_excel

(r"d:\python temp\分光光度計資料\pel_data\export\資料預處理.xls"

, index=false)

data = pd.

read_excel

(r"d:\python temp\分光光度計資料\pel_data\export\資料預處理.xls"

, skiprows=90)

data.

to_excel

(r"d:\python temp\分光光度計資料\pel_data\export\資料預處理2.xls"

, index=false)

solar_spectrum = pd.

read_excel

(r"d:\python temp\分光光度計資料\標準太陽光譜分布 jg-t 235-2014:iso 9485-1.xlsx"

, sheet_name=

"sheet1"

)near_infrared_spectrum = pd.

read_excel

(r"d:\python temp\分光光度計資料\標準太陽光譜分布 jg-t 235-2014:iso 9485-1.xlsx"

, sheet_name=

"sheet2"

)standard_whiteboard = pd.

read_excel

(r"d:\python temp\分光光度計資料\標準白板的反射率校準值.xlsx"

, sheet_name=

"sheet1"

)# print

(solar_spectrum)

# print

(standard_whiteboard)

# print

(data)

file1 = pd.

merge

(solar_spectrum, standard_whiteboard, on=

'λ/nm'

)solar_file = pd.

merge

(file1, data, left_on=

'λ/nm'

, right_on=

'#data'

)solar_file.

info()

solar_file.

to_excel

(r"d:\python temp\分光光度計資料\pel_data\export\太陽光反射比.xls"

, sheet_name=

"太陽光反射比"

, index=false)

solar_ratio =

((solar_file[

"ρ0"

]* solar_file[

"unnamed: 1"

]* solar_file[

"sλ*dλ"])

.sum()

)/(solar_file[

"sλ*dλ"].

sum())

/100

print

("太陽光反射比:"

, solar_ratio)

file2 = pd.

merge

(near_infrared_spectrum, standard_whiteboard, on=

'λ/nm'

)near_infrared_file = pd.

merge

(file2, data, left_on=

'λ/nm'

, right_on=

'#data'

)near_infrared_file.

info()

near_infrared_file.

to_excel

(r"d:\python temp\分光光度計資料\pel_data\export\近紅外反射比.xls"

, sheet_name=

"近紅外反射比"

, index=false)

nir=

((near_infrared_file[

"ρ0"

]* near_infrared_file[

"unnamed: 1"

]* near_infrared_file[

"sλ*dλ"])

.sum()

)/( near_infrared_file[

"sλ*dλ"].

sum())

/100

print

("近紅外反射比:"

,nir

)print

('done!'

)

雙向反射率BRDF的定義

自己看書一直沒有這麼明白這個定義,這個部落格解釋的確是很不錯,收藏了。物體表面對電磁波的反射有三種形式 鏡面反射 mirror reflection 反射能量集中在乙個方向,反射角 入射角 漫反射 diffuse reflection 整個表面都均勻地向各向反射入射光稱為漫反射 方向反射 direc...

輻射亮度與表觀反射率的轉換

其他資料 補充相關資料 就是影像灰度值 略大氣層頂的反射率 dn值轉輻射亮度 l gain dn offs etl lambda gain dn offset l ga in d n of fset 式中,l l lambda l 為衛星載荷通道入瞳處等效輻射亮度,單位為w m2 s r m w m...

WINDOW下ORACLE資料自動匯出並壓縮的指令碼

rem 如果檔案存在,那麼退出 rem if exists exp file goto end1 set store path f back yu yyyymmdd ppqq dmp.rar rem 匯出檔案 echo on exp system manager full n grants y ro...