public static JGrassRegion getJGrassRegionFromGridCoverage(GridCoverage2D gridCoverage2D)
throws InvalidGridGeometryException, TransformException {
Envelope2D env = gridCoverage2D.getEnvelope2D();
GridEnvelope2D worldToGrid = gridCoverage2D.getGridGeometry().worldToGrid(env);
double xRes = env.getWidth() / worldToGrid.getWidth();
double yRes = env.getHeight() / worldToGrid.getHeight();
JGrassRegion region =
new JGrassRegion(env.getMinX(), env.getMaxX(), env.getMinY(), env.getMaxY(), xRes, yRes);
return region;
}
public void testTest1() throws Exception {
double cellSize = 0.000249;
Point _pt = (new GeometryFactory()).createPoint(new Coordinate(80.333326615, 41.826389655));
Envelope2D _e = new Envelope2D(null, _pt.getX(), _pt.getY(), cellSize, cellSize);
System.out.println(_e.getX());
System.out.println(_e.getY());
System.out.println(_e.getWidth());
System.out.println(_e.getHeight());
String _wkt =
"PROJCS[\"Asia_North_Albers_Equal_Area_Conic\",GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137.0,298.257223563]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Albers\"],PARAMETER[\"False_Easting\",0.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",95.0],PARAMETER[\"Standard_Parallel_1\",15.0],PARAMETER[\"Standard_Parallel_2\",65.0],PARAMETER[\"Latitude_Of_Origin\",30.0],UNIT[\"Meter\",1.0]]";
CoordinateReferenceSystem _sourcePrj = CRS.decode("EPSG:4326", true);
// CoordinateReferenceSystem _targetPrj = CRS.decode("EPSG:21416");
CoordinateReferenceSystem _targetPrj = CRS.parseWKT(_wkt);
double _pixelArea =
(new ProjectTransformModel(_sourcePrj, _targetPrj))
.calculate(
Utilities.covertEnvelope2D(
new Envelope2D(null, _pt.getX(), _pt.getY(), cellSize, cellSize)))
.getArea();
double _pixelLen1 =
(new ProjectTransformModel(_sourcePrj, _targetPrj))
.calculate(
(new GeometryFactory())
.createLineString(
new Coordinate[] {
new Coordinate(80.333326615, 41.826389655),
new Coordinate(80.333326615 + cellSize, 41.826389655)
}))
.getLength();
double _pixelLen2 =
(new ProjectTransformModel(_sourcePrj, _targetPrj))
.calculate(
(new GeometryFactory())
.createLineString(
new Coordinate[] {
new Coordinate(80.333326615, 41.826389655),
new Coordinate(80.333326615, 41.826389655 + cellSize)
}))
.getLength();
// double _pixelArea = (new ProjectTransformModel("EPSG:4326",
// "EPSG:32644")).calculate(Utilities.covertEnvelope2D(new Envelope2D(null, _pt.getX(),
// _pt.getY(), cellSize, cellSize))).getArea();
System.out.println(_pixelArea);
System.out.println(_pixelLen1);
System.out.println(_pixelLen2);
// Polygon _box = Utilities.covertEnvelope2D(new Envelope2D(null, _pt.getX(), _pt.getY(),
// cellSize, cellSize));
Polygon _box =
(Polygon)
(new ProjectTransformModel(_sourcePrj, _targetPrj))
.calculate(
Utilities.covertEnvelope2D(
new Envelope2D(null, _pt.getX(), _pt.getY(), cellSize, cellSize)));
SimpleFeatureTypeBuilder _builder = new SimpleFeatureTypeBuilder();
_builder.setName("project1.shp");
_builder.add("Shape", Polygon.class);
_builder.add("Test", Double.class);
FeatureType _ft = _builder.buildFeatureType();
FeatureCollection _fc =
(CommonFactoryFinder.getFeatureCollections(GeoTools.getDefaultHints()).newCollection());
_fc.add(SimpleFeatureBuilder.build((SimpleFeatureType) _ft, new Object[] {_box, 10}, null));
ShapefileDataStore _fs =
new ShapefileDataStore(new File("D:\\Temp\\data3\\20081010\\shp\\test.shp").toURL());
_fs.createSchema((SimpleFeatureType) _ft);
FeatureStore _store = (FeatureStore) _fs.getFeatureSource();
_store.addFeatures(_fc);
}
