实验三 反演方法

实验三 反演方法

一、实验目的

利用统计经验方法，建立海洋光学的反演方法。

二、实验内容

1、利用70%的数据拟合所选定的函数，求出函数中的统计常数；

2、将70%数据代入得到的反演算法，计算反演参数，并与真值进行比较；

3、将剩余30%数据代入得到的反演算法，计算反演参数，并与真值进行比较。

4、分析结果

三、实验原理及方法

经验统计算法是基于一个实验数据集而建立的光学测量数据与水中组分浓度之间的关系。这个数据集中包含了现场测量的光谱反射率（辐亮度）以及同步获取的组分浓度的数据。在海洋光学反演中，最常用的关系是基于波段比值的方程式

其中，P 是待求的物理量，如叶绿素浓度、悬浮物浓度、可溶有机物及衰减系数等。Ri 表示第i 个通道的反射率或辐亮度，a0、a1和a2为通过回归都得到常数。

21102a R P a a R ??=+ ???

或者

Kd(490)=0.016

四、实验结果及分析

1.图1是利用70%的数据拟合所选定的函数，黑色点代表原始数据，红色曲线代表拟合曲线，用来拟合的函数为 Y=a1*x^2+b1*x+c1

求出a1=0.2979; b1=-1.3899; c1=0.1546;

计算出的反演参数与真实值的相关系数接近0.9。

)490(10)490(33,22,1,0,w x a x a x a a d K K j j j j +=+++23

100123log ((490)(490))d w K K a a x a x a x -=+++12()10()log ()rs rs R R x λλ=12()(490)(490)()(555)(665)

rs rs rs rs rs rs R R R R R R λλ=

或

图1

2.图二是将剩余30%数据代入得到的反演算法，计算反演参数，并与真值进行比较。蓝色点代表原始数据，红色曲线代表拟合曲线。计算出的反演参数与真实值间的相关误差大于0.9，均方差约为0.3，说明拟合曲线比较符合实际情况。

图2

附源程序代码：

clear all;

close all;

data=dlmread('team_1(coast)kd.txt','',1,0); Rrs443=data(:,1);

Rrs490=data(:,2);

Rrs555=data(:,3);

Rrs665=data(:,4);

kd490=data(:,5);

ratio1=Rrs490./Rrs555;

ratio2=Rrs443./Rrs555;

ratio3=Rrs555./Rrs665;

ratio4=Rrs490./Rrs665;

ratio5=Rrs443./Rrs490;

logratio1=log10(ratio1);

logratio2=log10(ratio2);

logratio3=log10(ratio3);

logratio4=log10(ratio4);

logratio5=log10(ratio5);

logkd490=log10(kd490);

n1=corrcoef(logratio1,logkd490);

n2=corrcoef(logratio2,logkd490);

n3=corrcoef(logratio3,logkd490);

n4=corrcoef(logratio4,logkd490);

n5=corrcoef(logratio5,logkd490);

figure(1);

loglog(ratio4(1:200),kd490(1:200),'dk');

xlim([0.2 10]);

ylim([0.01 100]);

xlabel('Rrs490/Rrs665');

ylabel('kd(490)');

newtype=fittype('a1*x^2+b1*x+c1');

m=fit(logratio4(1:200),logkd490(1:200),newtype);

a1=m.a1;

b1=m.b1;

c1=m.c1;

kd490t=10.^(a1*(logratio4(1:200)).^2+b1*logratio4(1:200)+c1);

hold on;

loglog(ratio4(1:200),kd490t(1:200),'r');

n6=corrcoef(kd490t,kd490(1:200));

hold off;

figure(2);

kd490tt=10.^(a1*(logratio4(201:279)).^2+b1*logratio4(201:279)+c1) ;

loglog(kd490(201:279),kd490tt,'*b');

xlabel('kd490-measured');

ylabel('kd490-derived');

hold on;

loglog(kd490(201:279),kd490(201:279),'r');

n7=corrcoef(kd490tt,kd490(201:279));

N=length(kd490(201:279));

e=((kd490tt-kd490(201:279))./kd490(201:279)).^2;

rmse=sqrt(sum(e)/(N-1));

数据处理的基本方法

第六节数据处理的基本方法 前面我们已经讨论了测量与误差的基本概念，测量结果的最佳值、误差和不确定度的计算。然而，我们进行实验的最终目的是为了通过数据的获得和处理，从中揭示出有关物理量的关系，或找出事物的内在规律性，或验证某种理论的正确性，或为以后的实验准备依据。因而，需要对所获得的数据进行正确的处理，数据处理贯穿于从获得原始数据到得出结论的整个实验过程。包括数据记录、整理、计算、作图、分析等方面涉及数据运算的处理方法。常用的数据处理方法有：列表法、图示法、图解法、逐差法和最小二乘线性拟合法等，下面分别予以简单讨论。 列表法是将实验所获得的数据用表格的形式进行排列的数据处理方法。列表法的作用有两种：一是记录实验数据，二是能显示出物理量间的对应关系。其优点是，能对大量的杂乱无章的数据进行归纳整理，使之既有条不紊，又简明醒目；既有助于表现物理量之间的关系，又便于及时地检查和发现实验数据是否合理，减少或避免测量错误；同时，也为作图法等处理数据奠定了基础。 用列表的方法记录和处理数据是一种良好的科学工作习惯，要设 计出一个栏目清楚、行列分明的表格，也需要在实验中不断训练，逐步掌握、熟练，并形成习惯。 一般来讲，在用列表法处理数据时，应遵从如下原则：

(1) 栏目条理清楚，简单明了，便于显示有关物理量的关系。 (2) 在栏目中，应给出有关物理量的符号，并标明单位(一般不重复写在每个数据的后面)。 (3) 填入表中的数字应是有效数字。 (4) 必要时需要加以注释说明。 例如，用螺旋测微计测量钢球直径的实验数据列表处理如下。 用螺旋测微计测量钢球直径的数据记录表 从表中，可计算出 D i D = n = 5.9967 ( mm)

测试技术与实验方法 柳昌庆 第十三章课后参考答案

13-4有一应变式测力传感器，弹性元件为实心圆柱，直径D=40mm 。在圆柱轴向 和周向各贴两片应变片（灵敏度系数s=2.），组成差动全桥电路，供桥电压为10v 。 设材料弹性模量E=2.1?1011 pa ，泊松比υ=0.3。试求测力传感器的灵敏度（该灵敏度用μv/kN 表示）。 解：设受压缩 F,轴向贴的应变片 R S R R ξ=?=?31 横向贴的应变片：SR R R μξ=?=?43 设原电阻 4321R R R R === ，则受力F 后： 11R R R ?-= , 33R R R ?-=, 22R R R ?-= , 44R R R ?-= 电桥输出电压变化： x U R R R R R R R R U ) )((43214 231+-= ?

x x x U R R R U R R R R R U R R R R R R 22)(2 1212 12 212 1?+?- ≈?+?-?+?- =?++?-?+?- = x U S U U 2 )1(ξ +- =?∴ E r F E r F A F 22πσξπσ?=?=?=?= ?∴ 代入上式 x U E r F S U U 22)1(π?+-=?∴ 测力传感器灵敏度 E r SU U F U K x 22)1(π+=??= ) (101.2)(02.0) (1022)3.01(11 22Pa m V ?????+= π 又因为： 218.91002.11m N Pa ??=- 所以： ) (8.9102.1101.202.0) (103.11112N V K ??????? =-π N V 10102.32103.1-???= N V 91089.41-?= KN V μ4189.0= 13-5在一受拉弯综合作用的构件上贴有四个电阻应变片。试分析各应变片感受的应变，将其值填写在应变表中。并分析如何组桥才能进行下述测试：(1) 只测弯矩，消除拉应力的影响；(2) 只测拉力，消除弯矩的影响。电桥输出各为多少？ 组桥如图。

数值分析实验三

贵州师范大学数学与计算机科学学院学生实验报告 课程名称： 数值分析 班级： 12级信本 实验日期： 2014年 11月5日 学 号： 120703010031 姓名： 李盼 指导教师： 实验成绩： 一、实验名称 实验三:数值积分 二、实验目的及要求 1. 让学生掌握复化梯形法, 复化Simpson 法和Romberg 公式以及变步长梯形法,变步长Simpson 法 2. 让学生能够用这些方法解决一些具体问题 三、实验环境 每人一台计算机，要求安装Windows XP 操作系统，Microsoft office2003、MATLAB6.5(或7.0). 四、实验内容 1 . 要求:分别用复合梯形法,复合Simpson 法和 Romberg 公式计算. 2．给定积分dx x ?3 11, 分别用下列方法计算积分值要求准确到510-, 并比较分析计算时间. 1) 变步长梯形法; 2) 变步长 Simpson 法. 五、算法描述及实验步骤 1复合梯形法 输入：被积函数数据点t,a ； 输出：积分值n T 。 复合Simpson 输入：被积函数)(x f ，积分区间[a,b]和n, 输出：复合Simpson 积分值n S 。 步1 :b a h n -?; (a)(b)n S f f ?-;x a ?. 步2：对k=1,2,…,n 执行 ;4();;2().22 n n n n h h x x S S f x x x S S f x ?+?-?+?-

步3：*6 n n h S S ?。 步4：输出n S 。 Romberg 公式 根据已知数据对其进行多项式拟合得出p(x);()()f x p x ?。 输入 ：被积函数)(x f ，积分区间端点a,b,允许误差ε， 输出：Romberg 积分值2n R 。 步1：a b h -?;));()((2 1b f a f h T +?.0;0;0;0111?===k R C S 步2：反复执行步3-步9 步3：.2 ;0h a x S +?? 步4：反复执行步5-步6 步5：.);(h x x x f S S +?+? 步6：若b x ≥,则退出本层循环。 步7：执行 ;3134;2212212T T S S h T T -?+?.63 16364;1511516122122C C R S S C -?-? 步8：执行 .1; ;; ;;2|;|2121212112+?????? -?k k R R C C S S T T h h R R e 步9：若ε≤e 且5≥k ，则退出循环。 步10：.22R R n ? 2：变步长梯形法 功能 求积分dx x ?3 1 1，允许误差为ε。 输入 被积函数)(x f ，a,b,ε. 步1：a b h -?。 步2：)).()((2 1b f a f h T +? 步3：反复执行步4--步10. 步4：.2;0h a x S +=?

地表温度反演实验报告

遥感原理与及应用 地表温度反演实验报告 专业：地理信息系统 班级： XXXXXXXX 姓名： XXX 学号： XXXXXX 成绩： 指导教师： XXX 2014年12月17日 一. 实验目的 1. 根据实际需要，学会在网上（如中国科学院遥

感与数字地球研究所数据共享网）下载研究区内的遥感数据； 2. 掌握在ENVI中实现简单的地表温度反演的原理与步骤。 二. 实验任务 1. 在中国科学院遥感与数字地球研究所数据共享网上订购并下载覆盖郫县的TM影像； 2. 在ENVI中实现简单的地表温度反演算法。 三. 实验数据 在中国科学院遥感与数字地球研究所数据共享网上下载的覆盖郫县地区的TM影像。

四. 实验原理 图1 TM 影像地表温度反演流程 1. 地表温度(Land Surface Temperature)反演公 式为: 2 1(1)K LST K In R ε=+， 其中，R m DN d =?+，2111607.76K W m sr m μ---=???，21260.56K K =。 2. 根据TM 辐射定标原理，热红外波段表观辐亮 度可以进一步写作: max min 6min 255L L R DN L -=?+， 其中LmaxBand6=15.303 , LminBand6=1.238。 3. 地表比辐射率ε为同温度下地表辐射能与黑体 辐射能的比率，其可以表示为: 1.0090.047(In )(0)NDVI NDVI ε=+>，

其中，4343 TM TM NDVI TM TM -=+，当0NDVI <=时(如水体)地表比辐射率取常数1。 五. 实验步骤 1. TM 数据下载 数 据查询和下载网址https://www.wendangku.net/doc/c114277375.html,/query.html ，界面如图2 所示。 图2 中国科学院遥感与数字地球研究所数据共享 网址界面

实验技术人员管理办法武汉理工大学教务处欢迎您

武汉理工大学实验技术人员管理办法 （经2005年第5次校长办公会定定通过） 为了加强对实验技术队伍的建设和管理，充分发挥实验技术人员的积极性、主动性和创造性，不断提高实验教学质量和实验技术人员的业务水平及科学研究水平，实现实验技术人员管理的规范化、科学化，特制定本办法。 第一条我校实验技术人员包括：在实验室工作的高级实验师（高级工程师），实验师（工程师），助理实验师（助理工程师），实验员（技术员），实验技师和技术工人，以及安排在实验室工作一年以上的教师。 第二条实验技术人员的职责： 1. 高级实验师（高级工程师） ①参与实验室建设规划的制定，并负责其中部分项目的执行。 ②每学年至少承担二门本实验室主要课程的实验教学工作。 ③每学年指导一位初级实验技术人员。 ④指导中级、初级实验技术人员对本实验室的大型、精密、贵重仪器设备进行维修保养。 ⑤主持本实验室新购大型仪器设备的验收以及大型仪器设备的功能开发。 2. 实验师（工程师） ①参与实验室建设规划的制定，并参与其中部分项目的执行。 ②每学年至少承担一门本实验室主要课程的实验教学工作。 ③负责对本实验室的大型、精密、贵重仪器设备的维修保养。 ④参与本实验室新购大型仪器设备的验收以及大型仪器设备的功能开发。 3. 助理实验师（助理工程师） ①参与实验室建设规划的制定，并协助参与其中部分项目的执行。 ②每学年至少承担一门本实验室主要课程的实验教学辅导工作。 ③负责对本实验室的仪器设备的维修保养。 ④协助实验师（工程师）的工作。 ⑤每学年参加一到两门本专业课程的学习，取得合格以上成绩。 4. 实验员（技术员） 协助实验师（工程师）完成实验教学工作，承担实验室实验教学的准备、实验室仪器设备的维护保养。 5. 实验技师和技术工人 ①实验室仪器设备的保养、维修及日常管理。 ②实验室环境的维护。 ③实验教学所需仪器设备的准备。 第三条各学院（部）每学年要检查聘约和合同执行情况，并结合工作岗位职责进行考核（具体考核实施办法见附件1）。考核情况应填写专门表格存入本人的技术档案，作为晋升职务的依据之一。考核连续两年不合格者，解除聘约或合同，被解除聘约或合同人员按校规定处理。 第四条实验技术人员实行坐班制，要按时上下班。要严格做好实验技术人员的考勤记录，每学期对实验技术人员进行一次日常工作考核（具体考核内容见附件2），报教务处存档。对不认真填写甚至弄虚作假的，给予批评教育。 第五条学校每两年进行一次评选先进实验室和先进实验技术人员的活动。以奖励在实验室工作中做出显着成绩的集体和个人。各院（部）根据实验技术人员所完成的工作量和业务考核成绩给予奖励，对工作认真、完成任务好、业务考核成绩优秀的实验技术人员给予精神及物质奖励。 第六条对工作不负责任、不遵守劳动纪律、出勤不出力者，对实验教学仪器设备特别是大型、精密、贵重仪器设备的使用、维修、保养不当，造成重大损失者，视情节轻重，给予批评教育、经济处罚、行政纪律处分。

基于ETM数据地表温度反演实验单通道算法操作文档

基于ETM 数据煤田火区地表温度反演的研究实验操作步骤与流程 算法：单通道算法，其公式为 Τs =γ ε?1 ψ1L sensor +ψ2 +ψ3 +δ（1） γ= c 2L sensor T sensor 2 λ4 c 1 L sensor +λ?1 ?1 （2） δ=?γL sensor +T sensor （3） L sensor =L min λ + L max λ ?L min λ Q DN Q max （4） T sensor = K 2 ln 1+K 1L λ （5） K 1=666.09 mW ?cm ?2?sr ?1?um ?1 , K 2=1282.71K ψ1=0.1471?ω2?0.1558ω+1.1234 （6） ψ2=?1.1836?ω2?0.3761ω?0.5289（7） ψ3=?0.0455?ω2+1.8719ω?0.3907 （8） ω=0.177e +0.339 （9） e =0.6108?exp 17.27 Τ0 ?273 237.3+Τ 0?273 ?RH （10） 先来说明单通道算法公式（1）中γ，L sensor ，δ，ψ1，ψ2，ψ3等这些参数的计算过程，地表比辐射率ε的计算过程稍后在说明。 （1）对于ψ1，ψ2，ψ3的计算，只要查阅资料得知相对湿度RH ，与温度Τ0后，就可以算出大气中水蒸汽的含量ω，进而可以根据公式算出ψ1，ψ2，ψ3。 （2）对于L sensor 的计算，也就是辐射校正的过程，主要目的在于把影像中像元的灰度值转化成辐亮度L sensor ，公式（4）中的L min λ ，L max λ ，Q max 在影像头文件中可以找到，Q DN 就是所要进行校正的影像。在ENVI 中的操作如下： Basic tools → band math ,然后点开出现如下左侧对话框： 对于ETM 数据热红外波段高增益就是L sensor =3.2+9.45?Q DN 255 然后点ok 出现如下右侧对话框：

大学物理实验数据处理基本方法

实验数据处理基本方法 实验必须采集大量数据，数据处理是指从获得数据开始到得出最后结 论的整个加工过程，它包括数据记录、整理、计算与分析等，从而寻找出 测量对象的内在规律，正确地给出实验结果。因此，数据处理是实验工作 不可缺少的一部分。数据处理涉及的内容很多，这里只介绍常用的四种方 法。 1列表法 对一个物理量进行多次测量，或者测量几个量之间的函数关系，往往 借助于列表法把实验数据列成表格。其优点是，使大量数据表达清晰醒目， 条理化，易于检查数据和发现问题，避免差错，同时有助于反映出物理量 之间的对应关系。所以，设计一个简明醒目、合理美观的数据表格，是每 一个同学都要掌握的基本技能。 列表没有统一的格式，但所设计的表格要能充分反映上述优点，应注意以下几点：1．各栏目均应注明所记录的物理量的名称(符号 )和单位； 2．栏目的顺序应充分注意数据间的联系和计算顺序，力求简明、齐全、有条理； 3．表中的原始测量数据应正确反映有效数字，数据不应随便涂改，确实要修改数据时， 应将原来数据画条杠以备随时查验； 4．对于函数关系的数据表格，应按自变量由小到大或由大到小的顺序排列，以便于判 断和处理。 2图解法 图线能够明显地表示出实验数据间的关系，并且通过它可以找出两个 量之间的数学关系，因此图解法是实验数据处理的重要方法之一。图解法 处理数据，首先要画出合乎规范的图线，其要点如下： 1.选择图纸作图纸有直角坐标纸 ( 即毫米方格纸 ) 、对数坐标纸和 极坐标纸等，根据 作图需要选择。在物理实验中比较常用的是毫米方格纸，其规格多为17 25 cm 。 2.曲线改直由于直线最易描绘 , 且直线方程的两个参数 ( 斜率和截距 ) 也较易算得。所以对于两个变量之间的函数关系是非线性的情形，在用图解法时 应尽可能通过变量代换 将非线性的函数曲线转变为线性函数的直线。下面为几种常用的变换方法。 ( 1) xy c ( c 为常数 ) 。 令 z 1，则 y cz，即 y 与 z 为线性关系。 x ( 2) x c y ( c 为常x2，y 1 z ，即 y 与为线性关系。

基于热红外波段的地表温度反演实验报告

遥感原理与应用 地表温度反演 实验报告 专业：地理信息系统 班级：XXXXXXXX 姓名：XXX 学号：XXXXXX 成绩： 指导教师：XXX 2014年12月17日

一. 实验目的 1. 根据实际需要，学会在网上（如中国科学院遥感与数字地球研究所数据共享网）下载研究区内的遥感数据； 2. 掌握在ENVI中实现简单的地表温度反演的原理与步骤。 二. 实验任务 1. 在中国科学院遥感与数字地球研究所数据共享网上订购并下载覆盖郫县的TM影像； 2. 在ENVI中实现简单的地表温度反演算法。 三. 实验数据 在中国科学院遥感与数字地球研究所数据共享网上下载的覆盖郫县地区的TM影像。 四. 实验原理 图1 TM影像地表温度反演流程

1. 地表温度(Land Surface Temperature)反演公式为: 2 1(1) K LST K In R ε= +， 其中，R m DN d =?+，2111607.76K W m sr m μ---=???，21260.56K K =。 2. 根据TM 辐射定标原理，热红外波段表观辐亮度可以进一步写作: max min 6min 255 L L R DN L -= ?+， 其中LmaxBand6=15.303 , LminBand6=1.238。 3. 地表比辐射率ε为同温度下地表辐射能与黑体辐射能的比率，其可以表示为: 1.0090.047(In ) (0)NDVI NDVI ε=+>， 其中，4343 TM TM NDVI TM TM -=+，当0NDVI <=时(如水体)地表比辐射率取常数1。 五. 实验步骤 1. TM 数据下载 数据查询和下载网址https://www.wendangku.net/doc/c114277375.html,/query .html ，界面如图2所示。 图2 中国科学院遥感与数字地球研究所数据共享网址界面

实验数据处理的几种方法

实验数据处理的几种方法 物理实验中测量得到的许多数据需要处理后才能表示测量的最终结果。对实验数据进行记录、整理、计算、分析、拟合等，从中获得实验结果和寻找物理量变化规律或经验公式的过程就是数据处理。它是实验方法的一个重要组成部分，是实验课的基本训练内容。本章主要介绍列表法、作图法、图解法、逐差法和最小二乘法。 1.4.1 列表法 列表法就是将一组实验数据和计算的中间数据依据一定的形式和顺序列成表格。列表法可以简单明确地表示出物理量之间的对应关系，便于分析和发现资料的规律性，也有助于检查和发现实验中的问题，这就是列表法的优点。设计记录表格时要做到：（1）表格设计要合理，以利于记录、检查、运算和分析。 （2）表格中涉及的各物理量，其符号、单位及量值的数量级均要表示清楚。但不要把单位写在数字后。 （3）表中数据要正确反映测量结果的有效数字和不确定度。列入表中的除原始数据外，计算过程中的一些中间结果和最后结果也可以列入表中。 （4）表格要加上必要的说明。实验室所给的数据或查得的单项数据应列在表格的上部，说明写在表格的下部。 1.4.2 作图法 作图法是在坐标纸上用图线表示物理量之间的关系，揭示物理量之间的联系。作图法既有简明、形象、直观、便于比较研究实验结果等优点，它是一种最常用的数据处理方法。 作图法的基本规则是： （1）根据函数关系选择适当的坐标纸（如直角坐标纸，单对数坐标纸，双对数坐标纸，极坐标纸等）和比例，画出坐标轴，标明物理量符号、单位和刻度值，并写明测试条件。 （2）坐标的原点不一定是变量的零点，可根据测试范围加以选择。，坐标分格最好使最低数字的一个单位可靠数与坐标最小分度相当。纵横坐标比例要恰当，以使图线居中。 （3）描点和连线。根据测量数据，用直尺和笔尖使其函数对应的实验点准确地落在相应的位置。一张图纸上画上几条实验曲线时，每条图线应用不同的标记如“+”、“×”、“·”、“Δ”等符号标出，以免混淆。连线时，要顾及到数据点，使曲线呈光滑曲线（含直线），并使数据点均匀分布在曲线（直线）的两侧，且尽量贴近曲线。个别偏离过大的点要重新审核，属过失误差的应剔去。 （4）标明图名，即做好实验图线后，应在图纸下方或空白的明显位置处，写上图的名称、作者和作图日期，有时还要附上简单的说明，如实验条件等，使读者一目了然。作图时，一般将纵轴代表的物理量写在前面，横轴代表的物理量写在后面，中间用“～”

Landsat8 TIRS 地表温度反演

热红外遥感（Infrared Remote Sensing）是指传感器工作波段限于红外波段范围之内的遥感。即利用星载或机载传感器收集、记录地物的热红外信息，并利用这种热红外信息来识别地物和反演地表参数如温度、湿度和热惯量等。目前有很多的卫星携带了热红外传感器，包括ASTER、AVHRR、MODIS、TM/ETM+/ TIRS等。 目前，地表温度反演算法主要有以下三种：大气校正法（也称为辐射传输方程：Radiative Transfer Equation——RTE）、单通道算法和分裂窗算法。 本实例是基于大气校正法，利用Landsat8 TIRS反演地表温度。 基本原理：首先估计大气对地表热辐射的影响, 然后把这部分大气影响从卫星传感器所观测到的热辐射总量中减去, 从而得到地表热辐射强度, 再把这一热辐射强度转化为相应的地表温度。 具体实现为：卫星传感器接收到的热红外辐射亮度值Lλ由三部分组成：大气向上辐射亮度L↑，地面的真实辐射亮度经过大气层之后到达卫星传感器的能量；大气向下辐射到达地面后反射的能量。卫星传感器接收到的热红外辐射亮度值Lλ的表达式可写为（辐射传输方程）： Lλ = [εB(T S) + (1-ε)L↓]τ+ L↑（1.1） 式中，ε为地表比辐射率，T S为地表真实温度(K)，B(T S)为黑体热辐射亮度，τ为大气在热红外波段的透过率。则温度为T的黑体在热红外波段的辐射亮度B(T S)为： B(T S) = [Lλ - L↑- τ(1-ε)L↓]/τε（1.2） T s可以用普朗克公式的函数获取。

T S = K2/ln(K1/ B(T S)+ 1) （1.3） 对于TM，K1 =607.76 W/(m2*μm*sr)，K2 =1260.56K。 对于ETM+，K1=666.09 W/(m2*μm*sr)，K2 =1282.71K。 对于TIRS Band10，K1= 774.89 W/(m2*μm*sr)，K2 = 1321.08K。 从上可知此类算法需要2个参数：大气剖面参数和地表比辐射率。大气剖面参数在NASA提供的网站(https://www.wendangku.net/doc/c114277375.html,/)中，输入成影时间以及中心经纬度可以获取大气剖面参数。适用于只有一个热红外波段的数据，如Landsat TM /ETM+/TIRS数据。 主要内容就是使用BandMath工具计算公式（1.2）和公式（1.3），处理流程如下图所示。

大学物理实验数据处理方法总结

有效数字 1、有效数字不同的数相加减时，以参加运算各量中有效数字最末一位位数最高的为准，最后结果与它对其，余下的尾数按舍入规则处理。 2、乘除法以参与运算的数值中有效位数最少的那个数为准，但当结果的第1位数较小，比如1、2、3时可以多保留一位（较小：结果的第一位数小于 有效数字最少的结果第一位数）！ 例如：n=tg56° θ=56° d θ=1° θθθθθ2cos d d d dtg dn == 为保留） （，带入848.156n 15605.018056cos 1cos 22=?=∴?=??=≈?=?= ?tg n θθπθθ 3、可以数字只出现在最末一位：对函数运算以不损失有效数字为准。 例如：20*lg63.4 可疑最小位变化0.1 Y=20lgx 01.04 .631.010ln 2010ln 20ln 10ln 20≈===x dx dx dx x d dy 04.364.63lg 20=∴ 4、原始数据记录、测量结果最后表示，严格按有效数字规定处理。（中间过程、结果多算几次） 5、4舍5入6凑偶 6、不估计不确定度时，有效数字按相应运算法则取位；计算不确定度时以不确定度的处理结果为准。 真值和误差 1、 误差=测量值-真值 ΔN=N-A 2、 误差既有大小、方向与政府。 3、 通常真值和误差都是未知的。 4、 相对约定真值，误差可以求出。 5、 用相对误差比较测量结果的准确度。 6、 ΔN/A ≈ΔN/N 7、 系统误差、随机误差、粗大误差 8、 随机误差：统计意义下的分布规律。粗大误差：测量错误 9、 系统误差和随机误差在一定条件下相互转化。 不确定度 1、P （x ）是概率密度函数 dx P dx x x P p )x (之间的概率是测量结果落在+当x 取遍所有可能的概率值为1. 2、正态分布且消除了系统误差，概率最大的位置是真值A 3、曲线“胖”精密度低“瘦”精密度高。 4、标准误差：无限次测量?∞∞-=-2 )()(dx X P A X x ）（σ 有限次测量且真值不知道标准偏

肿瘤学研究常用实验技术及方法

精心整理 肿瘤学研究常用实验技术及方法复习题 1.表观遗传的科学内涵 三个特点？ a)基因功能完全、稳定地丧失，有遗传性； b)细胞外环境长期作用的结果，有可逆性； c)分析方法灵敏，可检出少数细胞的变化 2. DHPLC 3. i. ii. 梯度中 iii. SDS- iv. 1.丙烯酰胺有神经毒性，可经皮肤、呼吸道等吸收，故操作时一定要注意防护。 2.制备聚丙烯酰胺凝胶时，小心防止凝胶渗漏。 3.蛋白加样量要合适，加样量太小，条带不清晰，加样量太大，则泳道超载，条带过宽而重 叠，甚至覆盖至相邻泳道。 4.对多种蛋白质而言，电流大则电泳条带清晰，但电流过大，玻璃板会因受热而破裂，故适 合的电流为玻璃板微热。 5.胶转膜后应在膜上标记好正反面及电泳方向。 6.电转移时应注意勿将滤膜和胶的位置放反，而且滤纸、滤膜和胶应等大，以免短路。

4.组织培养/基因的克隆、转化与表达 i.请列举至少3种流式细胞仪的用途? 1.定量检测细胞膜、细胞质和细胞核中的各种细胞成分 2.研究细胞的各种功能状态(细胞增殖，细胞凋亡、细胞分化、酶活性、细胞膜通透性、 氧化还原状态、吞噬性等)。 3.定量检测血清中的各种可溶性生物分子成分 临床中应用：1.白血病和淋巴瘤的免疫分型 1. 2. 3.细胞 磷酸酶（phosphatase） 激酶（Kinase） 连接酶（Ligase） 5.生物质谱简介 生物质谱在蛋白组学中的应用？ 1.ProteinsMWmeasurements； 2.ProteinIdentification；(ID) 3.Peptideandproteinprofiling； 4.ImagingMS； 5.ProteinPTMs；5.ProteinQuantification 6.《常用电泳技术基本原理及应用》课程复习题

landsat遥感影像地表温度反演教育教案(大气校正法)

基于辐射传输方程的Landsat数据地表温度反演教程 一、数据准备 Landsa 8遥感影像数据一景，本教程以重庆市2015年7月26日的=行列号为（128，049）影像（LC81280402016208LGN00）为例。 同时需提前查询影像的基本信息（详见下表） 注：基本信息在影像头文件中均可查询到，采集时间为格林尼治时间。 二、地表温度反演的总体流程 三、具体步骤

1、辐射定标 地表温度反演主要包括两部分，一是对热红外数据，二是多光谱数据进行辐射 定标。 （1）热红外数据辐射定标 选择Radiometric Correction/Radiometric Calibration。在File Selection 对话框中，选择数据LC81230322013132LGN02_MTL_Thermal，单击Spectral Subset选择Thermal Infrared1（10.9），打开Radiometric Calibration面板。 （2）多光谱数据辐射定标 Scale factor 不能改变，否则后续 计算会报错。保持默认1即可。

选择要校正的多光谱数据 “LC81230322013132LGN02_MTL_MultiSpectral” 进行辐射定标。 因为后续需要对多光谱数据进行大气校正，可直接单击Apply Flaash Settings，如下图。 注意与热红外数据辐射定标是的差 别，设置后Scale factor值为0.1。 2、大气校正 本教程选择Flaash 校正法。FLAASH Atmospheric Correction，双击此工具，打开辐射定标的数据，进行相关的参数设置进行大气校正。 注意：如果在多光谱数据辐射定标时Scale factor值忘记设置，可在本步骤中打开辐射定标数时设置single scale faceor 值为0.1，若已设置，则默认值为1即可。

实验技术岗位工作量计算办法

实验技术岗位工作量计算办法 为了加强实验队伍建设，规范我系的实验室管理，充分发挥实验教学的资源优势，进一步调动实验技术人员的工作积极性，逐步提高实验教学质量，结合我系情况，特制定本办法。 第一条内容及定额 实验人员工作量包括实验教学工作量、实验室管理工作量以及实验室建设奖励工作量等三大部分。 按照学校的有关规定，参照教师工作量定额，实验人员全年工作量按每年工作40周，每周5天，每天8小时计算，全年应完成工作量1600小时。为便于计算和考核，与教学工作量相一致，将此工作量折算为实验工作量，一个标准实验工作量约为7小时，一个实验人员全年工作量定额为228个。 第二条实验人员工作量的统计与考核 实验人员的考核办法参照《为了加强我系实验教学管理，客观评价实验技术人员的教学工作成绩和效果，依据《兰州城市学院教职工年度考核暂行办法》以及《兰州城市学院教师职务岗位考核实施方案》，结合我系具体情况，制定本考核办法。 一、考核组织安排 1、成立化学系实验技术人员考核领导小组具体负责每个学期实验技术人员的实验室工作考核。 2、实验技术人员考核，每年和教师同步进行，但每学期末须填写“实验教学任务统计表”和“实验技术人员工作量表”，一式叁份，一份本人自存，一份交教研室，一份交系办。 二、考核范围与内容 凡承担全校本、专科生、成教生(有学校计划正式下达教学任务的课程)教学任务的实验技术人员均应进行工作量考核。 (一)实验技术人员考核内容 1、实验教学：实验辅导、实验准备及实验室管理； 2、实验室建设：实验方法的改进、实验技术的研究、实验仪器设备的改进、实验教具、材料(图表、标本、组织切片等)制作； 3、实验室管理：实验材料的领用配制、整理、仪器设备的维修与管理，实验室清洁卫生的保持与环境优化等。 三、考核程序 实验技术人员考核由系实验技术人员考核领导小组组织实施，与全系教职工年度考核同步进行。具体考核程序 (一)实验技术人员考核 1、实验技术人员考核：本人应填写实验技术教学工作考核表，交实验室主任，由实验室主任召集相关专业的教研室进行工作质量测评与认定。 2、实验技术人员考核内容：思想政治及服务育人、工作态度及表现、实验任务完成情况、物资管理状况、实验室安全、仪器设备完好率、卫生状况、实验技能、实验室建设等。 四、考评结果的确定 1、系实验技术人员考核领导小组根据考核内容，对考核对象的全年工作进行综合评价，并写出综合评语。 2、根据综合考核情况，按优秀、合格、基本合格、不合格四个等级确定被考核对象的相应等级，填写《兰州城市学院化学系实验技术人员教学工作考核表》一

数值分析实验报告之常微分方程数值解

数学与计算科学学院实验报告 实验项目名称常微分方程数值解 所属课程名称数值方法B 实验类型验证 实验日期 2013.11.11 班级 学号 姓名 成绩

【实验过程】（实验步骤、记录、数据、分析） 注：以下图形是通过Excel 表格处理数据得出，并未通过MATLAB 编程序所得！ 1、1(0)1dy y x dx y ?=-++???=? 由题可知精确解为：x y x e -=+，当x=0时，y(x)=0。 h=0.1 表1 h=0.1时三个方法与精确值的真值表 图1 h=0.1时三个方法走势图 步长 Euler 法 预估校正法 经典四阶库 精确值 0.1 1.010000 1.005000 1.004838 1.249080 0.2 1.029000 1.019025 1.018731 1.055455 0.3 1.056100 1.041218 1.040818 1.091217 0.4 1.090490 1.070802 1.070320 1.131803 0.5 1.131441 1.107076 1.106531 1.176851 0.6 1.178297 1.149404 1.148812 1.226025 0.7 1.230467 1.197211 1.196586 1.279016 0.8 1.287421 1.249975 1.249329 1.335536 0.9 1.348678 1.307228 1.306570 1.395322 1.0 1.413811 1.368541 1.367880 1.458127

h=0.05（此时将源程序中i的围进行扩大，即for(i=0;i<20;i++)） 表2 h=0.05时三个方法与精确值的真值表步长Euler法预估校正法经典四阶库精确值 0.05 1.002500 1.001250 1.001229 1.011721 0.10 1.007375 1.004877 1.004837 1.024908 0.15 1.014506 1.010764 1.010708 1.039504 0.20 1.023781 1.018802 1.018731 1.055455 0.25 1.035092 1.028885 1.028801 1.072710 0.30 1.048337 1.040915 1.040818 1.091217 0.35 1.063421 1.054795 1.054688 1.110931 0.40 1.080250 1.070436 1.070320 1.131801 0.45 1.098737 1.087752 1.087628 1.153791 0.50 1.118800 1.106662 1.106531 1.176851 0.55 1.140360 1.127087 1.126950 1.200942 0.60 1.163342 1.148954 1.148812 1.226025 0.65 1.187675 1.172193 1.172046 1.252062 0.70 1.213291 1.196736 1.196585 1.279016 0.75 1.240127 1.222520 1.222367 1.306852 0.80 1.268121 1.249485 1.249329 1.335536 0.85 1.297215 1.277572 1.277415 1.365037 0.90 1.327354 1.306728 1.306570 1.395322 0.95 1.358486 1.336900 1.336741 1.426362 1.00 1.390562 1.368039 1.367880 1.458127 图2 h=0.05时三个方法走势图

erdas 北京地表温度反演_实习报告

Landsat TM6 地表温度反演实习报告 实习目的： 1、TM Level 1 数据的热红外波段辐射定标：学会阅读头文件，找出所需定标参数；利用定标参数将TM图像热红外波段DN值转换为辐射亮度； 2、运用单通道法，反演地表温度反演 实习步骤： 1.加载图像：import—>选择需要加载的图像 2.分部建模 2.1求算NDVI的建模如图所示

第三第四波段的辐射定标运算的增益和偏移均来自数据的头文件。 2.2第六波段辐射定标计算

说明：第六波段的辐射定标计算的增益和偏移不能再用头文件中的增益和偏移，否则误差会很大造成无法显示，因此必须在老师给的表格中查找。 然后再用老师给的公式进行计算. 从表格中找出L min 和L max 输入下面公式 255 G min max L L rescale -= min B L rescale = rescale cal rescale B Q G L +?=λ 即可求出增益和偏移，然后进行辐射定标运算即可求出所需结果。

2.3求解地表反射率（λε）的建模运算。 所使用的经验公式根据Van 的经验公式： )ln(047.00094.1NDVI +=λε 2.4求解)T (B s λ的建模。)T (B s λ为温度为s T 的黑体在热红外波段的辐射亮度。

使用的公式是 L L o o o s ↓ ↑ -- -=λ λ λ λλλλλεετε1L )T (B ）（，其中 L o ↓ λ 表示大气向下辐射亮度，模拟结果为1.68 Wm -2um -1Sr -1，L o ↑ λ表示大气向上辐射亮度， 模拟结果为1.74 Wm -2um -1Sr -1，λτo 为大气在热红外波段的透过率，模拟结果为0.77。λε为上一步求解的结果。 2.5反演温度的建模

技术试验及其方法(习题2)

技术试验及其方法（习题2） 1.阿什塔比拉河桥的垮塌事件提醒人们（C ） A.技术的发展离不开设计 B.设计是技术活动的核心 C.技术试验可使设计得以改进和完善，使风险和失误降到最低 D.人机关系是在技术设计中必须考虑的因素之一 2．陈晨同学动手制作了一个小板凳,他通过试验来检验小板凳承重力和稳定性，他的试验方法是( D ) A.优选试验法 B.模拟试验法 C.虚拟试验法 D.强化试验法 3.对安全帽进行超载试验是一种（ C ）试验。 A.预测 B.信息 C.性能 D.农业 4.以往电路图的设计都是手工绘制，繁琐且易出错。随着电子技术的发展，尤其是电子计算机技术的发展，现在用专门的电路图设计软件大大方便了电路设计，而且还可以进行仿真测试。从技术试验的角度看，“仿真测试”属于（ B ） A.优选试验法 B.模拟试验法 C.虚拟试验法 D.强化试验法 5.据报道，一名男子将镀铜的钢板当作纯铜欲出售给南京一家废品收购站，不料老板拿起“铜板”顺势在水泥地面上摩擦了两下，“铜板”边角处立刻变成灰白色，使骗子露馅。这位老板使用的试验方法是（ B ） A.优选试验法 B.强化试验法 C.模拟试验法 D.虚拟试验法 6．李宁自己动手制作了一张座椅，他对座椅的稳定性做了如下试验，你认为合理的一项是（ D ） A.亲自坐在椅子上感受一下 B.放一些重物在椅子上，看是否能承受 C.用电风扇对椅子吹风看能否吹倒 D.坐在椅子上，前后左右晃一晃看是否稳固 7.水利大坝从理论，设计到开工建设进行了许多次实验。一般采取下列的什么方法（B ） A.优选实验法 B.模拟实验法 C.虚拟实验法 D.强化实验法 8．安全带的研制是通过哪项技术试验法而得到检验的（ D ） A．虚拟试验法 B．强化试验法 C．优选试验法 D．模拟试验法 9.早在远古时代，人们就知道利用固体互相刻划来区分材料的软硬，并据此来选用材料。例如，皂石的硬度低，用于制作器皿和装饰品；炬石坚硬，用于制作工具和刀剑等。至今，硬度仍用来表示材料的软硬程度。硬度值的大小不仅取决于材料的成分和显微组织，而且还取决于测量方法,因此对于材料硬度的测试以下较为合适的是:( B ) Ａ、预测试验Ｂ、性能试验Ｃ、模拟试验Ｄ、强化试验 10．汽车的碰撞实验是（ B ） A.优选实验法 B.模拟实验法 C.移植实验法 D.强化实验法 11.某市铸管厂承接了一项供水水管的任务，为保证产品质量，需对所生产的水管进行质量检测，采用的是利用注水加压的方法检测水管的强度。这种试验方法属于（A ） A.性能试验 B.预测试验 C.优化试验 D.信息试验 12.对不同品种的水稻进行对比试验，这种试验方法是（B ） A.模拟试验法 B.优选试验法 C.移植试验法 D.强化试验法 15.航天员穿着舱外航天服在水下进行失重和出舱活动任务训练，这种试验方法是（ A ） A.模拟试验法 B.优选试验法 C.虚拟试验法 D.移植试验法 16.普通自行车设计过程中，以下试验不必要的是（A ） A. 与快速行驶的汽车进行撞击试验 B.承重试验

数值分析实验报告

学生实验报告实验课程名称 开课实验室 学院年级专业班 学生姓名学号 开课时间至学年学期

if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); %换行 end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); %消去end end x=zeros(length(b),1); %回代求解 x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end y=x; format short;%设置为默认格式显示,显示5位 （2）建立MATLAB界面 利用MA TLAB的GUI建立如下界面求解线性方程组： 详见程序。 五、计算实例、数据、结果、分析 下面我们对以上的结果进行测试，求解：

? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - - - 7 2 5 10 13 9 14 4 4 3 2 1 13 12 4 3 3 10 2 4 3 2 1 x x x x 输入数据后点击和，得到如下结果： 更改以上数据进行测试，求解如下方程组： 1 2 3 4 43211 34321 23431 12341 x x x x ?? ???? ?? ???? ?? ???? = ?? ???? - ?? ???? - ???? ?? 得到如下结果：

定量遥感_地表温度反演

遥感数字影像处理 作品名称：黄河三角洲地表温度反演 +学号： 小组成绩：

一、概述 1、作业背景： 地表温度是很多环境模型的一个重要参数，在大气与地表的能量与物质交换，天气预报，全球洋流循环，气候变化等研究领域有重要的应用。利用热红外遥感可以得到大围的地表温度面状信息，与传统的地表温度测量方式相比，具有快速、便捷、测量围大、信息连续等特点，因此利用热红外遥感数据反演地面温度得到了广泛的应用 2、作业意义： 黄河三角洲是黄河携带大量泥沙在渤海凹陷处沉积形成的冲积平原，位处黄河入处的黄河三角洲自然保护区正是以保护河口湿地生态系统和珍稀、濒危鸟类为主的湿地类型保护区。以利津为顶点，北到徒骇河口，南到小清河口，呈扇状三角形，面积5,450平方公里。地面平坦，在海拔10公尺以下。向东撒开的扇状地形，海拔高程低于15米，面积达5450平方公里。三角洲属，温带季风性气候。四季分明，光照充足，区自然资源丰富。 黄河口湿地生态旅游区占地23万亩，都处在黄河三角洲之，地貌以芦苇沼泽，湿地为主，其次为河口滩地，带翅碱蓬盐滩湿地，灌丛疏林湿地以及人工槐林湿地等。集自然景观与人文景观为一体，既有沧海桑田的神奇与壮阔，又有黄龙入海的壮观和长河落日的静美，是人们休闲、度假、观光科普的最佳场所。 二、数据介绍 数据来自地理空间数据云，Landsat 4-5 TM（陆地卫星4、5号，1982年发射后运行至今，携带有TM传感器）的相关遥感影像作为研究数据，研究黄河三角洲温度分布状况。 实验数据：2010年9月11号黄河三角洲图像（中心经度：118.8878w，中心纬度：37.4815n） 三、基本概念及技术流程图 3.1、基本概念：