Hi, Please help me. I have a mess in my head.
My input:
N lines with values:
-longitude, latittude, (a senter of a sector)
-radius (in meters 10...5000 meters)
-start angle (from "12 o'clock")
-angle width (in defrees)
My output:
Find sectors that intersects, output them to KML (or something else) and
visualize (I've chosen KML because it can be uploaded to google maps and
easilly seen)
Here is my code....
@Test(enabled = false)
public void visualizeRetailStoreMC7700001() throws ParserException,
IOException, FactoryException, TransformException {
List<GSMCellLocation> gsmCellLocations =
GeoDataSampleHelper.readGSMCellLocationFromMoscowRegion();
RetailStore retailStore =
GeoDataSampleHelper.readRetailStoreByCode(MC7700001);
GeoPoint storeGeoPoint =
retailStore.getStoreLocation().getGeoPoint();
List<Geometry> geometryList = new LinkedList<Geometry>();
GeometricShapeFactory gsFactory = new GeometricShapeFactory();
gsFactory.setCentre(transformToPlanar(new
Coordinate(storeGeoPoint.getLongitude(),
storeGeoPoint.getLatitude())));//storeGeoPoint.getLatitude(),
storeGeoPoint.getLongitude())));
gsFactory.setSize(200);
Geometry storeCircle = gsFactory.createEllipse();
geometryList.add(storeCircle);
for(GSMCellLocation cellLoc : gsmCellLocations){
gsFactory = new GeometricShapeFactory();
gsFactory.setCentre(transformToPlanar(new
Coordinate(cellLoc.getGeoPoint().getLongitude(),
cellLoc.getGeoPoint().getLatitude())));//cellLoc.getGeoPoint().getLatitude(),
cellLoc.getGeoPoint().getLongitude())));
gsFactory.setSize(cellLoc.getMidDist() * 2);
Geometry cellSector =
gsFactory.createArc(cellLoc.getStartAngle(), cellLoc.getAngleWidth());
if(cellSector.intersects(storeCircle)){
geometryList.add(cellSector);
}
}
LOG.info("geometryList.size {}",geometryList.size());
Encoder encoder = new Encoder(new KMLConfiguration());
encoder.setIndenting(true);
SimpleFeatureTypeBuilder typeBuilder = new
SimpleFeatureTypeBuilder();
typeBuilder.setName("poi");
typeBuilder.add("geometry", Polygon.class,
DefaultGeographicCRS.WGS84);//DefaultGeographicCRS.WGS84);
SimpleFeatureType TYPE = typeBuilder.buildFeatureType();
DefaultFeatureCollection features = new
DefaultFeatureCollection("test",TYPE);
for(int i=0; i<10; i++){
SimpleFeatureBuilder featureBuilder = new
SimpleFeatureBuilder(TYPE);
featureBuilder.add(
transformToGeo(geometryList.get(i)));//transformToGeo(
SimpleFeature simpleFeature =
featureBuilder.buildFeature(String.valueOf(i + 1));
features.add( simpleFeature );
}
encoder.encode(features, KML.kml, System.out);
}
What do I expect:
Correctly calculated intersections between sectors
correct kml output with intersected sectors
Ofcourse, I don't get anyting correctly... (((
Please help, I've spent more than a week trying to deal with geotools. Seems
like I completely don't understand the idea of reprojecting CRS...
--
View this message in context:
http://osgeo-org.1560.n6.nabble.com/WGS84-to-a-Cartesian-transformation-tp4318606p5034219.html
Sent from the geotools-gt2-users mailing list archive at Nabble.com.
------------------------------------------------------------------------------
Free Next-Gen Firewall Hardware Offer
Buy your Sophos next-gen firewall before the end March 2013
and get the hardware for free! Learn more.
http://p.sf.net/sfu/sophos-d2d-feb
_______________________________________________
GeoTools-GT2-Users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
