I actually tried deleting this post because it seems like the file without a readable character is not the problem. Sometimes my code works, and sometimes it does not, which makes me think it is rather a case of Snaphu failing to unwrap on certain scenes but I need to investigate this further.
For those interested, here is python code so far to automatically export to snaphu, unwrap using snaphu, import the results and convert to displacement using python and snappy. It includes code to ingest 2 files into one operator and running snaphu in subprocess. This has not been tested well yet but I essentially use gpf.create product to write a “snaphu” file during snaphu export, I then find the correct line to run in the snaphu.conf file that SNAP creates, I feed this to snaphu using subprocess, I then read in the unrwapped .hdr file and write it again as BEAM_DIMAP file (This step might be unnecessary), I read that file, together with the original interferogram file into snaphu import and Phase to displacment.
Product being the filepath + filename of the filtered interferogram .dim product, temp_path being the processing directory and outpath being the required output filepath + filename - paths as pathlib objects
parameters = HashMap()
parameters_snaphu = HashMap()
parameters_snaphu.put(“targetFolder”, str(temp_path))
inputfile = snappy.ProductIO.readProduct(str(product))
result_SNE = snappy.GPF.createProduct(“SnaphuExport”, parameters_snaphu, inputfile)
snappy.ProductIO.writeProduct(result_SNE, temp_path), “Snaphu”)
infile = temp_path / “snaphu.conf”
with open(str(infile)) as lines:
line = lines.readlines()[6]
snaphu_string = line[1:].strip()
snaphu_args = snaphu_string.split()
process = subprocess.Popen(snaphu_args, cwd=str(temp_path))
process.communicate()
process.wait()
unwrapped_list = glob.glob(str(temp_path)) + “/UnwPhase*.hdr”)
unwrapped_hdr = str(unwrapped_list[0])
unwrapped_read = snappy.ProductIO.readProduct(unwrapped_hdr)
snappy.ProductIO.writeProduct(unwrapped_read, str(temp_path / “unwrapped_read.dim”), “BEAM-DIMAP”)
unwrapped = snappy.ProductIO.readProduct(str(temp_path / “unwrapped_read.dim”))
snaphu_files = snappy.jpy.array(‘org.esa.snap.core.datamodel.Product’, 2)
snaphu_files[0] = inputfile
snaphu_files[1] = unwrapped
result_SI = snappy.GPF.createProduct(“SnaphuImport”, parameters , snaphu_files)
result_PD = snappy.GPF.createProduct(“PhaseToDisplacement”, parameters , result_SI)
snappy.ProductIO.writeProduct(result_PD, str(outpath), “BEAM-DIMAP”)