halcon 光度立体

* 读取4张图像

* strFiles = 'C:/Users/wuyk/Desktop/liujing/huashang2/10'
* set_system('image_dir',strFiles)
* read_image(Images,['e','n','w','s'])

read_image(Images_e,'e.bmp')
read_image(Images_n,'n.bmp')
read_image(Images_w,'w.bmp')
read_image(Images_s,'s.bmp')

 

concat_obj(Images_e, Images_n, Images_1)
concat_obj(Images_w, Images_s, Images_2)
concat_obj(Images_1, Images_2, Images)


* 初始化参数
Tilts :=[0,90,18,270]
Slants := [45,45,45,45]
ResultType := ['gradient','albedo'] 

* 4张原始图-->1张梯度图
photometric_stereo (Images, HeightField,
Gradient, Albedo, Slants, Tilts, ResultType, 'poisson', [], [])

* 1张梯度图-->立体图
derivate_vector_field (Gradient,
MeanCurvature, 1.0, 'mean_curvature')

* 显示反射率图
dev_display (Albedo)

 

* 灰度映射，ROI

scale_image_max (MeanCurvature,ImageScaleMax)
gen_rectangle1 (ROI_0, 1001.24, -118.652, 2653.7, 5590.4)
reduce_domain(ImageScaleMax, ROI_0, ImageReduced)
write_image( ImageReduced , 'bmp' , 0 ,'save_name')