RAPID批处理产生不同叶绿素含量的树

最新推荐文章于 2022-07-06 17:10:27 发布

B博士

最新推荐文章于 2022-07-06 17:10:27 发布

阅读量154

点赞数

分类专栏：遥感文章标签： rapid 叶绿素含量高光谱激光雷达波形模拟

本文链接：https://blog.csdn.net/qq_37970770/article/details/119764007

版权

遥感专栏收录该内容

24 篇文章 16 订阅

订阅专栏

实验中利用MATLAB的randperm函数为四株植物（三棵树和一株草）生成不同叶绿素含量，并尝试用RAPID进行61波段的高光谱激光雷达波形模拟。过程中遇到错误，提示可能仅支持32个波段模拟，同时在执行过程中出现‘非法字符在数值输入中’的FORTRAN运行时错误，导致程序无法继续。问题可能与文件读取、波段数量限制或数据格式有关。

摘要由CSDN通过智能技术生成

本实验中给定三棵树和一株草，树高分别为20m，10m，2m（树干1m），草高0.5m。
1.需要产生对应不同的叶绿素含量，使用matlab中的随机函数randperm(n,k)
p = randperm(n,k) 返回一行从1到n的整数中的k个，而且这k个数也是不相同的。

randperm(80,4)

ans =

41    34    69    80

将其记录并赋值给上述四株植物，并执行rapid高光谱激光雷达波形模拟。

模拟61个波段的回波信号时报错，可能是只支持32个波段模拟？

New Decidious Tree has polygon with 70 current total polygon is 4025

ans =

    1.5394    0.9624

New Decidious Tree has polygon with 80 current total polygon is 4105
scene done at 4105
scene LAI is 0.13195
tree number is 4
 ============================================== 
    Welcome to use RAPID3  
 ╔═╤═╤═╤═╤═╤═╤═╤═╤═╤═╤══╗ 
 ║　↓Ｒ↓　↓　↓Ｒ↓　↓　↓Ｒ↓　↓　↓Ｒ 
 ║　↓Ａ↓　↓　↓Ａ↓　↓　↓Ａ↓　↓　↓Ａ 
 ║　↓Ｐ↓　↓　↓Ｐ↓　↓　↓Ｐ↓　↓　↓Ｐ 
 ║　↓Ｉ↓　↓　↓Ｉ↓　↓　↓Ｉ↓　↓　↓Ｉ 
 ║　↓Ｄ↓　↓　↓Ｄ↓　↓　↓Ｄ↓　↓　↓Ｄ 
 ╚═╧═╧═╧═╧═╧═╧═╧═╧═╧═╧══╝  
    ref64.exe, 2019.12.12, Huaguo Huang 
 ============================================== 
   >>parsing command line 
    >>start to run the solar direct light code  
 ------------------------------------------ 
     >>use HIGH RESOLUTION  2000 
     >>Polygon number is  4105 
     >>Components=  6  BANDS=  61 
   >>allocate all done 
     ->memset all done 
     >>sun zenith  is   0.000000 -1.00000 
     >>sun azimuth is   0.000000 -1.00000 
   >>checking if cat.rt exists ... 
   >>cat.rt not found,establish it for each poly... 
   >>reading polygon vertices from "poly.in        " 
   >>number of polygons will be   4105 
   >>W A N R I N G: shoot.in missing 
   >>default values applied for all canopies 
   >>  3  cats established 
     >>leaf00     
     >>stem00     
     >>soil       
     >>the order, cat.name, rhos/taus & fsign          of each poly stored in cat.rt 
   >>W A N R I N G: leaf_lad.in missing 
   >>Sphere LAD distribution used 
     ->right after getpol in refhigh4.f 
     >>zspan = max z =   20.4070 
     >>xspan = max x =   30.0000  dx =  30.0000 
     >>yspan = max y =   30.0000  dy =  30.0000 
     >>rr=  23.5396 
     >>delx, dely =   47.0791  47.0791 
 ****************************************** 
     >>a,b,c,d=  42.4392  1000.00  42.4392  1000.00 
   >>will sort   16420  vertex 
   >>  11057  common points found 
 **************************** 
 Number of polygons is   4105 
 **************************** 
   >>finding common edges... 
   >>total   16420  edges 
   >>max zpoint =   16420 
     >>bits4 one  
     >>memset skips done  
   >>  7296  edges won't be drawn... 
 default soil high limit = 0.001 
 no. of soil cat & the order:  1  1 
 no. of leaf cat & the order:  1  3 
   >>a = c =   42.4392 
   >>b     =   1000.00 
   >>d     =   1000.00 
   >>Computing centers, normals and sizes... 
   >>pixel & real LAI is   3  0.289795  0.330212 
   >>LAD cannot be computed at this time. 
   >>XLAD cannot be computed at this time. 
     ->after polyinfo  42.4392  1000.00  42.4392  1000.00 
     >>reading texture.in and leaf_clump.in for HETOBJ 
     >>PAI =   0.127068 m2*m-2 
     >>Max number of sub-leafs =   379 
     >>Min Sub-Leaf size=  8.16251 cm > defined  0.200000 cm 
     >>reading pnorm.dat,pcenter.dat 
   >>slot is   2.923977E-02 
   >>slots assigned. 
   >>doing solar direct light projection 
   >>begin do_far 
     >>Saving projected Image s0.raw          ... 
     >>save done  
     >>saving projected picsub Image sub0.raw   ... 
     >>save done  
     >>Saving projected Zbuffer z0.raw          ... 
     >>save done  
     >>Saving projected component ID c0.raw          ... 
     >>save done  
   >>begin hist_for_irad 
     >>subcos done for hist_for_irad 
     >>iterate each polygon for darea 
     >>darea stored in "darea.dat" 
     >>deallocate crsubleaf done. 
   >>Specular part done. 
 restoring projected Image s0.raw     ... 
 restore done  
 restoring projected Image z0.raw     ... 
 restore done  
 s0.raw and z0.raw with size of   2000 
 memory used =   18341065  bytes 
 ------------------------------------------ 
    ...starting to solve lidar waveform... 
 ------------------------------------------ 
 example: -foot 2.0 -space 1.5 -FWHM 4 sgl 
 -foot 2.0 means radius of footprint is 35 meter 
 -space 5 means pulse spacing is 1.5 meter 
 -FWHM 4 means FWHM is 4 ns 
 sgl means only single scattering 
 reading polygon LAI,LAD from poly.in 
 W A N R I N G: leaf_lad.in missing 
 Sphere LAD distribution used 
 reading rhos/taus from optics.in 
 nb exceeds maxband, maxband used 
? FORTRAN Runtime Error: 
? Illegal character in numeric input 
? READ(UNIT=3,...
索引超出矩阵维度。

出错 runRAPID (line 30)
    fid = fopen(a(1).name,'rt');

出错 singlescene (line 115)
    runRAPID(strPath1,strPath2,strPath3,strPath4,SN);