Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Frame ID and TopsApp Processing #93

Closed
cmarshak opened this issue Jan 31, 2023 · 14 comments
Closed

Frame ID and TopsApp Processing #93

cmarshak opened this issue Jan 31, 2023 · 14 comments
Labels
bug Something isn't working

Comments

@cmarshak
Copy link
Collaborator

cmarshak commented Jan 31, 2023

Describe the bug
@dbekaert @sssangha @ehavazli @mgovorcin @gracebato - please let me know your thoughts.

The frame processing is taking an unreasonable amount of time particularly for topo step (see these topsapp steps). I want to make sure that I am doing things correctly with ISCE2 (or if this frame approach is possible as we currently have scoped it out 🙃).

I am specifying a region of interest that should be significantly smaller than the 3 slcs fed in. But my guess is that topo is still being performed for the 3 SLCs in the reference/secondary, which defeats the purpose of specifying the region of interest (see the image at the end of this issue ticket).

Here is the time for each of the above steps on our server leffe:

TopsApp Steps: 0%| | 0/24 [00:00<?, ?it/s]
TopsApp Steps: 4%|▍ | 1/24 [00:00<00:11, 2.00it/s]
TopsApp Steps: 8%|▊ | 2/24 [01:45<22:38, 61.77s/it]
TopsApp Steps: 12%|█▎ | 3/24 [01:48<12:15, 35.03s/it]
TopsApp Steps: 17%|█▋ | 4/24 [01:49<07:09, 21.49s/it]
TopsApp Steps: 21%|██ | 5/24 [16:43:58<114:30:38, 21696.76s/it]
TopsApp Steps: 25%|██▌ | 6/24 [16:44:02<71:36:20, 14321.11s/it]
TopsApp Steps: 29%|██▉ | 7/24 [16:44:02<45:31:10, 9639.42s/it]
TopsApp Steps: 33%|███▎ | 8/24 [16:44:03<29:12:14, 6570.90s/it]
TopsApp Steps: 38%|███▊ | 9/24 [16:44:03<18:49:14, 4516.94s/it]
TopsApp Steps: 42%|████▏ | 10/24 [16:44:04<12:08:37, 3122.67s/it]
TopsApp Steps: 46%|████▌ | 11/24 [16:44:05<7:49:32, 2167.14s/it]
TopsApp Steps: 50%|█████ | 12/24 [16:44:05<5:01:36, 1508.07s/it]
TopsApp Steps: 54%|█████▍ | 13/24 [17:20:23<5:13:41, 1711.06s/it]
TopsApp Steps: 58%|█████▊ | 14/24 [17:38:27<4:13:36, 1521.60s/it]
TopsApp Steps: 62%|██████▎ | 15/24 [17:38:30<2:39:34, 1063.82s/it]
TopsApp Steps: 67%|██████▋ | 16/24 [17:57:10<2:24:06, 1080.78s/it]
TopsApp Steps: 71%|███████ | 17/24 [17:59:14<1:32:31, 793.09s/it]
TopsApp Steps: 75%|███████▌ | 18/24 [17:59:35<56:05, 560.92s/it]
TopsApp Steps: 79%|███████▉ | 19/24 [18:08:23<45:55, 551.15s/it]
TopsApp Steps: 83%|████████▎ | 20/24 [18:08:24<25:43, 385.86s/it]
TopsApp Steps: 88%|████████▊ | 21/24 [18:23:38<27:12, 544.33s/it]
TopsApp Steps: 92%|█████████▏| 22/24 [18:23:38<12:42, 381.19s/it]
TopsApp Steps: 96%|█████████▌| 23/24 [18:23:39<04:27, 267.07s/it]
TopsApp Steps: 100%|██████████| 24/24 [18:23:40<00:00, 187.14s/it]
TopsApp Steps: 100%|██████████| 24/24 [18:23:40<00:00, 2759.19s/it]

I am running two jobs like this simultaneously on leffe (and both have the result of taking too long on topo).

To Reproduce

See these examples: https://github.com/ACCESS-Cloud-Based-InSAR/DockerizedTopsApp#using-fixed-frames-experimental

Here is the xml:

<topsApp>
  <component name="topsinsar">
    <component name="reference">
      <property name="orbit directory">orbits</property>
      <property name="auxiliary data directory">aux_cal</property>
      <property name="output directory">reference</property>
      <property name="safe">['S1A_IW_SLC__1SDV_20230113T140044_20230113T140111_046766_059B44_FBB8.zip', 'S1A_IW_SLC__1SDV_20230113T140019_20230113T140046_046766_059B44_A9C1.zip', 'S1A_IW_SLC__1SDV_20230113T135954_20230113T140021_046766_059B44_981B.zip']</property>
    </component>
    <component name="secondary">
      <property name="orbit directory">orbits</property>
      <property name="auxiliary data directory">aux_cal</property>
      <property name="output directory">secondary</property>
      <property name="safe">['S1A_IW_SLC__1SDV_20221208T140046_20221208T140113_046241_05897B_28A9.zip', 'S1A_IW_SLC__1SDV_20221208T140021_20221208T140048_046241_05897B_8EBC.zip', 'S1A_IW_SLC__1SDV_20221208T135956_20221208T140023_046241_05897B_86FA.zip']</property>
    </component>
    <property name="swaths">[1, 2, 3]</property>
    <property name="range looks">19</property>
    <property name="azimuth looks">7</property>
    <property name="filter strength">0.5</property>
    <property name="region of interest">[35.257855, 37.231403, -120.902529, -117.740922]</property>
    <property name="demFilename">full_res.dem.wgs84</property>
    <property name="geocodeDemFilename">low_res.dem.wgs84</property>
    <property name="do unwrap">True</property>
    <property name="unwrapper name">snaphu_mcf</property>
    <property name="do ionosphere correction">False</property>
    <property name="apply ionosphere correction">False</property>
    <property name="do ESD">False</property>
    <property name="ESD coherence threshold">-1</property>
    <property name="use virtual files">True</property>
    <property name="geocode list">['merged/phsig.cor', 'merged/filt_topophase.unw', 'merged/los.rdr', 'merged/topophase.flat', 'merged/filt_topophase.flat', 'merged/filt_topophase_2stage.unw', 'merged/topophase.cor', 'merged/filt_topophase.unw.conncomp']</property>
  </component>
</topsApp>

Below is the image of the SLCs (blue) and the frame in yellow.
image

@cmarshak cmarshak added the bug Something isn't working label Jan 31, 2023
@cmarshak
Copy link
Collaborator Author

My only other guess is to remove the integer truncation we do here.

The DEM might control the topo step. For the above example, the dem bounds are: BoundingBox(left=-122.0, bottom=35.0001012, right=-117.000002, top=38.0001)

@cmarshak
Copy link
Collaborator Author

Another option is to decrease the frame size per @mgovorcin's suggestion.

@ehavazli
Copy link
Collaborator

ehavazli commented Feb 1, 2023

Another option is to decrease the frame size per @mgovorcin's suggestion.

I second this idea. I remember topo step taking too long on my tests with virtual processing as well. I would stay make the yellow frame smaller enough so it doesn't touch the overlapping areas.

@cmarshak
Copy link
Collaborator Author

cmarshak commented Feb 1, 2023

It appears all the experts seem to be converging on decreasing the frame size. Thank you, @ehavazli.

Here are some relevant insights from @piyushrpt:

  1. Topo does not care so much about DEM bounds. It performs a search in 3D space for all points that project to a given 2D point - (azimuth, range). Since the dimensionality increases, it's slow. geo2rdr projects a 3D point to a 2D space - and is a lot faster due to reduction in dimensionality. If you want to solve for every point precisely, there is no option but to pay for this extra cost. If stacking you only do this once. If using GSLC you use the faster geo2rdr operation that reduces dimensionality and skip topo altogether. On the other hand, topo is trivially parallelizable and if a GPU is available - you will see the most speedup. Decreasing frame size will help as this will reduce the number of bursts you have to run topo for. If there is willingness to compromise a little on accuracy, you could skip every other pixel / 2nd pixel in each direction and reduce the number of topo operations. You can then interpolate on the fly for the pixels you skipped using bilinear / bicubic interpolation. This is what most people who don't have a lot of compute resources do. But that's something to experiment with and get approval from science users which I think will be the main challenge.
  2. topsApp will select all the bursts that intersect the region of interest at the very first preprocessing steps and will carry those bursts through the rest of the processing. topo will be run for all those bursts. You can throw as many SLCs as you want at it - only the bursts touching the region of interest will be used for analysis.
  3. You can get more efficient about cropping if you know the min/max altitudes in your region of interest already. https://github.com/isce-framework/isce2/blob/a5e662359151df52c8b76f02251837fc8aaff1b5/components/isceobj/Sensor/TOPS/BurstSLC.py#L338 Default is to use global min / max - which can cause a slightly larger extent to be chosen. Usually, this only matters when you are looking to subset by swath as well.

@cmarshak
Copy link
Collaborator Author

cmarshak commented Feb 1, 2023

We should do some rough counting to ensure frame size is suitable.

ESA Frames (images) come in packages of 9 bursts (I am grouping bursts across all 3 swaths as is done in the relative burst map). If our fixed frames are 10 bursts in size, we will get 3 images when the fixed frame covers all of an intersecting ESA image frame because an ESA image frame will have at least 1 burst overlap between its neighbors*. I would imagine similar edge cases for 9 burst frames because of issues identified below. I would say 8 burst size would prevent more than 2 SLCs from being identified at the outset.

*There are additional issues by our grouping of bursts:

  1. Our burst grouping may differ from how determines a given ESA image frame (?)
  2. We are using ASF metadata extents (not burst coverage) to determine intersections when determining relevant SLCs.
  3. We are feeding in bboxes to ISCE2 and so additional bursts will be captured during processing if included during enumeration of the IFG.

@piyushrpt
Copy link

Also, there are multiple regions of interest if I remember.

One specified in the sensor - crops things right up front and only processes the bursts that intersect the ROI.
One specified in topsApp is for geocoding only a small section of output but processes the whole thing. This was meant to geocode a small part of a large product possibly at a higher resolution etc.

Read up the docs on this. All of this is well documented in the UNAVCO classes. I haven't used this in more than 2 years.

@cmarshak
Copy link
Collaborator Author

cmarshak commented Feb 1, 2023

@piyushrpt
Copy link

I would also check topsApp.py source code. A quick glance reveals there are two variables ROI and REGION_OF_INTEREST that might be in conflict with each other. Hopefully, when you run with a number of SLCs as input - the screen log shows you that cropping is actually taking place.

@dbekaert
Copy link
Collaborator

dbekaert commented Feb 2, 2023

Im still a bit puzzled on the long runtime. From what i seen on the ISCE log file that @cmarshak showed me. The region of interest seems to work. We are extracting like 10-13 burst per subswath which is pretty much what we expect knowing we increased the overlap by 2 burst so about 4 extra than a regular SLC.

Perhaps the resources were capped, or multiple jobs were in competition? What was the number of threads/ memory availability at the time of the test?

What should not matter:

  • How we group burst ets and differences from an ESA frame etc should not be applicable to this. ISCE would always have extracted the common burst between reference and secondary and do the topo step on that. We successfully limit them based on the region of interest as confirmed by the log. The topo step then cycles through this list of extracted burst.
  • We are not using the virtual command that points to S3. The virtual processing with vrt to locally downloaded zip files had never a problem before or computational impact to topo.

@piyushrpt
Copy link

Check OMP_NUM_THREADS - the number of threads used should also be in the screen log

@cmarshak
Copy link
Collaborator Author

cmarshak commented Feb 2, 2023

Thank you again, Piyush, for all your help and insight to this.

@dbekaert - I proceeded to make the frames 8 bursts with 1 overlap and didn't really look back. As you noted, there are a number of issues that I may failed to include in the bug report:

  1. The number of people running intensive jobs on our research server when I reported this
  2. The fact that I was running two of these jobs simultaneously.

That said, I had 3 SLCs within my region of interest with 13 bursts (as I recall). Should not have been significantly more than the roughly 10-11 burst products we have seen in GUNWs.

The isce.log is an excellent place as David mentioned to confirm region_of_interest is doing what we expect. I will include a snippet from the first few steps of ISCE2 to illustrate (using the 8 frames; it appears that there are 1-2 bursts on each side of the image that are added because the region of interest is a box):

023-02-02 13:46:59,583 - isce.insar - INFO - ISCE VERSION = 2.6.1, RELEASE_SVN_REVISION = ,RELEASE_DATE = 20220811, CURRENT_SVN_REVISION = 
2023-02-02 13:46:59,612 - isce.insar - WARNING - Some filenames in insarApp.geocode_list configuration are different from those in InsarProc. Using names given to insarApp.
2023-02-02 13:47:00,013 - isce.insar - INFO - ISCE VERSION = 2.6.1, RELEASE_SVN_REVISION = ,RELEASE_DATE = 20220811, CURRENT_SVN_REVISION = 
2023-02-02 13:47:00,042 - isce.insar - WARNING - Some filenames in insarApp.geocode_list configuration are different from those in InsarProc. Using names given to insarApp.
2023-02-02 13:48:38,726 - isce.topsinsar.runPreprocessor - INFO - 
####################################################################################################
    runPreprocessor
----------------------------------------------------------------------------------------------------
common.Input list of swaths to process:  = [1, 2, 3]
reference.burstLength_1 = 1499
reference.burstLength_2 = 1511
reference.burstLength_3 = 1517
reference.burstWidth_1 = 21672
reference.burstWidth_2 = 25563
reference.burstWidth_3 = 24672
reference.numberOfBursts_1 = 12
reference.numberOfBursts_2 = 12
reference.numberOfBursts_3 = 11
reference.sensor.ascendingnodetime = 2023-01-25 13:20:59.914647
reference.sensor.family = topsswathslc
reference.sensor.folder = None
reference.sensor.imagingmode = TOPS
reference.sensor.mission = S1A
reference.sensor.name = topsswathslc_name
reference.sensor.numberofbursts = 11
reference.sensor.processingfacility = None
reference.sensor.processingsoftwareversion = 003.52
reference.sensor.processingsystem = None
reference.sensor.spacecraftname = None
secondary.burstLength_1 = 1499
secondary.burstLength_2 = 1511
secondary.burstLength_3 = 1517
secondary.burstWidth_1 = 21673
secondary.burstWidth_2 = 25563
secondary.burstWidth_3 = 24673
secondary.numberOfBursts_1 = 12
secondary.numberOfBursts_2 = 12
secondary.numberOfBursts_3 = 11
secondary.sensor.ascendingnodetime = 2022-12-20 13:21:01.322125
secondary.sensor.family = topsswathslc
secondary.sensor.folder = None
secondary.sensor.imagingmode = TOPS
secondary.sensor.mission = S1A
secondary.sensor.name = topsswathslc_name
secondary.sensor.numberofbursts = 11
secondary.sensor.processingfacility = None
secondary.sensor.processingsoftwareversion = 003.52
secondary.sensor.processingsystem = None
secondary.sensor.spacecraftname = None
####################################################################################################
2023-02-02 13:48:39,162 - isce.insar - INFO - ISCE VERSION = 2.6.1, RELEASE_SVN_REVISION = ,RELEASE_DATE = 20220811, CURRENT_SVN_REVISION = 
2023-02-02 13:48:39,191 - isce.insar - WARNING - Some filenames in insarApp.geocode_list configuration are different from those in InsarProc. Using names given to insarApp.
2023-02-02 13:48:41,119 - isce.topsinsar.runPreprocessor - INFO - 
####################################################################################################
    runComputeBaseline
----------------------------------------------------------------------------------------------------
baseline.IW-1 Bpar at midrange for first common burst = -36.04160300831371
baseline.IW-1 Bpar at midrange for last common burst = -39.4209536637765
baseline.IW-1 Bperp at midrange for first common burst = 47.352539109543734
baseline.IW-1 Bperp at midrange for last common burst = 53.76018676783484
baseline.IW-1 First common burst in reference = 0
baseline.IW-1 First common burst in secondary = 0
baseline.IW-1 Last common burst in reference = 12
baseline.IW-1 Last common burst in secondary = 12
baseline.IW-1 Number of bursts in reference = 12
baseline.IW-1 Number of bursts in secondary = 12
baseline.IW-1 Number of common bursts = 12
baseline.IW-2 Bpar at midrange for first common burst = -39.87334932713999
baseline.IW-2 Bpar at midrange for last common burst = -43.7524507360454
baseline.IW-2 Bperp at midrange for first common burst = 44.477662323577086
baseline.IW-2 Bperp at midrange for last common burst = 50.595975915712884
baseline.IW-2 First common burst in reference = 0
baseline.IW-2 First common burst in secondary = 0
baseline.IW-2 Last common burst in reference = 12
baseline.IW-2 Last common burst in secondary = 12
baseline.IW-2 Number of bursts in reference = 12
baseline.IW-2 Number of bursts in secondary = 12
baseline.IW-2 Number of common bursts = 12
baseline.IW-3 Bpar at midrange for first common burst = -43.33677388123395
baseline.IW-3 Bpar at midrange for last common burst = -47.22672798274666
baseline.IW-3 Bperp at midrange for first common burst = 42.378930259911115
baseline.IW-3 Bperp at midrange for last common burst = 47.688060313717685
baseline.IW-3 First common burst in reference = 0
baseline.IW-3 First common burst in secondary = 0
baseline.IW-3 Last common burst in reference = 11
baseline.IW-3 Last common burst in secondary = 11
baseline.IW-3 Number of bursts in reference = 11
baseline.IW-3 Number of bursts in secondary = 11
baseline.IW-3 Number of common bursts = 11
####################################################################################################
2023-02-02 13:48:41,579 - isce.insar - INFO - ISCE VERSION = 2.6.1, RELEASE_SVN_REVISION = ,RELEASE_DATE = 20220811, CURRENT_SVN_REVISION = 
2023-02-02 13:48:41,608 - isce.insar - WARNING - Some filenames in insarApp.geocode_list configuration are different from those in InsarProc. Using names given to insarApp.
2023-02-02 13:48:42,053 - isce.insar - INFO - ISCE VERSION = 2.6.1, RELEASE_SVN_REVISION = ,RELEASE_DATE = 20220811, CURRENT_SVN_REVISION = 
2023-02-02 13:48:42,082 - isce.insar - WARNING - Some filenames in insarApp.geocode_list configuration are different from those in InsarProc. Using names given to insarApp.
2023-02-02 13:48:42,389 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.909942770397513
2023-02-02 13:49:44,098 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9102072600692686
2023-02-02 13:50:45,551 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9104687140743533
2023-02-02 13:51:54,549 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.910727552870707
2023-02-02 13:53:11,934 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.910983990022114
2023-02-02 13:54:35,516 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9112378584517216
2023-02-02 13:55:44,284 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9114891893910424
2023-02-02 13:56:56,985 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9117380025334847
2023-02-02 13:58:04,875 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.911984322427526
2023-02-02 13:59:10,519 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9122281771497884
2023-02-02 14:00:16,693 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9124694133092315
2023-02-02 14:01:19,923 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.912708587790923
2023-02-02 14:02:18,437 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.910034451895691
2023-02-02 14:03:31,524 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.910297827328545
2023-02-02 14:04:56,033 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9105585699167977
2023-02-02 14:06:03,678 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9108167003216012
2023-02-02 14:07:18,625 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9110724336441547
2023-02-02 14:08:36,780 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9113252305515136
2023-02-02 14:09:33,497 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.911575686500642
2023-02-02 14:10:22,838 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.911823632817647
2023-02-02 14:11:23,921 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9120690949202923
2023-02-02 14:12:12,232 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9123119222192146
2023-02-02 14:12:41,154 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9125524981411526
2023-02-02 14:14:00,564 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9127910157299284
2023-02-02 14:14:37,975 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9103892496371646
2023-02-02 14:15:26,001 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.910648883935855
2023-02-02 14:16:03,239 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.910906300012814
2023-02-02 14:16:35,570 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.911160946308382
2023-02-02 14:17:01,611 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.911412859271944
2023-02-02 14:17:30,876 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9116624372657456
2023-02-02 14:18:02,479 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.911909697452096
2023-02-02 14:18:41,281 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.9121541167748024
2023-02-02 14:19:23,033 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.912396094254099
2023-02-02 14:20:09,724 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.912635826582291
2023-02-02 14:20:55,640 - isce.zerodop.topozero - WARNING - Default Peg heading set to: -2.912873507631279
2023-02-02 14:21:40,023 - isce.topsinsar.topo - INFO - 
####################################################################################################
    runTopo
----------------------------------------------------------------------------------------------------
topo.Bounding box for burst 0 - IW-1 = [37.25628223915817, 37.59052484462655, -118.71526398320441, -117.64319142278666]
topo.Bounding box for burst 0 - IW-2 = [37.340258864329584, 37.66913560671154, -119.67224247187424, -118.58355921173512]
topo.Bounding box for burst 0 - IW-3 = [37.24734595950941, 37.55560413731323, -120.55619916341831, -119.61166313202254]
topo.Bounding box for burst 1 - IW-1 = [37.08683623169173, 37.427344332452265, -118.7495593341275, -117.68893465451079]
topo.Bounding box for burst 1 - IW-2 = [37.17340134844576, 37.50246775959021, -119.69787166351654, -118.64074466929362]
topo.Bounding box for burst 1 - IW-3 = [37.080628949265275, 37.389973247677915, -120.58974070893065, -119.64053585339614]
topo.Bounding box for burst 10 - IW-1 = [35.59544096700589, 35.92717811836909, -119.0410802130512, -118.03522650998923]
topo.Bounding box for burst 10 - IW-2 = [35.673101617523685, 36.00714705358678, -120.00643534176227, -118.93810505035438]
topo.Bounding box for burst 10 - IW-3 = [35.58746658157326, 35.89809731662738, -120.89545293139564, -119.94971609461786]
topo.Bounding box for burst 11 - IW-1 = [35.429302169690736, 35.76114478056149, -119.07814490704337, -118.07547194026347]
topo.Bounding box for burst 11 - IW-2 = [35.50682871944659, 35.84110128167662, -120.04244170533889, -118.97651868443883]
topo.Bounding box for burst 2 - IW-1 = [36.92326572111398, 37.261120281203034, -118.78670632587546, -117.70204965302624]
topo.Bounding box for burst 2 - IW-2 = [37.006685262760605, 37.33504500387918, -119.72540204101465, -118.68659076691799]
topo.Bounding box for burst 2 - IW-3 = [36.914263730534614, 37.22405980850803, -120.62345341307636, -119.66642454749875]
topo.Bounding box for burst 3 - IW-1 = [36.753680598305955, 37.09547045301237, -118.80778572651253, -117.75552164393409]
topo.Bounding box for burst 3 - IW-2 = [36.839159757126225, 37.168753169353856, -119.75807955523909, -118.71242565044383]
topo.Bounding box for burst 3 - IW-3 = [36.74834059796568, 37.05829231286205, -120.65734315013349, -119.699737540962]
topo.Bounding box for burst 4 - IW-1 = [36.59363321030065, 36.926673121366775, -118.8475666264489, -117.77256406074748]
topo.Bounding box for burst 4 - IW-2 = [36.673558497557494, 37.002276059959335, -119.79318145275386, -118.74004341495014]
topo.Bounding box for burst 4 - IW-3 = [36.582558952941554, 36.8926645286163, -120.692509327384, -119.7353916561527]
topo.Bounding box for burst 5 - IW-1 = [36.42397171371149, 36.76146064853745, -118.87080653949117, -117.85142218098576]
topo.Bounding box for burst 5 - IW-2 = [36.50384778828086, 36.83647947105979, -119.82846459639585, -118.78827973727186]
topo.Bounding box for burst 5 - IW-3 = [36.416655306887904, 36.72704339571964, -120.73081944613695, -119.77085986252848]
topo.Bounding box for burst 6 - IW-1 = [36.258749045457456, 36.59336008485939, -118.89514378487048, -117.8663366704227]
topo.Bounding box for burst 6 - IW-2 = [36.33833475216675, 36.6705716429879, -119.86381689395351, -118.79578329354999]
topo.Bounding box for burst 6 - IW-3 = [36.25062120724144, 36.56152613434157, -120.76679595885832, -119.80635558401207]
topo.Bounding box for burst 7 - IW-1 = [36.092000397845695, 36.42584274562938, -118.93182712032636, -117.91025554731922]
topo.Bounding box for burst 7 - IW-2 = [36.17124773395679, 36.50467566470148, -119.8991774052434, -118.83380104367299]
topo.Bounding box for burst 7 - IW-3 = [36.08483333488086, 36.39664882127913, -120.79892064575351, -119.84179976544942]
topo.Bounding box for burst 8 - IW-1 = [35.92915587242853, 36.25974372278767, -118.9681170797399, -117.94721601193001]
topo.Bounding box for burst 8 - IW-2 = [36.004639374644526, 36.33877801546501, -119.93464948734017, -118.87066516114547]
topo.Bounding box for burst 8 - IW-3 = [35.919106996949424, 36.23037220455552, -120.82715978155541, -119.87723010264712]
topo.Bounding box for burst 9 - IW-1 = [35.76199677326822, 36.093191045542866, -119.00298159552801, -118.00868921868488]
topo.Bounding box for burst 9 - IW-2 = [35.83878142197146, 36.17302334814909, -119.97115187316041, -118.90166080177268]
topo.Bounding box for burst 9 - IW-3 = [35.753290905502006, 36.06437024131458, -120.86028597621406, -119.91297818185522]
topo.Dem Used = /mnt/leffe-data2/cmarshak/2023-frames-test/22438/full_res.dem.wgs84
topo.Number of common bursts IW-1 = 12
topo.Number of common bursts IW-2 = 12
topo.Number of common bursts IW-3 = 11
topo.Number of lines for burst 0 - IW-1 = 1499
topo.Number of lines for burst 0 - IW-2 = 1511
topo.Number of lines for burst 0 - IW-3 = 1517
topo.Number of lines for burst 1 - IW-1 = 1499
topo.Number of lines for burst 1 - IW-2 = 1511
topo.Number of lines for burst 1 - IW-3 = 1517
topo.Number of lines for burst 10 - IW-1 = 1499
topo.Number of lines for burst 10 - IW-2 = 1511
topo.Number of lines for burst 10 - IW-3 = 1517
topo.Number of lines for burst 11 - IW-1 = 1499
topo.Number of lines for burst 11 - IW-2 = 1511
topo.Number of lines for burst 2 - IW-1 = 1499
topo.Number of lines for burst 2 - IW-2 = 1511
topo.Number of lines for burst 2 - IW-3 = 1517
topo.Number of lines for burst 3 - IW-1 = 1499
topo.Number of lines for burst 3 - IW-2 = 1511
topo.Number of lines for burst 3 - IW-3 = 1517
topo.Number of lines for burst 4 - IW-1 = 1499
topo.Number of lines for burst 4 - IW-2 = 1511
topo.Number of lines for burst 4 - IW-3 = 1517
topo.Number of lines for burst 5 - IW-1 = 1499
topo.Number of lines for burst 5 - IW-2 = 1511
topo.Number of lines for burst 5 - IW-3 = 1517
topo.Number of lines for burst 6 - IW-1 = 1499
topo.Number of lines for burst 6 - IW-2 = 1511
topo.Number of lines for burst 6 - IW-3 = 1517
topo.Number of lines for burst 7 - IW-1 = 1499
topo.Number of lines for burst 7 - IW-2 = 1511
topo.Number of lines for burst 7 - IW-3 = 1517
topo.Number of lines for burst 8 - IW-1 = 1499
topo.Number of lines for burst 8 - IW-2 = 1511
topo.Number of lines for burst 8 - IW-3 = 1517
topo.Number of lines for burst 9 - IW-1 = 1499
topo.Number of lines for burst 9 - IW-2 = 1511
topo.Number of lines for burst 9 - IW-3 = 1517
topo.Number of pixels for bursts 0 - IW-1 = 21672
topo.Number of pixels for bursts 0 - IW-2 = 25563
topo.Number of pixels for bursts 0 - IW-3 = 24672
topo.Number of pixels for bursts 1 - IW-1 = 21672
topo.Number of pixels for bursts 1 - IW-2 = 25563
topo.Number of pixels for bursts 1 - IW-3 = 24672
topo.Number of pixels for bursts 10 - IW-1 = 21672
topo.Number of pixels for bursts 10 - IW-2 = 25563
topo.Number of pixels for bursts 10 - IW-3 = 24672
topo.Number of pixels for bursts 11 - IW-1 = 21672
topo.Number of pixels for bursts 11 - IW-2 = 25563
topo.Number of pixels for bursts 2 - IW-1 = 21672
topo.Number of pixels for bursts 2 - IW-2 = 25563
topo.Number of pixels for bursts 2 - IW-3 = 24672
topo.Number of pixels for bursts 3 - IW-1 = 21672
topo.Number of pixels for bursts 3 - IW-2 = 25563
topo.Number of pixels for bursts 3 - IW-3 = 24672
topo.Number of pixels for bursts 4 - IW-1 = 21672
topo.Number of pixels for bursts 4 - IW-2 = 25563
topo.Number of pixels for bursts 4 - IW-3 = 24672
topo.Number of pixels for bursts 5 - IW-1 = 21672
topo.Number of pixels for bursts 5 - IW-2 = 25563
topo.Number of pixels for bursts 5 - IW-3 = 24672
topo.Number of pixels for bursts 6 - IW-1 = 21672
topo.Number of pixels for bursts 6 - IW-2 = 25563
topo.Number of pixels for bursts 6 - IW-3 = 24672
topo.Number of pixels for bursts 7 - IW-1 = 21672
topo.Number of pixels for bursts 7 - IW-2 = 25563
topo.Number of pixels for bursts 7 - IW-3 = 24672
topo.Number of pixels for bursts 8 - IW-1 = 21672
topo.Number of pixels for bursts 8 - IW-2 = 25563
topo.Number of pixels for bursts 8 - IW-3 = 24672
topo.Number of pixels for bursts 9 - IW-1 = 21672
topo.Number of pixels for bursts 9 - IW-2 = 25563
topo.Number of pixels for bursts 9 - IW-3 = 24672
topo.Overall bounding box = [35.429302169690736, 37.66913560671154, -120.89545293139564, -117.64319142278666]

@cmarshak
Copy link
Collaborator Author

cmarshak commented Feb 2, 2023

I am going to close this issue because we are using 8 frames to ensure at most 2 SLCs in our interferogram. It might be wise to investigate this further should this come back up, but this will serve as a helpful reference for such issues.

@cmarshak cmarshak closed this as completed Feb 2, 2023
@dbekaert
Copy link
Collaborator

dbekaert commented Feb 3, 2023

@cmarshak Would still be instructive for your fame-code to know if there are any runtime issues. Once you have settled on the frames could you run it without competing resources, ensuring we set the number of treads correctly available etc. As you mentioned i would not expect much more than 10-15% increase in processing time as its only slightly larger.

We should do a back of the envelope calculation to ensure the unwrapping with snaphu has enough memory available (It was one of the reasons to drop to 90m and frame sized product at the time)

@cmarshak
Copy link
Collaborator Author

cmarshak commented Feb 3, 2023

So the topo issue was resolved in the latest runs. Again, I may have opened the issue ticket prematurely. I think there is a lot of valuable information in the exercise.

Can continue the frame discussion here (included some of the info and more there): #94

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

4 participants