Comments (5)
The PR #196 change worked and the geom_reference
output files were created, including the watermask download and reprojection. I was able to load the resulting UAVSAR stack into MintPy. Here is the view of the geometry data:
It looks like there is another much smaller problem with the incidenceAngle and azimuthAngle files that might not have the right data type set.
from isce2.
I added gdal.UseExceptions()
to the isce/applications/gdal2isce_xml.py
code and I get a different message now:
GDAL open (R): /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/merged/geom_reference/hgt.rdr.vrt
API open (WR): /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/merged/geom_reference/simamp.rdr
GDAL open (R): /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/merged/geom_reference/hgt.rdr.vrt
GDAL close: /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/merged/geom_reference/hgt.rdr.vrt
API close: /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/merged/geom_reference/simamp.rdr
Warning 1: Geotransform matrix has non rotational terms
C pointer already created. Finalize and recreate if image dimensions changed.
C pointer already created. Finalize and recreate if image dimensions changed.
--------------------------------------------------
generate multilooked geometry files with alks=16 and rlks=4 using gdal.Translate() ...
multilook /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/merged/geom_reference/hgt.rdr
.vrt
/u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/geom_reference/hgt.rdr
Traceback (most recent call last):
File "/u/pez0/fielding/tools/ISCE2_test/isce2/contrib/stack/stripmapStack/stripmapWrapper.py", line 155, in <module>
main(args.start,args.end)
File "/u/pez0/fielding/tools/ISCE2_test/isce2/contrib/stack/stripmapStack/stripmapWrapper.py", line 146, in main
cfgParser.runCmd()
File "/u/pez0/fielding/tools/ISCE2_test/isce2/contrib/stack/stripmapStack/stripmapWrapper.py", line 51, in runCmd
func_modules.main(self.funcParams[section])
File "/u/pez0/fielding/tools/ISCE2_test/isce2/contrib/stack/stripmapStack/topo.py", line 525, in main
runMultilook(in_dir=info.outdir, out_dir=out_dir, alks=inps.alks, rlks=inps.rlks)
File "/u/pez0/fielding/tools/ISCE2_test/isce2/contrib/stack/stripmapStack/topo.py", line 423, in runMultilook
gdal2isce_xml(out_file+'.vrt')
File "/u/pez0/fielding/tools/ISCE2_test/install/isce/applications/gdal2isce_xml.py", line 70, in gdal2isce_xml
ds = gdal.Open(fname, gdal.GA_ReadOnly)
File "/u/vento-r0/ericf/packages/Anaconda/anaconda3/envs/isce2/lib/python3.8/site-packages/osgeo/gdal.py", line 4192, in Open
return _gdal.Open(*args)
RuntimeError: /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s1/geom_reference/hgt.rdr.vrt: No such file or directory
It is creating the s1/geom_reference/hgt.rdr
file along with s1/geom_reference/hgt.hdr
and s1/geom_reference/hgt.rdr.aux.xml
files, but it is not creating the s1/geom_reference/hgt.rdr.vrt
file for some reason.
from isce2.
Thanks @yunjunz. I tried adding that and it did fix the issue for me. I sent a PR.
from isce2.
Hmm, do you see anything concerning here: https://github.com/OSGeo/gdal/blob/v3.1.0/gdal/MIGRATION_GUIDE.TXT#L1-L68
Maybe you could also try gdal.UseExceptions()
to see if it gives a more helpful message?
from isce2.
It seems like the *.vrt file is not created by the gdal.Translate()
call in topo.py#L418, adding the following two lines below L418 should fix it. Could you confirm @EJFielding?
ds = gdal.Open(out_file, gdal.GA_ReadOnly)
gdal.Translate(out_file+'.vrt', ds, options=gdal.TranslateOptions(format='VRT'))
from isce2.
Related Issues (20)
- Ion parameter file
- make_single_reference_isce; how to set the value of range_looks and azimuth_looks
- run_07_pairs_misreg - Segmentation fault (core dumped) HOT 2
- how to move stack to another path or drive HOT 3
- Issue in Stamps Step (1,1) HOT 3
- IndexError in Offset.py HOT 9
- How to use stackStripMap.py to process a stack of TerraSAR-X images HOT 1
- Inconsistency in the format of the parameter "-b/--bbox"
- fetchOrbit.py Error to download POD HOT 2
- alosStack create_cmds error : startingswath HOT 1
- alosStack estimated interferogram coherence larger than 1 HOT 1
- How can change coherence chip dimention
- topsApp.py error HOT 2
- encounter buffer overflow problem when extracting ALOS level-1 data from alaska archive HOT 3
- Masking interferograms before unwrapping
- Including brief description for downloading Sentienl-1 SLC data in topsStack documentation. HOT 1
- An import error when using the function of estCpxCoherence HOT 1
- run_07_pairs_misreg ,Bad match at level 1,how to deal it? HOT 7
- runESD error
- stripmapApp.py lutan.xml --------ImportError: /lib/libgdal.so.30: undefined symbol: GEOSMakeValidWithParams_r and core dumped:Segmentation Fault HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from isce2.