Halcon photometric stereoscopic

*  Read 4 Zhang image 

* 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)


*  Initialize parameters 
Tilts :=[0,90,18,270]
Slants := [45,45,45,45]
ResultType := ['gradient','albedo'] 

* 4 Original picture -->1 Gradient graph 
photometric_stereo (Images, HeightField,
Gradient, Albedo, Slants, Tilts, ResultType, 'poisson', [], [])

* 1 Gradient graph --> Stereogram 
derivate_vector_field (Gradient,
MeanCurvature, 1.0, 'mean_curvature')

*  Show the reflectance map 
dev_display (Albedo)

 

*  Grayscale mapping ,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')