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
╔═╤═╤═╤═╤═╤═╤═╤═╤═╤═╤══╗
║ ↓R↓ ↓ ↓R↓ ↓ ↓R↓ ↓ ↓R
║ ↓A↓ ↓ ↓A↓ ↓ ↓A↓ ↓ ↓A
║ ↓P↓ ↓ ↓P↓ ↓ ↓P↓ ↓ ↓P
║ ↓I↓ ↓ ↓I↓ ↓ ↓I↓ ↓ ↓I
║ ↓D↓ ↓ ↓D↓ ↓ ↓D↓ ↓ ↓D
╚═╧═╧═╧═╧═╧═╧═╧═╧═╧═╧══╝
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.