Skip to content

Commit

Permalink
Commiting recent changes
Browse files Browse the repository at this point in the history
  • Loading branch information
@mchristoffersen
mchristoffersen committed Jan 26, 2024
1 parent a7d1853 commit 385b341
Show file tree
Hide file tree
Showing 12 changed files with 163 additions and 31 deletions.
36 changes: 36 additions & 0 deletions config/sharad_ctx.ini
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
[paths]
outPath = /zippy/MARS/targ/modl/simc/SHARAD_CTX_SIM/ctDist45k_at30_ct
demPath =
navPath =
sigPath =
[navigation]
navsys = +proj=longlat +a=3396190 +no_defs
navfunc = GetNav_FPBgeom_PDS

[simParams]
speedLight = 299792458
dt = 37.5e-9
traceSamples = 3600
body = mars
demBump = False
demInterp = False

[facetParams]
atDist = 400
ctDist = 45000
atStep = 40
ctStep = 5

[outputs]
showNadir = False
showFret = False
combinedAdj = True
combined = True
left = True
right = True
binary = True
echomap = True
echomapAdj = True
red = False
nadir = True
fret = True
10 changes: 5 additions & 5 deletions config/sharad_fpb.ini
View file
Original file line number Diff line number Diff line change
@@ -1,12 +1,12 @@
[paths]
outPath = /home/mchristo/proj/simc/out/fpb/
outPath = /zippy/MARS/targ/modl/simc/SHARAD_MOLA_SIM
demPath =
navPath =
sigPath =

[navigation]
navsys = +proj=longlat +a=3396190 +no_defs
navfunc = GetNav_FPBgeom
navfunc = GetNav_FPBgeom_PDS

[simParams]
speedLight = 299792458
Expand All @@ -18,12 +18,12 @@ demInterp = False

[facetParams]
atDist = 400
ctDist = 45000
ctDist = 24000
atStep = 200
ctStep = 30

[outputs]
showNadir = True
showNadir =True
showFret = True
combinedAdj = True
combined = True
Expand All @@ -34,4 +34,4 @@ echomap = True
echomapAdj = True
red = True
nadir = True
fret = True
fret = true
36 changes: 36 additions & 0 deletions config/sharad_hrsc_dt4.ini
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
[paths]
outPath = /zippy/MARS/targ/modl/simc/SHARAD_HRSC_DT4_SIM
demPath =
navPath =
sigPath =
[navigation]
navsys = +proj=longlat +a=3396190 +no_defs
navfunc = GetNav_FPBgeom_PDS

[simParams]
speedLight = 299792458
dt = 37.5e-9
traceSamples = 3600
body = mars
demBump = False
demInterp = False

[facetParams]
atDist = 200
ctDist = 45000
atStep = 25
ctStep = 15

[outputs]
showNadir = True
showFret = True
combinedAdj = True
combined = True
left = True
right = True
binary = True
echomap = True
echomapAdj = True
red = False
nadir = False
fret = False
14 changes: 7 additions & 7 deletions config/sharad_hrsc_roberto.ini
View file
Original file line number Diff line number Diff line change
Expand Up @@ -18,19 +18,19 @@ demInterp = False
[facetParams]
atDist = 400
ctDist = 45000
atStep = 50
ctStep = 30
atStep = 100
ctStep = 15

[outputs]
showNadir = True
showFret = True
showNadir = False
showFret = False
combinedAdj = True
combined = True
left = True
right = True
binary = True
echomap = True
echomapAdj = True
red = False
nadir = False
fret = False
red = True
nadir = True
fret = True
29 changes: 29 additions & 0 deletions example_ctx.sh
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,29 @@
#!/bin/bash

#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00894601_GEOM.TAB -d ../pds/MOLA_SHARAD_128ppd_radius_mthellasmask.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00894601_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0680XX/S_06802001_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00891301_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0125XX/S_01252601_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0106XX/S_01061401_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_ctx_dt4.ini -n S_00894601_GEOM.TAB -d ../dem/ctx/h1887_0000_dt4.img
#python main.py ./config/sharad_ctx_dt4.ini -n S_01061401_GEOM.TAB -d ../dem/ctx/h1887_0000_dt4.img
#python main.py ./config/sharad_ctx_dt4.ini -n S_00891301_GEOM.TAB -d ../dem/ctx/h1887_0000_dt4.img
#python main.py ./config/sharad_ctx_dt4.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0680XX/S_06802001_GEOM.TAB -d ../dem/ctx/h1423_0001_dt4.img
#python main.py ./config/sharad_ctx_dt4.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00891301_GEOM.TAB -d ../dem/ctx/h1887_0000_dt4.img
#python main.py ./config/sharad_ctx.ini -n S_00891301_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0125XX/S_01252601_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00894601_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0106XX/S_01061401_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n GEOM_LDA/S_00915701_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n GEOM_LDA/S_01255901_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n GEOM_LDA/S_01252601_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif


#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00891301_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00894601_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0091XX/S_00915701_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0095XX/S_00951301_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0106XX/S_01061401_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0125XX/S_01252601_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
#python main.py ./config/sharad_ctx.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0125XX/S_01255901_GEOM.TAB -d /zippy/MARS/targ/modl/MRO/CTX_DEM/blended_r12/blended_r12-tile-0.tif
20 changes: 20 additions & 0 deletions example_hrsc.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4,4 +4,24 @@
#python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00894601_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0680XX/S_06802001_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00891301_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0125XX/S_01252601_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0106XX/S_01061401_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_dt4.ini -n S_00894601_GEOM.TAB -d ../dem/hrsc/h1887_0000_dt4.img
#python main.py ./config/sharad_hrsc_dt4.ini -n S_01061401_GEOM.TAB -d ../dem/hrsc/h1887_0000_dt4.img
#python main.py ./config/sharad_hrsc_dt4.ini -n S_00891301_GEOM.TAB -d ../dem/hrsc/h1887_0000_dt4.img
#python main.py ./config/sharad_hrsc_dt4.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0680XX/S_06802001_GEOM.TAB -d ../dem/hrsc/h1423_0001_dt4.img
#python main.py ./config/sharad_hrsc_dt4.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00891301_GEOM.TAB -d ../dem/hrsc/h1887_0000_dt4.img
#python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00891301_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0126XX/S_01266602_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n S_00891301_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n GEOM_LDA/S_00915701_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n GEOM_LDA/S_01255901_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
#python main.py ./config/sharad_hrsc_roberto.ini -n GEOM_LDA/S_01252601_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif

python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00891301_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0089XX/S_00894601_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0091XX/S_00915701_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0095XX/S_00951301_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0106XX/S_01061401_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0125XX/S_01252601_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
python main.py ./config/sharad_hrsc_roberto.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0125XX/S_01255901_GEOM.TAB -d ../dem/Mars_HRSC_MOLA_BlendDEM_Global_200mp_v2.tif
7 changes: 7 additions & 0 deletions example_mola.sh
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
#!/bin/bash

#python main.py /path/to/config.ini -n /path/to/nav.tab -d /path/to/dem.tif

#python main.py ./config/sharad_fpb.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0126XX/S_01266602_GEOM.TAB -d ../pds/MOLA_SHARAD_128ppd_radius_mthellasmask.tif
#python main.py ./config/sharad_fpb.ini -n S_00891301_GEOM.TAB -d ../pds/MOLA_SHARAD_128ppd_radius_mthellasmask.tif
python main.py ./config/sharad_fpb.ini -n /zippy/MARS/orig/supl/SHARAD/PDS_v2/DATA/GEOM/S_0106XX/S_01061401_GEOM.TAB -d ../pds/MOLA_SHARAD_128ppd_radius_mthellasmask.tif
2 changes: 1 addition & 1 deletion ingest.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -191,7 +191,7 @@ def readConfig(argDict):
}

lleD = {
"mars": "+proj=longlat +a=3396190 +b=3396190 +no_defs",
"mars": "+proj=longlat +R=3396190 +no_defs",
#"mars": "+proj=longlat +a=3396190 +b=3376200 +no_defs",
"moon": "+proj=longlat +a=1737400 +b=1737400 +no_defs",
"earth": "+proj=longlat +a=6378140 +b=6356750 +no_defs",
Expand Down
7 changes: 4 additions & 3 deletions main.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -43,13 +43,17 @@ def main():
confDict["navigation"]["xyzsys"], demCrs
)

print("-----------------------------------------------")
print("xform {}".format(xform))
print(confDict["navigation"]["xyzsys"])
print("-----------------------------------------------")

nav, oDict, inv = prep.prep(confDict, dem, nav)

bounds = prep.calcBounds(
confDict,
dem,
demCrs,
nav,
confDict["navigation"]["xyzsys"],
confDict["facetParams"]["atdist"],
Expand All @@ -69,9 +73,6 @@ def main():
win = rio.windows.Window.from_slices(rowSub, colSub)
print("win {}".format(win))
demData = dem.read(1, window=win)
print("dem shape {}".format(dem.read(1).shape))
print("demData shape {}".format(demData.shape))
#demData = dem.read(1)

with open(confDict["paths"]["logpath"], "a") as fd:
fd.write("Simulating %d traces\n" % len(nav))
Expand Down
22 changes: 11 additions & 11 deletions output.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -60,10 +60,10 @@ def build(confDict, oDict, fcalc, nav, i, oi):
frFacets = cti[cbin == cbin.min()]
oDict["frmap"][frFacets, j] = 1

oDict["pwr"][j,:] = fcalc[:,0]
oDict["twtt"][j,:] = fcalc[:,1]
oDict["theta"][j, :] = fcalc[:,9]
oDict["phi"][j,:] = fcalc[:,10]
#oDict["pwr"][j,:] = fcalc[:,0]
#oDict["twtt"][j,:] = fcalc[:,1]
#oDict["theta"][j, :] = fcalc[:,9]
#oDict["phi"][j,:] = fcalc[:,10]

'''
mlon, mlat, melev = pyproj.transform(
Expand Down Expand Up @@ -142,9 +142,11 @@ def save(confDict, oDict, nav, dem, demData, demCrs, win):
nx, ny, nz = pyproj.transform(
demCrs, confDict["navigation"]["xyzsys"], gx, gy, demz
)

## CHECK THIS LINE <------------------------------------------------------
nlon, nlat, nelev = pyproj.transform(
#dem.crs, confDict["navigation"]["llesys"], gx, gy, demz
confDict["navigation"]["xyzsys"], demCrs, nx, ny, demz # from pds
demCrs, confDict["navigation"]["llesys"], gx, gy, demz
#confDict["navigation"]["xyzsys"], demCrs, nx, ny, demz # from pds
)

nr = np.sqrt((nx - x) ** 2 + (ny - y) ** 2 + (nz - z) ** 2)
Expand Down Expand Up @@ -186,8 +188,8 @@ def save(confDict, oDict, nav, dem, demData, demCrs, win):

flon, flat, felev = pyproj.transform(
confDict["navigation"]["xyzsys"],
"+proj=longlat +R=3396190 +no_defs", #from pds
#confDict["navigation"]["llesys"],
#"+proj=longlat +R=3396190 +no_defs", #from pds
confDict["navigation"]["llesys"],
fret[:, 0],
fret[:, 1],
fret[:, 2],
Expand Down Expand Up @@ -305,7 +307,7 @@ def save(confDict, oDict, nav, dem, demData, demCrs, win):
).mean()
ySquish = confDict["facetParams"]["ctstep"] / postSpace
yDim = np.floor(emap.shape[0] * ySquish).astype(np.int32) + 1

print(yDim)
idx = np.arange(0, emap.shape[0])
nidx = np.floor(idx * ySquish).astype(np.int32)
egram = np.zeros((yDim, emap.shape[1]))
Expand Down Expand Up @@ -333,15 +335,13 @@ def save(confDict, oDict, nav, dem, demData, demCrs, win):
egram_color1 = np.copy(egram)
egram_color2 = np.copy(egram)
egram_color3 = np.copy(egram)
print("----------------------TTTTTTTTTTTTTTT--------------")
print(egram.shape)
print(egram_angles.shape)
x = (egram_angles / 3) % 2 < 1
print(x.shape)
print(egram_angles)
angle = 5
shift = 2.5
print("----------------------TTTTTTTTTTTTTTT--------------")
egram_color1[(((egram_angles / 12) + shift) / angle) % 2 < 1] = egram[(((egram_angles / 12) + shift) / angle) % 2 < 1] * 0.7
#egram_color2[(((egram_angles / 12) + shift) / angle) % 2 >= 1] = egram[(((egram_angles / 12) + shift) / angle) % 2 >= 1] * 0.7
egram_color3 = egram * 0.7#[(((egram_angles / 12) + shift) / angle) % 2 >= 1] = egram[(((egram_angles / 12) + shift) / angle) % 2 >= 1] * 0.6
Expand Down
6 changes: 5 additions & 1 deletion parseNav.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -241,9 +241,13 @@ def GetNav_FPBgeom_PDS(navfile, navsys, xyzsys):
)
df["z"] = (df["elev"] * 1000) * np.sin(np.radians(df["lat"]))

# CHECK THIS LINE <------------------------------------------------------------------------------
#df["datum"] = (1e3*(df["elev"] - df["marsRad"])*2.0/c) - (1800.0*37.5e-9)
df["datum"] = (1e3*(df["elev"] - 3396000.0)*2.0/c) - (1800.0*37.5e-9) # modifiying this line from the PDS

df["areoid"] = 1e3*df["marsRad"] * 0#np.zeros_like(df["x"]) #zval +3396000
# CHECK THIS LINE <------------------------------------------------------------------------------
#df["areoid"] = 1e3*df["marsRad"] * 0#np.zeros_like(df["x"]) #zval +3396000
df["areoid"] = np.zeros_like(df["x"]) #zval +3396000

return df[["x", "y", "z", "datum", "time", "areoid"]]

Expand Down
5 changes: 2 additions & 3 deletions prep.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -105,7 +105,7 @@ def prep(confDict, dem, nav):
return nav, oDict, inv


def calcBounds(confDict, dem, nav, xyzsys, atDist, ctDist):
def calcBounds(confDict, dem, demCrs, nav, xyzsys, atDist, ctDist):
corners = np.zeros((len(nav) * 9, 3))

for i in range(len(nav)):
Expand All @@ -114,8 +114,7 @@ def calcBounds(confDict, dem, nav, xyzsys, atDist, ctDist):
corners[i * 9 : (i * 9) + 9, :] = np.stack((gx, gy, gz), axis=1)

demX, demY, demZ = pyproj.transform(
#xyzsys, dem.crs, corners[:, 0], corners[:, 1], corners[:, 2]
xyzsys, "+proj=longlat +R=3396190 +no_defs", corners[:, 0], corners[:, 1], corners[:, 2]
xyzsys, demCrs, corners[:, 0], corners[:, 1], corners[:, 2]
)
gt = ~dem.transform
print("gt {}".format(gt))
Expand Down

0 comments on commit 385b341

Please sign in to comment.