226 |
|
} |
227 |
|
|
228 |
|
|
229 |
< |
void |
229 |
> |
int |
230 |
|
viewloc( /* find image location for point */ |
231 |
|
FVECT ip, |
232 |
|
VIEW *v, |
249 |
|
ip[2] = -ip[2]; |
250 |
|
d = -d; |
251 |
|
} |
252 |
< |
if (d > FTINY) { |
253 |
< |
d = 1.0/d; |
254 |
< |
disp[0] *= d; |
255 |
< |
disp[1] *= d; |
256 |
< |
disp[2] *= d; |
257 |
< |
} |
252 |
> |
if (d <= FTINY) |
253 |
> |
return(0); /* at infinite edge */ |
254 |
> |
d = 1.0/d; |
255 |
> |
disp[0] *= d; |
256 |
> |
disp[1] *= d; |
257 |
> |
disp[2] *= d; |
258 |
|
ip[2] *= (1.0 - v->vfore*d); |
259 |
|
break; |
260 |
|
case VT_HEM: /* hemispherical fisheye */ |
269 |
|
d = DOT(disp,v->hvec); |
270 |
|
d2 = DOT(disp,v->vdir); |
271 |
|
ip[0] = 180.0/PI * atan2(d,d2) / v->horiz + 0.5 - v->hoff; |
272 |
< |
d = 1.0/sqrt(d*d + d2*d2); |
272 |
> |
d = d*d + d2*d2; |
273 |
> |
if (d <= FTINY*FTINY) |
274 |
> |
return(0); /* at pole */ |
275 |
> |
d = 1.0/sqrt(d); |
276 |
|
ip[1] = DOT(disp,v->vvec)*d/v->vn2 + 0.5 - v->voff; |
277 |
|
ip[2] = VLEN(disp); |
278 |
|
ip[2] *= (1.0 - v->vfore*d); |
279 |
< |
return; |
279 |
> |
return(1); |
280 |
|
case VT_ANG: /* angular fisheye */ |
281 |
|
ip[0] = 0.5 - v->hoff; |
282 |
|
ip[1] = 0.5 - v->voff; |
283 |
|
ip[2] = normalize(disp) - v->vfore; |
284 |
|
d = DOT(disp,v->vdir); |
285 |
|
if (d >= 1.0-FTINY) |
286 |
< |
return; |
286 |
> |
return(1); |
287 |
|
if (d <= -(1.0-FTINY)) { |
288 |
|
ip[0] += 180.0/v->horiz; |
289 |
< |
return; |
289 |
> |
return(1); |
290 |
|
} |
291 |
|
d = (180.0/PI)*acos(d) / sqrt(1.0 - d*d); |
292 |
|
ip[0] += DOT(disp,v->hvec)*d/v->horiz; |
293 |
|
ip[1] += DOT(disp,v->vvec)*d/v->vert; |
294 |
< |
return; |
294 |
> |
return(1); |
295 |
|
case VT_PLS: /* planispheric fisheye */ |
296 |
|
ip[0] = 0.5 - v->hoff; |
297 |
|
ip[1] = 0.5 - v->voff; |
298 |
|
ip[2] = normalize(disp) - v->vfore; |
299 |
|
d = DOT(disp,v->vdir); |
300 |
|
if (d >= 1.0-FTINY) |
301 |
< |
return; |
301 |
> |
return(1); |
302 |
|
if (d <= -(1.0-FTINY)) |
303 |
< |
return; /* really an error */ |
303 |
> |
return(0); |
304 |
|
ip[0] += DOT(disp,v->hvec)/((1. + d)*sqrt(v->hn2)); |
305 |
|
ip[1] += DOT(disp,v->vvec)/((1. + d)*sqrt(v->vn2)); |
306 |
< |
return; |
306 |
> |
return(1); |
307 |
|
} |
308 |
|
ip[0] = DOT(disp,v->hvec)/v->hn2 + 0.5 - v->hoff; |
309 |
|
ip[1] = DOT(disp,v->vvec)/v->vn2 + 0.5 - v->voff; |
310 |
+ |
return(1); |
311 |
|
} |
312 |
|
|
313 |
|
|