IDL绘制黑体辐射曲线

普朗克定律是热红外遥感中常常使用的三大定律之一，描述了黑体辐射能量的情况。绝对黑体的辐射光谱对于研究一切物体的辐射规律具有根本的意义。1900年普朗克引进量子概念，将辐射当做不连续的量子发射，成功地从理论上得出了与实验精确符合的绝对黑体辐射出射度随波长的分布函数。试着用IDL绘出了黑体辐射曲线，利于以后计算黑体辐射出射度，以备不时之需！

普朗克公式：

代码：

; :Author: caozhigang

; :Copyright:UCAS

; :blog:blog.sina.com.cn/ahnucao

PRO PLANK_SPECTRUMCURVES

;生成普朗克光谱曲线

colors=['r','g','b','c','m','y','k']

wavelengths=FINDGEN(1000,start = 1,increment = 0.1);1_100 um

radiance = CALC_OMMITRADIANCE(wavelengths,250)

plots = PLOT(wavelengths,radiance,Name = STRING(250)+'K',colors[0],$

AXIS_STYLE = 1,xtitle ='wavelength(um)',ytitle = '辐射出射度/(m^2*um)')

PLOTS.TITLE = '黑体辐射曲线'

lege = legend(target = plots,$

/DATA, /AUTO_TEXT_COLOR)

c_index = 1

FOR T= 270,340,20 DO BEGIN

;计算黑体辐射出射度

radiance = CALC_OMMITRADIANCE(wavelengths,T)

; ;随机产生一个颜色

; colorNames = TAG_NAMES(!color)

; index = ROUND(RANDOMU(seed,1)*N_ELEMENTS(colorNames)+0)

; color = colorNames[index]

plots = PLOT(wavelengths,radiance,Name = STRING(T)+'K',colors[c_index],$

AXIS_STYLE = 1,xtitle ='wavelength(um)',ytitle = '辐射出射度/(m^2*um)',/OVERPLOT)

lege.ADD,plots

c_index = c_index+1

ENDFOR

END

FUNCTION CALC_OMMITRADIANCE,wavelengths,T

;

;wavelengths=FINDGEN(1000,start = 1,increment = 0.1);1_100 um

spectrum= FINDGEN(1000);保存强度

colors = !color

c = 3*10.0^14; m/s

h = 6.6256*10.0^(-34);Plank statics num

k = 1.38*10.0^(-23);S_P statics num

FOR i = 0,N_ELEMENTS(wavelengths)-1 DO BEGIN

;

first_factor = 2*!pi*c*c*h/(wavelengths[i]^5.0)

sencond_factor = 1.0/(EXP(c*h/(k*wavelengths[i]*T))-1)

spectrum[i] = first_factor * sencond_factor*10^12.0;这里乘上10^12次方，是为了将um2转换为转化为m2的单位

ENDFOR

RETURN,spectrum

END

效果图：

转载于:https://www.cnblogs.com/ahnucao/p/4938752.html

IDL绘制黑体辐射曲线相关推荐

R语言构建logistic回归模型并评估模型：模型预测结果抽样、可视化模型分类预测的概率分布情况、使用WVPlots包绘制ROC曲线并计算AUC值
R语言构建logistic回归模型并评估模型:模型预测结果抽样.可视化模型分类预测的概率分布情况.使用WVPlots包绘制ROC曲线并计算AUC值目录
R语言使用pROC包绘制ROC曲线、获取最优阈值（threshold）及最优阈值对应的置信区间
R语言使用pROC包绘制ROC曲线并获取最佳阈值(threshold)及最佳阈值对应的置信区间 #ROC曲线 ROC(receiver operating characteristic curve)接 ...
R语言使用pROC包绘制ROC曲线并使用smooth函数绘制平滑的ROC曲线（方法包括：binormal、density、fitdistr、logcondens、logcondens.smooth）
R语言使用pROC包绘制ROC曲线并使用smooth函数绘制平滑的ROC曲线(方法包括:binormal.density.fitdistr.logcondens.logcondens.smooth) ...
R语言使用pROC包绘制ROC曲线实战：roc函数计算AUC值、plot.roc函数绘制ROC曲线、添加置信区间、为回归模型中的每个因子绘制ROC曲线并在同一个图中显示出来
R语言使用pROC包绘制ROC曲线实战:roc函数计算AUC值.plot.roc函数绘制ROC曲线.添加置信区间.为回归模型中的每个因子绘制ROC曲线并在同一个图中显示出来目录
R使用pROC和ggplot2包绘制ROC曲线
R使用ggplot2包绘制ROC曲线目录 R使用ggplot2包绘制ROC曲线 logistic回归模型构建使用pROC包
使用R构建Xgboost模型并绘制ROC曲线
使用R构建Xgboost模型并绘制ROC曲线 xgboost算法论文全称为<XGBoost: A Scalable Tree Boosting System>,由陈天奇于2016年发表的, ...
验证曲线（ validation curve）是什么？如何绘制验证曲线（ validation curve）？验证曲线（ validation curve）详解及实践
验证曲线( validation curve)是什么?如何绘制验证曲线( validation curve)?验证曲线( validation curve)详解及实践验证曲线( validation ...
通过交叉验证（Cross Validation）KFold绘制ROC曲线并选出最优模型进行模型评估、测试、包含分类指标、校准曲线、混淆矩阵等
通过交叉验证(Cross Validation,CV)KFold绘制ROC曲线并选出最优模型进行模型评估.测试.包含分类指标.校准曲线.混淆矩阵等 Cross Validation cross val ...
校准曲线(calibration curve)是什么？如何绘制校准曲线(calibration curve)？如何通过过校准曲线进行分析？什么是高估？什么是低估？
校准曲线(calibration curve)是什么? 如何绘制校准曲线(calibration curve)? 如何通过过校准曲线进行分析? 什么是高估? 什么是低估? 注意:横纵坐标的差异,但是本 ...

IDL绘制黑体辐射曲线

IDL绘制黑体辐射曲线相关推荐

最新文章

热门文章