本发明属于图形图像处理领域，特别是涉及一种三维模型与照片融合方法。
背景技术：
一直以来，在建设工程领域，都是利用效果图对建设工程建行效果展示，使规划管理者、民众对建设工程的效果和影响有了初步的认识。但效果图存在工程形态不够精细、工程周边仅为示意性、不真实、效果图视角单一等问题。近年，三维仿真技术逐渐在建设工程方案展示中得到应用，通过建立工程方案的三维模型，并补充建设周边现状模型，实现对建设工程及周边的模拟展示。但随着近年城乡规划等管理工作的精细化管理要求，手工建设的三维模型存在的问题开始突显，主要包括三维模型与现实世界存在色彩差异较大、工程的纹理、材质表现与现实不一致、地面细节环境不真实等问题，从而导致建设工程在三维仿真场景中进行展示时，不能完全反映建设工程建成后的效果，导致工程实施效果达不到真正管理的目的。此外，三维仿真模型还存在两个问题，一是用于建设周边现状城市建模的费用也较高，一般可达到10余万元每平方公里，从而导致利用三维仿真开展建设工程方案审查成本高；二是三维仿真对计算机的计算性能、图形展示性能要求较大，难以在互联网，尤其是当前最为流行的移动互联网设备上进行展示，极大限制了其应用范围。
全景图像是利用一个地点覆盖周边360度及头顶、脚底的多张照片拼接形成的，可模拟一个观察者在该位置看到的前后左右及上下的真实世界，全景照片因制作成本低、效果真实、可在浏览器和移动端展示等特点，在城市展示、旅游宣传等领域得到了广泛的应用。全景图像利用照片制作，具有场景真实、色彩真实、细节丰富等特点，这些特点正好补足了三维仿真技术的问题。利用无人机等辅助工具，全景图像也能从地面、空间等多视角对城市进行展示，可解决建设工程效果图视角单一的问题。但照片和全景地图目前针对的都是现实场景，与建设工程的模型进行融合难度较大，因此难以在建设工程方案设计、审查阶段开展应用。
可见，利用现实照片的优势，可以很好的补充三维仿真技术中存在的问题，各自存在优缺点，将两者结合成为必须趋势。当前，将虚拟世界与现实世界进行整合的技术，一般叫做现实增强(augmentedreality，简称ar)，可以实现将各种虚拟物体叠加到视频、照片中。由于通过摄像头、照片等获取的图像，图像中的各类物体不具备空间位置，因此，ar仅能做到将图像进行覆盖叠加，无法体现前、后关系，一般用于游戏、场景展示等工作中，难以应用到建设工程审查等需要进行精确空间关系表达的工程建设领域。也有部分效果图通过将建设工程的效果与现实照片进行整合，通过photoshop等图像处理软件进行人工合成，但一方面只能从图像中手工提取建设工程轮廓，容易造成轮廓部分缺失(过度剔除)、产生多余的非建设工程部分(过度保留)或产生多余空洞(误删除)等情况，工作量大，且影响最终合成效果；另一方面，由于没有空间定位，只能依靠人工判断现实中的物体与待建工程的空间位置关系，不仅工作量大，容易存在前后关系错乱的问题，并且现实物体和建设工程在体量(长、宽、高)上的关系也与实际情况不一致，
因此应用效果较差。
技术实现要素：
有鉴于现有技术的上述缺陷，本发明所要解决的技术问题是提供一种效果更好的三维模型与照片融合方法。
为实现上述目的，本发明提供了一种三维模型与照片融合方法，其特征是按以下步骤进行：
步骤s1、记录建筑环境照片的相机参数；
步骤s2、制作建筑工程三维模型，并加载到三维场景中；
步骤s3、设置三维场景中的虚拟相机参数，输出三维模型图像；
步骤s4、混合三维模型图像与照片，得到混合图像。
较佳的，所述步骤s1按以下步骤进行：
步骤s11、选择m个照片拍摄地点；所述m为自然数且m≥1；
步骤s12、在各拍摄地点获取n个预设方位
步骤s13、对拍摄的照片，记录每张照片的相机参数cameraparam{x，y，z，ori，pitch，roll，f}
所述x、y、z表示相机坐标；ori为相机方位，且0≤ori<360；pitch为相机俯仰角，且-90≤pitch≤90；roll为相机翻滚角，且-90≤roll≤90；f为镜头焦距，为大于0的自然数。
较佳的，所述步骤s2按以下步骤进行：
步骤s21、制作建筑三维模型；
步骤s22、加载所述建筑三维模型。
较佳的，所示步骤s3包括：
s31、依次对每个照片拍摄地点m，选择一张建设工程待建位置居中的照片pc；
s32、利用照片pc对应的相机参数，设置三维场景中虚拟相机的参数，包括位置x、y、z、相机方位ori、相机俯仰角pitch、相机翻滚角roll、相机焦距f；设置三维场景出图参数，使图像长度、宽度与照片pc尺寸一致；
s33、输出建设工程三维模型图像imgt；
s34、逐像素对图像imgt的颜色进行统计，找出图像imgt不包含的颜色cn：
s341、构建一个数组colorused{c0,c1,……,cn}，用于保存获取每一像素的rgb颜色值，其中，c为自然数，n为图像的像素数；
s342、依次获取图像的每个像素i的rgb值{ri,gi,bi}，计算其对应的自然数颜色ci，存入colorused中；
ci＝ri+gi*256+bi*65536；
所述ri、gi、bi为8位rgb颜色的红、绿、蓝分量，0≤ri≤255,0≤gi≤255，0≤bi≤255，0≤ci≤16777215；
s343、对colorused中的自然数，按从小到大进行排序；
s344、从0开始，寻找数组colorused中第一个不存在的自然数m，0≤m≤16777215，并计算图像不包含的颜色cn的rgb值cn{rn,gn,bn}，
s35、设置三维场景的背景颜色为cn，再次输出建设工程三维模型的图像imgc；
s36、选择与照片pc位置上相邻的照片pa，依照步骤s32至s35，利用照片pa对应的相机参数设置相机位置和出图大小后，输出建设工程三维模型图像imga。
较佳的，如果数组colorused包含了所有m的可能值，依次统计colorused中每个值重复的次数，选取重复最少的一个值作为不使用的颜色m，再计算cn；
较佳的，所述步骤s4按以下步骤执行：
s41、加载图像imgc，提取图像中颜色值不为cn的区域，该区域即为建设工程图像的轮廓sketch；
s42、利用轮廓sektch提取pc中的部分图像pcp；
s43、提取图像pcp中位于建设工程轮廓sketch前方的部分pcp_front：
s431、依照步骤s41、s42，从图像imga中提取出建设工程图像轮廓sketch_a，利用sketch_a从照片pa中提取出部分图像pap；
s432、针对图像imgc、imga利用shif特征点匹配算法，进行特征点匹配，寻找两张照片中的同名点，得到同名点对fti{(tc1,ta1),(tc2,ta2),……,(tcn,tan)}；其中，(tcn,tan)为第n对同名点，表示为(tcn(xcn,ycn)，tan(xan,yan))，x、y为同名点在图像中的坐标；同时按相同方法计算pcp、pap中的同名点对ftp{(tp1,tpa1),(tp2,tpa2),……,(tpn,tpan)}；
s433、根据imgc、imga特征点匹配结果fti，计算建设工程模型中最小同名点偏移量mvmin，计算
s434、根据pcp、pap特征点匹配结果ftp，计算图像pcp中各同名点的偏移量，并在图像pcp上标注各同名点偏移量；
s435、根据pcp上各同名点偏移量，针对图像上同名点偏移量小于mvmin的物体，将该物体轮廓提取出来，将所有提取的部分图像保存为pcp_front；
s44、按从下到上的次序依次叠加pc、sketch和pcp_front，并合并得到混合后的照片pf。
本发明一方面利用现实世界的照片来反映待建项目周边的真实环境，细节丰富、真实，尤其是建筑、景观的形态和颜色，与现实世界保持一致。另一方面，利用三维模型输出的图像来反映待建项目的形态和色彩。最后，通过图像融合，使模型能与真实世界的环境融为一体，并从多个视角真实反映项目建成后在现实世界中的真实效果。
较佳的，本发明还包括混合三维模型图像与照片得到混合图像的步骤：
s51、依次对每个照片拍摄地点m，利用混合后的照片pf替换原照片pc；
s52、加载所有照片，生成全景图像；
s53、发布全景图像。
本发明的有益效果是：本发明通过查找建筑模型图像不使用的颜色作为背景色，可以自动、完全、精确提取出建设工程的形状轮廓；通过使用与真实相机一致的虚拟相机参数输出建设工程模型图像，使现实照片与建设工程模型在图像中的体量关系正确、透视一致；基于特征点匹配技术，比较相邻两张照片中同名点的位移关系，可以较好的判断现实照片中各物体与建设工程的相对关系，解决手工操作空间关系不准确问题；利用全景图像技术，不仅可以从地面上某个视点对建设工程进行观察，还可以从空中等多个角度，全方面反映建设工程效果；最后，成果还可以全面应用在互联网、移动端等展示渠道中。本发明通过三维模型和现实照片的整合，能够开创性的应用在城市重大项目规划方案审查、建设工程未来形象展示、城乡规划和管理、地产商业宣传等领域。
附图说明
图1是本发明一具体实施方式的流程示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明：
如图1所示，一种三维模型与照片融合方法，按以下步骤进行：
步骤s1、现场拍摄建筑周边环境，并记录照片相机参数；
步骤s2、制作建筑工程三维模型，并加载到三维场景中；
步骤s3、设置三维场景中的虚拟相机参数，输出三维模型图像；
步骤s4、混合三维模型图像与照片，得到混合图像；
步骤s5、利用混合图像与其他照片，生产全景图像。
本实施例中，所述步骤s1按以下步骤进行：
步骤s11、根据建设工程效果展示需要，选择m个照片拍摄地点；所述m为自然数且m≥1；
为更好展示建设工程效果，所述拍摄地点，有可能在地面，也有可能在空中。一般情况下，地面拍摄地点应选择具有代表性的位置，如主干道、城市广场、观景台等。无论从地面或空中拍摄，拍摄地点应根据建设工程的占地范围、高度等，选择适当的拍摄距离或相机镜头焦距，以保证后续融合后，建设工程模型完全包含在照片中。
步骤s12、在各拍摄地点获取n个预设方位
步骤s13、对拍摄的照片，记录每张照片的相机参数cameraparam；
cameraparam{x，y，z，ori，pitch，roll，f}
所述x、y、z表示相机坐标；ori为相机方位，且0≤ori<360；pitch为相机俯仰角，且-90≤pitch≤90；roll为相机翻滚角，且-90≤roll≤90；f为镜头焦距，为大于0的自然数。
本实施例中，所述步骤s2按以下步骤进行：
步骤s21、通过3dsmax制作建筑三维模型，当然在其他实施例中也可以采用其它三维模型制作软件；
步骤s22、利用三维模型浏览软件，加载所述建筑三维模型。
本实施例中，所示步骤s3包括：
s31、依次对每个照片拍摄地点m，选择一张建设工程待建位置居中的照片pc；
s32、利用照片pc对应的相机参数，将三维场景中虚拟相机的参数设置成与照片一致，虚拟相机的参数包括位置x、y、z、相机方位ori、相机俯仰角pitch、相机翻滚角roll、相机焦距f；设置三维场景出图参数，使其图像长度、宽度与照片pc尺寸一致；
s33、输出1张建设工程三维模型图像imgt；
s34、逐像素对图像imgt的颜色进行统计，找出图像imgt不包含的颜色cn：
s341、构建一个数组colorused{c0,c1,……,cn}，用于保存获取每一像素的rgb颜色值，其中，c为自然数，n为图像的像素数；
s342、依次获取图像的每个像素i的rgb值{ri,gi,bi}，计算其对应的自然数颜色ci，存入colorused中；
ci＝ri+gi*256+bi*65536；
所述ri、gi、bi为8位rgb颜色的红、绿、蓝分量，0≤ri≤255,0≤gi≤255，0≤bi≤255，0≤ci≤16777215；
s343、对colorused中的自然数，按从小到大进行排序；
s344、从0开始，寻找数组colorused中第一个不存在的自然数m，0≤m≤16777215，并计算图像不包含的颜色cn的rgb值cn{rn,gn,bn}，
所述(int)表示对数值向下取整；
s345、如果数组colorused包含了所有m的可能值，依次统计colorused中每个值重复的次数，选取重复最少的一个值作为不使用的颜色m，再依据上述公式计算cn；
s35、设置三维场景的背景颜色为cn，再次输出1张建设工程三维模型的图像imgc；
s36、选择与照片pc位置上相邻的照片pa，pa与pc为不同位置，依照步骤s32至s35，利用照片pa对应的相机参数设置相机位置和出图大小后，输出一张建设工程三维模型图像imga。
本实施例中，所述步骤s4按以下步骤执行：
s41、在图像处理软件中，加载图像imgc，提取图像中颜色值不为cn的区域，该区域即为建设工程图像的轮廓sketch；
s42、利用轮廓sektch提取pc中的部分图像pcp；
s43、提取图像pcp中位于建设工程轮廓sketch前方的部分pcp_front：
s431、依照步骤s41、s42，从图像imga中提取出建设工程图像轮廓sketch_a，利用sketch_a从照片pa中提取出部分图像pap；
s432、针对图像imgc、imga利用shif特征点匹配算法，进行特征点匹配，寻找两张照片中的同名点(即shif的匹配点)，得到同名点对fti{(tc1,ta1),(tc2,ta2),……,(tcn,tan)}；其中，(tcn,tan)为第n对同名点，表示为(tcn(xcn,ycn)，tan(xan,yan))，x、y为同名点在图像中的坐标；同时按相同方法计算pcp、pap中的同名点对ftp{(tp1,tpa1),(tp2,tpa2),……,(tpn,tpan)}；
s433、根据imgc、imga特征点匹配结果fti，计算建设工程模型中最小同名点偏移量mvmin，计算
s434、根据pcp、pap特征点匹配结果ftp，计算图像pcp中各同名点的偏移量，并在图像pcp上标注各同名点偏移量；
s435、根据pcp上各同名点偏移量，针对图像上同名点偏移量小于mvmin的物体，利用图像处理软件的选择工具，将该物体轮廓提取出来，将所有提取的部分图像保存为pcp_front。在其他实施例中，为了保证需要保留在前景的物体都能找到同名点，可以适当放宽shif特征点匹配工程同的匹配精度，以得到更多同名点。
s44、按从下到上的次序依次叠加pc、sketch和pcp_front，并合并得到混合后的照片pf。
本实施例中，还包括混合三维模型图像与照片得到混合图像的步骤：
s51、依次对每个照片拍摄地点m，利用混合后的照片pf替换原照片pc；
s52、利用全景图像生成软件，加载所有照片，生成全景图像；
s53、利用现有的全景图像发布系统发布全景图像，或利用webgl技术开发全景发布系统进行全景图像发布。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。