Skip to content
Projects
Groups
Snippets
Help
Loading...
Help
Support
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
M
MariaDB
Project overview
Project overview
Details
Activity
Releases
Repository
Repository
Files
Commits
Branches
Tags
Contributors
Graph
Compare
Issues
0
Issues
0
List
Boards
Labels
Milestones
Merge Requests
0
Merge Requests
0
CI / CD
CI / CD
Pipelines
Jobs
Schedules
Analytics
Analytics
CI / CD
Repository
Value Stream
Wiki
Wiki
Snippets
Snippets
Members
Members
Collapse sidebar
Close sidebar
Activity
Graph
Create a new issue
Jobs
Commits
Issue Boards
Open sidebar
nexedi
MariaDB
Commits
5150bd8f
Commit
5150bd8f
authored
Oct 31, 2005
by
hf@deer.(none)
Browse files
Options
Browse Files
Download
Plain Diff
Merge bk@192.168.21.1:/usr/home/bk/mysql-4.1
into deer.(none):/home/hf/work/mysql-4.1.12839
parents
b919b29d
16c66d49
Changes
6
Hide whitespace changes
Inline
Side-by-side
Showing
6 changed files
with
353 additions
and
79 deletions
+353
-79
mysql-test/r/gis.result
mysql-test/r/gis.result
+6
-0
mysql-test/t/gis.test
mysql-test/t/gis.test
+3
-0
sql/field.cc
sql/field.cc
+3
-5
sql/item_geofunc.cc
sql/item_geofunc.cc
+28
-49
sql/spatial.cc
sql/spatial.cc
+297
-22
sql/spatial.h
sql/spatial.h
+16
-3
No files found.
mysql-test/r/gis.result
View file @
5150bd8f
...
...
@@ -655,3 +655,9 @@ t1 where object_id=85984;
object_id geometrytype(geo) ISSIMPLE(GEO) ASTEXT(centroid(geo))
85984 MULTIPOLYGON 0 POINT(-114.87787186923 36.33101763469)
drop table t1;
select (asWKT(geomfromwkb((0x000000000140240000000000004024000000000000))));
(asWKT(geomfromwkb((0x000000000140240000000000004024000000000000))))
POINT(10 10)
select (asWKT(geomfromwkb((0x010100000000000000000024400000000000002440))));
(asWKT(geomfromwkb((0x010100000000000000000024400000000000002440))))
POINT(10 10)
mysql-test/t/gis.test
View file @
5150bd8f
...
...
@@ -360,4 +360,7 @@ t1 where object_id=85984;
drop
table
t1
;
select
(
asWKT
(
geomfromwkb
((
0x000000000140240000000000004024000000000000
))));
select
(
asWKT
(
geomfromwkb
((
0x010100000000000000000024400000000000002440
))));
# End of 4.1 tests
sql/field.cc
View file @
5150bd8f
...
...
@@ -5745,8 +5745,7 @@ void Field_blob::get_key_image(char *buff,uint length,
return
;
}
get_ptr
(
&
blob
);
gobj
=
Geometry
::
create_from_wkb
(
&
buffer
,
blob
+
SRID_SIZE
,
blob_length
-
SRID_SIZE
);
gobj
=
Geometry
::
construct
(
&
buffer
,
blob
,
blob_length
);
if
(
gobj
->
get_mbr
(
&
mbr
,
&
dummy
))
bzero
(
buff
,
SIZEOF_STORED_DOUBLE
*
4
);
else
...
...
@@ -6039,8 +6038,7 @@ void Field_geom::get_key_image(char *buff, uint length, CHARSET_INFO *cs,
return
;
}
get_ptr
(
&
blob
);
gobj
=
Geometry
::
create_from_wkb
(
&
buffer
,
blob
+
SRID_SIZE
,
blob_length
-
SRID_SIZE
);
gobj
=
Geometry
::
construct
(
&
buffer
,
blob
,
blob_length
);
if
(
gobj
->
get_mbr
(
&
mbr
,
&
dummy
))
bzero
(
buff
,
SIZEOF_STORED_DOUBLE
*
4
);
else
...
...
@@ -6100,7 +6098,7 @@ int Field_geom::store(const char *from, uint length, CHARSET_INFO *cs)
uint32
wkb_type
;
if
(
length
<
SRID_SIZE
+
WKB_HEADER_SIZE
+
SIZEOF_STORED_DOUBLE
*
2
)
goto
err
;
wkb_type
=
uint4korr
(
from
+
WKB_HEADER_SIZE
);
wkb_type
=
uint4korr
(
from
+
SRID_SIZE
+
1
);
if
(
wkb_type
<
(
uint32
)
Geometry
::
wkb_point
||
wkb_type
>
(
uint32
)
Geometry
::
wkb_end
)
return
-
1
;
...
...
sql/item_geofunc.cc
View file @
5150bd8f
...
...
@@ -78,8 +78,7 @@ String *Item_func_geometry_from_wkb::val_str(String *str)
str
->
q_append
(
srid
);
if
((
null_value
=
(
args
[
0
]
->
null_value
||
!
Geometry
::
create_from_wkb
(
&
buffer
,
wkb
->
ptr
(),
wkb
->
length
())
||
str
->
append
(
*
wkb
))))
!
Geometry
::
create_from_wkb
(
&
buffer
,
wkb
->
ptr
(),
wkb
->
length
(),
str
))))
return
0
;
return
str
;
}
...
...
@@ -96,8 +95,7 @@ String *Item_func_as_wkt::val_str(String *str)
if
((
null_value
=
(
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
)))))
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
())))))
return
0
;
str
->
length
(
0
);
...
...
@@ -123,8 +121,7 @@ String *Item_func_as_wkb::val_str(String *str)
if
((
null_value
=
(
args
[
0
]
->
null_value
||
!
(
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
)))))
!
(
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
())))))
return
0
;
str
->
copy
(
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
,
...
...
@@ -142,8 +139,7 @@ String *Item_func_geometry_type::val_str(String *str)
if
((
null_value
=
(
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
)))))
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
())))))
return
0
;
/* String will not move */
str
->
copy
(
geom
->
get_class_info
()
->
m_name
.
str
,
...
...
@@ -164,8 +160,7 @@ String *Item_func_envelope::val_str(String *str)
if
((
null_value
=
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))))
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))))
return
0
;
srid
=
uint4korr
(
swkb
->
ptr
());
...
...
@@ -188,8 +183,7 @@ String *Item_func_centroid::val_str(String *str)
uint32
srid
;
if
((
null_value
=
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))))
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))))
return
0
;
str
->
set_charset
(
&
my_charset_bin
);
...
...
@@ -218,8 +212,7 @@ String *Item_func_spatial_decomp::val_str(String *str)
if
((
null_value
=
(
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
)))))
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
())))))
return
0
;
srid
=
uint4korr
(
swkb
->
ptr
());
...
...
@@ -267,8 +260,7 @@ String *Item_func_spatial_decomp_n::val_str(String *str)
if
((
null_value
=
(
args
[
0
]
->
null_value
||
args
[
1
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
)))))
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
())))))
return
0
;
str
->
set_charset
(
&
my_charset_bin
);
...
...
@@ -475,10 +467,8 @@ longlong Item_func_spatial_rel::val_int()
if
((
null_value
=
(
args
[
0
]
->
null_value
||
args
[
1
]
->
null_value
||
!
(
g1
=
Geometry
::
create_from_wkb
(
&
buffer1
,
res1
->
ptr
()
+
SRID_SIZE
,
res1
->
length
()
-
SRID_SIZE
))
||
!
(
g2
=
Geometry
::
create_from_wkb
(
&
buffer2
,
res2
->
ptr
()
+
SRID_SIZE
,
res2
->
length
()
-
SRID_SIZE
))
||
!
(
g1
=
Geometry
::
construct
(
&
buffer1
,
res1
->
ptr
(),
res1
->
length
()))
||
!
(
g2
=
Geometry
::
construct
(
&
buffer2
,
res2
->
ptr
(),
res2
->
length
()))
||
g1
->
get_mbr
(
&
mbr1
,
&
dummy
)
||
g2
->
get_mbr
(
&
mbr2
,
&
dummy
))))
return
0
;
...
...
@@ -543,8 +533,7 @@ longlong Item_func_isclosed::val_int()
null_value
=
(
!
swkb
||
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
is_closed
(
&
isclosed
));
return
(
longlong
)
isclosed
;
...
...
@@ -566,9 +555,7 @@ longlong Item_func_dimension::val_int()
null_value
=
(
!
swkb
||
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
dimension
(
&
dim
,
&
dummy
));
return
(
longlong
)
dim
;
}
...
...
@@ -583,9 +570,8 @@ longlong Item_func_numinteriorring::val_int()
Geometry
*
geom
;
null_value
=
(
!
swkb
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
num_interior_ring
(
&
num
));
return
(
longlong
)
num
;
}
...
...
@@ -600,9 +586,8 @@ longlong Item_func_numgeometries::val_int()
Geometry
*
geom
;
null_value
=
(
!
swkb
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
num_geometries
(
&
num
));
return
(
longlong
)
num
;
}
...
...
@@ -618,9 +603,8 @@ longlong Item_func_numpoints::val_int()
null_value
=
(
!
swkb
||
args
[
0
]
->
null_value
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
num_points
(
&
num
));
return
(
longlong
)
num
;
}
...
...
@@ -635,9 +619,8 @@ double Item_func_x::val()
Geometry
*
geom
;
null_value
=
(
!
swkb
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
get_x
(
&
res
));
return
res
;
}
...
...
@@ -652,9 +635,8 @@ double Item_func_y::val()
Geometry
*
geom
;
null_value
=
(
!
swkb
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
get_y
(
&
res
));
return
res
;
}
...
...
@@ -670,9 +652,8 @@ double Item_func_area::val()
const
char
*
dummy
;
null_value
=
(
!
swkb
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
area
(
&
res
,
&
dummy
));
return
res
;
}
...
...
@@ -686,9 +667,8 @@ double Item_func_glength::val()
Geometry
*
geom
;
null_value
=
(
!
swkb
||
!
(
geom
=
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
))
||
!
(
geom
=
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()))
||
geom
->
length
(
&
res
));
return
res
;
}
...
...
@@ -700,9 +680,8 @@ longlong Item_func_srid::val_int()
Geometry_buffer
buffer
;
null_value
=
(
!
swkb
||
!
Geometry
::
create_from_wkb
(
&
buffer
,
swkb
->
ptr
()
+
SRID_SIZE
,
swkb
->
length
()
-
SRID_SIZE
));
!
Geometry
::
construct
(
&
buffer
,
swkb
->
ptr
(),
swkb
->
length
()));
if
(
null_value
)
return
0
;
...
...
sql/spatial.cc
View file @
5150bd8f
...
...
@@ -119,25 +119,23 @@ Geometry::Class_info *Geometry::find_class(const char *name, uint32 len)
return
0
;
}
Geometry
*
Geometry
::
create_from_wkb
(
Geometry_buffer
*
buffer
,
const
char
*
data
,
uint32
data_len
)
Geometry
*
Geometry
::
construct
(
Geometry_buffer
*
buffer
,
const
char
*
data
,
uint32
data_len
)
{
uint32
geom_type
;
Geometry
*
result
;
char
byte_order
;
if
(
data_len
<
1
+
4
)
if
(
data_len
<
SRID_SIZE
+
1
+
4
)
return
NULL
;
data
++
;
/*
FIXME: check byte ordering
Also check if we could replace this with one byte
*/
geom_type
=
uint4korr
(
data
);
data
+=
4
;
byte_order
=
data
[
SRID_SIZE
];
geom_type
=
uint4korr
(
data
+
SRID_SIZE
+
1
);
data
+=
SRID_SIZE
+
WKB_HEADER_SIZE
;
if
(
!
(
result
=
create_by_typeid
(
buffer
,
(
int
)
geom_type
)))
return
NULL
;
result
->
m_data
=
data
;
result
->
m_data_end
=
data
+
data_len
;
result
->
m_data_end
=
data
+
(
data_len
-
(
SRID_SIZE
+
WKB_HEADER_SIZE
))
;
return
result
;
}
...
...
@@ -168,13 +166,71 @@ Geometry *Geometry::create_from_wkt(Geometry_buffer *buffer,
return
NULL
;
if
(
init_stream
)
{
result
->
init_from_wkb
(
wkt
->
ptr
(),
wkt
->
length
());
result
->
set_data_ptr
(
wkt
->
ptr
(),
wkt
->
length
());
result
->
shift_wkb_header
();
}
return
result
;
}
static
double
wkb_get_double
(
const
char
*
ptr
,
Geometry
::
wkbByteOrder
bo
)
{
double
res
;
if
(
bo
!=
Geometry
::
wkb_xdr
)
float8get
(
res
,
ptr
);
else
{
char
inv_array
[
8
];
inv_array
[
0
]
=
ptr
[
7
];
inv_array
[
1
]
=
ptr
[
6
];
inv_array
[
2
]
=
ptr
[
5
];
inv_array
[
3
]
=
ptr
[
4
];
inv_array
[
4
]
=
ptr
[
3
];
inv_array
[
5
]
=
ptr
[
2
];
inv_array
[
6
]
=
ptr
[
1
];
inv_array
[
7
]
=
ptr
[
0
];
float8get
(
res
,
inv_array
);
}
return
res
;
}
static
uint32
wkb_get_uint
(
const
char
*
ptr
,
Geometry
::
wkbByteOrder
bo
)
{
if
(
bo
!=
Geometry
::
wkb_xdr
)
return
uint4korr
(
ptr
);
/* else */
{
char
inv_array
[
4
];
inv_array
[
0
]
=
ptr
[
3
];
inv_array
[
1
]
=
ptr
[
2
];
inv_array
[
2
]
=
ptr
[
1
];
inv_array
[
3
]
=
ptr
[
0
];
return
uint4korr
(
inv_array
);
}
}
int
Geometry
::
create_from_wkb
(
Geometry_buffer
*
buffer
,
const
char
*
wkb
,
uint32
len
,
String
*
res
)
{
uint32
geom_type
;
Geometry
*
geom
;
if
(
len
<
WKB_HEADER_SIZE
)
return
1
;
geom_type
=
wkb_get_uint
(
wkb
+
1
,
(
wkbByteOrder
)
wkb
[
0
]);
if
(
!
(
geom
=
create_by_typeid
(
buffer
,
(
int
)
geom_type
))
||
res
->
reserve
(
WKB_HEADER_SIZE
,
512
))
return
1
;
res
->
q_append
((
char
)
wkb_ndr
);
res
->
q_append
(
geom_type
);
return
geom
->
init_from_wkb
(
wkb
+
WKB_HEADER_SIZE
,
len
-
WKB_HEADER_SIZE
,
(
wkbByteOrder
)
wkb
[
0
],
res
);
}
bool
Geometry
::
envelope
(
String
*
result
)
const
{
MBR
mbr
;
...
...
@@ -344,6 +400,20 @@ bool Gis_point::init_from_wkt(Gis_read_stream *trs, String *wkb)
}
uint
Gis_point
::
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
{
double
x
,
y
;
if
(
len
<
POINT_DATA_SIZE
||
res
->
reserve
(
POINT_DATA_SIZE
))
return
0
;
x
=
wkb_get_double
(
wkb
,
bo
);
y
=
wkb_get_double
(
wkb
+
SIZEOF_STORED_DOUBLE
,
bo
);
res
->
q_append
(
x
);
res
->
q_append
(
y
);
return
POINT_DATA_SIZE
;
}
bool
Gis_point
::
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
{
double
x
,
y
;
...
...
@@ -413,6 +483,33 @@ bool Gis_line_string::init_from_wkt(Gis_read_stream *trs, String *wkb)
}
uint
Gis_line_string
::
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
{
uint32
n_points
,
proper_length
;
const
char
*
wkb_end
;
Gis_point
p
;
if
(
len
<
4
)
return
0
;
n_points
=
wkb_get_uint
(
wkb
,
bo
);
proper_length
=
4
+
n_points
*
POINT_DATA_SIZE
;
if
(
len
<
proper_length
||
res
->
reserve
(
proper_length
))
return
0
;
res
->
q_append
(
n_points
);
wkb_end
=
wkb
+
proper_length
;
for
(
wkb
+=
4
;
wkb
<
wkb_end
;
wkb
+=
POINT_DATA_SIZE
)
{
if
(
!
p
.
init_from_wkb
(
wkb
,
POINT_DATA_SIZE
,
bo
,
res
))
return
0
;
}
return
proper_length
;
}
bool
Gis_line_string
::
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
{
uint32
n_points
;
...
...
@@ -592,7 +689,7 @@ bool Gis_polygon::init_from_wkt(Gis_read_stream *trs, String *wkb)
trs
->
check_next_symbol
(
')'
))
return
1
;
ls
.
init_from_wkb
(
wkb
->
ptr
()
+
ls_pos
,
wkb
->
length
()
-
ls_pos
);
ls
.
set_data_ptr
(
wkb
->
ptr
()
+
ls_pos
,
wkb
->
length
()
-
ls_pos
);
if
(
ls
.
is_closed
(
&
closed
)
||
!
closed
)
{
trs
->
set_error_msg
(
"POLYGON's linear ring isn't closed"
);
...
...
@@ -607,6 +704,43 @@ bool Gis_polygon::init_from_wkt(Gis_read_stream *trs, String *wkb)
}
uint
Gis_polygon
::
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
{
uint32
n_linear_rings
;
const
char
*
wkb_orig
=
wkb
;
if
(
len
<
4
)
return
0
;
n_linear_rings
=
wkb_get_uint
(
wkb
,
bo
);
if
(
res
->
reserve
(
4
,
512
))
return
0
;
wkb
+=
4
;
len
-=
4
;
res
->
q_append
(
n_linear_rings
);
while
(
n_linear_rings
--
)
{
Gis_line_string
ls
;
uint32
ls_pos
=
res
->
length
();
int
ls_len
;
int
closed
;
if
(
!
(
ls_len
=
ls
.
init_from_wkb
(
wkb
,
len
,
bo
,
res
)))
return
0
;
ls
.
set_data_ptr
(
res
->
ptr
()
+
ls_pos
,
res
->
length
()
-
ls_pos
);
if
(
ls
.
is_closed
(
&
closed
)
||
!
closed
)
return
0
;
wkb
+=
ls_len
;
}
return
wkb
-
wkb_orig
;
}
bool
Gis_polygon
::
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
{
uint32
n_linear_rings
;
...
...
@@ -895,6 +1029,36 @@ bool Gis_multi_point::init_from_wkt(Gis_read_stream *trs, String *wkb)
}
uint
Gis_multi_point
::
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
{
uint32
n_points
;
uint
proper_size
;
Gis_point
p
;
const
char
*
wkb_end
;
if
(
len
<
4
)
return
0
;
n_points
=
wkb_get_uint
(
wkb
,
bo
);
proper_size
=
4
+
n_points
*
(
WKB_HEADER_SIZE
+
POINT_DATA_SIZE
);
if
(
len
<
proper_size
||
res
->
reserve
(
proper_size
))
return
0
;
res
->
q_append
(
n_points
);
wkb_end
=
wkb
+
proper_size
;
for
(
wkb
+=
4
;
wkb
<
wkb_end
;
wkb
+=
(
WKB_HEADER_SIZE
+
POINT_DATA_SIZE
))
{
res
->
q_append
((
char
)
wkb_ndr
);
res
->
q_append
((
uint32
)
wkb_point
);
if
(
!
p
.
init_from_wkb
(
wkb
+
WKB_HEADER_SIZE
,
POINT_DATA_SIZE
,
(
wkbByteOrder
)
wkb
[
0
],
res
))
return
0
;
}
return
proper_size
;
}
bool
Gis_multi_point
::
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
{
uint32
n_points
;
...
...
@@ -1004,6 +1168,42 @@ bool Gis_multi_line_string::init_from_wkt(Gis_read_stream *trs, String *wkb)
}
uint
Gis_multi_line_string
::
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
{
uint32
n_line_strings
;
const
char
*
wkb_orig
=
wkb
;
if
(
len
<
4
)
return
0
;
n_line_strings
=
wkb_get_uint
(
wkb
,
bo
);
if
(
res
->
reserve
(
4
,
512
))
return
0
;
res
->
q_append
(
n_line_strings
);
for
(
wkb
+=
4
;
n_line_strings
;
n_line_strings
--
)
{
Gis_line_string
ls
;
int
ls_len
;
if
((
len
<
WKB_HEADER_SIZE
)
||
res
->
reserve
(
WKB_HEADER_SIZE
,
512
))
return
0
;
res
->
q_append
((
char
)
wkb_ndr
);
res
->
q_append
((
uint32
)
wkb_linestring
);
if
(
!
(
ls_len
=
ls
.
init_from_wkb
(
wkb
+
WKB_HEADER_SIZE
,
len
,
(
wkbByteOrder
)
wkb
[
0
],
res
)))
return
0
;
wkb
+=
(
ls_len
+
WKB_HEADER_SIZE
);
len
-=
(
ls_len
+
WKB_HEADER_SIZE
);
}
return
wkb
-
wkb_orig
;
}
bool
Gis_multi_line_string
::
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
{
...
...
@@ -1109,7 +1309,7 @@ int Gis_multi_line_string::length(double *len) const
double
ls_len
;
Gis_line_string
ls
;
data
+=
WKB_HEADER_SIZE
;
ls
.
init_from_wkb
(
data
,
(
uint32
)
(
m_data_end
-
data
));
ls
.
set_data_ptr
(
data
,
(
uint32
)
(
m_data_end
-
data
));
if
(
ls
.
length
(
&
ls_len
))
return
1
;
*
len
+=
ls_len
;
...
...
@@ -1138,7 +1338,7 @@ int Gis_multi_line_string::is_closed(int *closed) const
Gis_line_string
ls
;
if
(
no_data
(
data
,
0
))
return
1
;
ls
.
init_from_wkb
(
data
,
(
uint32
)
(
m_data_end
-
data
));
ls
.
set_data_ptr
(
data
,
(
uint32
)
(
m_data_end
-
data
));
if
(
ls
.
is_closed
(
closed
))
return
1
;
if
(
!*
closed
)
...
...
@@ -1220,6 +1420,41 @@ bool Gis_multi_polygon::init_from_wkt(Gis_read_stream *trs, String *wkb)
}
uint
Gis_multi_polygon
::
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
{
uint32
n_poly
;
const
char
*
wkb_orig
=
wkb
;
if
(
len
<
4
)
return
0
;
n_poly
=
wkb_get_uint
(
wkb
,
bo
);
if
(
res
->
reserve
(
4
,
512
))
return
0
;
res
->
q_append
(
n_poly
);
for
(
wkb
+=
4
;
n_poly
;
n_poly
--
)
{
Gis_polygon
p
;
int
p_len
;
if
(
len
<
WKB_HEADER_SIZE
||
res
->
reserve
(
WKB_HEADER_SIZE
,
512
))
return
0
;
res
->
q_append
((
char
)
wkb_ndr
);
res
->
q_append
((
uint32
)
wkb_polygon
);
if
(
!
(
p_len
=
p
.
init_from_wkb
(
wkb
+
WKB_HEADER_SIZE
,
len
,
(
wkbByteOrder
)
wkb
[
0
],
res
)))
return
0
;
wkb
+=
(
p_len
+
WKB_HEADER_SIZE
);
len
-=
(
p_len
+
WKB_HEADER_SIZE
);
}
return
wkb
-
wkb_orig
;
}
bool
Gis_multi_polygon
::
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
{
uint32
n_polygons
;
...
...
@@ -1356,7 +1591,7 @@ int Gis_multi_polygon::area(double *ar, const char **end_of_data) const
Gis_polygon
p
;
data
+=
WKB_HEADER_SIZE
;
p
.
init_from_wkb
(
data
,
(
uint32
)
(
m_data_end
-
data
));
p
.
set_data_ptr
(
data
,
(
uint32
)
(
m_data_end
-
data
));
if
(
p
.
area
(
&
p_area
,
&
data
))
return
1
;
result
+=
p_area
;
...
...
@@ -1388,7 +1623,7 @@ int Gis_multi_polygon::centroid(String *result) const
while
(
n_polygons
--
)
{
data
+=
WKB_HEADER_SIZE
;
p
.
init_from_wkb
(
data
,
(
uint32
)
(
m_data_end
-
data
));
p
.
set_data_ptr
(
data
,
(
uint32
)
(
m_data_end
-
data
));
if
(
p
.
area
(
&
cur_area
,
&
data
)
||
p
.
centroid_xy
(
&
cur_cx
,
&
cur_cy
))
return
1
;
...
...
@@ -1442,7 +1677,7 @@ uint32 Gis_geometry_collection::get_data_size() const
if
(
!
(
geom
=
create_by_typeid
(
&
buffer
,
wkb_type
)))
return
GET_SIZE_ERROR
;
geom
->
init_from_wkb
(
data
,
(
uint
)
(
m_data_end
-
data
));
geom
->
set_data_ptr
(
data
,
(
uint
)
(
m_data_end
-
data
));
if
((
object_size
=
geom
->
get_data_size
())
==
GET_SIZE_ERROR
)
return
GET_SIZE_ERROR
;
data
+=
object_size
;
...
...
@@ -1482,6 +1717,46 @@ bool Gis_geometry_collection::init_from_wkt(Gis_read_stream *trs, String *wkb)
}
uint
Gis_geometry_collection
::
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
{
uint32
n_geom
;
const
char
*
wkb_orig
=
wkb
;
if
(
len
<
4
)
return
0
;
n_geom
=
wkb_get_uint
(
wkb
,
bo
);
if
(
res
->
reserve
(
4
,
512
))
return
0
;
res
->
q_append
(
n_geom
);
for
(
wkb
+=
4
;
n_geom
;
n_geom
--
)
{
Geometry_buffer
buffer
;
Geometry
*
geom
;
int
g_len
;
uint32
wkb_type
;
if
(
len
<
WKB_HEADER_SIZE
||
res
->
reserve
(
WKB_HEADER_SIZE
,
512
))
return
0
;
res
->
q_append
((
char
)
wkb_ndr
);
wkb_type
=
wkb_get_uint
(
wkb
+
1
,
(
wkbByteOrder
)
wkb
[
0
]);
res
->
q_append
(
wkb_type
);
if
(
!
(
geom
=
create_by_typeid
(
&
buffer
,
wkb_type
))
||
!
(
g_len
=
geom
->
init_from_wkb
(
wkb
+
WKB_HEADER_SIZE
,
len
,
(
wkbByteOrder
)
wkb
[
0
],
res
)))
return
0
;
wkb
+=
(
g_len
+
WKB_HEADER_SIZE
);
len
-=
(
g_len
+
WKB_HEADER_SIZE
);
}
return
wkb
-
wkb_orig
;
}
bool
Gis_geometry_collection
::
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
{
...
...
@@ -1506,7 +1781,7 @@ bool Gis_geometry_collection::get_data_as_wkt(String *txt,
if
(
!
(
geom
=
create_by_typeid
(
&
buffer
,
wkb_type
)))
return
1
;
geom
->
init_from_wkb
(
data
,
(
uint
)
(
m_data_end
-
data
));
geom
->
set_data_ptr
(
data
,
(
uint
)
(
m_data_end
-
data
));
if
(
geom
->
as_wkt
(
txt
,
&
data
))
return
1
;
if
(
txt
->
append
(
","
,
1
,
512
))
...
...
@@ -1541,7 +1816,7 @@ bool Gis_geometry_collection::get_mbr(MBR *mbr, const char **end) const
if
(
!
(
geom
=
create_by_typeid
(
&
buffer
,
wkb_type
)))
return
1
;
geom
->
init_from_wkb
(
data
,
(
uint32
)
(
m_data_end
-
data
));
geom
->
set_data_ptr
(
data
,
(
uint32
)
(
m_data_end
-
data
));
if
(
geom
->
get_mbr
(
mbr
,
&
data
))
return
1
;
}
...
...
@@ -1582,7 +1857,7 @@ int Gis_geometry_collection::geometry_n(uint32 num, String *result) const
if
(
!
(
geom
=
create_by_typeid
(
&
buffer
,
wkb_type
)))
return
1
;
geom
->
init_from_wkb
(
data
,
(
uint
)
(
m_data_end
-
data
));
geom
->
set_data_ptr
(
data
,
(
uint
)
(
m_data_end
-
data
));
if
((
length
=
geom
->
get_data_size
())
==
GET_SIZE_ERROR
)
return
1
;
data
+=
length
;
...
...
@@ -1635,7 +1910,7 @@ bool Gis_geometry_collection::dimension(uint32 *res_dim, const char **end) const
data
+=
WKB_HEADER_SIZE
;
if
(
!
(
geom
=
create_by_typeid
(
&
buffer
,
wkb_type
)))
return
1
;
geom
->
init_from_wkb
(
data
,
(
uint32
)
(
m_data_end
-
data
));
geom
->
set_data_ptr
(
data
,
(
uint32
)
(
m_data_end
-
data
));
if
(
geom
->
dimension
(
&
dim
,
&
end_data
))
return
1
;
set_if_bigger
(
*
res_dim
,
dim
);
...
...
sql/spatial.h
View file @
5150bd8f
...
...
@@ -202,6 +202,10 @@ class Geometry
virtual
const
Class_info
*
get_class_info
()
const
=
0
;
virtual
uint32
get_data_size
()
const
=
0
;
virtual
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
)
=
0
;
/* returns the length of the wkb that was read */
virtual
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
)
=
0
;
virtual
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
=
0
;
virtual
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
=
0
;
virtual
bool
dimension
(
uint32
*
dim
,
const
char
**
end
)
const
=
0
;
...
...
@@ -231,11 +235,13 @@ class Geometry
return
my_reinterpret_cast
(
Geometry
*
)(
buffer
);
}
static
Geometry
*
c
reate_from_wkb
(
Geometry_buffer
*
buffer
,
const
char
*
data
,
uint32
data_len
);
static
Geometry
*
c
onstruct
(
Geometry_buffer
*
buffer
,
const
char
*
data
,
uint32
data_len
);
static
Geometry
*
create_from_wkt
(
Geometry_buffer
*
buffer
,
Gis_read_stream
*
trs
,
String
*
wkt
,
bool
init_stream
=
1
);
static
int
Geometry
::
create_from_wkb
(
Geometry_buffer
*
buffer
,
const
char
*
wkb
,
uint32
len
,
String
*
res
);
int
as_wkt
(
String
*
wkt
,
const
char
**
end
)
{
uint32
len
=
get_class_info
()
->
m_name
.
length
;
...
...
@@ -249,7 +255,7 @@ class Geometry
return
0
;
}
inline
void
init_from_wkb
(
const
char
*
data
,
uint32
data_len
)
inline
void
set_data_ptr
(
const
char
*
data
,
uint32
data_len
)
{
m_data
=
data
;
m_data_end
=
data
+
data_len
;
...
...
@@ -293,6 +299,7 @@ class Gis_point: public Geometry
public:
uint32
get_data_size
()
const
;
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
);
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
);
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
;
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
;
...
...
@@ -339,6 +346,7 @@ class Gis_line_string: public Geometry
public:
uint32
get_data_size
()
const
;
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
);
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
);
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
;
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
;
int
length
(
double
*
len
)
const
;
...
...
@@ -364,6 +372,7 @@ class Gis_polygon: public Geometry
public:
uint32
get_data_size
()
const
;
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
);
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
);
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
;
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
;
int
area
(
double
*
ar
,
const
char
**
end
)
const
;
...
...
@@ -389,6 +398,7 @@ class Gis_multi_point: public Geometry
public:
uint32
get_data_size
()
const
;
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
);
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
);
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
;
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
;
int
num_geometries
(
uint32
*
num
)
const
;
...
...
@@ -410,6 +420,7 @@ class Gis_multi_line_string: public Geometry
public:
uint32
get_data_size
()
const
;
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
);
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
);
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
;
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
;
int
num_geometries
(
uint32
*
num
)
const
;
...
...
@@ -433,6 +444,7 @@ class Gis_multi_polygon: public Geometry
public:
uint32
get_data_size
()
const
;
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
);
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
);
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
;
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
;
int
num_geometries
(
uint32
*
num
)
const
;
...
...
@@ -456,6 +468,7 @@ class Gis_geometry_collection: public Geometry
public:
uint32
get_data_size
()
const
;
bool
init_from_wkt
(
Gis_read_stream
*
trs
,
String
*
wkb
);
uint
init_from_wkb
(
const
char
*
wkb
,
uint
len
,
wkbByteOrder
bo
,
String
*
res
);
bool
get_data_as_wkt
(
String
*
txt
,
const
char
**
end
)
const
;
bool
get_mbr
(
MBR
*
mbr
,
const
char
**
end
)
const
;
int
num_geometries
(
uint32
*
num
)
const
;
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment