search

AdrianKriger / osm_LoD1_3DCityModel

LoD1 3D City Model from OpenStreetMap contributions and raster DEM
https://adriankriger.github.io/osm_LoD1_3DCityModel/
MIT License
8 stars 1 forks source link

Districts -> 4 dem rasters #7

Open AdrianKriger opened 2 years ago

AdrianKriger commented 2 years ago

Perhaps it would be more useful to mosaic, clip, re-project -> extract -> TIN -> a minimum of 4 adjacent dem raster's. I'm thinking specifically of Pretoria CBD. 2528 CA, CB, CC, CD.

No more than that.

AdrianKriger commented 2 years ago

try:



'in_raster': "./raster/LO29_050M_2528CA.tif ./raster/LO29_050M_2528CB.tif ./raster/LO29_050M_2528CC.tif ./raster/LO29_050M_2528CD.tif"

a_string = jparams['in_raster']
imgs = a_string.split()

if len(imgs) > 1:
    imgs = '", "'.join(a_string.split())
    new = '"' + imgs + '"'

    OutTile = gdal.Warp(jparams['projClip_raster'], 
                            list(new),
                            resampleAlg='bilinear',
                            dstSRS=jparams['crs'],
                            #-- outputBounds=[minX, minY, maxX, maxY]
                            outputBounds = [extent[0], extent[1], extent[2], extent[3]])
     OutTile = None