应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » matlab 如何批量计算envi标准格式影像的NDVI值
303/3返回列表上一页123
查看: 3682  |  回复: 29
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

1987勇往直前

新虫 (小有名气)
引用回帖:
19楼: Originally posted by weiyinan at 2014-10-12 16:25:54
我用这个算了，会出现这个：Warning: Divide by zero，得到的图像是一幅全黑的图。但是用ENVI计算出来的结果却不是这样。我想知道原因，是不是因为我是读图以后，直接计算，没有对数据进行什么处理啊？...

首先，你先确定下你读出来的数据，各个波段是正常的矩阵，是有值的，变量查看那块点开看就行；然后，看看你的数据类型，是不是float型，如果是整型的有可能保留为0；再然后，如果是矩阵的除法，应该加个点，比如band1./band2，希望可以帮到你，有问题再交流
赞一下 回复此楼
21楼2014-10-13 10:16:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyinan

新虫 (初入文坛)
引用回帖:
21楼: Originally posted by 1987勇往直前 at 2014-10-13 10:16:42
首先，你先确定下你读出来的数据，各个波段是正常的矩阵，是有值的，变量查看那块点开看就行；然后，看看你的数据类型，是不是float型，如果是整型的有可能保留为0；再然后，如果是矩阵的除法，应该加个点，比如ba ...

我是新手，有很多地方不懂。我是这样操作的：
>> a=imread('D:\B3.tif');
>> b=imread('D:\B4.tif');
>> a1=imcrop(a,[3000 3000 200 200]);
>> b1=imcrop(b,[3000 3000 200 200]);
>> NDVI=(b1-a1)./(b1+a1);
>> imshow(NDVI)
由于图像太大，我对图像进行了裁剪。数据不是float型，是uint8型，最后的图像任然是一幅全黑色的图。这个操作过程没有错吧？还用是不是因为数据类型的问题才导致这个结果？
赞一下 回复此楼
22楼2014-10-13 19:20:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyinan

新虫 (初入文坛)
我今天对之前的两幅图做了减法，结果是全黑的图，说明两幅图是一样的吧？然后我有换了两幅图，用了相同的方法，把数据格式换了，操作是这样的：
>> a1=imread('D:\3B3.tif');
>> b1=imread('D:\4B4.tif');
>> a2=imcrop(a1,[3704 5119 200 200]);
>> b2=imcrop(b1,[3704 5119 200 200]);
>> a=im2double(a2);
>> b=im2double(b2);
>> NDVI=(b-a)./(b+a);
>> imshow(NDVI)

得到了图像。应该是NDVI的图像了吧，和ENVI做出来的NDVI进行了比较，还有点小差距。

总算是求出来了，谢谢你了。非常感谢。有问题还会请教你的。
赞一下 回复此楼
23楼2014-10-14 09:56:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1987勇往直前

新虫 (小有名气)
引用回帖:
23楼: Originally posted by weiyinan at 2014-10-14 09:56:56
我今天对之前的两幅图做了减法，结果是全黑的图，说明两幅图是一样的吧？然后我有换了两幅图，用了相同的方法，把数据格式换了，操作是这样的：
>> a1=imread('D:\3B3.tif');
>> b1=imread('D:\4B4.t ...

不必客气，我也不专业，互相学习撒~
赞一下 回复此楼
24楼2014-10-14 20:37:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyinan

新虫 (初入文坛)
有问题想请教：怎么用MATLAB读取hdf格式文件呢？我知道有hdftool,可以得到相关数据，然后怎么把数据呈现为图像呢？
赞一下 回复此楼
25楼2014-10-23 08:15:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1987勇往直前

新虫 (小有名气)
引用回帖:
25楼: Originally posted by weiyinan at 2014-10-23 08:15:26
有问题想请教：怎么用MATLAB读取hdf格式文件呢？我知道有hdftool,可以得到相关数据，然后怎么把数据呈现为图像呢？

我没直接读过hdf的，是之前网上下载过读envi标准格式的，就先envi转了，再这么读。你说的要直接读hdf的，百度一下，我试了，好几个看起来靠谱的答案，多试试
赞一下(1人) 回复此楼
26楼2014-10-23 10:03:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyinan

新虫 (初入文坛)
引用回帖:
26楼: Originally posted by 1987勇往直前 at 2014-10-23 10:03:09
我没直接读过hdf的，是之前网上下载过读envi标准格式的，就先envi转了，再这么读。你说的要直接读hdf的，百度一下，我试了，好几个看起来靠谱的答案，多试试...

恩恩，我自己再研究研究。谢谢
赞一下(1人) 回复此楼
27楼2014-10-23 15:03:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Shallow8

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by 1987勇往直前 at 2012-12-06 14:14:12
ndvi.m(835B)
http://kuai.xunlei.com/d/GAOEXGTOIQZF?p=130497
      希望可以帮到有同样问题的筒子们

真心感谢楼主的帖子！正开始学用matlab批量作图，走了太多弯路，谢谢楼主
赞一下 回复此楼
28楼2018-03-27 21:36:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Shallow8

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by 1987勇往直前 at 2012-12-06 14:14:12
ndvi.m(835B)
http://kuai.xunlei.com/d/GAOEXGTOIQZF?p=130497
      希望可以帮到有同样问题的筒子们

楼主如何查看你上传的数据啊？
赞一下 回复此楼
29楼2018-03-27 21:47:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1987勇往直前

新虫 (小有名气)
引用回帖:
29楼: Originally posted by Shallow8 at 2018-03-27 21:47:50
楼主如何查看你上传的数据啊？...

历时久远 我忘了问题是啥了 你再把你的问题详细说一下
赞一下 回复此楼
30楼2018-03-28 09:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1987勇往直前 的主题更新
303/3返回列表上一页123
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[基金申请] fileCode有新解读？ +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于代码变化问题，想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭