RAPID模型导入场景注意事项

1.RAPID中方格内导入并显示场景时场景文件必须为.map后缀，in不可以。


2.RGMDISP.exe中导入场景时需要导入poly.in，后缀为in的文件，方可显示。

3.需要运行rapid模型时，.map文件为非必须文件，in文件为必须文件。
3.1在软件中手动运行模型时，需要先创建.map文件，然后创建对应的in文件，设定一下太阳和观测角度，，然后需要设置camera.in文件，里面需要填写激光雷达的x y，z位置，扫描的水平扫描间隔、竖直扫描间隔、视场角等，
20 20 50
0.8 0.8 0.4
0
设置完成后即可开始运行Go—64bit模型。

3.2使用代码创建场景并运行模型时，先使用代码直接创建poly.in文件，将省略手动创建.map文件的过程，然后代码中设定太阳和观测角度，然后代码执行rapid的ref64和lidar64，得到waveform_sgl_1000_1000.dat高光谱激光雷达回波波形文件。
其中，代码中没有设置高光谱激光雷达的位置，也没有产生camera.in文件，
lidar64 sgl -FWHM 4 -foot 35 -size 2000 ，
源代码默认光斑直径为35米，图像分辨率为2000，即如果以水平切割场景的话，切割层展开到二维平面上时最多有2000个像元。

以下为一次运行过后，matlab输出显示结果：

Welcome to use Cylinder Scene Generation Program...
input from:D:\0\poly\stand1.dat
result in:D:\0\poly\poly_stand1.in
Start running...

ans =

   95.0332   87.8759

...


ans =

   95.0332   87.8759

New Decidious Tree has polygon with 247 current total polygon is 6070
scene done at 6070
scene LAI is 1.0559
tree number is 10
 ============================================== 
    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  6070 
     >>Components=  3  BANDS=  4 
   >>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   6070 
   >>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.4089 
     >>xspan = max x =   32.2731  dx =  32.2731 
     >>yspan = max y =   32.0651  dy =  32.0651 
     >>rr=  24.3007 
     >>delx, dely =   48.6014  48.6014 
 ****************************************** 
     >>a,b,c,d=  41.1099  1000.00  41.1099  1000.00 
   >>will sort   24280  vertex 
   >>  12369  common points found 
 **************************** 
 Number of polygons is   6070 
 **************************** 
   >>finding common edges... 
   >>total   24280  edges 
   >>max zpoint =   24280 
     >>bits4 one  
     >>memset skips done  
   >>  8010  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 =   41.1099 
   >>b     =   1000.00 
   >>d     =   1000.00 
   >>Computing centers, normals and sizes... 
   >>pixel & real LAI is   3  2.59075  2.62877 
   >>LAD cannot be computed at this time. 
   >>XLAD cannot be computed at this time. 
     ->after polyinfo  41.1099  1000.00  41.1099  1000.00 
     >>reading texture.in and leaf_clump.in for HETOBJ 
     >>PAI =   0.755007 m2*m-2 
     >>Max number of sub-leafs =   355 
     >>Min Sub-Leaf size=  8.42643 cm > defined  0.200000 cm 
     >>reading pnorm.dat,pcenter.dat 
   >>slot is   1.977261E-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 =   27473110  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 
 done 
 number of directions is   40 
 before allocated   306 
 before memset xsol  0.000000 
 allocate done 
   
  laser pulse shape is:  
             ..... 
             ...... 
             ......... 
             ........... 
             .............. 
             ................ 
             .................. 
             ................... 
             .................... 
             ................... 
             .................. 
             ................ 
             .............. 
             ........... 
             ......... 
             ...... 
             ..... 
    
 reading from "darea.dat" 
 before solve 
 penetrating position   1000  1000 
 doing single scattering ... 
 Solution found!, Saving results.... 
 begin write 
 return pulses stored in "waveform_sgl_*.dat" 
Solving done.