private static void applyTiePointGeoCoding(
TiffFileInfo info, double[] tiePoints, Product product) {
final SortedSet<Double> xSet = new TreeSet<>();
final SortedSet<Double> ySet = new TreeSet<>();
for (int i = 0; i < tiePoints.length; i += 6) {
xSet.add(tiePoints[i]);
ySet.add(tiePoints[i + 1]);
}
final double xMin = xSet.first();
final double xMax = xSet.last();
final double xDiff = (xMax - xMin) / (xSet.size() - 1);
final double yMin = ySet.first();
final double yMax = ySet.last();
final double yDiff = (yMax - yMin) / (ySet.size() - 1);
final int width = xSet.size();
final int height = ySet.size();
int idx = 0;
final Map<Double, Integer> xIdx = new HashMap<>();
for (Double val : xSet) {
xIdx.put(val, idx);
idx++;
}
idx = 0;
final Map<Double, Integer> yIdx = new HashMap<>();
for (Double val : ySet) {
yIdx.put(val, idx);
idx++;
}
final float[] lats = new float[width * height];
final float[] lons = new float[width * height];
for (int i = 0; i < tiePoints.length; i += 6) {
final int idxX = xIdx.get(tiePoints[i + 0]);
final int idxY = yIdx.get(tiePoints[i + 1]);
final int arrayIdx = idxY * width + idxX;
lons[arrayIdx] = (float) tiePoints[i + 3];
lats[arrayIdx] = (float) tiePoints[i + 4];
}
String[] names = Utils.findSuitableLatLonNames(product);
final TiePointGrid latGrid =
new TiePointGrid(names[0], width, height, xMin, yMin, xDiff, yDiff, lats);
final TiePointGrid lonGrid =
new TiePointGrid(names[1], width, height, xMin, yMin, xDiff, yDiff, lons);
product.addTiePointGrid(latGrid);
product.addTiePointGrid(lonGrid);
final SortedMap<Integer, GeoKeyEntry> geoKeyEntries = info.getGeoKeyEntries();
final Datum datum = getDatum(geoKeyEntries);
product.setGeoCoding(new TiePointGeoCoding(latGrid, lonGrid, datum));
}
private static void applyGcpGeoCoding(
final TiffFileInfo info, final double[] tiePoints, final Product product) {
int numTiePoints = tiePoints.length / 6;
final GcpGeoCoding.Method method;
if (numTiePoints >= GcpGeoCoding.Method.POLYNOMIAL3.getTermCountP()) {
method = GcpGeoCoding.Method.POLYNOMIAL3;
} else if (numTiePoints >= GcpGeoCoding.Method.POLYNOMIAL2.getTermCountP()) {
method = GcpGeoCoding.Method.POLYNOMIAL2;
} else if (numTiePoints >= GcpGeoCoding.Method.POLYNOMIAL1.getTermCountP()) {
method = GcpGeoCoding.Method.POLYNOMIAL1;
} else {
return; // not able to apply GCP geo coding; not enough tie points
}
final int width = product.getSceneRasterWidth();
final int height = product.getSceneRasterHeight();
final GcpDescriptor gcpDescriptor = GcpDescriptor.getInstance();
final ProductNodeGroup<Placemark> gcpGroup = product.getGcpGroup();
for (int i = 0; i < numTiePoints; i++) {
final int offset = i * 6;
final float x = (float) tiePoints[offset + 0];
final float y = (float) tiePoints[offset + 1];
final float lon = (float) tiePoints[offset + 3];
final float lat = (float) tiePoints[offset + 4];
if (Double.isNaN(x) || Double.isNaN(y) || Double.isNaN(lon) || Double.isNaN(lat)) {
continue;
}
final PixelPos pixelPos = new PixelPos(x, y);
final GeoPos geoPos = new GeoPos(lat, lon);
final Placemark gcp =
Placemark.createPointPlacemark(
gcpDescriptor, "gcp_" + i, "GCP_" + i, "", pixelPos, geoPos, product.getGeoCoding());
gcpGroup.add(gcp);
}
final Placemark[] gcps = gcpGroup.toArray(new Placemark[gcpGroup.getNodeCount()]);
final SortedMap<Integer, GeoKeyEntry> geoKeyEntries = info.getGeoKeyEntries();
final Datum datum = getDatum(geoKeyEntries);
product.setGeoCoding(new GcpGeoCoding(method, gcps, width, height, datum));
}
