1 |
3044
|
perry
|
/**
|
2 |
|
|
* '$RCSfile$'
|
3 |
|
|
* Copyright: 2000 Regents of the University of California and the
|
4 |
|
|
* National Center for Ecological Analysis and Synthesis
|
5 |
|
|
*
|
6 |
|
|
* Author: Matthew Perry
|
7 |
|
|
* '$Date$'
|
8 |
|
|
* '$Revision$'
|
9 |
|
|
*
|
10 |
|
|
* This program is free software; you can redistribute it and/or modify
|
11 |
|
|
* it under the terms of the GNU General Public License as published by
|
12 |
|
|
* the Free Software Foundation; either version 2 of the License, or
|
13 |
|
|
* (at your option) any later version.
|
14 |
|
|
*
|
15 |
|
|
* This program is distributed in the hope that it will be useful,
|
16 |
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
17 |
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
18 |
|
|
* GNU General Public License for more details.
|
19 |
|
|
*
|
20 |
|
|
* You should have received a copy of the GNU General Public License
|
21 |
|
|
* along with this program; if not, write to the Free Software
|
22 |
|
|
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
|
23 |
|
|
*/
|
24 |
|
|
|
25 |
|
|
package edu.ucsb.nceas.metacat.spatial;
|
26 |
|
|
|
27 |
|
|
import org.geotools.data.shapefile.ShapefileDataStore;
|
28 |
|
|
import org.geotools.data.FeatureStore;
|
29 |
|
|
import org.geotools.data.FeatureSource;
|
30 |
|
|
import org.geotools.data.DefaultTransaction;
|
31 |
|
|
import org.geotools.data.Transaction;
|
32 |
|
|
import org.geotools.filter.Filter;
|
33 |
|
|
import org.geotools.filter.AbstractFilter;
|
34 |
|
|
import org.geotools.filter.Expression;
|
35 |
|
|
import org.geotools.filter.CompareFilter;
|
36 |
|
|
import org.geotools.filter.GeometryFilter;
|
37 |
|
|
import org.geotools.filter.FilterFactory;
|
38 |
|
|
import org.geotools.filter.FilterFactoryFinder;
|
39 |
|
|
import org.geotools.filter.IllegalFilterException;
|
40 |
|
|
import org.geotools.feature.FeatureCollection;
|
41 |
|
|
import org.geotools.feature.FeatureCollections;
|
42 |
5829
|
leinfelder
|
import org.opengis.feature.simple.SimpleFeature;
|
43 |
|
|
import org.opengis.feature.simple.SimpleFeatureType;
|
44 |
|
|
|
45 |
3044
|
perry
|
import com.vividsolutions.jts.geom.Envelope;
|
46 |
3048
|
perry
|
import com.vividsolutions.jts.geom.Geometry;
|
47 |
|
|
import com.vividsolutions.jts.geom.GeometryFactory;
|
48 |
|
|
import com.vividsolutions.jts.geom.Polygon;
|
49 |
|
|
import com.vividsolutions.jts.geom.Coordinate;
|
50 |
3044
|
perry
|
|
51 |
|
|
import java.io.File;
|
52 |
|
|
import java.net.URI;
|
53 |
|
|
import java.net.URL;
|
54 |
|
|
import java.net.MalformedURLException;
|
55 |
|
|
import java.sql.ResultSet;
|
56 |
|
|
import java.sql.PreparedStatement;
|
57 |
|
|
import java.util.Vector;
|
58 |
|
|
import java.util.Iterator;
|
59 |
|
|
import java.io.IOException;
|
60 |
|
|
|
61 |
|
|
import org.apache.log4j.Logger;
|
62 |
|
|
|
63 |
4699
|
daigle
|
import edu.ucsb.nceas.metacat.util.MetacatUtil;
|
64 |
3044
|
perry
|
|
65 |
|
|
/**
|
66 |
|
|
* Class to query the persistent spatial cache
|
67 |
|
|
* and returns docids matching spatial constraints
|
68 |
|
|
*/
|
69 |
|
|
public class SpatialQuery {
|
70 |
|
|
|
71 |
|
|
private static Logger log = Logger.getLogger(SpatialQuery.class.getName());
|
72 |
|
|
|
73 |
|
|
/**
|
74 |
|
|
* empty constructor to initialize spatial query
|
75 |
|
|
*/
|
76 |
|
|
public SpatialQuery() { }
|
77 |
|
|
|
78 |
|
|
/**
|
79 |
|
|
* Querys all features in the spatial cache
|
80 |
|
|
* and filters based on bouding coordinates.
|
81 |
|
|
* Returns Vector of docids.
|
82 |
|
|
*
|
83 |
|
|
* @param w West bounding coordinate
|
84 |
|
|
* @param s South bounding coordinate
|
85 |
|
|
* @param e East bounding coordinate
|
86 |
|
|
* @param n North bounding coordinate
|
87 |
|
|
*
|
88 |
|
|
*/
|
89 |
4080
|
daigle
|
public Vector<String> filterByBbox( float w, float s, float e, float n ) {
|
90 |
|
|
Vector<String> docids = new Vector<String>();
|
91 |
3044
|
perry
|
SpatialFeatureSchema featureSchema = new SpatialFeatureSchema();
|
92 |
|
|
|
93 |
|
|
ShapefileDataStore store = null;
|
94 |
|
|
FeatureSource features = null;
|
95 |
|
|
FeatureCollection collection = null;
|
96 |
|
|
FilterFactory filterFactory = FilterFactoryFinder.createFilterFactory();
|
97 |
|
|
|
98 |
|
|
try {
|
99 |
3045
|
perry
|
// read the spatial cache (polygons)
|
100 |
3044
|
perry
|
store = new ShapefileDataStore( (new File( featureSchema.polygonShpUri )).toURL() );
|
101 |
|
|
features = store.getFeatureSource(store.getTypeNames()[0]);
|
102 |
|
|
|
103 |
3045
|
perry
|
// Construct bounding box
|
104 |
3044
|
perry
|
Envelope envelope = new Envelope( w, e, s, n );
|
105 |
|
|
Expression bbox = filterFactory.createBBoxExpression( envelope );
|
106 |
|
|
|
107 |
3048
|
perry
|
// Construct the bbox as an actual geometry
|
108 |
|
|
Coordinate[] linestringCoordinates = new Coordinate[5];
|
109 |
|
|
linestringCoordinates[0] = new Coordinate( w, s );
|
110 |
|
|
linestringCoordinates[1] = new Coordinate( w, n );
|
111 |
|
|
linestringCoordinates[2] = new Coordinate( e, n );
|
112 |
|
|
linestringCoordinates[3] = new Coordinate( e, s );
|
113 |
|
|
linestringCoordinates[4] = new Coordinate( w, s );
|
114 |
|
|
|
115 |
|
|
GeometryFactory geomFac = new GeometryFactory();
|
116 |
|
|
Polygon bboxGeom = geomFac.createPolygon( geomFac.createLinearRing(linestringCoordinates), null);
|
117 |
|
|
|
118 |
3045
|
perry
|
// Set up geometry filter based on bbox
|
119 |
5829
|
leinfelder
|
SimpleFeatureType featureType = store.getSchema( store.getTypeNames()[0] );
|
120 |
|
|
Expression geometry = filterFactory.createAttributeExpression( featureType.getGeometryDescriptor().getName().toString());
|
121 |
3044
|
perry
|
GeometryFilter bboxFilter = filterFactory.createGeometryFilter(AbstractFilter.GEOMETRY_BBOX);
|
122 |
|
|
bboxFilter.addLeftGeometry( geometry );
|
123 |
|
|
bboxFilter.addRightGeometry( bbox );
|
124 |
|
|
|
125 |
3045
|
perry
|
// Iterate through the filtered feature collection
|
126 |
|
|
// and add matches to the docid Vector
|
127 |
3044
|
perry
|
collection = features.getFeatures(bboxFilter);
|
128 |
|
|
Iterator iterator = collection.iterator();
|
129 |
|
|
try {
|
130 |
|
|
for( Iterator i=collection.iterator(); i.hasNext(); ) {
|
131 |
5829
|
leinfelder
|
SimpleFeature feature = (SimpleFeature) i.next();
|
132 |
3048
|
perry
|
|
133 |
|
|
Geometry geom = (Geometry)feature.getAttribute(0);
|
134 |
3492
|
leinfelder
|
if ( geom.within( bboxGeom ) ) {
|
135 |
3048
|
perry
|
// assumes docid is attribute number 1
|
136 |
|
|
// in a zero-based index of dbf columns
|
137 |
|
|
docids.add( (String) feature.getAttribute(1) );
|
138 |
|
|
}
|
139 |
3044
|
perry
|
}
|
140 |
|
|
} finally {
|
141 |
|
|
collection.close( iterator );
|
142 |
|
|
}
|
143 |
|
|
|
144 |
3045
|
perry
|
/*
|
145 |
|
|
* Also query the point cache since there may be point-only documents
|
146 |
|
|
* Filter by the bbox AND check against docids Vector so that
|
147 |
|
|
* docids already in the Vector don't get duplicated.
|
148 |
|
|
*/
|
149 |
|
|
// read the spatial cache (points)
|
150 |
|
|
store = new ShapefileDataStore( (new File( featureSchema.pointShpUri )).toURL() );
|
151 |
|
|
features = store.getFeatureSource(store.getTypeNames()[0]);
|
152 |
3044
|
perry
|
|
153 |
3045
|
perry
|
// Set up geometry filter based on bbox
|
154 |
|
|
featureType = store.getSchema( store.getTypeNames()[0] );
|
155 |
5829
|
leinfelder
|
geometry = filterFactory.createAttributeExpression( featureType.getGeometryDescriptor().getName().toString() );
|
156 |
3045
|
perry
|
bboxFilter = filterFactory.createGeometryFilter(AbstractFilter.GEOMETRY_BBOX);
|
157 |
|
|
bboxFilter.addLeftGeometry( geometry );
|
158 |
|
|
bboxFilter.addRightGeometry( bbox );
|
159 |
3044
|
perry
|
|
160 |
3045
|
perry
|
// Iterate through the filtered feature collection
|
161 |
|
|
// and add matches to the docid Vector IF
|
162 |
|
|
// they aren't already present
|
163 |
|
|
collection = features.getFeatures(bboxFilter);
|
164 |
|
|
iterator = collection.iterator();
|
165 |
|
|
String docid = null;
|
166 |
|
|
try {
|
167 |
|
|
for( Iterator i=collection.iterator(); i.hasNext(); ) {
|
168 |
5829
|
leinfelder
|
SimpleFeature feature = (SimpleFeature) i.next();
|
169 |
3048
|
perry
|
Geometry geom = (Geometry)feature.getAttribute(0);
|
170 |
|
|
if ( geom.intersects( bboxGeom ) ) {
|
171 |
|
|
// assumes docid is attribute number 1
|
172 |
|
|
// in a zero-based index of dbf columns
|
173 |
|
|
docid = (String) feature.getAttribute(1);
|
174 |
|
|
if( !docids.contains( docid ) ) {
|
175 |
|
|
docids.add( docid );
|
176 |
|
|
}
|
177 |
3045
|
perry
|
}
|
178 |
|
|
}
|
179 |
|
|
} finally {
|
180 |
|
|
collection.close( iterator );
|
181 |
|
|
}
|
182 |
|
|
|
183 |
|
|
} catch (MalformedURLException ex) {
|
184 |
|
|
ex.printStackTrace();
|
185 |
|
|
} catch (IOException ex) {
|
186 |
|
|
ex.printStackTrace();
|
187 |
|
|
} catch (IllegalFilterException ex) {
|
188 |
|
|
ex.printStackTrace();
|
189 |
|
|
}
|
190 |
|
|
|
191 |
3044
|
perry
|
return docids;
|
192 |
|
|
|
193 |
|
|
}
|
194 |
|
|
|
195 |
|
|
}
|