c11308e85106b9ea84f44e04f2e93958b8f7f4b8
[blender-staging.git] / source / blender / blenlib / intern / uvproject.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * ***** END GPL LICENSE BLOCK *****
19  */
20
21 /** \file blender/blenlib/intern/uvproject.c
22  *  \ingroup bli
23  */
24
25
26 #include <math.h>
27
28 #include "MEM_guardedalloc.h"
29
30 #include "DNA_camera_types.h"
31 #include "DNA_object_types.h"
32
33 #include "BLI_math.h"
34 #include "BLI_uvproject.h"
35
36 typedef struct UvCameraInfo {
37         float camangle;
38         float camsize;
39         float xasp, yasp;
40         float shiftx, shifty;
41         float rotmat[4][4];
42         float caminv[4][4];
43         short do_persp, do_pano, do_rotmat;
44 } UvCameraInfo;
45
46 void project_from_camera(float target[2], float source[3], UvCameraInfo *uci)
47 {
48         float pv4[4];
49
50         copy_v3_v3(pv4, source);
51         pv4[3]= 1.0;
52
53         /* rotmat is the object matrix in this case */
54         if(uci->do_rotmat)
55                 mul_m4_v4(uci->rotmat, pv4);
56
57         /* caminv is the inverse camera matrix */
58         mul_m4_v4(uci->caminv, pv4);
59
60         if(uci->do_pano) {
61                 float angle= atan2f(pv4[0], -pv4[2]) / ((float)M_PI * 2.0f); /* angle around the camera */
62                 if (uci->do_persp==0) {
63                         target[0]= angle; /* no correct method here, just map to  0-1 */
64                         target[1]= pv4[1] / uci->camsize;
65                 }
66                 else {
67                         float vec2d[2]; /* 2D position from the camera */
68                         vec2d[0]= pv4[0];
69                         vec2d[1]= pv4[2];
70                         target[0]= angle * ((float)M_PI / uci->camangle);
71                         target[1]= pv4[1] / (len_v2(vec2d) * (uci->camsize * 2.0f));
72                 }
73         }
74         else {
75                 if (pv4[2]==0.0f) pv4[2]= 0.00001f; /* don't allow div by 0 */
76
77                 if (uci->do_persp==0) {
78                         target[0]= (pv4[0]/uci->camsize);
79                         target[1]= (pv4[1]/uci->camsize);
80                 }
81                 else {
82                         target[0]= (-pv4[0]*((1.0f/uci->camsize)/pv4[2])) / 2.0f;
83                         target[1]= (-pv4[1]*((1.0f/uci->camsize)/pv4[2])) / 2.0f;
84                 }
85         }
86
87         target[0] *= uci->xasp;
88         target[1] *= uci->yasp;
89         
90         /* adds camera shift + 0.5 */
91         target[0] += uci->shiftx;
92         target[1] += uci->shifty;
93 }
94
95 /* could rv3d->persmat */
96 void project_from_view(float target[2], float source[3], float persmat[4][4], float rotmat[4][4], float winx, float winy)
97 {
98         float pv[3], pv4[4], x= 0.0, y= 0.0;
99
100         mul_v3_m4v3(pv, rotmat, source);
101
102         copy_v3_v3(pv4, source);
103         pv4[3]= 1.0;
104
105         /* rotmat is the object matrix in this case */
106         mul_m4_v4(rotmat, pv4); 
107
108         /* almost project_short */
109         mul_m4_v4(persmat, pv4);
110         if(fabsf(pv4[3]) > 0.00001f) { /* avoid division by zero */
111                 target[0] = winx/2.0f + (winx/2.0f) * pv4[0] / pv4[3];
112                 target[1] = winy/2.0f + (winy/2.0f) * pv4[1] / pv4[3];
113         }
114         else {
115                 /* scaling is lost but give a valid result */
116                 target[0] = winx/2.0f + (winx/2.0f) * pv4[0];
117                 target[1] = winy/2.0f + (winy/2.0f) * pv4[1];
118         }
119
120         /* v3d->persmat seems to do this funky scaling */ 
121         if(winx > winy) {
122                 y= (winx - winy)/2.0f;
123                 winy = winx;
124         }
125         else {
126                 x= (winy - winx)/2.0f;
127                 winx = winy;
128         }
129
130         target[0]= (x + target[0]) / winx;
131         target[1]= (y + target[1]) / winy;
132 }
133
134 /* 'rotmat' can be obedit->obmat when uv project is used.
135  * 'winx' and 'winy' can be from scene->r.xsch/ysch */ 
136 UvCameraInfo *project_camera_info(Object *ob, float (*rotmat)[4], float winx, float winy)
137 {
138         UvCameraInfo uci;
139         Camera *camera= ob->data;
140
141         uci.do_pano = (camera->flag & CAM_PANORAMA);
142         uci.do_persp = (camera->type==CAM_PERSP);
143
144         uci.camangle= focallength_to_hfov(camera->lens, camera->sensor_x) / 2.0f;
145         uci.camsize= uci.do_persp ? tanf(uci.camangle) : camera->ortho_scale;
146
147         /* account for scaled cameras */
148         copy_m4_m4(uci.caminv, ob->obmat);
149         normalize_m4(uci.caminv);
150
151         if (invert_m4(uci.caminv)) {
152                 UvCameraInfo *uci_pt;
153
154                 /* normal projection */
155                 if(rotmat) {
156                         copy_m4_m4(uci.rotmat, rotmat);
157                         uci.do_rotmat= 1;
158                 }
159                 else {
160                         uci.do_rotmat= 0;
161                 }
162
163                 /* also make aspect ratio adjustment factors */
164                 if (winx > winy) {
165                         uci.xasp= 1.0f;
166                         uci.yasp= winx / winy;
167                 }
168                 else {
169                         uci.xasp= winy / winx;
170                         uci.yasp= 1.0f;
171                 }
172                 
173                 /* include 0.5f here to move the UVs into the center */
174                 uci.shiftx = 0.5f - (camera->shiftx * uci.xasp);
175                 uci.shifty = 0.5f - (camera->shifty * uci.yasp);
176                 
177                 uci_pt= MEM_mallocN(sizeof(UvCameraInfo), "UvCameraInfo");
178                 *uci_pt= uci;
179                 return uci_pt;
180         }
181
182         return NULL;
183 }
184
185 void project_from_view_ortho(float target[2], float source[3], float rotmat[4][4])
186 {
187         float pv[3];
188
189         mul_v3_m4v3(pv, rotmat, source);
190
191         /* ortho projection */
192         target[0] = -pv[0];
193         target[1] = pv[2];
194 }
195
196
197 void project_camera_info_scale(UvCameraInfo *uci, float scale_x, float scale_y)
198 {
199         uci->xasp *= scale_x;
200         uci->yasp *= scale_y;
201 }