Mercurial > gemma

comparison pkg/imports/sr.go @ 4768:a2f16bbcc846 direct-diff

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Morph differences: Directly raster A and subtract B as a raster.
author Sascha L. Teichmann <sascha.teichmann@intevation.de>
date Mon, 21 Oct 2019 02:01:56 +0200
parents d786c37b02c0
children 0feddad5bb6d
comparison
equal deleted inserted replaced
4767:f66e5b2fa894 4768:a2f16bbcc846
891 defer func() { 891 defer func() {
892 feedback.Info("Generating iso areas took %s.", 892 feedback.Info("Generating iso areas took %s.",
893 time.Since(total)) 893 time.Since(total))
894 }() 894 }()
895 895
896 areas := octree.TraceAreas(heights, isoCellSize, tree.Min(), tree.Max(), tree.Value) 896 box := octree.Box2D{
897 X1: tree.Min().X,
898 Y1: tree.Min().Y,
899 X2: tree.Max().X,
900 Y2: tree.Max().Y,
901 }
902
903 raster := octree.NewRaster(box, isoCellSize)
904 raster.Rasterize(tree.Value)
905 areas := raster.Trace(heights)
897 906
898 return storeAreas( 907 return storeAreas(
899 ctx, tx, feedback, 908 ctx, tx, feedback,
900 areas, tree.EPSG(), heights, id) 909 areas, tree.EPSG(), heights, id)
901 } 910 }