Hi,
I have a problem when I realize the querys with CQL. The result point is
a little bit out of the street.
I put the code that do it the job, the rest of the method are not so
important. I' am doing something like this:
Point point = null;
StringBuffer query = new StringBuffer();
query.append("ST_NAME = 'Av. Santa Fe'");
query.append("and L_REFADDR <=");
query.append(address.getDoor());
query.append(" and L_NREFADDR >=");
query.append(address.getDoor());
log.info("query "+ query);
Filter filter = CQL.toFilter(query.toString());
//get the shape from a folder on the disk
FeatureSource featureSource =
(FeatureSource)featuresSources.get("cap_fed_ejes"); //cap_fed_ejes.shp
FeatureCollection feactures = featureSource.getFeatures(filter);
CoordinateReferenceSystem targetCRS1 = null;
try {
targetCRS1 = CRS.decode("EPSG:4326");//CRS.parseWKT(WGS84);
} catch (FactoryException e1) {
e1.printStackTrace();
}
MathTransform trans1 = null;
try {
CoordinateReferenceSystem sourceCRS = targetCRS1;
log.info("Sistema de coordenadas ordiginal:
"+sourceCRS.toString());
trans1 = CRS.findMathTransform(sourceCRS, targetCRS1);
} catch (NoSuchAuthorityCodeException e) {
e.printStackTrace();
} catch (FactoryException e) {
e.printStackTrace();
}
FeatureIterator iterator = feactures.features();
if (feactures.size() != 1) new GeoCodeException("Se encontro más de
un valor de geocodificación");
try {
SimpleFeature feature = (SimpleFeature)iterator.next();
log.info("Puntos del segmento: ");
log.info(" inicio impar x " + feature.getAttribute("AII"));
log.info("; fin impar y " + feature.getAttribute("AFI"));
log.info("; inicio par x " + feature.getAttribute("AIP"));
log.info("; fin par y " + feature.getAttribute("AFP"));
log.info("---------------------------------------------------");
log.info(" getMinX " + feature.getBounds().getMinX());
log.info("; getMinY " + feature.getBounds().getMinY());
log.info("; getMaxX " + feature.getBounds().getMaxX());
log.info("; getMaxY " + feature.getBounds().getMaxY());
log.info
("----------------------------------------------------");
Map points = new Hashtable();
points.put("x1", String.valueOf(feature.getBounds().getMinX()));
points.put("y1", String.valueOf(feature.getBounds().getMinY()));
DirectPosition p1 = new GeneralDirectPosition(new
Double(feature.getBounds().getMinX()).doubleValue(),new
Double(feature.getBounds().getMinY()).doubleValue());
log.info
("------------------------------------------------------------");
log.info("Punto original p1 : " + p1);
points.put("x2", String.valueOf(feature.getBounds().getMaxX()));
points.put("y2", String.valueOf(feature.getBounds().getMaxY()));
DirectPosition p2 = new GeneralDirectPosition(new
Double(feature.getBounds().getMaxX()).doubleValue(),new
Double(feature.getBounds().getMaxY()).doubleValue());
log.info
("------------------------------------------------------------");
log.info("Punto original p2 : " + p2);
point = UtilInterpolation.interpolate(points,
address.getDoor());
log.info("------------------------------------------------");
log.info("Punto interpolado - x: " + point.getX() + " y: " +
point.getY());
} catch (MismatchedDimensionException e) {
e.printStackTrace();
}
This is the log.
11:04:10.980 INFO [main] com.ases.geo.services.impl.GeoService:
---------------------------------------------------
11:04:10.980 INFO [main] com.ases.geo.services.impl.GeoService: getMinX
-58.416011810302734
11:04:10.980 INFO [main] com.ases.geo.services.impl.GeoService: ; getMinY
-34.585758209228516
11:04:10.980 INFO [main] com.ases.geo.services.impl.GeoService: ; getMaxX
-58.41497039794922
11:04:10.980 INFO [main] com.ases.geo.services.impl.GeoService: ; getMaxY
-34.585121154785156
11:04:10.980 INFO [main] com.ases.geo.services.impl.GeoService:
----------------------------------------------------
11:04:10.981 INFO [main] com.ases.geo.services.impl.GeoService:
------------------------------------------------------------
11:04:10.981 INFO [main] com.ases.geo.services.impl.GeoService: Punto
original p1 : GeneralDirectPosition[-58.416011810302734,
-34.585758209228516]
11:04:10.981 INFO [main] com.ases.geo.services.impl.GeoService:
------------------------------------------------------------
11:04:10.981 INFO [main] com.ases.geo.services.impl.GeoService: Punto
original p2 : GeneralDirectPosition[-58.41497039794922, -34.585121154785156]
11:04:10.982 INFO [main] com.ases.geo.services.impl.GeoService:
------------------------------------------------
11:04:10.982 INFO [main] com.ases.geo.services.impl.GeoService: Punto
interpolado - x: -58.41507453918457 y: -34.58518486022949
11:04:10.982 INFO [main] com.ases.geo.services.impl.GeoService: Fin del
proceso de geo, Duracion: 1s:675ms
11:04:10.982 INFO [main] com.ases.geo.GeoServiceTest: Y,X =
-34.58518486022949,-58.41507453918457
When I search the point Y,X = -34.58518486022949,-58.41507453918457 into
google maps, it's a little bit out of the street. To see it please go to
http://maps.google.com/ and put this point (just copy like this
-34.58518486022949,-58.41507453918457). If you see the arrow, you will see
in the middle of the block.
Anybody
------------------------------------------------------------------------------
Come build with us! The BlackBerry® Developer Conference in SF, CA
is the only developer event you need to attend this year. Jumpstart your
developing skills, take BlackBerry mobile applications to market and stay
ahead of the curve. Join us from November 9-12, 2009. Register now!
http://p.sf.net/sfu/devconf
_______________________________________________
Geotools-gt2-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
