Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Support GeoJSON uploads for creating opportunity datasets. #902

Open
wants to merge 2 commits into
base: dev
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Original file line number Diff line number Diff line change
Expand Up @@ -358,6 +358,9 @@ private OpportunityDatasetUploadStatus createOpportunityDataset(Request req, Res
} else if (uploadFormat == FileStorageFormat.SHP) {
LOG.info("Detected opportunity dataset stored as ESRI shapefile.");
pointsets.addAll(createGridsFromShapefile(fileItems, zoom, status));
} else if (uploadFormat == FileStorageFormat.GEOJSON) {
LOG.info("Detected opportunity dataset stored as GeoJSON.");
pointsets.addAll(Grid.fromGeoJson(fileItems.get(0).getInputStream(), zoom, status));
} else if (uploadFormat == FileStorageFormat.CSV) {
LOG.info("Detected opportunity dataset stored as CSV");
// Create a grid even when user has requested a freeform pointset so we have something to visualize.
Expand Down
67 changes: 67 additions & 0 deletions src/main/java/com/conveyal/r5/analyst/Grid.java
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,9 @@
import org.geotools.data.FileDataStore;
import org.geotools.data.FileDataStoreFinder;
import org.geotools.data.Transaction;
import org.geotools.data.geojson.GeoJSONReader;
import org.geotools.data.simple.SimpleFeatureCollection;
import org.geotools.data.simple.SimpleFeatureIterator;
import org.geotools.feature.simple.SimpleFeatureTypeBuilder;
import org.geotools.gce.geotiff.GeoTiffFormat;
import org.geotools.gce.geotiff.GeoTiffWriteParams;
Expand Down Expand Up @@ -753,6 +756,70 @@ public static List<Grid> fromShapefile (File shapefile, int zoom, ProgressListen
return new ArrayList<>(grids.values());
}

/**
* Take an `InputStream` containing GeoJson Features and turn it into an opportunity grid.
*/
public static List<Grid> fromGeoJson (InputStream geoJsonInputStream, int zoom, ProgressListener progressListener)
throws IOException {
GeoJSONReader reader = new GeoJSONReader(geoJsonInputStream);
SimpleFeatureCollection features = reader.getFeatures();
Envelope envelope = features.getBounds();

checkWgsEnvelopeSize(envelope, "GeoJSON file");
WebMercatorExtents extents = WebMercatorExtents.forWgsEnvelope(envelope, zoom);

int total = features.size();
if (progressListener != null) {
progressListener.setTotalItems(total);
}

AtomicInteger count = new AtomicInteger(0);
HashMap<String, Grid> grids = new HashMap<>();

SimpleFeatureIterator featureIterator = features.features();
while (featureIterator.hasNext()) {
SimpleFeature feature = featureIterator.next();
Geometry geom = (Geometry) feature.getDefaultGeometry();

for (var p : feature.getProperties()) {
var val = p.getValue();

if (!(val instanceof Number)) continue;
double numericVal = ((Number) val).doubleValue();
if (numericVal == 0) continue;

String attributeName = p.getName().getLocalPart();

Grid grid = grids.get(attributeName);
if (grid == null) {
grid = new Grid(extents);
grid.name = attributeName;
grids.put(attributeName, grid);
}

if (geom instanceof Point) {
Point point = (Point) geom;
// already in WGS 84
grid.incrementPoint(point.getY(), point.getX(), numericVal);
} else if (geom instanceof Polygon || geom instanceof MultiPolygon) {
grid.rasterize(geom, numericVal);
} else {
throw new IllegalArgumentException("Unsupported geometry type: " + geom);
}
}

int currentCount = count.incrementAndGet();
if (progressListener != null) {
progressListener.setCompletedItems(currentCount);
}
if (currentCount % 10000 == 0) {
LOG.info("{} / {} features read", human(currentCount), human(total));
}
}
reader.close();
return new ArrayList<>(grids.values());
}

@Override
public double sumTotalOpportunities() {
double totalOpportunities = 0;
Expand Down
Loading