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
100
101
102
103
104
|
/*
* rcoords.c
*/
#include "csf.h"
#include "csfimpl.h"
/* compute true world co-ordinate of a pixel
* RrowCol2Coords computes the true world co-ordinate from a
* row, column index.
* The row, column co-ordinate
* don't have to be on the map. They are just relative to upper left position.
* For example (row,col) = (-1,0) computes the (x,y) co-ordinate of
* the pixel that is right above upper left pixel.
*
* returns
* 0 if the co-ordinate is outside the map.
* 1 if inside.
* -1 in case of an error.
*
* Merrno
* ILL_CELLSIZE
*/
int RgetCoords(
const MAP *m, /* map handle */
int inCellPos, /* nonzero if you want the co-ordinate
* at the centre of the cell, 0 if you
* want the upper left co-ordinate of the cell
*/
size_t row, /* Row number (relates to y position). */
size_t col, /* Column number (relates to x position). */
double *x, /* write-only. Returns x of true co-ordinate */
double *y) /* write-only. Returns y of true co-ordinate */
{
return RrowCol2Coords(m,
(double)row+(inCellPos ? 0.5 : 0.0),
(double)col+(inCellPos ? 0.5 : 0.0),
x,y);
}
/* compute true world co-ordinate from row, column index
* RrowCol2Coords computes the true world co-ordinate from a
* row, column index. The row,column co-ordinate can be fractions.
* For example (row,col) = (0.5,0.5) computes the (x,y) co-ordinate of
* the centre of the upper left pixel. Secondly, the row and column co-ordinate
* don't have to be on the map. They are just relative to upper left position.
* For example (row,col) = (-0.5,0.5) computes the (x,y) co-ordinate of
* the centre of the pixel that is right above upper left pixel.
*/
void RasterRowCol2Coords(
const CSF_RASTER_LOCATION_ATTRIBUTES *m, /* raster handle */
double row, /* Row number (relates to y position). */
double col, /* Column number (relates to x position). */
double *x, /* write-only. x co-ordinate */
double *y) /* write-only. y co-ordinate */
{
double cs = m->cellSize;
double c = m->angleCos;
double s = m->angleSin;
double yRow = cs * row;
double xCol = cs * col;
double xCol_t = xCol * c - yRow * s;
double yRow_t = xCol * s + yRow * c;
*x = m->xUL + xCol_t;
if (m->projection == PT_YINCT2B)
*y = m->yUL + yRow_t;
else /* all other projections */
*y = m->yUL - yRow_t;
}
/* compute true world co-ordinate from row, column index
* RasterRowCol2Coords computes the true world co-ordinate from a
* row, column index. The row,column co-ordinate can be fractions.
* For example (row,col) = (0.5,0.5) computes the (x,y) co-ordinate of
* the centre of the upper left pixel. Secondly, the row and column co-ordinate
* don't have to be on the map. They are just relative to upper left position.
* For example (row,col) = (-0.5,0.5) computes the (x,y) co-ordinate of
* the centre of the pixel that is right above upper left pixel.
*
* returns
* 0 if the co-ordinate is outside the map.
* 1 if inside.
* -1 in case of an error.
*
* Merrno
* ILL_CELLSIZE
*/
int RrowCol2Coords(const MAP *m, /* map handle */
double row, /* Row number (relates to y position). */
double col, /* Column number (relates to x position). */
double *x, /* write-only. x co-ordinate */
double *y) /* write-only. y co-ordinate */
{
if (m->raster.cellSize <= 0
|| m->raster.cellSize != m->raster.cellSizeDupl )
{
M_ERROR(ILL_CELLSIZE);
goto error;
}
RasterRowCol2Coords(&(m->raster),row,col,x,y);
return( (m->raster.nrRows > row) && (m->raster.nrCols > col) &&
(row >= 0) && (col >= 0));
error: return(-1);
}
|
