def _build_standard_products(idaho_id, bbox, proj):
wkt = box(*bbox).wkt
dem = ipe.GeospatialCrop(ipe.IdahoRead(bucketName="idaho-dems", imageId=idaho_id, objectStore="S3"), geospatialWKT=str(wkt))
if proj is not "EPSG:4326":
dem = ipe.Reproject(dem, **reproject_params(proj))
return {
"dem": dem
}
评论列表
文章目录