1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
/*
****************************************************************************
*
* MODULE: v.transform
* AUTHOR(S): See other files as well...
* Eric G. Miller <egm2@jps.net>
* PURPOSE: To transform a vector layer's coordinates via a set of tie
* points.
* COPYRIGHT: (C) 2002 by the GRASS Development Team
*
* This program is free software under the GNU General Public
* License (>=v2). Read the file COPYING that comes with GRASS
* for details.
*
*****************************************************************************/
/*
* Functions in this file:
* setup_transform ( n_points)
* returns: ALL_OK, POINTS_NOT_SPREAD, NEED_MORE_POINTS
* init_transform_arrays()
* print_transform_error(stat)
*
* Written by the GRASS Team, 02/16/90, -mh .
*/
#include <stdlib.h>
#include <stdio.h>
#include "glocale.h"
#include "gis.h"
#include "trans.h"
int
setup_transform (int n_points)
{
int status ;
/* compute_transformation_coef() returns:
* -2, Not enough points
* 1, everything is okay
* -1, points weren't spread out enough
*/
/* if there are enough points registered compute residuals */
if(n_points >= MIN_COOR)
status = compute_transformation_coef (ax, ay, bx, by, use, MAX_COOR) ;
else
return(NEED_MORE_POINTS) ;
if(status != ALL_OK)
return(POINTS_NOT_SPREAD) ;
residuals_a_predicts_b (ax, ay, bx, by, use, MAX_COOR, residuals, &rms) ;
return (ALL_OK) ;
} /* setup_transform() */
int
init_transform_arrays (void)
{
int i ;
/* initiliaze use[], no points valid */
for (i=0 ; i<MAX_COOR ; ++i)
{
*(use+i) = 0 ; *(bx+i) = 0 ; *(by+i) = 0 ;
*(residuals+i) = 0 ;
}
reg_cnt = 0 ;
return 0;
}
int
print_transform_error (int stat)
{
char buff[128];
switch(stat)
{
case POINTS_NOT_SPREAD:
G_fatal_error ( _("The points weren't spread out enough.") );
break ;
case NEED_MORE_POINTS:
G_fatal_error ( _("You need to enter at least %d points."), MIN_COOR) ;
break ;
default:
G_fatal_error ( "Your calling print_transform_error() with no error." );
break ;
}
return 0;
}
|
