24小时热门版块排行榜    

查看: 2686  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zzahkj

木虫 (著名写手)

[求助] 单窗算法用Landsat TM反演地表温度出错? 已有1人参与

我采用单窗算法反演地表温度 得到的亮温范围为203.371307-216.949997,这个应该是不对的 可我仔细检查了好几遍,没发现什么错误?我的步骤为:
1.求L6=(15.303-1.238)*b1/255.0+1.238(b1为第六波段像元灰度值);(在这里需要做个说明,当然也有个疑问:有的文献写:L6=0.1238+(0.005632156*b1),我看结果大概相差十倍,可能是单位问题吧,这个我也不懂,但如果用这个公式得出的L6代入下式得出的结果平均亮度温度150多,这个更加不对
2.T6=K2/ln(1+K1/L6)  其中:K1=607.76,K2=1260.56
这样的得出的像元亮度温度值范围203.371307-216.949997,平均210多 这个肯定不多吧?
如果继续做的话Ts={a(1-C-D)+[b(1-C-D)+C+D]T6-DTa}/C得出的地表温度值大多都是-60多 我想应该是亮度温度的那一步错了!
求救各位虫友,我已经检查多遍,我反演的是青岛地区2006年10月27号的TM影像,非常感谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

不要迷恋哥,哥只是一个过客!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzahkj

木虫 (著名写手)

引用回帖:
3楼: Originally posted by thangrila at 2013-11-29 11:40:23
我想请教一下楼主,您的地表发射率是怎么计算的啊?我也看了单窗算法,但是道理上我懂,具体在envi里怎么计算我就不知道了。把地表分为水体、城镇和自然表面,但是怎么才能对这三种地表进行分别计算其地表发射率啊? ...

具体可以参考http://blog.sina.com.cn/s/blog_764b1e9d010129lt.html
大哥,不是应助不要点应助回帖
你让我咋办?
不要迷恋哥,哥只是一个过客!
4楼2013-11-29 18:56:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

zzahkj

木虫 (著名写手)

我太粗心了 第一步代入的b1已经是辐射定标后的图像,应该是原图像,我等于对于定标后的图像又进行了一次定标!版主把我的帖子关了吧!
不要迷恋哥,哥只是一个过客!
2楼2012-08-21 15:49:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thangrila

新虫 (初入文坛)

【答案】应助回帖

我想请教一下楼主,您的地表发射率是怎么计算的啊?我也看了单窗算法,但是道理上我懂,具体在envi里怎么计算我就不知道了。把地表分为水体、城镇和自然表面,但是怎么才能对这三种地表进行分别计算其地表发射率啊?
谢谢
3楼2013-11-29 11:40:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thangrila

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by zzahkj at 2013-11-29 18:56:45
具体可以参考http://blog.sina.com.cn/s/blog_764b1e9d010129lt.html
大哥,不是应助不要点应助回帖
你让我咋办?...

啊?不好意思不好意思哈,当时没太注意,抱歉!
你是按照博客里那种方法计算地表发射率的吗?可我怎么总觉得那样计算是不正确的啊。,这个其实我之前也看过,(b1 le 0)*0.995+(b1 gt 0 and b1 lt 0.7)*(0.9589 + 0.086*b2 - 0.0671*b2^2)+(b1 ge 0.7)*(0.9625 + 0.0614*b2 - 0.0461*b2^2),但是他这样计算,也就是认为NDVI<0的一定是水体,0<NDVI<0.7的一定是城镇,而NDVI>0.7的一定是自然表面,这实际上是一种简便的估算,我认为这个不够准确。
5楼2013-12-02 16:12:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 5/250 2026-08-12 16:19 by 云上清扬
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 35/1750 2026-08-12 14:38 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见