def DigitalGlobe2Grid_main(aNameFile, targetSyst, altiMin, altiMax, nbLayers):
inputSyst = "+proj=longlat +datum=WGS84 "
#input syst proj4
refineCoef = ""
binaire = True
stepPixel = 100.
stepCarto = 50.
# Reading Inverse RPC, computing Direct RPC and setting up RPC object
aRPC = RPC()
aRPC.ReadRPB(aNameFile)
print("RPB File read")
# Generating a 50*50*50 grid on the normalized space with random normalized heights
aGridSz = Pt3di(50, 50, 50)
aGridGeoNorm = aRPC.GenerateNormGrid(aGridSz)
print(len(aGridGeoNorm))
# Converting the points to image space
aGridImNorm = []
for i in range(len(aGridGeoNorm)):
aGridImNorm.append(aRPC.InverseRPCNorm(aGridGeoNorm[i]))
print(len(aGridImNorm))
